Closed khufkens closed 2 years ago
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) )
For @lauramarques, checks coverage of GPP values on HH p-model formatted output