In inbo/dhcurve: Automated Modelling of Diameter Height Curves for Trees

library(ggplot2)
library(assertthat)

Boomsoort: `r unique(SelectieModel$BMS)`; Domein: `r unique(SelectieModel$BOS_BHI)`

DOMEIN_ID: r unique(SelectieModel$DOMEIN_ID)

r unique(SelectieModel$Reden)

aantal bomen = r unique(SelectieModel$nBomen); aantal bomen met omtrek > 0,5 m = r unique(SelectieModel$nBomenOmtrek05); aantal bomen in bruikbaar interval = r unique(SelectieModel$nBomenInterval); aantal bomen in bruikbaar interval met omtrek > 0,5 m = r unique(SelectieModel$nBomenIntervalOmtrek05); aantal bomen boven bruikbaar interval: r unique(SelectieModel$nExtra)

Model bruikbaar tussen de omtrekklassen met klassenmiddens r unique(SelectieModel$Q5k) - 0.3 en r unique(SelectieModel$Q95k) + 0.2 m

DiffMediaan: r round(unique(SelectieModel$DiffMediaan), 2); DiffMin: r round(unique(SelectieModel$DiffMin), 2); DiffMax: r round(unique(SelectieModel$DiffMax), 2)

SelectieModel <- SelectieModel %>%
  mutate(
    VoorModelFit = ifelse(is.na(VoorModelFit), Omtrek < 2, VoorModelFit)
  )

Grafiek <-
  ggplot(SelectieModel,
         aes(x = Omtrek, y = HOOGTE, colour = VoorModelFit)) +
  scale_colour_manual(values = alpha(c("blue", "black"))) +
  geom_jitter(width = 0.03, height = 0) +
  geom_line(aes(y = H_D_finaal, linetype = "Domein"),
            size = 1) +
  geom_point(aes(y = H_D_finaal),
             shape = 21, fill = "white", size = 2, stroke = 1) +
  scale_linetype_manual(name = "Model",
                        values = c("Vlaams" = "dashed", "Domein" = "solid")) +
  # geom_vline(xintercept = SelectieModel$Q5) +
  # geom_vline(xintercept = SelectieModel$Q95) +
  scale_x_continuous(
    minor_breaks = seq(0.05, unique(SelectieModel$Q95k) + 0.2, 0.1)
  )#+
  # geom_vline(xintercept = SelectieModel$Omtrek_Buigpunt, colour = "red") +
  # geom_vline(xintercept = SelectieModel$Omtrek_Extr_Hoogte, colour = "red") #nolint

if (has_name(SelectieModel, "H_VL_finaal")) {
  Grafiek <- Grafiek +
    geom_line(aes(y = H_VL_finaal, linetype = "Vlaams"),
              colour = "black", size = 1)
}

print(Grafiek)