Hi All,

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

Thanks

[<img alt="image: plot" src="http://www.springerimages.com/img/Images/BMC/VOL=2011.3/ISU=1/ART=4/MediaObjects/MEDIUM_13148_2011_4_Fig4_HTML.jpg" style="height:448px; width:535px" />][1]

 [1]: http://www.springerimages.com/img/Images/BMC/VOL=2011.3/ISU=1/ART=4/MediaObjects/MEDIUM_13148_2011_4_Fig4_HTML.jpg
There might be a package that will produce those exactly, but if not you can do it in ggplot2:

```r
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.