John Logsdon j.logsdon at quantex-research.com writes:
I am trying to extract the solution from a simple lme calculation.
y-c(-1.118,-.5,.5,1.118,10)
gp-factor(c(rep('one',4),'two'))
res-lme(y~1,rand=~1|gp)
...
Random effects:
Formula: ~1 | gp
(Intercept) Residual
I am trying to extract the solution from a simple lme calculation.
For example (the first 4 have a mean 0, sd 1):
y-c(-1.118,-.5,.5,1.118,10)
gp-factor(c(rep('one',4),'two'))
res-lme(y~1,rand=~1|gp)
Linear mixed-effects model fit by REML
Data: NULL
Log-restricted-likelihood: -8.67141