Hledání společných témat
Balíček chipenrich nabízí funkci chipenrich(), která identifikuje skupiny genů, jež jsou s ChIP-seq peaky spojeny častěji, než by se dalo očekávat náhodou. Klíčové je přitom rozhodnout, jak geny do skupin sdružit. V tomto cvičení se podíváš na Hallmark genové sady definované na Broad Institute.
Zpravidla bys analýzu omezil/a na diferenciálně vázané peaky, aby lépe vynikly molekulární procesy odlišující obě skupiny vzorků. Kvůli malé velikosti datové sady tu ale jednoduše zahrneš peaky se silnějším signálem ve vzorcích nádorů rezistentních na léčbu.
Toto cvičení je součástí kurzu
ChIP-seq s Bioconductor v R
Pokyny k cvičení
- Vyber všechny peaky s vyšší intenzitou ve vzorcích rezistentních nádorů než ve vzorcích primárních nádorů.
- Spusť analýzu obohacení.
- Vypiš výsledky analýzy.
Interaktivní cvičení na vyzkoušení si v praxi
Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.
# Select all peaks with higher intensity in treatment resistant samples
turp_peaks <- peaks_binding[, "GSM1598218"] + peaks_binding[, "GSM1598219"] < ___[, "GSM1598223"] + ___[, "GSM1598225"]
# Run enrichment analysis
enrich_turp <- ___(peaks_comb[turp_peaks, ], genome="hg19",
genesets = "hallmark", out_name = NULL,
locusdef = "nearest_tss", qc_plots=FALSE)
# Print the results of the analysis
___(enrich_turp$results)