# #——风格,回声= FALSE,结果= '飞机 '--------------------------------- 库(BiocStyle ) ## -------------------------------------------------------------------------- 库(希尔达)inputFile < -执行(“extdata / esophageal.mp.txt.gz”、包=“希尔达”)G < - hildaReadMPFile (inputFile numBases = 5, trDir = TRUE ) ## -------------------------------------------------------------------------- 负载(执行(“extdata /样品。rdata”,包= "希尔达”)(G)类 ) ## -------------------------------------------------------------------------- set.seed (123) hildaGlobal < - hildaTest (inputG = G, numSig = 3, localTest = FALSE, refGroup = 1:4,硝石= 1000)hildaLocal < hildaTest (inputG = G, numSig = 3, localTest = TRUE, refGroup = 1:4,硝石= 1000 ) ## -------------------------------------------------------------------------- 参数< - pmgetSignature (G,K = 3 ) ## -------------------------------------------------------------------------- set.seed (123) hildaGlobal < - hildaTest (inputG = G, numSig = 3, useInits =参数,localTest = TRUE, refGroup = 1:4,硝石= 1000)hildaLocal < hildaTest (inputG = G, numSig = 3, useInits =参数,localTest = TRUE, refGroup = 1:4,硝石= 1000 ) ## -------------------------------------------------------------------------- hildaRhat hildaGlobal hildaRhat (hildaLocal ) ## -------------------------------------------------------------------------- pmPlots < - pmBarplot (G, Param refGroup = 1:4, sigOrder = c(1、3、2))cowplot:: plot_grid (pmPlots sigPlot美元,pmPlots propPlot美元,rel_widths = c (1, 3 )) ## -------------------------------------------------------------------------- hildaPlots < - hildaBarplot (G, hildaLocal refGroup = 1:4,sigOrder = c(1、3、2))cowplot:: plot_grid (pmPlots sigPlot美元,pmPlots propPlot美元,rel_widths = c (1, 3 )) ## -------------------------------------------------------------------------- hildaDiffPlots < - hildaDiffPlot (G, hildaLocal, sigOrder = c(1、3、2))cowplot:: plot_grid (hildaDiffPlots sigPlot美元,hildaDiffPlots diffPlot美元,rel_widths = c(1、3))