SimpleRleList to GRanges

I am using 'rtracklayer' package in R to import a ".wig" file. Its a Chip-Seq data.

After importing I get the SimpleRleList which I don't know how to use in "ChIPseekr". I am hoping converting to GRanges might help. Any help how to do it will be great.

here are my codes

temp <- import("GSM1382416_H2AX_MCF7_wz_E2_treat_afterfiting_all.wig.gz", trackline =FALSE)> newpeak = coverage(temp)

class(newpeaks) 
[1] "SimpleRleList" attr(,"package") 
[1] "IRanges"

chip-seq rtracklayer

• 1,855 views • link

written 7.8 years ago by prabin.dm

There was a simple solution. I got it. Thanks

• 1 views • link

written 7.8 years ago by prabin.dm

post the solution if you have found it.

• 1 views • link

written 7.8 years ago by Mark

0 answers

No answers yet.

More posts like this

GRange: Average Metadata Value from 3 Replicates
written 3.3 years ago by RelativeOptics •

I have CPM-normalized genome coverage data from three replicate experiments in the form of bigwig files. I read these into R using rtracklayer and now …
How to resolve the error which I get when I try to make objects for GRanges.
written 4.4 years ago by s.singh •

Hi, I am trying to use ChIPseeker to annotate the ChIP peaks which I called using MACS2. I don't know why I get this error: …
conumee - adding bed file instead of sample detailed regions
written 4.6 years ago by esimonova.me •

I am analyzing the Methylation Array (EPIC) with conumee package. I am trying to import a list of regions in which I am interested to …
Error when trying to import GTF files using rtracklayer's import function
written 5.5 years ago by Marcel •

I'm trying to use rtracklayer's import function to import a GTF file. I downloaded the current comprehensive genome annotation for human from GENCODE, gunzipped the …
ROSE Algorithm: Error when importing upper and join from string
written 6.1 years ago by mropri

Hello, I am trying to run the ROSE algorithm created by the young lab, url here: http://younglab.wi.mit.edu/super_enhancer_code.html Specifically, I am running the ROSE_main.py script: http://younglab.wi.mit.edu/super_enhancer_code.html …
converting Bed to Wig
written 7.3 years ago by rahel14350 •

Dear all, I want to convert a bed file to wig one for a better visualization on UCSC but I have an error (below). Do …
locateVariants missing some mutation site
written 8.8 years ago by noeD

Hello! I am using Variant Annotation package from Bioconductor to annotate my variant but it missed some variants in the output. For example, in `var` …
RTrackLayer import bed without name and strand orientation
written 8.9 years ago by Seigfried

I have a bed file like this ChrIII_A_nidulans_FGSC_A4 0 10 0.83700 ChrIII_A_nidulans_FGSC_A4 10 43 1.11600 ChrIII_A_nidulans_FGSC_A4 43 59 1.39500 ChrIII_A_nidulans_FGSC_A4 59 65 1.67400 ChrIII_A_nidulans_FGSC_A4 65 …
findOverlapsOfPeaks() of package ChIPpeakAnno in R
written 9.3 years ago by salamandra

When running the following code with ChIPpeakAnno package in R: gr3 = GRanges(seqnames = c("chr1"), strand = c("+","+","+"), ranges = IRanges(start = c(1,3,5), width = …
How to convert the data from .wig file to .table or .csv?
written 10.3 years ago by wql_zcq •

Hi, all I read a data from wig file with import.wig function of rtracklayer package. The class() of this data showed :"UCSCdata". I am failed …