extractRange: extractRange: Extract subset from object where records lie in...

Description Usage Arguments Value Author(s) Examples

Description

Searches in object for data which lie inside the given range and returns an object of same type containing extracted data.

Usage

1

Arguments

object

gapSites or cRanges (or derived). Object inside which the data is searched for.

seqid

character. Reference sequence identifier.

start

numeric. Start position of given range.

end

numeric. End position of given range.

Value

Same type as object

Author(s)

Wolfgang Kaisers

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
# A) Read gapSites
bam<-system.file("extdata","rna_fem.bam",package="spliceSites")
reader<-bamReader(bam,idx=TRUE)
ga<-alignGapList(reader)
bamClose(reader)
# B) Load refGenome
ucf<-system.file("extdata","uc_small.RData",package="spliceSites")
uc<-loadGenome(ucf)
# C) For gapSites
extractRange(ga,seqid="chr1",start=14000,end=30000)
# D) For cRanges
lj<-lJunc(ga,featlen=3,gaplen=6,strand='+')
extractRange(lj,seqid="chr1",start=14000,end=30000)

spliceSites documentation built on May 6, 2019, 3:05 a.m.