WGCNA解析による遺伝子共発ネットワークの構築
データの読み込みとクラスタリング
library(WGCNA)
# ディレクトリ名の確認
dir()
# データの読み込み
load('DEG_TB_LTBI_step13.Rdata')
# 行がサンプル名、列が遺伝子名
################################ サンプルクラスタリング ####################
datExpr = t(dataset_TB_LTBI_DEG)
# 初期クラスタリング
sampleTree = hclust(dist(datExpr), method = "av ...
6月4日 22:01 投稿