針對 3 個群組進行差異表現檢定
現在你已經指定了設計矩陣與對比矩陣,可以開始進行差異表現檢定。
本練習屬於課程
使用 R 中的 limma 進行差異表達分析
練習說明
含有缺氧資料的 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)