Model DESeq2 - mengeksplorasi dispersi
CATATAN: Memuat latihan ini mungkin memerlukan sedikit waktu lebih lama.
Setelah menyesuaikan model pada latihan sebelumnya, mari kita telaah kecocokan data smoc2 kita terhadap model binomial negatif dengan memplot estimasi dispersi menggunakan fungsi plotDispEsts(). Ingat bahwa estimasi dispersi digunakan untuk memodelkan hitungan mentah; jika dispersi tidak mengikuti asumsi yang dibuat oleh DESeq2, maka ragam dalam data bisa diperkirakan kurang tepat dan hasil DE bisa menjadi kurang akurat.
Asumsi yang dibuat DESeq2 adalah bahwa dispersi umumnya menurun seiring meningkatnya rerata dan kurang lebih mengikuti garis hasil pemfitingan.
Latihan ini adalah bagian dari kursus
RNA-Seq dengan Bioconductor di R
Petunjuk latihan
- Plot estimasi dispersi untuk data
smoc2menggunakan fungsiplotDispEsts(). Asumsikan semua langkah sebelumnya telah dijalankan, termasuk pembuatan objek DESeq2,dds_smoc2, dan menjalankan fungsiDESeq().
Latihan interaktif praktis
Cobalah latihan ini dengan menyelesaikan kode contoh berikut.
# Plot dispersions
___(___)