library(PRcalc)
data(jp_upper_2019)
pr_obj1 <- prcalc(jp_upper_2019, m = 50, method = "dt")
pr_obj2 <- prcalc(jp_upper_2019, m = 50, method = "dt", threshold = 0.025)
pr_obj3 <- prcalc(jp_upper_2019, m = 50, method = "dt", threshold = 0.05)
list("t = 0%" = pr_obj1,
"t = 2.5%" = pr_obj2,
"t = 5%" = pr_obj3) |>
compare() |>
plot()