# #设置,回声= FALSE --------------------------------------------------- suppressPackageStartupMessages({库(UseBioconductor)库(GenomicRanges)}) stopifnot (BiocInstaller: biocVersion() = = 3.1” ") ## ---- 风格,回声= FALSE,结果= ' asis '------------------------------- BiocStyle:减价()knitr: opts_chunk设置(整洁=假美元 ) ## ------------------------------------------------------------------------ 库(GenomicRanges) showClass(“GRangesList ") ## ------------------------------------------------------------------------ 库(TxDb.Dmelanogaster.UCSC.dm3.ensGene)txdb < - TxDb.Dmelanogaster.UCSC.dm3。ensGene tx_by_gn <- transcriptsBy(txdb, by="gene") unlisted <- unlist(tx_by_gn) TSS <- ifelse(strand(not listed) == "+"， start(not listed)， end(not listed)) TSS <- GRanges(seqnames(not listed)， IRanges(TSS, width=1)， strand(not listed)) TSS_by_gn <- relist(TSS, tx_by_gn) mcls (TSS) <- mcls (not listed) TSS_by_gn <- relist(TSS, tx_by_gn)