Displaying 3 results from an estimated 3 matches for "originaldata".
2012 Oct 29
2
naming datasubsets in a loop
...reality I have got
about 70 subsets I want to realize this in a loop
Somehow like this:
names=c("nam1",
"nam2","nam3","nam4","nam5","nam6","nam7","nam8","nam9","nam10")
for (i in 1:10){
subset(originaldata, originaldata[ ,130]<=time[i]) #This subset should
get the name safed in "names" on position i
}
I don't know how to realized it. None of my attempts worked,- so I
would be grateful for help.
Best regards
Claudia
2011 May 17
1
simprof test using jaccard distance
...ind(hclust.results$merge, matrix(data = NA,
nrow = nrow(hclust.results$merge), ncol = 1))
simprof.results <- simprof.body(rawdata = rawdata, num.expected = num.expected,
num.simulated = num.simulated, method.cluster = method.cluster,
method.distance = method.distance, originaldata = rawdata,
alpha = alpha, clust.order = hclust.results$merge, startrow = nrow(hclust.results$merge),
pMatrix = pMatrix, const = const, silent = silent, increment = increment)
results <- list()
results[["numgroups"]] <- length(simprof.results$samples)
resu...
2012 Mar 09
0
pdMat class in LME to mimic SAS proc mixed group option? Group-specific random slopes
...the 5 subgroups (as well as the residual R-side variances).
If I'm willing to assume that each subject's random slope and intercepts is independent (a dangerous assumption), then I can get close to the model using the following code:
#R CODE MODEL1
G.Data<-groupedData(y~1|Subject,data=OriginalData)
G.Data$subgroup<-as.factor(G.Data$subgroup)
pd<-pdDiag(~subgroup+subgroup:time-1)
# subgroup is the subgroup-specific intercept and subgroup:time is the subgroup-specific intercept.
model1<-lme(fixed=y~x*time,
random=pd,
weights=varIdent(form=~1|subgroup),...