# # - eval = FALSE -------------------------------------------------------------- # 如果(!requireNamespace(“BiocManager”,悄悄地= TRUE)) # install.packages (BiocManager) # BiocManager::安装(“biscuiteer ") ## ---- eval = FALSE -------------------------------------------------------------- # 如果(!requireNamespace(“BiocManager”,悄悄地= TRUE)) # install.packages (BiocManager) # BiocManager::安装(“trichelab / biscuiteerData”)# BiocManager::安装(“trichelab / biscuiteer ") ## ----------------------------------------------------------------------------- 库(biscuiteer) orig_bed < -系统。file("extdata"， "MCF7_Cunha_chr11p15.bed.gz"， package="饼干")文件(“extdata”、“MCF7_Cunha_header_only.vcf.gz”、包=“biscuiteer”)包钢< readBiscuit (BEDfile = orig_bed VCFfile = orig_vcf合并= FALSE ) ## ----------------------------------------------------------------------------- 包钢biscuitMetadata ( ) ## ----------------------------------------------------------------------------- 元数据(包钢)vcfHeader美元 ## ----------------------------------------------------------------------------- shuf_bed < -系统。file("extdata"， " mcf7_cunha_chr11p15_shuffed .bed.gz"， package="饼干")shuf_vcf <- system. file(" mcf7_cunha_chr11p15_shuffle .bed.gz")文件(“extdata”、“MCF7_Cunha_shuffled_header_only.vcf.gz”、包=“biscuiteer”)bisc2 < - readBiscuit (BEDfile = shuf_bed VCFfile = shuf_vcf合并= FALSE)梳< -成立工会(包钢,bisc2 ) ## ----------------------------------------------------------------------------- epibed。系统。文件("extdata"， "hct116. name .epiread.gz"， package="饼干")Bsseq <- system。file("extdata"， "hct116.bsseq.epiread.gz"， package="饼干")epibed. name .gr <- readEpibed(epibed = epibed. gz)省,是。nome = TRUE, genome = "hg19"， chr = "chr1") epibed.bsseq.gr <- readEpibed(epibed = epibed. gr)bsseq基因组=“hg19”,对应= " chr1 ") ## ----------------------------------------------------------------------------- reg < -农庄(seqnames =代表(“chr11”,5),链=代表(“*”,5),范围= IRanges(开始= c (e7 e6 0, 2.8, 1.17, 1.38 e7, 1.69 e7),结束= c (1.17 2.8 e6、e7, 1.38 e7, 1.69 e7, 2.2 e7)))压裂< - getLogitFracMeth(包钢,minSamp = 1, r =注册)压裂 ## ----------------------------------------------------------------------------- 年龄< - WGBSage(梳子,“霍”)