[R] A question on REML in R

From: Shuangshuang Jin ^_^ <charmyss_at_hotmail.com>
Date: Thu 04 Jan 2007 - 22:06:40 GMT


Hello, everyone, I'm using R to deal with a REML problem. I found "lmer" is the right function for this. But I got stuck because I couldn't interpret the result. I'm attaching a short example of my executing log. Please have a look and give me some advice on it. Thanks a lot!

  Plot Block Treatment Data

    1     1        2         7.8
    2     1        1         5.9
    3     1        3        10.3
    4     2        3        10.9
    5     2        2         8.9
    6     2        1         7.2
    7     3        2        11.1
    8     3        3        12.8
    9     3        1         9.1
   10     4        1         9.8
   11     4        2        12.2
   12     4        3        14.0

>anova(lm(Data~as.factor(Treatment)+as.factor(Block)))
Analysis of Variance Table

Response: Data

                     Df Sum Sq Mean Sq F value    Pr(>F)
as.factor(Treatment)  2 32.000  16.000  282.35 1.162e-06 ***
as.factor(Block)      3 30.000  10.000  176.47 3.066e-06 ***
Residuals             6  0.340   0.057
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1


>lmer(Data~Treatment+(Treatment|Block))
Linear mixed-effects model fit by REML Formula: Data ~ Treatment + (Treatment | Block) AIC BIC logLik MLdeviance REMLdeviance 28.68 31.1 -9.339 16.88 18.68 Random effects: Groups Name Variance Std.Dev. Corr Block (Intercept) 3.3116e+00 1.8198e+00 Treatment 2.4303e-11 4.9298e-06 0.000 Residual 4.8606e-02 2.2047e-01 number of obs: 12, groups: Block, 4 Fixed effects: Estimate Std. Error t value (Intercept) 6.00000 0.92533 6.484 Treatment 2.00000 0.07795 25.658 Correlation of Fixed Effects: (Intr) Treatment -0.168 Warning message: Estimated variance-covariance for factor 'Block' is singular in: `LMEoptimize<-`(`*tmp*`, value = list(maxIter = 200, tolerance = 1.49011611938477e-08, Charmy _________________________________________________________________ Communicate instantly! Use your Hotmail address to sign into Windows Live Messenger now. http://get.live.com/messenger/overview ______________________________________________ R-help@stat.math.ethz.ch 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 Fri Jan 05 10:35:48 2007

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Fri 05 Jan 2007 - 00:30:25 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.