empdist: compute the empirical cumulative density.

Description Usage Arguments Details Value References Examples

Description

compute the empirical cumulative density given two random variables, which will be the input for the mixture copula model

Usage

1
empdist(x1, x2)

Arguments

x1

a vector of data for replicate 1, i.g., sequencing counts for each genomic locus.

x2

a vector of data for replicate 2. It must have the same length with x1.

Details

The data was first ranked with a ties method equals to "random". The cumulative density function (cdf) is then obtained using the ecdf function. To avoid infinity, a factor is multiplied to the cdf. The factor is the length(x1)/[length(x1)+1].

Value

The empdist function returns a matrix that has two columns, each stores the empirical cumulative density of one replicate.

References

Q. Li, J. B. Brown, H. Huang and P. J. Bickel. (2011) Measuring reproducibility of high-throughput experiments. Annals of Applied Statistics, Vol. 5, No. 3, 1752-1779.

Examples

1
2
3
4
data(Chipseq_TF)
x1 <- Chipseq_TF[,1]
x2 <- Chipseq_TF[,2]
empdist(x1, x2)

MonkeyLB/TDR documentation built on May 7, 2019, 4:59 p.m.