How to calculate eigengenes on arbitrary expression data, not inside WGCNA
1
1
Entering edit mode
4.4 years ago
tospo ▴ 50

I would like to compare the results from different algorithms for clustering data from an RNAseq experiment. One methodology we used is WGCNA, which represents each cluster/module by a module eigengene.

What I am now looking for is a way to calculate such eigengenes on arbitrary lists of genes and their expression profiles. I want to use that to reduce the clusters generated by other methods to a set of "representative genes" and compare that to the WGCNA output.

Unfortunately, I'm not clear about how to get from a list of genes and their expression profiles to an eigengene. The functions inside the WGCNA R package are deeply tied into the WGCNA analysis and I can't see how to use them on arbitrary data frames of gene expression data. Any hints would be highly welcome!

WGCNA eigengene rnaseq rna-seq • 5.6k views
9
Entering edit mode
4.4 years ago

Eigengenes were first defined in this paper. Singular value decomposition (SVD) is what you're looking for. The eigengenes are the right singular vectors of the SVD of the expression matrix. If X is your data with genes as rows and samples as columns, the SVD of X is X=USV' and the eigengenes are defined as the vectors in V. In R,

eigengenes <- svd(X)\$v


0
Entering edit mode

thank you, I will try that!

0
Entering edit mode

BTW: the link to "this paper" doesn't work

1
Entering edit mode

0
Entering edit mode

Thank you very much for this!!! I was looking for a similar solution for months and somehow just found this one line of code that solved my problem!