2x2 因子実験での差次的発現の検定
コントラストの数はこれまでの例より多くなりますが、limma で差次的発現を検定するためのパイプライン自体は同じです。
この演習はコースの一部です
R での limma を用いた Differential Expression 解析
演習の手順
Populus データを含む ExpressionSet オブジェクト eset、デザイン行列（design）、コントラスト行列（cm）はすでにワークスペースに読み込まれています。
lmFitでモデル係数を当てはめます。contrasts.fitでコントラストを当てはめます。eBayesで t 統計量を計算します。decideTestsで結果を要約します。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Load package
library(limma)
# Fit the model
fit <- ___(eset, ___)
# Fit the contrasts
fit2 <- ___(fit, contrasts = ___)
# Calculate the t-statistics for the contrasts
fit2 <- ___(fit2)
# Summarize results
results <- ___(fit2)
summary(results)