Lequal: Hypothesis matrix for equality of factor level effects

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

Usage

Lequal(fit, pat)

Arguments

fit
pat

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, pat) 
{
    L <- Lmat(fit, pat)
    nam <- rownames(L)
    n <- nrow(L)
    if (n < 2) 
        return(L)
    plus <- unlist(apply(rbind(2:n), 2, seq, n))
    minus <- rep(1:(n - 1), (n - 1):1)
    Lp <- L[plus, ]
    Lm <- L[minus, ]
    Lret <- rbind(Lp - Lm)
    rn <- paste(nam[plus], nam[minus], sep = " - ")
    rownames(Lret) <- rn
    Lret
  }

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