Dear all, I have problems with the lmer output (see below)(2.6.2, WindowsXP), it is incomplete (no DFs, no p-values)and also generating an anova-table doesn?t work out. What is the fault? Does anybody know? This example is about plant height within different treatments (each with 5 replicates) measured several times. I?m interested in the treatment effect.> pap<-groupedData(heig~meas|box,outer=~treat,pap) > mod<-lmer(heig~treat+(meas|box),pap) > summary(mod)Linear mixed-effects model fit by REML Formula: heig ~ treat + (meas | box) Data: pap AIC BIC logLik MLdeviance REMLdeviance 134.1 156.6 -58.05 97.14 116.1 Random effects: Groups Name Variance Std.Dev. Corr box (Intercept) 0.621046 0.78806 meas 0.790570 0.88914 -1.000 Residual 0.037141 0.19272 number of obs: 90, groups: box, 30 Fixed effects: Estimate Std. Error t value (Intercept) 1.79473 0.08353 21.487 treatB1 0.23334 0.11813 1.975 treatB2 0.62677 0.11813 5.306 treatD1 0.07911 0.11813 0.670 treatD2 0.43942 0.11813 3.720 treatU -0.11933 0.11813 -1.010 Correlation of Fixed Effects: (Intr) tretB1 tretB2 tretD1 tretD2 treatB1 -0.707 treatB2 -0.707 0.500 treatD1 -0.707 0.500 0.500 treatD2 -0.707 0.500 0.500 0.500 treatU -0.707 0.500 0.500 0.500 0.500> anova(mod)Analysis of Variance Table Df Sum Sq Mean Sq treat 5 2.11106 0.42221 Thanks a lot in advance Maria Meyen Department of Conservation Biology Faculty of Biology Philipps-University of Marburg Karl-von-Frisch-Str.8 35043 Marburg __________________________________________________________ Der Lieblings-Mailbox der Welt.
Maria, you might get an answer in this thread, from Douglas Bates "lmer, p-values and all that", dated May 19 2006 HTH David Maria Meyen <maria_meyen at yahoo.de> ha escrito:> Dear all, > > I have problems with the lmer output (see > below)(2.6.2, WindowsXP), it is incomplete (no DFs, no > p-values)and also generating an anova-table doesn?t > work out. What is the fault? Does anybody know? > > This example is about plant height within different > treatments (each with 5 replicates) measured several > times. I?m interested in the treatment effect. > >> pap<-groupedData(heig~meas|box,outer=~treat,pap) >> mod<-lmer(heig~treat+(meas|box),pap) >> summary(mod) > Linear mixed-effects model fit by REML > Formula: heig ~ treat + (meas | box) > Data: pap > AIC BIC logLik MLdeviance REMLdeviance > 134.1 156.6 -58.05 97.14 116.1 > Random effects: > Groups Name Variance Std.Dev. Corr > box (Intercept) 0.621046 0.78806 > meas 0.790570 0.88914 -1.000 > Residual 0.037141 0.19272 > number of obs: 90, groups: box, 30 > > Fixed effects: > Estimate Std. Error t value > (Intercept) 1.79473 0.08353 21.487 > treatB1 0.23334 0.11813 1.975 > treatB2 0.62677 0.11813 5.306 > treatD1 0.07911 0.11813 0.670 > treatD2 0.43942 0.11813 3.720 > treatU -0.11933 0.11813 -1.010 > > Correlation of Fixed Effects: > (Intr) tretB1 tretB2 tretD1 tretD2 > treatB1 -0.707 > treatB2 -0.707 0.500 > treatD1 -0.707 0.500 0.500 > treatD2 -0.707 0.500 0.500 0.500 > treatU -0.707 0.500 0.500 0.500 0.500 > >> anova(mod) > Analysis of Variance Table > Df Sum Sq Mean Sq > treat 5 2.11106 0.42221 > > Thanks a lot in advance > > > > Maria Meyen > Department of Conservation Biology > Faculty of Biology > Philipps-University of Marburg > Karl-von-Frisch-Str.8 > 35043 Marburg > > > > __________________________________________________________ > > Der Lieblings-Mailbox der Welt. > > ______________________________________________ > R-help at 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. > >
Hi Maria, the output is not incomplete. See the FAQ here: http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-are-p_002dvalues-not-displayed-when-using-lmer_0028_0029_003f and the lengthy discussions about p-values and lmer (try Google). Andrew On Tue, Apr 15, 2008 at 10:51:46AM +0200, Maria Meyen wrote:> Dear all, > > I have problems with the lmer output (see > below)(2.6.2, WindowsXP), it is incomplete (no DFs, no > p-values)and also generating an anova-table doesn?t > work out. What is the fault? Does anybody know? > > This example is about plant height within different > treatments (each with 5 replicates) measured several > times. I?m interested in the treatment effect. > > > pap<-groupedData(heig~meas|box,outer=~treat,pap) > > mod<-lmer(heig~treat+(meas|box),pap) > > summary(mod) > Linear mixed-effects model fit by REML > Formula: heig ~ treat + (meas | box) > Data: pap > AIC BIC logLik MLdeviance REMLdeviance > 134.1 156.6 -58.05 97.14 116.1 > Random effects: > Groups Name Variance Std.Dev. Corr > box (Intercept) 0.621046 0.78806 > meas 0.790570 0.88914 -1.000 > Residual 0.037141 0.19272 > number of obs: 90, groups: box, 30 > > Fixed effects: > Estimate Std. Error t value > (Intercept) 1.79473 0.08353 21.487 > treatB1 0.23334 0.11813 1.975 > treatB2 0.62677 0.11813 5.306 > treatD1 0.07911 0.11813 0.670 > treatD2 0.43942 0.11813 3.720 > treatU -0.11933 0.11813 -1.010 > > Correlation of Fixed Effects: > (Intr) tretB1 tretB2 tretD1 tretD2 > treatB1 -0.707 > treatB2 -0.707 0.500 > treatD1 -0.707 0.500 0.500 > treatD2 -0.707 0.500 0.500 0.500 > treatU -0.707 0.500 0.500 0.500 0.500 > > > anova(mod) > Analysis of Variance Table > Df Sum Sq Mean Sq > treat 5 2.11106 0.42221 > > Thanks a lot in advance > > > > Maria Meyen > Department of Conservation Biology > Faculty of Biology > Philipps-University of Marburg > Karl-von-Frisch-Str.8 > 35043 Marburg > > > > __________________________________________________________ > > Der Lieblings-Mailbox der Welt. > > ______________________________________________ > R-help at 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.-- Andrew Robinson Department of Mathematics and Statistics Tel: +61-3-8344-6410 University of Melbourne, VIC 3010 Australia Fax: +61-3-8344-4599 http://www.ms.unimelb.edu.au/~andrewpr http://blogs.mbs.edu/fishing-in-the-bay/