View source: R/findThreshold.R
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 <firstname.lastname@example.org>
Liliana Lopez Kleine <email@example.com>
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.
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") cor_pearson
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.