การเชื่อมโยง Peaks กับยีน
ตอนนี้เกือบพร้อมแล้วที่จะสำรวจความสัมพันธ์ระหว่าง peaks, ยีน และหน้าที่ของยีนอย่างละเอียด แต่ก่อนอื่น มาดูให้ชัดขึ้นก่อนว่าตำแหน่งของ peaks สัมพันธ์กับยีนอย่างไร แพ็กเกจ chipenrich มีฟังก์ชันสำหรับสร้างกราฟที่เป็นประโยชน์สำหรับจุดประสงค์นี้ ออบเจกต์ peaks เก็บพิกัดของ peaks ที่ใช้ไปก่อนหน้านี้ ฟังก์ชัน plot_dist_to_tss() สร้างกราฟแสดงการกระจายของระยะทางระหว่างตำแหน่ง peaks กับ transcription start site (TSS) ของยีน
ฟังก์ชัน plot_chipenrich_spline() ช่วยให้ประเมินความถี่ของ peaks เทียบกับความยาวของยีนได้ นอกจากความถี่ที่สังเกตได้แล้ว ยังแสดงการกระจายที่คาดหวังสำหรับ Fisher's exact test, โมเดล binomial และโมเดลที่ใช้โดยฟังก์ชัน chipenrich() อีกด้วย
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
ChIP-seq with Bioconductor in R
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Plot distribution of distances between peaks and transcription start sites
___(___, genome = "hg19")