还是老习惯,给出官网教程,至于你是看还是不看,它就在那里,等着你的深入研究~
https://horvath.genetics.ucla.edu/html/CoexpressionNetwork/Rpackages/WGCNA/Tutorials/
往期教程:
WGCNA分析图文详解专题中要解释的第三张图,所有结果中最主要的结果图恐怕就是这张了:
官方注释:
Figure 1: Module-trait associations. Each row corresponds to a module eigengene, column to a trait. Each cell contains the corresponding correlation and p-value. The table is color-coded by correlation according to the color legend.
4,每个性状与模块之间的相关性计算是否独立的:即表型放在一起分析和分开单独分析是否有不同?
小面我们来一一解读。如有理解错误,还请各位大侠批评指正。
这里小编自己的理解就是这个模块有134个样本,n个基因,然后用pca主成分分析对这n个基因进行降维取其第一主成分作为这个模块的特征。
这里我们得到横坐标的表型性状向量和纵坐标的每个模块的特征代表eigengene值向量,然后就可以计算每个模块和性状之间的相关性了。
嗯,用的cor函数。默认的method=”pearson”
那么这里图中的小格子中的数值就代表每个性状和每个模块的特征值之间两两计算的相关性值以及对应的pvalue。
重点,敲黑板!!!
这个图说难也不难,但是,却是整个WGCNA分析的最主要的结果了。有个地方很关键,就是将临床表型转换为数值的trait的时候,直接涉及到了后面相关性的计算。如果转换不合理很可能得到不好的结果哟。
今天就说到这里,欢迎大家留言讨论。我们下期再见~
参考资料:
1,https://horvath.genetics.ucla.edu/html/CoexpressionNetwork/Rpackages/WGCNA/Tutorials/