库（DiffBind）Diffbind < - system.file（包=“diffbind”，“extra”）dir（diffbind）csv < - read.csv（file.path（diffbind，“tamoxifen.csv”））vls < - dir（File.Path（DiffBind，“Peaks”），Full = True）名称（FL）< - 子（“Bed.gz”，“。，BaseName（FLS））长度（FLS）库（Rtracklayer）Peak1 < - 导入（FLS [1]）Peak1 hist（宽度（峰值1））峰值< - lapply（vls，import）所有< -​​ refile（do.call（c，Unname（峰）））占用< - sapply（峰值，倒数，查询=全部）HEAD（占用）COLSUMS（占用）绘图（表（Rowsums（占用！= 0）））库（TXDB.HSAPIENS.CUCSC.HG18.knowngene）GN < - 基因（TXDB.HSAPIENS.CUCSC.HG18.knowngene）CHR18 < - GN [SEQNAMES（GN）==“CHR18”]总和（倒数（全部，CHR18）！= 0）##与SUM相同（倒数（倒数（全部，CHR18，TYPE =“内部”）！= 0）##与SUM相同（％CHR18内）上游<侧面（CHR18,5000）和（全部％up％上游）下游< - 侧面（CHR18,5000，start = false）总和（均超过百分比下游）附近< - 调整大小（CHR18，宽度（CHR18）+ 10000，FIX =“CENTER”）总和（靠近近乎％））dist < - distancetonearest（全部[全百分比上游），开始）图（密度（mcols（dist）$距离））mcols（全部）$ is_near_gene < - ally int％附近百分比