Skip to content

Instantly share code, notes, and snippets.

Embed
What would you like to do?
mydir <- "/data2/minion/R7/Ebola/PRJEB10571_ERR1014225/014370/pass/"
meta <- read.meta.info(dir=mydir, path.t="/Analyses/Basecall_2D_000/", path.c="Analyses/Basecall_2D_000/")
yield <- plot.cumulative.yield(meta)
plot.length.histogram(meta)
meta.s <- summarise.by.channel(meta)
# the number of times the channel appears in any context
plot.channel.summary(meta.s)
# cumulative 2D length
plot.channel.summary(meta.s, report.col="l2d")
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment