momega {tatest}

Null rho for testing multiple null hypotheses

Description

Momega is a function that is used to estimate omega for testing multiple null hypotheses. Omege is a null \rho that is used as a threshold for \rho from real data. Since this omega is majorly applied to micoarray, we named this function as momega.

Usage

momega(dat, na, nb, distr = "norm")

Arguments

dat

a matrix dataset containing string columns and two-condition numeric columns

na

numeric value, size of sample 1 or number of replicates in condition A.

nb

numeric value, size of sample 1 or number of replicates in condition A.

distr

string, data distribution.

Details

This function is to use null data to calculate omega value with rho = 1.

Value

return a numeric value

Author(s)

Yuan-De Tan [e-mail]

References

Yuan-De Tan Anita M. Chandler, Arindam Chaudhury, and Joel R. Neilson (2015) A Powerful Statistical Approach for Large-scale Differential Transcription Analysis. Plos One. 2015 DOI: 10.1371/journal.pone.0123658.

See Also

omega,rhov

Examples

data(Ecadherin)
Ecad<-as.data.frame(Ecadherin)
Ecad1<-Ecad[,-1]
w<-momega(dat=Ecad1[1:2000,], na=3, nb=3, distr = "norm")