Error with a function in ATAC-seq analysis
Hi ads,
I got an error when trying to apply code from a paper: https://epigeneticsandchromatin.biomedcentral.com/articles/10.1186/s13072-020-00342-y#MOESM1
library(GenomicRanges)
library(csaw)
diseased1.peaks <- GRanges(diseased1.peaks)
diseased.peaks <- union(diseased1.peaks, diseased2.peaks)
Error in as.vector(x) : no method for coercing this S4 class to a vector
Would you suggest me what wrong here? This error happens in other cases as well, not only with union(). Thank you so much!
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
ATAC seq Normalization and analysis pipeline problems
written by bridgetaylward13 •
-
How can I use a vcf file in trackViewer? Issues with TabixFile - coercing S4 class to a vector
written by cristinadias.gen •
-
Can't merge gene names with accession numbers
written by mintdrink •
-
R+Gviz: Fake AlignmentsTrack() to visualize Enhancer-Promoter connections
written by Matthias Zepper
-
Error in as.vector(x) : no method for coercing this S4 class to a vector
written by pr2009 •
-
how to query biomaRt from a GRanges object
written by sebastianzeki0
-
Filtering Bam file based on specific regions using Rsamtools
written by laraib85 •
-
Problem With Iranges And Findoverlap
written by e.karasmani •
sounds like
diseased1.peaksis not a vector. What do you get when you runstr(diseased1.peaks)?Yes, it is not a vector. Formal class 'GRanges' Because GRanges() function made it become a Granges class which is a data frame before that.
run
unionon the data framesSo you skip the GRanges()? Yes, no error anymore but I am not sure if that causes changing in downstream analysis.