findThreshold: Find the threshold value to create a co-expression network

View source: R/findThreshold.R

findThresholdR Documentation

Find the threshold value to create a co-expression network


Finds the threshold value to establish the cutoff in the process to define the edges in the co-expression network final from two steps. In the first one, obtains the subtraction from clustering coefficient values of the real and random networks created from the possible threshold values in the correlation matrix. In the second one, a Kolmogorov-Smirnov test is made to evaluate the degree distribution respect normality.


findThreshold(expData, method, plotting = FALSE)



A whole expression matrix or the expression matrix to differentially expressed genes, it may be stored in a SummarizedExperiment object.


The method name to create the correlation matrix, this can be "correlation" to obtain the Pearson Correlation Coefficient. On the other hand, can be "mutual information" to obtain the correlation values from an entropy-based method.


The option to show the result in a plot. By default FALSE.


The best threshold value found using the two criteria and a plot showing the result.


Juan David Henao Sanchez <>

Liliana Lopez Kleine <>


Elo, L. L., Jarvenpaa, H., Oresic, M., Lahesmaa, R., & Aittokallio, T. (2007). Systematic construction of gene coexpression networks with applications to human T helper cell differentiation process. Bioinformatics, 23(16), 2096-2103.

Leal, L. G., Lopez, C., & Lopez-Kleine, L. (2014). Construction and comparison of gene co-expression networks shows complex plant immune responses. PeerJ, 2, e610.

See Also

difExprs to find the differentially expressed genes matrix.


# Loading data

pathfile <- system.file("extdata","expression_example.txt",package = "coexnet")
data <- read.table(pathfile,stringsAsFactors = FALSE)

# Finding threshold value

cor_pearson <- findThreshold(expData = data,method = "correlation")

jdhenaos/coexnet documentation built on April 22, 2022, 12:43 a.m.