## ---- echo = FALSE------------------------------------------------------- options(width = 100) ## ------------------------------------------------------------------------ library(iCOBRA) data(cobradata_example) class(cobradata_example) cobradata_example ## ------------------------------------------------------------------------ cobradata_example <- calculate_adjp(cobradata_example) ## ------------------------------------------------------------------------ cobraperf <- calculate_performance(cobradata_example, binary_truth = "status", cont_truth = "logFC", splv = "none", maxsplit = 4) slotNames(cobraperf) ## ------------------------------------------------------------------------ cobraplot <- prepare_data_for_plot(cobraperf, colorscheme = "Dark2", facetted = TRUE) ## ---- fig.width = 7, fig.height = 4, fig.cap = "", warning = FALSE------- plot_tpr(cobraplot) ## ---- fig.width = 7, fig.height = 4, fig.cap = "", warning = FALSE------- plot_fdrtprcurve(cobraplot) ## ---- fig.width = 7, fig.height = 4, fig.cap = "", warning = FALSE------- plot_overlap(cobraplot) ## ---- fig.width = 7, fig.height = 5, warning = FALSE--------------------- cobraperf <- calculate_performance(cobradata_example, binary_truth = "status", cont_truth = "status", splv = "expr_cat") cobraplot <- prepare_data_for_plot(cobraperf, colorscheme = "Dark2", facetted = TRUE) plot_tpr(cobraplot) ## ---- fig.width = 7, fig.height = 5, warning = FALSE--------------------- library(ggplot2) pp <- plot_tpr(cobraplot, stripsize = 7.5, pointsize = 3) pp + theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust = 1, size = 10), axis.text.y = element_text(size = 10), axis.title.x = element_text(size = 10), axis.title.y = element_text(size = 10)) ## ---- fig.width = 7, fig.height = 4, warning = FALSE--------------------- pp + theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust = 1, size = 10), axis.text.y = element_text(size = 10), axis.title.x = element_text(size = 10), axis.title.y = element_text(size = 10), legend.position = "bottom") + facet_wrap(~splitval, nrow = 1) ## ---- fig.width = 7, fig.height = 5, warning = FALSE--------------------- plot_overlap(cobraplot) plot_overlap(cobraplot, cex = c(1, 0.7, 0.7)) ## ---- fig.width = 7, fig.height = 5, warning = FALSE--------------------- cobraplot <- prepare_data_for_plot(cobraperf, colorscheme = c("blue", "green", "pink"), facetted = TRUE) pp <- plot_tpr(cobraplot, stripsize = 7.5, pointsize = 3) pp + theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust = 1, size = 10), axis.text.y = element_text(size = 10), axis.title.x = element_text(size = 10), axis.title.y = element_text(size = 10)) ## ---- eval = FALSE------------------------------------------------------- ## COBRAapp(cobradata_example) ## COBRAapp() ## ------------------------------------------------------------------------ COBRAData_to_text(cobradata = cobradata_example, truth_file = "cobradata_truth.txt", result_files = "cobradata_results.txt", feature_id = "feature") ## ------------------------------------------------------------------------ cobra <- COBRAData_from_text(truth_file = "cobradata_truth.txt", result_files = "cobradata_results.txt", feature_id = "feature") cobra