## ----消息= false -----------------------------------------------------------------------------------------------------------库（minfi）minfidatadir <-system.file（“ extdata”，package =“ minfidata”）目标< - read.metharray.sheet.sheet（minfidatadir，verbose = =false）rgset < - read.metharray.exp（targets = targets）## ---消息= false ---------------------------------------------------------------------------------------------------------------------------------------------------------库（Flowsorted.blood.450k）库（MCSEA）数据（MCSEADATA）Cellcounts= estionAteCellCounts（rgset）print（cellcounts）## ----消息= false，结果='hide'----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------库（MCSEA）数据（MCSEADATA）## -------------------------------------------------------------------------------------------------------------------------------------------------- 头（betatest，3）print（表词）## --------------------------------------------------------------------------------------------------------------------------------------------------------------------- myRank <- rankProbes(betaTest, phenoTest,refgroup =“控制”）## ---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------警告= false---------------------------------------------------------------------------------------------- myResults <- mCSEATest(myRank, betaTest, phenoTest, regionsTypes = "promoters", platform = "EPIC") ## ----------------------------------------------------------------------------- ls(myResults) ## ----------------------------------------------------------------------------- head(myResults[["promoters"]][,-7]) ## ----------------------------------------------------------------------------- head(myResults[["promoters_association"]], 3) ## ----------------------------------------------------------------------------- head(assocGenes450k, 3) ## ---- message = FALSE, results='hide'----------------------------------------- mCSEAPlot(myResults, regionType = "promoters", dmrName = "CLIC6", transcriptAnnotation = "symbol", makePDF = FALSE) ## ----------------------------------------------------------------------------- mCSEAPlotGSEA(myRank, myResults, regionType = "promoters", dmrName = "CLIC6") ## ----------------------------------------------------------------------------- # Explore expression data head(exprTest, 3) # Run mCSEAIntegrate function resultsInt <- mCSEAIntegrate(myResults, exprTest, "promoters", "ENSEMBL") resultsInt ## ----sessionInfo, echo=FALSE-------------------------------------------------- sessionInfo()