[This article was first published on Getting Genetics Done, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.
A couple of months ago I posted about how to visualize exome coverage with bedtools and R. But if you’re looking to get a basic handle on genome arithmetic, take a look at Aaron Quinlan’s bedtools tutorials from the 2013 CSHL course. The tutorial uses data from the Maurano et al exploration of DnaseI hypersensitivity sites in hundreds of primary tissue types (Science 337:1190-1195).Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.
The tutorial provides examples with pictures and code to do things like:
Intersections to find features that overlap (or do NOT overlap):
Merging features like {ChIP,MEDIP,DNAse,*}-seq peaks:
Examining coverage:
Advanced usage using the Jaccard statistic to measure similarity of all 20 tissue types to all other 20 20 tissues, and visualizing this with principal components analysis and ggplot2 in R:
See the full bedtools documentation for more.
2013 CSHL bedtools tutorial: http://quinlanlab.org/tutorials/cshl2013/bedtools.html
Getting Genetics Done by Stephen Turner is licensed under a Creative Commons Attribution-ShareAlike 3.0 Unported License.
To leave a comment for the author, please follow the link and comment on their blog: Getting Genetics Done.
R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.