## ----message=FALSE, warning=FALSE, cache=TRUE, eval=FALSE--------------------- # if (!requireNamespace("BiocManager", quietly = TRUE)) { # install.packages("BiocManager") # } # if (!requireNamespace("CytoMethIC", quietly = TRUE)) { # BiocManager::install("CytoMethIC") # } ## ----cyto2, message=FALSE, warning=FALSE-------------------------------------- library(CytoMethIC) library(ExperimentHub) library(sesame) sesameDataCache() ## ----cyto3, result="asis", echo=FALSE----------------------------------------- library(knitr) kable(cmi_models[, c("ModelID", "PredictionLabelDescription")], caption = "CytoMethIC supported models" ) ## ----cyto4, message=FALSE----------------------------------------------------- cmi_predict(sesameDataGet("HM450.1.TCGA.PAAD")$betas, ExperimentHub()[["EH8395"]], lift_over=TRUE) cmi_predict(sesameDataGet("HM450.1.TCGA.PAAD")$betas, ExperimentHub()[["EH8396"]], lift_over=TRUE) ## ----cyto5, message=FALSE----------------------------------------------------- cmi_predict(sesameDataGet("HM450.1.TCGA.PAAD")$betas, ExperimentHub()[["EH8422"]]) ## ----cyto6, message=FALSE----------------------------------------------------- cmi_predict(sesameDataGet("HM450.1.TCGA.PAAD")$betas, ExperimentHub()[["EH8421"]]) ## ----cyto7, message=FALSE----------------------------------------------------- cmi_predict(sesameDataGet("HM450.1.TCGA.PAAD")$betas, ExperimentHub()[["EH8423"]]) ## ----cyto8-------------------------------------------------------------------- base_url <- "https://github.com/zhou-lab/CytoMethIC_models/raw/main/models" cmi_model <- readRDS(url(sprintf("%s/Race3_rfcTCGA_InfHum3.rds", base_url))) betas <- openSesame(sesameDataGet("EPICv2.8.SigDF")[[1]]) cmi_predict(betas, cmi_model, lift_over=TRUE) ## ----------------------------------------------------------------------------- sessionInfo()