Peaks Genen zuordnen
Du bist fast so weit, die Beziehung zwischen Peaks, Genen und ihrer Funktion genauer zu untersuchen. Aber zuerst schauen wir uns an, wie sich Peak-Positionen auf Gene beziehen. Das Paket chipenrich stellt dafür hilfreiche Plotfunktionen bereit. Das Objekt peaks enthält die Koordinaten der Peaks, die du zuvor verwendet hast. Die Funktion plot_dist_to_tss() kann ein Diagramm der Verteilung der Abstände zwischen Peak-Positionen und dem Transkriptionsstartpunkt (TSS) von Genen erstellen.
Mit der Funktion plot_chipenrich_spline() kannst du die Häufigkeit von Peaks in Relation zur Genlänge beurteilen. Zusätzlich zu den beobachteten Häufigkeiten werden auch die erwarteten Verteilungen für den exakten Test nach Fisher, ein Binomialmodell und das von der Funktion chipenrich() verwendete Modell geplottet.
Diese Übung ist Teil des Kurses
<Kurs>ChIP-seq mit Bioconductor in R</Kurs>Interaktive praktische Übung
Versuche dich an dieser Übung, indem du diesen Beispielcode vervollständigst.
# Plot distribution of distances between peaks and transcription start sites
___(___, genome = "hg19")