Displaying 2 results from an estimated 2 matches for "cloneid".
Did you mean:
cloned
2003 May 14
1
lme speedup question
...00){
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,level=0.95,which="fixed")
cp1[j,2] <- tmpInt$fixed[2,2]*2
cp1[j,3] <- tmpInt$fixed[2,1]
cp1[j,4] <- tmpInt$fixed[2,3]
}
}
An Aside: I am using lme for a mixed model with treating rp as random and nested in the fixed effect vg. I am col...
2007 Nov 13
0
Difficult to set a quiet formula in maanova
...t;,sep="\t")
> todos.norm <- vsn2(todos)
> write.table(exprs(todos.norm),"todos.norm.txt",sep="\t")
> library(maanova)
> fabiana.raw <- read.madata("todos.norm.vsn2.maanova.txt", designfile="design.txt", header=TRUE, spotflag=FALSE,CloneID=1,metarow=2, metacol=3, pmt=4)
> fabiana <- createData(fabiana.raw, n.rep=2, avgreps=1, log.trans=FALSE)
> model.full.mix <- makeModel(data=fabiana, formula=~Var+Trat+Time+Sample+Var:Trat+Var:Time+Trat:Time+Var:Trat:Time, random=~Sample)
> summary(model.full.mix)
Model Summa...