linest {doBy} | R Documentation |
Compute linear estimates for a range of models. One example of linear estimates is population means (also known as LSMEANS).
linest(object, K=NULL, level=0.95, ...)
object |
Model object |
K |
Either |
level |
The level of the (asymptotic) confidence interval. |
... |
Additional arguments; currently not used. |
A dataframe with results from computing the contrasts.
Søren Højsgaard, sorenh@math.aau.dk
## Make balanced dataset dat.bal <- expand.grid(list(AA=factor(1:2), BB=factor(1:3), CC=factor(1:3))) dat.bal$y <- rnorm(nrow(dat.bal)) ## Make unbalanced dataset # 'BB' is nested within 'CC' so BB=1 is only found when CC=1 # and BB=2,3 are found in each CC=2,3,4 dat.nst <- dat.bal dat.nst$CC <-factor(c(1,1,2,2,2,2,1,1,3,3,3,3,1,1,4,4,4,4)) mod.bal <- lm(y ~ AA + BB*CC, data=dat.bal) mod.nst <- lm(y ~ AA + BB : CC, data=dat.nst) K <- LSmatrix(mod.nst, effect=c("BB","CC")) linest( mod.nst, K )