row.oneway.anova: Perform one-way ANOVA for many variables.

Description Usage Arguments Details Value Author(s) References Examples

View source: R/row.oneway.anova.R

Description

For each row of Y, use one-way ANOVA to compare means across groups defined by grplbl.

Usage

1
row.oneway.anova(Y, grplbl)

Arguments

Y

data matrix with variables in rows and subjects in columns

grplbl

vector of group labels for the subjects

Details

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

Value

A data.frame with three columns:

stat

a vector with the ANOVA F-statistic for each row of Y

pval

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

ebp

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

Author(s)

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

References

Chambers, J. M. and Hastie, T. J. (1992) Statistical Models in S, Wadsworth & Brooks/Cole.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
  
####################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] 
# ANOVA test
row.oneway.anova(brain.express.set,brain.pheno.data[,1])
 

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