Displaying 1 result from an estimated 1 matches for "mlma".
Did you mean:
lma
2004 Jan 22
0
problem fitting linear mixed models
...1.5320000000000000,11.7100000000000010,11.3519999999999990,11.0000000000000000,10.6300000000000010,
10.8720000000000020,10.6133333333333350)
Gene<-data.frame(Subject=as.factor(as.character(Subject)),Tissue=as.factor(as.character(Tissue)),Intensity=IntensityA)
# fitted model (does not work!)
mlmA<-lme(Intensity ~
Tissue-1,weights=varIdent(form=~1|Tissue),correlation=corSymm(),data=Gene,random=~1|Subject);
summary(mlmA)
#sample means
lapply(split(Gene$Intensity,Gene$Tissue),mean)
#fitted model with rounded data (it works)
Gene$Intensity<-round(Gene$Intensity,4)
mlmA<-lme(I...