Displaying 1 result from an estimated 1 matches for "exprmat".
Did you mean:
exprmap
2007 Dec 26
1
nlme package
...the syntax of the NLME package that I am not aware.
This is the R script:
library(nlme)
treat=as.factor(c(1,2,1,2,1,2,1,2))
mouse=as.factor(c(1,1,2,2,3,3,4,4))
time=as.factor(c(1,1,1,1,2,2,2,2))
mouse.anova=function(vec){lme(vec~treat+time+treat*time,random=~1|mouse/time)}
lme.out=apply(exprmat,1,mouse.anova)
This is the ERROR message:
Error in lme.formula(vec ~ treat + time + treat * time, random = ~1 | :
nlminb problem, convergence error code = 1
message = singular convergence (7)
I will appreciate your help in pointing mistakes in this R script.
Thank you very much.
Roger...