MBttest (version 1.0.0)

oddratio: Calculation of Zeta

Description

Zeta is used to measure homogeneity intensity of two subdatasets. If zeta >1, these two subdatasets have good homogeneity; otherwise, zeta

Usage

oddratio(XX, na, nb)

Arguments

XX

count data of RNA reads generated by next generation sequencing.

na

number of replicate libraries in condition A.

nb

number of replicate libraries in condition B.

Value

• oddrat
  list of zeta values

Details

Zeta is defined as $$zeta=log(1+(mean(XA,XB)*var(XA,XB))/(mean(XA)*var(XA)+mean(XB)*var(XB)))$$ zeta is different from psi. If two subdatasets have big a gap and good homogeneity,then seta value has much larger than 1.

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

pratio, mbetattest.

Examples

XX<-matrix(NA,2,8)
XX[1,]<-c(112,122, 108,127,302, 314, 322, 328)
XX[2,]<-c(511, 230, 754, 335,771, 842, 1014,798)
#XX
#     [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8]
#[1,]  112  122  108  127  302  314  322  328
#[2,]  511  230  754  335  771  842 1014  798
oddratio(XX=XX,na=4,nb=4)

#[1] 3.9432676 0.8762017	
	
# see example in mbetattest