corfdrci.inv {GeneNT}R Documentation

Inverse screening procedure based on Pearson correlation coefficient

Description

The inverse screening procedure is presented in the the paper below. It allows the FDR p-value of a gene pair's strength of association measured using Pearson correlation coefficient to be computed

Usage

corfdrci.inv(cormin)

Arguments

cormin Pearson correlation threshold

Value

This function returns a list of fdr adjusted p-values.

Author(s)

Dongxiao Zhu (http://www-personal.umich.edu/~zhud)

References

Zhu, D., Hero, A.O., Qin, Z.S. and Swaroop, A. High throughput screening of co-expressed gene pairs with controlled False Discovery Rate (FDR) and Minimum Acceptable Strength (MAS). Submitted.

Fisher, R.A. (1921). On the 'probable error' of a coefficient of correlation deduced from a small sample. Metron, 1, 1–32.

See Also

pcorfdrci.inv

Examples


# load GeneNT and GeneTS library
library(GeneTS)
library(GeneNT)

#EITHER use the example dataset
data(dat) 
#OR use the following if you want to import external data 
#dat <- read.table("gal.txt", h = T, row.names = 1) 
#Note, data matrix name has to be "dat"
#Subset data to run faster
dat <- dat[1:10,]
#Set MAS level cormin = 0.8, calculate fdr adjusted p-values. 
#fdrp8 <- corfdrci.inv(0.8)                                   

[Package GeneNT version 1.2 Index]