clusterGenes: Cluster Genes

Description Usage Arguments Value Note Author(s) References

Description

Compute topological overlap matrix and gene dendrogram given a data frame of normalized gene expression data.

Usage

1
clusterGenes(expressionData, beta = NULL, RsquaredCut = 0.9)

Arguments

expressionData

A normalized matrix of gene expression data, given as a data frame with rows as samples and columns as genes/probes.

beta

Soft thresholding exponent, typically 5-6. Setting to NULL (the default) causes the function to select the optimal beta between 1 and 12.

RsquaredCut

the threshold for linearity of the 'scale-free' (log(p(k)) vs. log(k)) plot after for accepting a candidate value of beta. The function returns the minimum beta meeting this criterion. Note: If no value of beta causes the R-squared value of the regression to exceed the given threshold, then the threshold is decremented by 0.05 until a solution is found.

Value

beta

The 'soft threshold' beta, or NULL (the default) if the function is to determine what value to use.

sftStatistics

A data frame having the statistics collected while determining beta, as produced by WGCNA::pickSoftThreshold, or NULL if beta is passed in.

dichotCor

|cc|^beta (elementwise, where 'cc' is the matrix of inter-gene Pearson correlation coefficients) but with zero diagonal.

tomDist

The toplogical overlap matrix, in the form of a distance metric.

geneTree

The gene dendrogram created from WGCNA::flashClust using 'tomDist' as the distance metric.

Note

This is both the most computationally intensive and the most rote part of coexpression analysis. Therefore it is broken out as an independent function. Variations on the downstream steps can then be run without repeating this step.

Author(s)

Bruce Hoff

References

http://www.genetics.ucla.edu/labs/horvath/CoexpressionNetwork/

Zhang, B. & Horvath, S. (2005) Statistical Applications in Genetics and Molecular Biology 4, Article 17.

http://www.genetics.ucla.edu/labs/horvath/CoexpressionNetwork/Rpackages/WGCNA/


Sage-Bionetworks/SageBionetworksCoex documentation built on May 9, 2019, 12:11 p.m.