Hej,

Calculation of  confidence intervals for  means
based on a model fitted with lmer 
using the package multcomp

- yields results for  calpha=adjusted_calpha
- NA's for            calpha=univariate_calpha


Example:
library(lme4)
library(multcomp)
###  Generate data
set.seed(8)
d<-expand.grid(treat=1:2,block=1:3)
e<-rnorm(3)
names(e)<-1:3
d$y<-rnorm(nrow(d)) + e[d$block] 
d<-transform(d,treat=factor(treat),block=factor(block))
##### lmer fit 
Mod<-lmer(y~treat+ (1|block), data=d)
### estimate  treatment means
L<-cbind(c(1,0),c(0,1))
s<-glht(Mod,linfct=L)
## confidence intervals
confint(s,calpha=adjusted_calpha())
#produces NA's for the confidence limits
confint(s,calpha=univariate_calpha())

#for  models fitted with lm the problem does not occur
G<-lm(y~treat+ block, data=d)
L<-matrix( c(1,0,1/3,1/3,1,1,1/3,1/3),2,4,byrow=TRUE)
s<-glht(G,linfct=L)
confint(s,calpha=adjusted_calpha())
confint(s,calpha=univariate_calpha())



multcomp  version 1.2-7
R:platform       i386-pc-mingw32                             
version.string R version 2.13.1 Patched (2011-08-19 r56767)


Regards

Ulrich Halekoh
Department of Molecular Biology and Genetics,
Aarhus University, Denmark
ulrich.hale...@agrsci.dk

______________________________________________
R-help@r-project.org 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 via email to