Check coverage on HH output #6
Description
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)
)