[R-sig-ME] lme4 and confidence intervals

Ben Bolker bbolker at gmail.com
Wed Oct 23 21:04:28 CEST 2013


   This works for me with the development version of lme4 (although I do
get some warnings about cov2cor() that have been referenced on another
recent thread, and the confidence intervals for the Insitution RE are
not getting calculated).  I strongly suspect it will also work with the
release version, although I haven't tested it.

  Can you start with a clean session and give the results of
sessionInfo() ??


On 13-10-23 01:03 PM, Nicholas Mitsakakis wrote:
> Hi Ben,
> 
> I am attaching a dataset, brazil. The code I used is
> 
> I <- which(as.numeric(brazil$state) == 1 & as.numeric(brazil$drug) == 1)
> modref2 <- lmer(log.Std.Unit.Price ~ 1 + (1|Institution) + (1|supplier)
> + log(Quantity) + as.numeric(Date), brazil, subset = I)
> confint(modref2)
> 
> 
> Please let me know if you are getting the same problem.
> 
> Thanks a lot,
> Nicholas
> 
> 
> 
> 
> On Wed, Oct 23, 2013 at 12:51 PM, Ben Bolker <bbolker at gmail.com
> <mailto:bbolker at gmail.com>> wrote:
> 
>      [cc'ing to r-sig-mixed models]
> 
>       This doesn't happen for me: I would need a reproducible example in
>     order to be able to see what's going on.
> 
>        sincerely
>        Ben Bolker
> 
> 
>     On 13-10-23 12 <tel:13-10-23%2012>:31 PM, Nicholas Mitsakakis wrote:
>     > Hello Bob,
>     >
>     > I want to get the (or some) confidence intervals for the coefficients
>     > from a mixed effects model fitted by the lmer function. I noticed
>     in the
>     > lme4 documentation that confint function can do this but it does not
>     > work for me giving errors messages:
>     >
>     > Error: $ operator not defined for this S4 class
>     >
>     > I have installed the most recent version of lme4 from CRAN.
>     >
>     > Am I missing something?
>     >
>     > I would appreciate it if you could let me know.
>     >
>     > Regards,
>     > nicholas
>     >
>     > --
>     > Nicholas Mitsakakis, MSc, PhD
>     > Biostatistician
>     > Toronto Health Economics and Technology Assessment (THETA)
>     Collaborative
>     > Assistant Professor
>     > Leslie Dan Faculty of Pharmacy, University of Toronto
>     > 6th Floor, Room 641
>     > 144 College Street
>     > Toronto ON, M5S 3M2
>     > e: n.mitsakakis at theta.utoronto.ca
>     <mailto:n.mitsakakis at theta.utoronto.ca>
>     <mailto:n.mitsakakis at theta.utoronto.ca
>     <mailto:n.mitsakakis at theta.utoronto.ca>>
>     > t: (416) 946 - 3700 <tel:%28416%29%20946%20-%203700>
>     > f: (416) 946 - 3719 <tel:%28416%29%20946%20-%203719>
>     > w: www.theta.utoronto.ca <http://www.theta.utoronto.ca>
>     <http://www.theta.utoronto.ca>
> 
> 
> 
> 
> -- 
> Nicholas Mitsakakis, MSc, PhD
> Biostatistician
> Toronto Health Economics and Technology Assessment (THETA) Collaborative
> Assistant Professor
> Leslie Dan Faculty of Pharmacy, University of Toronto
> 6th Floor, Room 641
> 144 College Street
> Toronto ON, M5S 3M2
> e: n.mitsakakis at theta.utoronto.ca <mailto:n.mitsakakis at theta.utoronto.ca>
> t: (416) 946 - 3700
> f: (416) 946 - 3719
> w: www.theta.utoronto.ca <http://www.theta.utoronto.ca>



More information about the R-sig-mixed-models mailing list