
## S3 method for class 'RNAseq':
findIslands(obj, max.gap = integer(1), min.cov = 1L,
min.length = integer(1), plot = TRUE, ...)
RNAseq
RNAseq
object with the readIsland slot set with a
RangedData containing the selected islands and the readCount slot
actualized with a list containing the count table per island.hist
## NOTE that this function might need to be actualized
obj <- new('RNAseq',
organismName="Dmelanogaster",
readLength=36L,
chrSize=as.list(seqlengths(Dmelanogaster))
)
obj <- fetchCoverage(
obj,
format="bam",
filename=system.file(
"extdata",
"ACACTG.bam",
package="RnaSeqTutorial")
)
obj <- findIslands(
obj,
max.gap=10L,
min.cov=10L,
min.length=200L)
Run the code above in your browser using DataLab