MWUTest: MWU test for RNA-Seq Data

View source: R/MWUTest.R

MWUTestR Documentation

MWU test for RNA-Seq Data

Description

Computes p-values using a Mann-Whitney U Test for the null hypothesis H_0: β_j = 0.

Usage

MWUTest(dge, design)

Arguments

dge

A dgeList object, created with edgeR, containing the normalization factors as computed by edgeR.

design

A model matrix; The first column should be all 1s. The second column should have two unique values, corresponding to the groups

Author(s)

Jakob Walter

Examples

Y <- rnbinom(20*10, mu = 10, size = 1/0.2)
Y <- data.frame(array(Y, dim = c(20, 10)))
X1 <- as.factor(rep(c("A", "B"), each = 20/2))
design <- model.matrix(~X1, contrasts.arg = list(X1 = "contr.sum"))
dge <- edgeR::DGEList(counts = t(Y), group = X1)
dge <- edgeR::calcNormFactors(dge)
MWUTest(dge, design)

jakobwalter/RNAPerm documentation built on Sept. 17, 2022, 3:57 p.m.