Christoph Buser wrote:
Hi all
I tried to reproduce an example with lme and used the Orthodont
dataset.
library(nlme)
fm2a.1 - lme(distance ~ age + Sex, data = Orthodont, random = ~ 1 | Subject)
anova(fm2a.1)
numDF denDF F-value p-value
(Intercept) 180 4123.156 .0001
age
Douglas Bates wrote:
Christoph Buser wrote:
Hi all
I tried to reproduce an example with lme and used the Orthodont
dataset.
library(nlme)
fm2a.1 - lme(distance ~ age + Sex, data = Orthodont, random = ~ 1 |
Subject)
anova(fm2a.1)
...
Regards,
Christoph Buser
No. The calculation of denominator
Hi, Christoph:
As documented in Pinheiro and Bates (2000) Mixed-Effects Models
for S and S-Plus (Springer), the nlme package includes a function
simulate.lme. They used that function to study the performance of the
likelihood ratio statistic under the null hypothesis of no effect. The