## ----setup, include=FALSE----------------------------------------------------- knitr::opts_chunk$set(dpi = 300) knitr::opts_chunk$set(cache=FALSE) ## ---- echo = FALSE, hide=TRUE, message=FALSE, warning=FALSE------------------- library('oppti') ## ---- eval = TRUE------------------------------------------------------------- set.seed(1) cohort1.proteomes = as.data.frame(matrix(abs(rnorm(100*30)), 100, 30)) rownames(cohort1.proteomes) = paste('marker', 1:100, sep = '') colnames(cohort1.proteomes) = paste('cohort1.sample', 1:30, sep = '') ## ---- eval = TRUE------------------------------------------------------------- library(oppti) result = oppti(cohort1.proteomes) ## ---- eval = TRUE------------------------------------------------------------- cohort1.outlier.scores = result[[1]] ## ---- eval = TRUE, echo = FALSE, size = 8------------------------------------- knitr::kable(cohort1.outlier.scores[1:10,1:4], digits = 2, caption = "Example matrix of the outlier scores, displayed for the first 10 proteins (rows) and the first 4 samples (columns)", row.names = TRUE) ## ---- eval = TRUE------------------------------------------------------------- cohort1.normal.states = result[[2]] ## ---- eval = TRUE, echo = FALSE, size = 8------------------------------------- knitr::kable(cohort1.normal.states[1:10,1:4], digits = 2, caption = "Example matrix of the normal states", row.names = TRUE) ## ---- eval = TRUE------------------------------------------------------------- cohort1.markers.tests = result[[3]] ## ---- eval = TRUE, echo = FALSE, size = 8------------------------------------- knitr::kable(cohort1.markers.tests[1:10,], digits = 4, caption = "Statistical significance of outlying markers", row.names = TRUE) ## ---- eval = TRUE------------------------------------------------------------- cohort2.proteomes = as.data.frame(matrix(abs(rnorm(80*20)), 80, 20)) rownames(cohort2.proteomes) = paste('marker', 51:130, sep = '') colnames(cohort2.proteomes) = paste('cohort2.sample', 31:50, sep = '') ## ---- eval = TRUE------------------------------------------------------------- result = oppti(list(cohort1.proteomes,cohort2.proteomes)) ## ---- eval = TRUE------------------------------------------------------------- outlier.scores = result[[1]] ## ---- eval = TRUE------------------------------------------------------------- cohort1.outlier.scores = outlier.scores[[1]] ## ---- eval = TRUE, echo = FALSE, size = 8------------------------------------- knitr::kable(cohort1.outlier.scores[1:10,1:4], digits = 2, caption = "Example outlier scores in cohort1", row.names = TRUE) ## ---- eval = TRUE------------------------------------------------------------- cohort2.outlier.scores = outlier.scores[[2]] ## ---- eval = TRUE, echo = FALSE, size = 8------------------------------------- knitr::kable(cohort2.outlier.scores[1:10,1:4], digits = 2, caption = "Example outlier scores in cohort2", row.names = TRUE) ## ---- eval = FALSE------------------------------------------------------------ # result = oppti(list(cohort1.proteomes,cohort2.proteomes), draw.sc.plots = TRUE, # panel.markers = rownames(cohort1.proteomes)[46:55]) ## ---- eval = FALSE------------------------------------------------------------ # result = oppti(list(cohort1.proteomes,cohort2.proteomes), draw.ou.plots = TRUE, # panel.markers = rownames(cohort1.proteomes)[46:55]) ## ---- eval = FALSE------------------------------------------------------------ # result = oppti(list(cohort1.proteomes,cohort2.proteomes), # draw.ou.markers = c('marker50', 'marker55'), # panel.markers = rownames(cohort1.proteomes)[46:55])