Problem with ci.lmer() in package:gmodels

Previous Topic Next Topic
 
classic Classic list List threaded Threaded
2 messages Options
Reply | Threaded
Open this post in threaded view
|

Problem with ci.lmer() in package:gmodels

Michael Kubovy
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/

______________________________________________
[hidden email] 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.
Reply | Threaded
Open this post in threaded view
|

Re: 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/
[hidden email] 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.

______________________________________________
[hidden email] 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.