Hello,

 

I’m trying to perform a REML analysis using the lmer() function (lme4
package). Well, it seems to work well, except that I’m not getting any
p-value (see example below). Can someone tell me what I did wrong?

 

Thanks for your help,

 

Amélie

 

> library(gdata)

> dive <- read.xls("C:/Documents and Settings/Amelie/My Documents/Postdoc/CE
2005-2006/divebydive.xls", perl="C:/perl/bin/perl.exe")

> library(lme4)

Loading required package: Matrix

Loading required package: lattice

> reml.res <- lmer(UNDS~SUCCESSMN+(1|BIRD), dive)

> summary(reml.res)

Linear mixed-effects model fit by REML 

Formula: UNDS ~ SUCCESSMN + (1 | BIRD) 

          Data: dive 

      AIC     BIC    logLik MLdeviance REMLdeviance

 60032.37 60053.8 -30013.19    60031.9     60026.37

Random effects:

 Groups   Name        Variance Std.Dev.

 BIRD     (Intercept)  4.4504  2.1096  

 Residual             36.4240  6.0352  

number of obs: 9324, groups: BIRD, 12

 

Fixed effects:

            Estimate Std. Error t value

(Intercept) 13.39764    0.63887 20.9707

SUCCESSMN    4.22197    4.11527  1.0259

 

Correlation of Fixed Effects:

          (Intr)

SUCCESSMN -0.276

> anova(reml.res)

Analysis of Variance Table

          Df Sum Sq Mean Sq

SUCCESSMN  1 38.337  38.337


	[[alternative HTML version deleted]]


