Lrm: Construct hypothesis matrix to test repeated measures factor...

In gmonette/spida15: Collection of miscellaneous functions for mixed models etc. prepared for SPIDA 2009+ (development version)

Description

Construct hypothesis matrix to test repeated measures factor effects

Usage

Lrm(fit, nam, vals = 1:nrow(L.mu))

Arguments

fit
nam
vals

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (fit, nam, vals = 1:nrow(L.mu)) 
{
    L.mu <- Lmu(fit, nam)
    pp <- cbind(1, Poly(vals, nrow(L.mu) - 1))
    ortho <- Q(pp)[, -1]
    ortho <- ortho[, -1]
    maxp <- max(5, nrow(L.mu))
    colnames(ortho) <- c("linear", "quadratic", "cubic", paste("^", 
        4:maxp, sep = ""))[1:ncol(ortho)]
    L <- t(ortho) %*% L.mu
    L
  }

gmonette/spida15 documentation built on May 17, 2019, 7:26 a.m.