## ----------------------------------------------------------------------------- suppressMessages(库(MultiAssayExperiment))稍后通知< - readRDS(执行(“extdata”、“training.tba.toydata。rds”,包= " AffiXcan ")) # #这是MultiAssayExperiment对象名称(化验(稍后通知 )) ## ----------------------------------------------------------------------------- 化验(稍后通知)$ ENSG00000139269.2 [1:5, 1:4 ] ## ----------------------------------------------------------------------------- suppressMessages(库(SummarizedExperiment))负载(“. . /数据/ exprMatrix.RData”)化验(exprMatrix)美元值[1:5 ] ## -----------------------------------------------------------------------------负载(“. . /数据/ regionAssoc.RData”)regionAssoc [1:3 ,] ## ----------------------------------------------------------------------------- 负载(“. . /数据/ trainingCovariates.RData”)头(trainingCovariates ) ## ----------------------------------------------------------------------------- suppressMessages(库(“AffiXcan”))trainingTbaPaths < -执行(“extdata”、“training.tba.toydata。rds"， package="AffiXcan") data(exprMatrix) data(regionAssoc) data(trainingCovariates) assay <- "values" training <- affiXcanTrain(exprMatrix=exprMatrix, assay=assay, tbappaths = trainingtbappaths, regionAssoc=regionAssoc, cov=trainingCovariates, varexained =80, scale=TRUE) testingtbappaths <- system.file("extdata"，"testing.tba.toydata. "rds”,package="AffiXcan") exprmatrix <- affiXcanImpute(tbaPaths=testingTbaPaths, affiXcanTraining=training, scale=TRUE) grexMatrix <- assays(exprmatrix)$GReX as.data.frame(grexMatrix)[,1:5] ## ----------------------------------------------------------------------------- trainingTbaPaths <- system.file("extdata","training.tba.toydata.rds", package="AffiXcan") data(exprMatrix) data(regionAssoc) data(trainingCovariates) assay <- "values" training <- affiXcanTrain(exprMatrix=exprMatrix, assay=assay, tbaPaths=trainingTbaPaths, regionAssoc=regionAssoc, cov=trainingCovariates, varExplained=80, scale=TRUE, kfold=5)