24 lines
706 B
R
24 lines
706 B
R
|
# calibration plot
|
||
|
|
|
||
|
|
|
||
|
|
ds <- openxlsx2::read_xlsx("1 PA Decline/Fra DDV/calibration_imp.xlsx")
|
||
|
|
|
||
|
|
p <- ds |> split(ds$model) |>
|
||
|
|
purrr::imap(\(.x,.i){
|
||
|
|
.x |> predtools::calibration_plot(obs="y",pred="pred")|>
|
||
|
|
purrr::pluck("calibration_plot")+
|
||
|
|
ggplot2::ggtitle(.i)+
|
||
|
|
ggplot2::scale_x_continuous(breaks=seq(0,1,.25),limits=c(0,1))+
|
||
|
|
ggplot2::scale_y_continuous(breaks=seq(0,1,.25),limits=c(-.1,1.1))
|
||
|
|
}) |>
|
||
|
|
patchwork::wrap_plots(ncol=2)
|
||
|
|
|
||
|
|
|
||
|
|
ggplot2::ggsave(filename = "1 PA Decline/calibration_imp.pdf",
|
||
|
|
plot=p,
|
||
|
|
units = "mm",
|
||
|
|
width = 200,
|
||
|
|
height = 100,
|
||
|
|
# pointsize = 30,
|
||
|
|
dpi = 1200)
|
||
|
|
|