Repitools_vignette的源代码。Rnw“# # #编码:utf - 8 ################################################### ### 代码块1号:loadPkg ################################################### 选项(提示= " ",继续= " ")set.seed(4)库(Repitools ) ################################################### ### 代码块2号:findClusters ################################################### 统计数据。表<- data.frame(chr = "chr1", pos = 1:50 00, t.stat = rnorm(500))表[100:104,"t.stat"] <- rnorm(5,5) stats. >表[200:204,"t.stat"] <- rnorm(5, -5) stats. >集群< - findClusters(统计数据。表,得分。Col = 3, w.size = 5, n.med = 2, n.consec = 3, cut。samps = seq(-2, -10, -2), maxFDR = 0.05, trend = "down", n.perm = 10)。1 <——(统计数据。集群$cluster == 1) stats.clustered[集群。1 , ] ################################################### ### 代码块3号:genomeBlocks ################################################### 空空的< - c(50000、10000)的名字(杆)< - c(“chr1”、“chr2”)genomeBlocks(科、宽度= 5000)