From 19da0ed95d01694bd7577be01908db02b7e0138f Mon Sep 17 00:00:00 2001
From: ppoyk <138587496+ppoyk@users.noreply.github.com>
Date: Wed, 12 Jun 2024 14:15:46 +0300
Subject: [PATCH] add 95% vlines to coord uncertainty effect plots

---
 00-base.Rmd | 16 ++++++++++++++++
 1 file changed, 16 insertions(+)

diff --git a/00-base.Rmd b/00-base.Rmd
index 1c2fa41..1674e6b 100644
--- a/00-base.Rmd
+++ b/00-base.Rmd
@@ -903,6 +903,10 @@ points(density(shifted_dtm_elev_maxerr$dtm_elev - kaikki_paikat$dtm_elev,
 text(x=0.4, y=3, paste("n =",
        length(na.omit(shifted_dtm_elev_maxerr$dtm_elev - kaikki_paikat$dtm_elev))))
 abline(v = 0, col="gray", lty=2)
+abline(v = Quantile(kaikki_paikat$dtm_elev - shifted_dtm_elev_maxerr$dtm_elev,
+           probs = c(.025,.975), na.rm = T), lty=3, lwd=1.5, col="black")
+abline(v = Quantile(kaikki_paikat$km2_elev - shifted_km2_elev_maxerr$elev,
+           probs = c(.025,.975), na.rm = T), lty=3, lwd=1.5, col="dodgerblue")
 legend(0.5, 15, legend=c("Lidar-DTM","KM2"), col=c("black","dodgerblue"), lty=1,lwd=2)
 dev.off()
 
@@ -938,6 +942,10 @@ points(density(shifted_dtm_elev_meanerr$dtm_elev - kaikki_paikat$dtm_elev,
 text(x=0.4, y=4, paste("n =",
        length(na.omit(shifted_dtm_elev_meanerr$dtm_elev - kaikki_paikat$dtm_elev))))
 abline(v = 0, col="gray", lty=2)
+abline(v = Quantile(kaikki_paikat$dtm_elev - shifted_dtm_elev_meanerr$dtm_elev,
+           probs = c(.025,.975), na.rm = T), lty=3, lwd=1.5, col="black")
+abline(v = Quantile(kaikki_paikat$km2_elev - shifted_km2_elev_meanerr$elev,
+           probs = c(.025,.975), na.rm = T), lty=3, lwd=1.5, col="dodgerblue")
 legend(0.5, 25, legend=c("Lidar-DTM","KM2"), col=c("black","dodgerblue"), lty=1,lwd=2)
 dev.off()
 
@@ -973,6 +981,10 @@ points(density(shifted_dtm_elev_mederr$dtm_elev - kaikki_paikat$dtm_elev,
 text(x=0.35, y=9, paste("n =",
        length(na.omit(shifted_dtm_elev_mederr$dtm_elev - kaikki_paikat$dtm_elev))))
 abline(v = 0, col="gray", lty=2)
+abline(v = Quantile(kaikki_paikat$dtm_elev - shifted_dtm_elev_mederr$dtm_elev,
+           probs = c(.025,.975), na.rm = T), lty=3, lwd=1.5, col="black")
+abline(v = Quantile(kaikki_paikat$km2_elev - shifted_km2_elev_mederr$elev,
+           probs = c(.025,.975), na.rm = T), lty=3, lwd=1.5, col="dodgerblue")
 legend(0.5, 25, legend=c("Lidar-DTM","KM2"), col=c("black","dodgerblue"), lty=1,lwd=2)
 dev.off()
 
@@ -1008,6 +1020,10 @@ points(density(shifted_dtm_elev_dblmaxerr$dtm_elev - kaikki_paikat$dtm_elev,
 text(x=0.7, y=2.5, paste("n =",
        length(na.omit(shifted_dtm_elev_dblmaxerr$dtm_elev - kaikki_paikat$dtm_elev))))
 abline(v = 0, col="gray", lty=2)
+abline(v = Quantile(kaikki_paikat$dtm_elev - shifted_dtm_elev_dblmaxerr$dtm_elev,
+           probs = c(.025,.975), na.rm = T), lty=3, lwd=1.5, col="black")
+abline(v = Quantile(kaikki_paikat$km2_elev - shifted_km2_elev_dblmaxerr$elev,
+           probs = c(.025,.975), na.rm = T), lty=3, lwd=1.5, col="dodgerblue")
 legend(0.6, 5, legend=c("Lidar-DTM","KM2"), col=c("black","dodgerblue"), lty=1,lwd=2)
 dev.off()