การเพิ่ม Annotations

ดูกราฟนี้:

กราฟนี้คล้ายกับที่สร้างในแบบฝึกหัดก่อนหน้ามาก แต่มี track เพิ่มเติมชื่อ Genes ซึ่งแสดง gene annotations ควรเรียงลำดับการรันโค้ดส่วนต่าง ๆ ต่อไปนี้อย่างไร เพื่อเพิ่ม track นี้ลงในกราฟ?

ส่วนที่ 1

plotTracks(list(ideogram, cover_track,peak_track, tx, 
           GenomeAxisTrack()), chromosome="chr20", 
           from = start_pos, to=end_pos)

ส่วนที่ 2

tx <- AnnotationTrack(TxDb.Hsapiens.UCSC.hg19.knownGene, name="Genes")

ส่วนที่ 3

tx <- GeneRegionTrack(TxDb.Hsapiens.UCSC.hg19.knownGene, name="Genes")

ส่วนที่ 4

library(TxDb.Hsapiens.UCSC.hg19.knownGene)