p 值直方圖
一個實用的診斷圖是 p 值的直方圖。若低 p 值的密度很高,代表有許多差異表現基因;相反地,若直方圖接近均勻分佈,則表示很少。請為白血病研究建立一個 p 值直方圖。
本練習屬於課程
使用 R 中的 limma 進行差異表達分析
練習說明
第 2 章白血病研究的擬合模型物件 fit2 已載入到你的工作空間。limma 套件也已載入。
使用
topTable取得每個基因的摘要統計。若要取得所有基因的結果,將引數
number設為fit2的列數。若要停用依顯著性排序結果,將引數
sort.by設為"none"。使用
hist繪製 p 值的直方圖。
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
# Obtain the summary statistics for every gene
stats <- ___(fit2, number = ___, sort.by = ___)
# Plot a histogram of the p-values
___(stats[, ___])