# # - eval = TRUE,消息= FALSE ------------------------------------------------ 库(DMCHMM) fn < list.files(系统。file("extdata"，package = "DMCHMM")) fn。f <- list.files(system.file("extdata"，package="DMCHMM")， full.names=TRUE) OBJ <- readBismark(fn。f, fn mc.cores = 2) cdOBJ < - DataFrame(=细胞因子(c(“公元前”、“TC”、“Mono”),标签= c(“公元前”、“TC”、“Mono”)),row.names = c(“BCU1568”、“BCU173”,“BCU551”))colData (OBJ) <——cdOBJ OBJ # # - eval = TRUE,消息= FALSE ------------------------------------------------ nr < - 150;数控< - 8 metht < -矩阵(as.integer (runif (nr *数控,0,20)),nr) methc < -矩阵(rbinom (n = nr *数控,c (metht),概率= runif (nr * nc)), nr, nc) r1 < -农庄(代表(“chr1”,nr) IRanges (1: nr,宽度= 1),链=“*”)的名字(r1) < - 1: nr cd1 < - DataFrame(组=代表(c (G1, G2),每个= nc / 2), row.names =信[1:数控])其中OBJ1 <——cBSData (rowRanges = r1 methReads = methc totalReads = metht colData = cd1)其中OBJ1 # # - eval = TRUE,消息= FALSE ------------------------------------------------ methoda < - methHMEM(其中OBJ1MaxK = 2) methoda # # - eval = TRUE,消息= FALSE ------------------------------------------------ OBJ3 < - methHMMCMC (methoda) OBJ3 # # - eval = TRUE,消息= FALSE ------------------------------------------------ OBJ4 < - findDMCs (OBJ3)头(元数据(OBJ4) DMCHMM美元 ) ## ----------------------------------------------------------------------------- sessionInfo ()