## ----包括= false -------------------------------------------------------------------------------------------------尼特尔:: opts_chunk $ set（collapse = true，注释=“＃>”）选项（width = 70）## ----设置----------------------------------------------------------------------------------------------------------------------------------------------库（Eisar）＃＃--------------------------------------------------------------------------------------------------- gtf <-system.file（“ extdata/gencode.v28.annotation.sub.gtf”，package =“ eisar”）grl <-getFeatUreranges（gtf = gtf，farteuretype = c（“剪接”，“内含子”），intonType =“ and syte”，fanklength = 50L，joinoverlappingIntrons = false，valse，drebose = true = true）## ------------------------------------------------------------------------------------------------------------------------------- grl ## ---------------------------------------------------------------------------------------------------------------- lapply(S4Vectors::metadata(grl)$featurelist, head) ## -----------------------------------------------------------------------------------------------------头（S4Vectors :: Metadata（Grl）$ Corrgene）## -------------------------------------------------------------------------------------------------------- suppressPackageStartupMessages({ library(BSgenome) library(BSgenome.Hsapiens.ucsc.hg38）}）seqs < - genomicFeatures :: extractTranscriptSeqs（x = bsgenome.hsapiens.ucsc.hg38，tresscripts = grl）------------------------------------------------------------------------------------exportTogtf（grl，filepath = file.path（tempdir（），“ exported.gtf”））## ------------------------------------------------------------------------------------------------------------------ df <-gettx2gene（grl）头（DF））尾巴（DF）## ----------------------------------------------------------------------------------------------------------------------- sessionInfo（）