Displaying 1 result from an estimated 1 matches for "philmaanova".
2003 May 14
1
lme speedup question
...1,2,1,1,2,1,2,2,1,2,1,1,2,1,2,2,1,2,1)
ag <- c(1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4,5,5,5,5,6,6,6,6)
dy <- c(1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2)
rp <- c(1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3)
for (perm in 1:216){
vg <- vgAll[perm,]
for(j in 1:15000){
ge <- c(philmaanova.rloess$adjdata[j,1:24])
dat <- data.frame(ge,vg=factor(vg),ag=factor(ag),dy=factor(dy),rp=factor(rp))
dat$vgrp <- getGroups(dat, form = ~ 1|vg/rp, level = 2)
ge.lme <- lme(fixed=ge~vg+ag+dy, data=dat, random=~1|vgrp)
cp1[j,1] <- philmaanova.rloess$cloneid[j]
tmpInt <- i...