search for: llgp1

...t;- data.frame(y=c(rnorm(5, 1, 2), rnorm(5,3,5)), gp=factor(rep(1:2, each=5))) fit22 <- lme(y~gp, lmePblm, random=~1|gp, method='ML') fitGp1 <- lm(y~1, lmePblm[lmePblm$gp==1, ]) fitGp2 <- lm(y~1, lmePblm[lmePblm$gp==2, ]) (ll22 <- logLik(fit22)) (llGp1 <- logLik(fitGp1)) (llGp2 <- logLik(fitGp2)) # Why isn't (ll22 = llGp1+llGp2)? (ll22 - llGp1-llGp2) # And secondarily, how can I get the residual standard deviations # within each gp from fit22? Thanks, Spencer Graves