Puncak mana yang berbeda?
Mari kita lihat seberapa berbeda sampel-sampel tersebut dalam hal pemanggilan puncak (peak calls). Dengan melihat puncak gabungan yang sudah Anda gunakan sebelumnya, berapa banyak puncak yang dibagi oleh berbagai sampel?
Salah satu plot yang dapat menjawab pertanyaan ini adalah diagram Venn. Diagram ini dapat dibuat langsung dari objek DBA menggunakan fungsi dba.plotVenn(). Alternatifnya adalah menggunakan plot UpSet, yang dapat Anda buat dengan paket UpSetR. Plot ini menggantikan elips pada diagram Venn dengan serangkaian diagram batang untuk menunjukkan ukuran berbagai himpunan bagian. Fungsi upset() memerlukan Anda menyediakan sebuah data frame dengan satu kolom untuk setiap himpunan, yaitu setiap sampel, yang menunjukkan keanggotaan himpunan untuk setiap puncak. Data yang diperlukan terdapat pada kolom called dari objek ar_diff.
Latihan ini adalah bagian dari kursus
ChIP-seq dengan Bioconductor di R
Latihan interaktif praktis
Cobalah latihan ini dengan menyelesaikan kode contoh berikut.
# Create Venn diagram
___(ar_diff, mask=1:4)