Displaying 7 results from an estimated 7 matches for "rm1".
Did you mean:
rc1
2012 Jul 19
1
Change log(J) to log(J+1) to stop log(0) from occurring in harModel
...as wondering if any one could tell me how
how I can achieve this as I am very naive with R still.
# Get the matrix for estimation of linear model
maxp = max(periods,periodsJ); #max number of aggregation levels
if(!is.null(leverage)){ maxp = max(maxp,leverage) }
n = length(RM1); #Number of Days
# Aggregate RV:
RVmatrix1 = aggRV(RM1,periods);
if( nest==2 ){ RVmatrix2 = aggRV(RM2,periods); } # In case a jumprobust
estimator is supplied
# Aggregate and subselect y:
y = aggY(RM1,h,maxp);
# Only keep useful parts:
x1 = RVmatrix1[(maxp:(n-h)),];...
2005 Oct 06
2
data.frame error using sem package
I keep getting this error when I try to use the sem package. I and
another person who has successfully used the sem package for similar
analysis (fMRI effective connectivity) cannot figure out what is
wrong with my code. I would appreciate any suggestions.
The error message:
Error in data.frame(object$coeff, se, z, 2 * (1 - pnorm(abs(z))),
par.code) :
arguments imply differing
2009 Apr 04
1
comparing columns in a dataframe
...; a file with 130 rows
Below is the code that I have started with.
Ideally I would like to identify the actual row numbers where there
is difference in the variables (either pairwise or between 3 variables).
x <- read.csv("c://rec_compare.csv",header=T, as.is=TRUE)
u <- union(x$rm1, x$redc1)
write.csv(u,"c:/union_test.csv")
i <- intersect(x$rm1, x$redc1)
write.csv(i,"c:/intersect_test.csv")
sd <- setdiff(x$rm1, x$redc1)
write.csv(sd,"c:/setdiff_test.csv")
Any suggestions are appreciated.
regards
Bob
2005 Apr 12
1
R Package: mmlcr and/or flexmix
...BetweenStays+ CROStays+ DaysSinceLastStay+
wghtmean_median_age_pop100+ perc_num_white+ perc_num_hispanic+
perc_CROStays+ numMonthsActive+ WEBStays+ property_loyalty+ perc_NoZone+
rho+ PCR_Dummy_class+ ltgold1+ ltgold3+ p_hat_PCR, data = data, k = 2, model
= FLXglm(family = "poisson"))
rm1<-refit(m1)
summary(rm1)
Scott Geller
Advanced Analytics,
Decision Sciences Department,
InterContinental Hotels Group
770-604-5149
2012 Jul 19
1
Switching log(J) to log(J+1) to avoid log(0) in HAR-RVJ model
...uot;);
return( model )
}#End HAR-RV-J if cond
and this
if( type == "HARRVCJ" ){
# Are the jumps significant? if not set to zero:
if( jumptest=="ABDJumptest" ){
TQ = apply.daily(data, TQfun);
J = J[,1];
teststats = ABDJumptest(RV=RM1,BPV=RM2,TQ=TQ );
}else{ jtest = match.fun(jumptest); teststats = jtest(data,...) }
Jindicators = teststats > qnorm(1-alpha);
*J[!Jindicators] = 0;*
to
if( type == "HARRVCJ" ){
# Are the jumps significant? if not set to zero:
if( jumptest=="ABDJumptest...
2008 Jul 23
5
Histogram
Hi,
how can I treat data organised in classes and frequencies?
Ex.
class frequency
20-23 9
23-25 7
26-28 5
29-31 5
32-34 3
Thanks
Angelo Scozzarella
2012 Oct 23
10
How to pick colums from a ragged array?
I have a large dataset (~1 million rows) of three variables: ID (patient's name), DATE (of appointment) and DIAGNOSIS (given on that date).
Patients may have been assigned more than one diagnosis at any one appointment - leading to two rows, same ID and DATE but different DIAGNOSIS.
The diagnoses may change between appointments.
I want to subset the data in two ways:
- define groups