row.kruskal.wallis: Apply the Kruskal-Wallis test many times

Description Usage Arguments Details Value Author(s) References Examples

View source: R/row.kruskal.wallis.R

Description

For each row of Y, use the Kruskal-Wallis test to compare medians across groups defined in grplbl

Usage

1
row.kruskal.wallis(Y, grplbl)

Arguments

Y

data matrix with variables in rows and observations (subjects) in columns

grplbl

vector of group labels

Details

The alternative hypothesis is that, for each gene, there are at least two groups of different median. The null hypothesis is that all groups have the same median for each gene studied.

Value

A data.frame with three columns

stat

a vector with the test-statistic for each row of Y

pval

a vector with the p-value for each row of Y

ebp

a vector with the empirical Bayes probability of equal medians for each row of Y

Author(s)

Stan Pounds <stanley.pounds@stjude.org>; Demba Fofana <demba.fofana@stjude.org>

References

Kruskal, W.H. and W.A. Wallis (1952) Use of ranks in one-criterion variance analysis. J. Amer. Stat. Assoc. 47:583-621.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
  
####################Three group comparison###################
# load data
data(GroupComp.data)
# Read the expression values   
brain.express.set <- exprs(GroupComp.data)
head(brain.express.set)
# Read the phenotype
brain.pheno.data <- pData(GroupComp.data)
brain.pheno.data[,1] 
row.kruskal.wallis(brain.express.set,brain.pheno.data[,1])
 

HybridMTest documentation built on Nov. 8, 2020, 8:29 p.m.