[R] lme4 model with no fixed effects?
Frank Samuelson
gmane0806.13.cudgel at neverbox.com
Tue Nov 21 20:22:51 CET 2006
Crossed random effects:
> lmer( y ~ (1 | i1) + (1|i2) ,data=dta)
Linear mixed-effects model fit by REML
Formula: y ~ (1 | i1) + (1 | i2)
Data: dta
AIC BIC logLik MLdeviance REMLdeviance
91.18 94.84 -42.59 85.2 85.18
Random effects:
Groups Name Variance Std.Dev.
i1 (Intercept) 0.068224 0.26120
i2 (Intercept) 0.463112 0.68052
Residual 1.469250 1.21213
number of obs: 25, groups: i1, 5; i2, 5
Fixed effects:
Estimate Std. Error t value
(Intercept) 0.09806 0.40625 0.2414
Crossed random effects with no fixed effects?
> lmer( y ~ -1 + (1 | i1) + (1|i2) ,data=dta)
Error in lmer(y ~ -1 + (1 | i1) + (1 | i2), data = dta) :
Xp must be a real matrix
If this isn't possible, it is easy enough to throw
together a slow ML estimator myself.
Thanks for any help.
-Frank
More information about the R-help
mailing list