Read counts laden
Zoals je in de video hebt gezien, moet je een set consensus-peak calls maken voordat je differentiële binding kunt testen. Dit kan met de volgende regel R-code:
ar_counts <- dba.count(ar_peaks, summits=200)
Bekijk de volgende uitspraken.
- In
ar_countshebben alle samples read counts voor dezelfde set peak calls. - Sommige read counts kunnen 0 zijn.
- Alle pieken in
ar_countszijn 200 bp breed. - Alle pieken in
ar_countszijn 400 bp breed.
Welke van deze uitspraken zijn waar?
Deze oefening maakt deel uit van de cursus
ChIP-seq met Bioconductor in R
Interactieve oefening met praktijkervaring
Zet theorie om in actie met een van onze interactieve oefeningen
Begin oefening