[R-sig-ME] models with no fixed effects
farewelld at Cardiff.ac.uk
Thu Sep 11 15:03:18 CEST 2008
I'm running into an error when using lmer to fit models with no fixed effects terms.
For example, generating some data with
df$y <- with(df <- data.frame(i = gl(5, 5), b = rep(rnorm(5), each = 5)), b + rnorm(25))
and fitting like this
fit1 <- lmer(y ~ 1 + (1 | i), df)
works fine. But fitting like this
fit0 <- lmer(y ~ 0 + (1 | i), df)
gives the following error:
CHOLMOD error: Pl?
Error in mer_finalize(ans) :
Cholmod error `invalid xtype' at file:../Cholesky/cholmod_solve.c, line 971
Am I missing something obvious?
Here's my sessionInfo(), in case it's useful:
R version 2.7.1 (2008-06-23)
attached base packages:
 stats graphics grDevices utils datasets methods base
other attached packages:
 lme4_0.999375-26 Matrix_0.999375-15 lattice_0.17-8
loaded via a namespace (and not attached):
More information about the R-sig-mixed-models