plot a bedgraph
2
2
Entering edit mode
6.6 years ago
sacha ★ 2.4k

Hi ,

I would like to plot a bed file with numeric values ( = bedgraph). My input data are : chr pos start end value ( value between -1 and 1)

exemple :

chr19   50906809        50907015          0.8173737057494912
chr19   50909381        50909562          0.2639404477044525
chr19   50909499        50909702          0.5421707627971464
chr19   50909697        50909863          0.7663584319162735

I would like something like this :

enter image description here

Do you have an easy way ? R package or matplotlib solution ?

bed bedgraph plot RNA-Seq • 8.4k views
ADD COMMENT
2
Entering edit mode

Maybe something interesting in : https://bernatgel.github.io/karyoploter_tutorial/

ADD REPLY
1
Entering edit mode

Haven't tried, maybe something in https://github.com/deeptools/pyGenomeTracks

ADD REPLY
0
Entering edit mode

Python !! This is awesome ! Thanks !

ADD REPLY
0
Entering edit mode

pyGenomeTracks can't plot multiple chromosomes at once

ADD REPLY
2
Entering edit mode
6.6 years ago
sacha ★ 2.4k

I just discovered the Sushi package which do the job pretty well !
https://bioconductor.org/packages/release/bioc/html/Sushi.html
See plotbedgraph

ADD COMMENT
0
Entering edit mode

Hi there, the plot you have in your original post has the whole genome and uses points intead of bars. Would this possible to do in Sushi? I went through the pdf manual but I could not find the exact usage. Best

ADD REPLY
0
Entering edit mode
6.6 years ago

qqman: an R package for visualizing GWAS results using Q-Q and manhattan plots

https://www.biorxiv.org/content/early/2014/05/14/005165

enter image description here

ADD COMMENT
0
Entering edit mode

Thanks pierre. I will try. I forgot to mention that my bed refer to amplicon.

ADD REPLY
0
Entering edit mode

yeah but for ~100bp will you really be able to see a width of a scaled interval ? should be <1 pixel for a human genome.

ADD REPLY

Login before adding your answer.

Traffic: 881 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6