Load Exon Count Data

Description

This function is used to load (sub-)exon count data. Exon count data can be got by the Python script count_in_exons.py.

Usage

1
loadExonCountData(case.files, control.files)

Arguments

case.files

a character vector containing the exon count file names for case samples

control.files

a character vector containing the exon count file names for control samples

Details

You may need the Python script count_in_exons.py (released with this package) to generate your exon count files from read mapping results (say BAM files). The detailed usage can be obtained by simply typing python \path\to\count_in_exons.py. Users can also use other scripts or software for exon read counting.

The format of the exon count file is:

1
2
3
4
5
6
GeneName1:001[tab]Count11
GeneName1:002[tab]Count12
...
GeneName1:00N[tab]Count1N
GeneName2:001[tab]Count21
...

Value

This function returns a ReadCountSet object.

Author(s)

Xi Wang, xi.wang@newcastle.edu.au

See Also

newReadCountSet, ReadCountSet-class

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
library(SeqGSEA)
dat.dir = system.file("extdata", package="SeqGSEA", mustWork=TRUE)
case.pattern <- "^SC"
ctrl.pattern <- "^SN"
case.files <- dir(dat.dir, pattern=case.pattern, full.names = TRUE)
control.files <- dir(dat.dir, pattern=ctrl.pattern, full.names = TRUE)

## Not run: 

RCS <- loadExonCountData(case.files, control.files)
RCS 

## End (Not run)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.