cal_cdf: This function calculates the weighted-pvalues.

In mqzhanglab/wHC: Goodness of fit Statistics Calculation based on weighted p-values

This function calculates the weighted-pvalues.

Description

It calculates the CDF(Cumulative Distribution Function) of weighted pvalues from the original pvalues and their weights and assign new pvalues based on this CDF.

Usage

cal_cdf(pm, w = 1)

Arguments

pm	is a statistics matrix of P-values from nrow=n genes(independent tests),ncol=d. Pm are not encouraged to have only 1 rows, if that happend, warning massage will produced and this function returns sqrt(1/pm-1), rather than the double side statistics from linear regression.Thus, we have S(t)~Binomial(n,p).
w	is the weight, if not specify, w=1, if specify w must have the same length as nrow Pm

Value

numeric vector with each elements is a Higher criticism values calculated from each colum of the Pm

References

Genovese, C. R., Roeder, K., & Wasserman, L. (2006). False Discovery Control with p-Value Weighting. Biometrika, 93(3), 509–524.

Examples

pval=matrix(runif(200,0,1),ncol=4,nrow=50)
w0=seq(0.5,1.5,length=50)
pwval=cal_cdf(pval,w=w0)

mqzhanglab/wHC documentation built on April 1, 2022, 6:28 p.m.