Biostar Beta. Not for public use.
plot a bedgraph
1
Entering edit mode
15 months ago
sacha ♦ 1.7k
France

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 ?

ADD COMMENTlink
2
Entering edit mode

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

ADD REPLYlink
1
Entering edit mode

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

ADD REPLYlink
0
Entering edit mode

Python !! This is awesome ! Thanks !

ADD REPLYlink
1
Entering edit mode
15 months ago
sacha ♦ 1.7k
France

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

ADD COMMENTlink
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 REPLYlink
0
Entering edit mode
16 months ago
France/Nantes/Institut du Thorax - INSE…

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 COMMENTlink
0
Entering edit mode

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

ADD REPLYlink
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 REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.3.1