Skip to content

Instantly share code, notes, and snippets.

@al2na
Created October 5, 2012 17:50
Show Gist options
  • Save al2na/3841293 to your computer and use it in GitHub Desktop.
Save al2na/3841293 to your computer and use it in GitHub Desktop.
reading BSMAP ratio files
library(methylKit)
obj=read("bsmap_output.txt",sample.id="test",assembly="mm9",header=TRUE, context="CpG", resolution="base",
pipeline=list(fraction=TRUE,chr.col=1,start.col=2,end.col=2,
coverage.col=6,strand.col=3,freqC.col=5 )
)
obj
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment