View source: R/BioTIP_update_04202022.R
getReadthrough | R Documentation |
The getReadthrough
() function is used to find long transcripts that cover more
than two coding regions for gene regions of interst.
getReadthrough(gr, cod_gr)
gr |
A GRanges object that shows the start and end loci on genome. |
cod_gr |
A GRanges object contaning coding regions. |
For details of findOverlaps, type.partialOverlap, type.50Overlap type.toPlot, queryhits, readthrough and subjecthits see, GenomicRanges https://www.bioconductor.org/packages/release/bioc/html/GenomicRanges.html, IRanges https://www.bioconductor.org/packages/release/bioc/html/IRanges.html, and BiocManager http://bioconductor.org/install/index.html.
A GRanges object that returns overlapping regions of the classified transcript biotypes.
Replace the path_file when loading data locally to the data directory.
Zhezhen Wang and Biniam Feleke
Reference GRCh37 genome https://www.gencodegenes.org/human/release_25lift37.html. For details on gtf format visit ensemble https://useast.ensembl.org/info/website/upload/gff.html.
Wang, Z. Z., J. M. Cunningham and X. H. Yang (2018).'CisPi: a transcriptomic score for disclosing cis-acting disease-associated lincRNAs.' Bioinformatics34(17): 664-670'
#First Load datasets and libraries
library(GenomicRanges)
data("gencode")
data("ILEF")
data("cod")
# Assigning datasets a GRanges object
gencode_gr = GRanges(gencode)
ILEF_gr = GRanges(ILEF)
cod_gr = GRanges(cod)
getReadthrough(ILEF_gr, cod_gr)
## Not run: getReadthrough(cod_gr)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.