mixedroastlmerlist: Multiple Contrast Rotation Tests based on Marginal Mixed...

Description Usage Arguments Details Value Author(s) References

Description

Adaption of the roast function in package limma to perform a maxT-test for fixed effects in multiple linear mixed models, controlling FWER for multiple contrasts.

Usage

1
mixedroastlmerlist(lmerlist, K, nrot = 9999, adjusted = TRUE)

Arguments

lmerlist

A list with multiple lmer modl objects of package lme4.

K

A contrast matrix (different contrasts in rows), with the same number of columns as number of parameters in the model

nrot

Number of Monte Carlo simulation steps

adjusted

logical, if TRUE single-step maxT FWER control is performed; if FALSE the comparison-wise error rate is controlled

Details

Most of the code was taken from function roast in package limma (Bioconductor).

Value

A list with

Y

The matrix of response variables

X

The design matrix

statistic

A matrix of test statistics

pvalue

A matrix of (adjusted) p-values

B

A matrix with contrast parameter estimates

mrot

A vector with simulated maxT statistics under H0

Author(s)

Daniel Gerhard

References

Wu, D, Lim, E, Francois Vaillant, F, Asselin-Labat, M-L, Visvader, JE, and Smyth, GK (2010). ROAST: rotation gene set tests for complex microarray experiments. Bioinformatics 26, 2176-2182.

Langsrud, O (2005). Rotation tests. Statistics and Computing 15, 53-60.


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