Extract (or "get") components from a KRmodcomp
object,
which is the result of the KRmodcomp
function.
getKR(
object,
name = c("ndf", "ddf", "Fstat", "p.value", "F.scaling", "FstatU", "p.valueU", "aux")
)
getSAT(object, name = c("ndf", "ddf", "Fstat", "p.value"))
Arguments
- object
A KRmodcomp
object, which is the result of the
KRmodcomp
function
- name
The available slots. If name
is missing or NULL
then everything is returned.
References
Ulrich Halekoh, Søren Højsgaard (2014)., A Kenward-Roger
Approximation and Parametric Bootstrap Methods for Tests in Linear Mixed
Models - The R Package pbkrtest., Journal of Statistical Software,
58(10), 1-30., https://www.jstatsoft.org/v59/i09/
Examples
data(beets, package='pbkrtest')
lg <- lmer(sugpct ~ block + sow + harvest + (1|block:harvest),
data=beets, REML=FALSE)
sm <- update(lg, .~. - harvest)
modcomp <- KRmodcomp(lg, sm)
getKR(modcomp, "ddf") # get denominator degrees of freedom.
#> [1] 2