library("RnaSeqTutorial")
library(BSgenome.Dmelanogaster.UCSC.dm3)
## creating a count table from 4 bam files
count.table <- easyRNASeq(filesDirectory=
system.file(
"extdata",
package="RnaSeqTutorial"),
pattern="[A,C,T,G]{6}\\.bam$",
format="bam",
readLength=36L,
organism="Dmelanogaster",
chr.sizes=as.list(seqlengths(Dmelanogaster)),
annotationMethod="rda",
annotationFile=system.file(
"data",
"gAnnot.rda",
package="RnaSeqTutorial"),
count="exons")
Run the code above in your browser using DataLab