search for: adjdata

...,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 <- intervals(ge.lme...