Displaying 1 result from an estimated 1 matches for "proclme".
Did you mean:
problme
2007 Nov 12
1
R - lme
....639
But am not able to get results from the lme model as below. (The model also seems to be wrong with lot of NaNs produced)
How do I frame the model such that I can get only ~1 | Batch %in% Supplier
Or something else is wrong
Please tell me where I have gone wrong
library(nlme)
proclme=lme(Purity~Supplier,random = ~1|Supplier/Batch,process)
summary(proclme)
VarCorr(proclme)
Linear mixed-effects model fit by REML
Data: process
AIC BIC logLik
154.8440 163.823 -71.42198
Random effects:
Formula: ~1 | Supplier
(Intercept)
StdDev: 1.250514
Formula: ~1...