[prev in list] [next in list] [prev in thread] [next in thread] 

List:       r-help
Subject:    Re: [R] variance estimates in lme biased?
From:       "Pascal A. Niklaus" <Pascal.Niklaus () unibas ! ch>
Date:       2003-12-17 10:18:20
[Download RAW message or body]

Running lme on your data set results exactly in what you expect - or do 
you expect something different?

Pascal

 > L1<-factor(F1f)
 > L2<-factor(F2f)
 > L3<-factor(F3f)
 > lme(value ~ 1,random = ~ 1 | L1/L2/L3)
Linear mixed-effects model fit by REML
  Data: NULL
  Log-restricted-likelihood: 438.9476
  Fixed: value ~ 1
(Intercept)
  0.2955631

Random effects:
 Formula: ~1 | L1
        (Intercept)
StdDev:  0.02472988              <== level F1 which is 0

 Formula: ~1 | L2 %in% L1
        (Intercept)
StdDev:    1.140782                <== level F2 which is 1

 Formula: ~1 | L3 %in% L2 %in% L1
         (Intercept)  Residual
StdDev: 0.0005512791 0.1020479   <== F3 which is 0 , and F4 which is 0.1

Number of Observations: 625
Number of Groups:
                L1         L2 %in% L1 L3 %in% L2 %in% L1
                 5                 25                125

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help

[prev in list] [next in list] [prev in thread] [next in thread]