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

## Usage

 1 Ldiff(fit, pat, levnames = c(reflevel, substring(rownames(L), cut + 1)), reflevel = "", cut = nchar(pat), verbose = F)

## Arguments

 fit pat levnames reflevel cut verbose

## Examples

 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 ##---- 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 = "", cut = nchar(pat), verbose = F) { L <- Lmat(fit, pat) nam <- rownames(L) n <- nrow(L) zm <- matrix(1:n, nrow = n, ncol = n) plus <- zm[col(zm) < row(zm)] minus <- rep(1:(n - 1), (n - 1):1) Lp <- L[plus, ] Lm <- L[minus, ] Lret <- rbind(L, Lp - Lm) pnames <- levnames[c(1:n, plus) + 1] mnames <- levnames[c(rep(0, n), minus) + 1] if (verbose) { print(levnames) print(plus) print(minus) print(Lp) print(Lm) print(L) print(Lret) print(pnames) print(mnames) } 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 14, 2017, 1:02 p.m.