mumotest: Simultaneous Inference for Parameters in Multiple Models

Description Usage Arguments Value Author(s) Examples

Description

Resampling-based testing of general linear hypotheses for parameters of multiple models

Usage

1
mumotest(mlist, K, B, margin=0)

Arguments

mlist

List of model objects

K

A contrast matrix (Number of columns = Number of coefficients in all models)

B

Number of iterations

margin

Test margin

Value

An object of class mumo

Author(s)

Daniel Gerhard

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
### data example from Hasler & Hothorn (2011), 
### A Dunnett-Type Procedure for Multiple Endpoints,
### The International Journal of Biostatistics: Vol. 7: Iss. 1, Article 3.
### DOI: 10.2202/1557-4679.1258
# but using two-sided inference

### see ?coagulation
data("coagulation", package = "SimComp")

### Marginal Models for each endpoint
(m1 <- lm(Thromb.count ~ Group, data = coagulation))
(m2 <- lm(ADP ~ Group, data = coagulation))
(m3 <- lm(TRAP ~ Group, data = coagulation))

### Dunnett contrast for comparisons to a control for each endpoint
K <- rbind("Thromb: B - S"=c(0,1,0,0,0,0,0,0,0),
           "Thromb: H - S"=c(0,0,1,0,0,0,0,0,0),
           "ADP: B - S"   =c(0,0,0,0,1,0,0,0,0),
           "ADP: H - S"   =c(0,0,0,0,0,1,0,0,0),
           "TRAP: B - S"  =c(0,0,0,0,0,0,0,1,0),
           "TRAP: H - S"  =c(0,0,0,0,0,0,0,0,1))

### Resampling of contrast test statistics
mm <- mumotest(list(m1, m2, m3), K, B=10000)

### Adjusted p-values
summary(mm)

### Simultaneous confidence intervals
confint(mm)

daniel-gerhard/mumocomp documentation built on May 14, 2019, 3:38 p.m.