doBy (version 4.5-10)

linest: Compute linear estimates

Description

Compute linear estimates for a range of models. One example of linear estimates is population means (also known as LSMEANS).

Usage

linest(object, K, level=0.95, ...)

Arguments

object

Model object

K

A contrast matrix. If given then the specifications in effect and at are ignored. Otherwise, the matrix is generated according to the specifications in effect and at.

level

The level of the (asymptotic) confidence interval.

...

Additional arguments; currently not used.

Value

• A dataframe with results from computing the contrasts.

See Also

LSmeans LSmatrix

Examples

## 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 )