Check coverage on HH output

[khufkens]

For @lauramarques, checks coverage of GPP values on HH p-model formatted output

# quality control on coverage of site level HH data

library(tidyverse)

df <- readRDS("site_based_drivers_HH.rds")

stats <- df %>%
  group_by(sitename, date) %>%
  summarize(
    na_values = length(which(is.na(gpp)))/length(gpp)
  ) %>%
  group_by(sitename) %>%
  summarize(
    coverage = 1 - mean(na_values)
  )