# # # R代码从装饰图案的fCCAC来源。Rnw“# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块1号:风格# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # BiocStyle::乳胶()# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块2号:设置# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #选项(宽度= 200)选项(继续= " ")选项(提示=“R >”) # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块3号:fCCACexample # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #选项(宽度= 50);如果(.Platform $ OS。类型= =“windows”) {print (“…rtracklayer无法读取位格式文件在Windows…”)}如果(.Platform $ OS。类型= = {# # hg19“unix”)。chr21:40000000 - 48129895 H3K4me3数据虽然et al . (2015) owd <——setwd (tempdir ());库(fCCAC) bigwig1 <——“chr21_H3K4me3_1。bw chr21_H3K4me3_2 bigwig2 < -。bw chr21_H3K4me3_3 bigwig3 < -。bw chr21_merged_ACT_K4 peakFile < -。床”标签< - c (“H3K4me3”、“H3K4me3”,“H3K4me3”) r1 < -系统。文件(“extdata bigwig1包=“fCCAC mustWork = TRUE) r2 < -系统。文件(“extdata bigwig2包=“fCCAC mustWork = TRUE) r3 < -系统。文件(“extdata bigwig3包=“fCCAC mustWork = TRUE) r4 < -系统。file("extdata", peakFile, package="fCCAC",mustWork = TRUE) ti <- "H3K4me3 peaks (chr21)" fc <- fccac(bar=NULL, main=ti, peaks=r4, bigwigs=c(r1,r2,r3), labels=labels, splines=15, nbins=100, ncan=15) head(fc) setwd(owd) } ################################################### ### code chunk number 4: heatmapPlot ################################################### options(width=50) if (.Platform$OS.type == "windows") { print("...rtracklayer is unable to read bigWig format files in Windows...") } if (.Platform$OS.type == "unix" ){ heatmapfCCAC(fc) } ################################################### ### code chunk number 5: sessionInfo ################################################### sessionInfo()