Ldiff: Hypothesis matrix for ...

Usage Arguments Examples

View source: R/wald.R

Usage

1
Ldiff(fit, pat, levnames = c(reflevel, substring(rownames(L), cut + 1)), reflevel = "<ref>", 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 = "<ref>", 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.