R: bisulfite squencing PCR plot
3
2
Entering edit mode
8.8 years ago
Shicheng Guo ★ 9.5k

Hi All,

Do you have R code for the bisulfite sequencing PCR result plotting script.

Thanks

image: plot

BSP R • 2.5k views
ADD COMMENT
2
Entering edit mode
8.8 years ago

There might be a package that will produce those exactly, but if not you can do it in ggplot2:

d <- data.frame(CpG=rep(1:10,4), read=rep(1:4, each=10), value=sample(c(1,16), 40, replace=T))
ggplot() + scale_shape_identity() + geom_point(data=d, aes(x=CpG, y=read, shape=value), size=5)

Note that you'll need to tweak the size, remove the background, etc., but that should serve as a start.

ADD COMMENT
3
Entering edit mode
8.8 years ago
Ying W ★ 4.2k

I made a plotter that will start with PCR data and will also do nucleosome footprinting (NOMe-seq)

https://github.com/ying-w/bioinformatics-figures/tree/master/methylcircleplot

image:fig from paper

ADD COMMENT
1
Entering edit mode
8.8 years ago
Shicheng Guo ★ 9.5k

Thanks. Devon, It would be more perfect with the suitable setting: coord_cartesian

library("ggplot2")
d <- data.frame(CpG=rep(1:10,4), read=rep(1:4, each=10), value=sample(c(1,16), 40, replace=T))
ggplot() + scale_shape_identity() + geom_point(data=d, aes(x=CpG, y=read, shape=value), size=9)+coord_cartesian(ylim = c(0, 5))
ADD COMMENT

Login before adding your answer.

Traffic: 1486 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