From: Douglas Bates (bates@stat.wisc.edu)
Date: Thu 29 Apr 2004 - 04:16:04 EST
Message-id: <6roepcq87v.fsf@bates4.stat.wisc.edu>
Deepayan Sarkar and I have a source package of release candidate 1 of
the 0.6 series of the lme4 package available at
http://www.stat.wisc.edu/~bates/lme4_0.6-0-1.tar.gz
This package requires Matrix_0.8-6 which has been uploaded to CRAN and
should be available in a few days. A copy of the source package is
available as
http://www.stat.wisc.edu/~bates/Matrix_0.8-6.tar.gz
Although this version of lme4 passes "R CMD check" on our GNU/Linux
systems we have not uploaded it to CRAN because it still lacks
capabilities that are available in lme4_0.5-2, which is currently on
CRAN. As soon as we have all the capabilities of the 0.5 series
available in the 0.6 series we will release lme4_0.6-1.tar.gz to CRAN.
This version of lme4 is a complete rewrite of the data structures and
algorithms for fitting linear mixed models. An incomplete draft
version of a paper describing the methods is available as a vignette.
Subsequent releases will contain a more polished version of this
paper.
The big change relative to earlier versions is that you can fit models
with crossed random effects quickly and easily. For example, using
the data on Scottish secondary school students achievement scores
(from http://multilevel.ioe.ac.uk/softrev/) we can fit a model with
random effects for both the secondary and the primary school attended
as
> library(lme4)
This package is in development. For production work use
lme from package nlme or glmmPQL from package MASS.
> data(ScotsSec)
> fm1 = lme(attain ~ verbal*sex, ScotsSec, random=list(primary=~1,second=~1))
> gc();system.time(lme(attain ~ verbal*sex, ScotsSec, random=list(primary=~1,second=~1)))
used (Mb) gc trigger (Mb)
Ncells 701438 18.8 1166886 31.2
Vcells 267929 2.1 786432 6.0
[1] 0.1 0.0 0.1 0.0 0.0
> summary(fm1)
Linear mixed-effects model fit by REML
Fixed: formula
AIC BIC logLik
14882.32 14925.32 -7434.162
Random effects:
Groups Name Variance Std.Dev.
primary (Intercept) 0.275458 0.52484
second (Intercept) 0.014748 0.12144
Residual 4.2531 2.0623
Fixed effects:
Estimate Std. Error DF t value Pr(>|t|)
(Intercept) 5.9147e+00 7.6795e-02 3431 77.0197 < 2e-16 ***
verbal 1.5836e-01 3.7872e-03 3431 41.8136 < 2e-16 ***
sexF 1.2155e-01 7.2413e-02 3431 1.6786 0.09332 .
verbal:sexF 2.5929e-03 5.3885e-03 3431 0.4812 0.63041
--- Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1Correlation of Fixed Effects: (Intr) verbal sexF verbal 0.177 sexF -0.482 -0.178 verbal:sexF -0.122 -0.680 0.161
Number of Observations: 3435 Number of Groups: primary second 148 19
There are other examples in the tests subdirectory.
The lme function behaves as previously *with one exception*. In the model specification there is no longer any distinction between crossed or nested or partially crossed random effects. This means that for nested random effects you must ensure that every inner grouping corresponds to a unique level of the inner grouping factor. For example, in the Pixel data there are two grouping factors, Dog and Side with Side nested within Dog. You must create a new grouping factor, say DS, with unique levels for each Dog/Side combination to be able to specify a model of "Side within Dog".
_______________________________________________ R-packages mailing list R-packages@stat.math.ethz.ch https://www.stat.math.ethz.ch/mailman/listinfo/r-packages
This archive was generated by hypermail 2.1.3 : Wed 27 Oct 2004 - 11:05:51 EST