[R] Problem with ci.lmer() in package:gmodels

[Gregory Warnes]

Hello Michael,

I checked the source code, and the lower and upper confidence interval endpoints were simply reversed.  The calculations themselves are correct.

I have uploaded a new version of the gmodels package to the CRAN repository, and it should show up in a couple of days.  In the mean time, you can either simply reverse them manually, or modify the source file est.lmer to swap the two endpoints.   

The new package also fixes the two places where my email address wasn't properly updated.

-Greg


From: Michael Kubovy <kubovy_at_virginia.edu>
Date: Fri 09 Mar 2007 - 10:08:28 GMT


Dear Friends,

Please note that in the following CI lower > CI higher:

> require(lmer)
> require(gmodels)
> fm2 <- lmer(Reaction ~ Days + (1|Subject) + (0+Days|Subject),
sleepstudy)
> ci(fm2)

              Estimate  CI lower   CI upper Std. Error p-value
(Intercept) 251.66693 266.06895 238.630280   7.056447       0
Days         10.52773  13.63372   7.389946   1.646900       0



_____________________________


Professor Michael Kubovy
University of Virginia
Department of Psychology

USPS:     P.O.Box 400400    Charlottesville, VA 22904-4400
Parcels:    Room 102        Gilmer Hall
         McCormick Road    Charlottesville, VA 22903
Office:    B011    +1-434-982-4729
Lab:        B019    +1-434-982-4751
Fax:        +1-434-982-4766


WWW: http://www.people.virginia.edu/~mk9y/
R-help at stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.