# #设置,包括= FALSE ----------------------------------------------------- knitr: opts_chunk设置(echo = TRUE) #美元避免加载消息库(NanoMethViz预加载 ) ## ----------------------------------------------------------------------------- 库(NanoMethViz ) ## ----------------------------------------------------------------------------- methy_calls < -系统。file(package = "NanoMethViz"， c(" sample1_nanopolsh .tsv.gz"， " sample2_nanopolsh .tsv.gz")) #看一下methy_data methy_calls_example <- read. gz的前10行。表(methy_calls [1], 9 = " \ t”,标题= TRUE, nrows = 6) methy_calls_example # # = F——消息 --------------------------------------------------------------- methy_tabix < - file.path (tempdir(),“methy_data.bgz”)样本< - c(“sample1”、“sample2”)#您应该看到消息当运行这个自己create_tabix_file (methy_calls methy_tabix,样本)与实际数据# #不这样做,我们必须用gzfile告诉R,我们有一个gzip压缩文件methy_data < -阅读。Table (gzfile(methy_tabix)， col.names = methy_col_names()， nrows = 6