chromVAR icon indicating copy to clipboard operation
chromVAR copied to clipboard

no output

Open rojinsafavi opened this issue 4 years ago • 0 comments

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)

rojinsafavi avatar Apr 10 '20 16:04 rojinsafavi