開始使用免費開始

探索 DESeq2 結果

NOTE: 載入此練習可能需要多一點時間。

為了降低回傳的差異表現(DE)基因數量,並減少這些 DE 基因在生物學上不具意義的可能性,我們將使用較小的 log2 倍數變化閾值來判定 DE 基因。

本練習屬於課程

使用 R 與 Bioconductor 進行 RNA-Seq

檢視課程

練習說明

  • 使用 results() 函式擷取 smoc2 的結果,作法與先前相同,設定 alpha 為 0.05,且將 condition 的基準層級設為 normal。但這次請使用 0.32 的 log2 倍數變化閾值。假設先前步驟都已完成,包括建立 DESeq2 物件 dds_smoc2,以及執行 DESeq() 函式。

  • 使用 lfcShrink() 函式對 log2 倍數變化進行縮減。

動手互動練習

試著完成這個範例程式碼,體驗一下這個練習。

# Explore the results() function
?results

# Extract results
smoc2_res <- ___(___, 
                contrast = ___, 
                alpha = ___, 
                lfcThreshold = ___)

# Shrink the log2 fold changes
smoc2_res <- ___(___, 
                    ___, 
                    res = ___)
編輯並執行程式碼