--- a +++ b/overview/calibration-plots.R @@ -0,0 +1,24 @@ +source('../lib/handymedical.R', chdir = TRUE) +requirePlus('ggplot2', 'cowplot') + +cox.calibration <- + read.csv('../../output/cox-bigdata-varsellogrank-01-calibration-table.csv') +rf.calibration <- + read.csv('../../output/rfsrc-cv-nsplit-try3-calibration-table.csv') + +cox.calibration.plot <- calibrationPlot(cox.calibration, max.points = 2000) +rf.calibration.plot <- calibrationPlot(rf.calibration, max.points = 2000) + +plot_grid( + cox.calibration.plot, rf.calibration.plot, + labels = c("C", ""), + align = "v", ncol = 2 +) + +ggsave( + '../../output/models-calibration-eg.pdf', + width = 16, + height = 8, + units = 'cm', + useDingbats = FALSE +) \ No newline at end of file