Re: [R] why incomplete lmer output?

From: DAVID ARTETA GARCIA <darteta001_at_ikasle.ehu.es>
Date: Tue, 15 Apr 2008 11:03:53 +0200

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.
>
>



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. Received on Tue 15 Apr 2008 - 09:06:56 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 15 Apr 2008 - 09:30:30 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive