chromVAR
chromVAR copied to clipboard
no output
I am not getting any output for counts_filtered <- filterSamples(example_counts, min_depth = 1500,min_in_peaks = 0.15)
any idea why this is happening?
library(chromVAR)
library(motifmatchr)
library(BSgenome.Hsapiens.UCSC.hg19)
### Using example counts from package ------------------------------------------
data(example_counts, package = "chromVAR")
example_counts <- addGCBias(example_counts,
genome = BSgenome.Hsapiens.UCSC.hg19)
counts_filtered <- filterSamples(example_counts, min_depth = 1500,
min_in_peaks = 0.15)
counts_filtered <- filterPeaks(counts_filtered)
motifs <- getJasparMotifs()
motif_ix <- matchMotifs(motifs, counts_filtered,
genome = BSgenome.Hsapiens.UCSC.hg19)
# computing deviations
dev <- computeDeviations(object = counts_filtered,
annotations = motif_ix)