# #设置,包括= FALSE ---------------------------------------------- knitr: opts_chunk美元集(崩溃= TRUE,发表评论 = "#>" ) ## ---- showDataExample、消息= FALSE警告= FALSE ----------------------- 库(PowerExplorer)数据(“exampleProteomicsData”)头(exampleProteomicsData dataMatrix美元 ) ## ------------------------------------------------------------------------ 显示(exampleProteomicsData groupVec美元 ) ## ------------------------------------------------------------------------ colnames (exampleProteomicsData dataMatrix美元)# #——运行,呼应= TRUE, fig.keep =‘没有’,消息= FALSE,warning=FALSE---- library(PowerExplorer) data("exampleProteomicsData") res <- estimatePower(inputObject = exampleProteomicsData$dataMatrix, groupVec = exampleProteomicsData$groupVec, islogtransform =FALSE, dataType = "Proteomics"， minLFC = 0.5, enableROTS = TRUE, paraROTS = list(B = 1000, K = NULL)， alpha = 0.05, ST = 50, seed = 345, showProcess =FALSE，saveResultData = FALSE ) ## ---- 情节估计的力量 ------------------------------------------------ plotEstPwr (res ) ## ------------------------------------------------------------------------ res ## ------------------------------------------------------------------------ 数据(exampleObject) listEstPwr (exampleObject减少= TRUE, = 10 ) ## ------------------------------------------------------------------------ listEstPwr (exampleObject选择= c(“ENSMUSG00000000303”,"ENSMUSG00000087272"， "ENSMUSG00000089921")) ## ----预测run2, message=FALSE, warning=FALSE----------------------- data("exampleProteomicsData") res <- predictPower(inputObject = res, groupVec = exampleProteomicsData$groupVec, islogtransform =FALSE, dataType = "Proteomics"， minLFC = 0.5, rangeSimNumRep = c(5,10,15,20)， enableROTS = TRUE, paraROTS = list(B = 1000, K = NULL)， alpha = 0.05, ST = 50，种子= 345)# #——LinePlot ------------------------------------------------------------ plotPredPwr (res LFCscale = 0.5 ) ## ------------------------------------------------------------------------ res ## ------------------------------------------------------------------------ listPredPwr (exampleObject减少= TRUE, = 10 ) ## ------------------------------------------------------------------------ listPredPwr (exampleObject选择= c(“ENSMUSG00000000303”、“ENSMUSG00000087272”,“ENSMUSG00000089921”))