Ldiff.old: Older version of Ldiff

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

Usage

Ldiff.old(fit, pat, levnames = c(reflevel, substring(rownames(L), cut + 1)), reflevel = "<ref>", cut = nchar(pat))

Arguments

fit
pat
levnames
reflevel
cut

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, levnames = c(reflevel, substring(rownames(L), 
    cut + 1)), reflevel = "<ref>", cut = nchar(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(L, Lp - Lm)
    rn <- paste(levnames[c(1:n, plus) + 1], levnames[c(rep(0, 
        n), minus) + 1], sep = " - ")
    rownames(Lret) <- rn
    Lret
  }

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