[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