Ldiff.rdc: Version of Ldiff used in RDC

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

Usage

Ldiff.rdc(fit, nam, ref = "no longer used")

Arguments

fit
nam
ref

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, ref = "no longer used") 
{
    Lm <- Lmu(fit, nam)
    levs <- rownames(Lm)
    n <- nrow(Lm)
    if (n < 2) 
        return(Lm)
    plus <- unlist(apply(rbind(2:n), 2, seq, n))
    minus <- rep(1:(n - 1), (n - 1):1)
    Lret <- Lm[plus, ] - Lm[minus, ]
    rn <- paste(levs[plus], levs[minus], sep = " - ")
    rownames(Lret) <- rn
    Lret
  }

gmonette/spida documentation built on May 17, 2019, 7:25 a.m.