parseEvent: Parse ATR events from genome annotation

Description Usage Arguments Value References

View source: R/prepare.R

Description

This function parses ATR events (including AS, ATI, and APA) from genome annotation. It also identifies location features for each event. The latter task could be computationally demanding for unfiltered (raw) genome annotation.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
parseEvent(
  anno.file,
  anno.format = tools::file_ext(anno.file),
  cores = max(1, detectCores() - 2),
  min.event.length = 6,
  location.feature = TRUE,
  depth.exon = 100,
  depth.intron = 300,
  remove.duplicate = TRUE,
  verbose = FALSE
)

Arguments

anno.file

character, directory to genome annotation file.

anno.format

character, the format of the annotation file. The format can be inferred from anno.file automatically, unless it is not implicit by the extension.

cores

integer, number of computing cores to use.

min.event.length

numeric (positive), minimum length of a valid event.

location.feature

logical, whether (default to TRUE) to add location features for each event. This usually takes the longest time in annotation parsing procedure.

depth.exon

integer, extended depth into exon, default 50 nt.

depth.intron

integer, extended depth into gene's flanks, default to 300 nt.

remove.duplicate

logical, whether (default to TRUE) to remove identical event duplicates (by keeping one).

verbose

logical, whether (default to TRUE) to print out progress.

Value

a surf object with ATR event annotation and updated genePartsList slot.

References

Chen, F., & Keles, S. (2020). SURF: integrative analysis of a compendium of RNA-seq and CLIP-seq datasets highlights complex governing of alternative transcriptional regulation by RNA-binding proteins. Genome Biology, 21(1), 1-23.


fchen365/surf documentation built on June 18, 2021, 12:02 p.m.