tests/testthat/test_gwlZypTest.R

In bcgov/bcgroundwater: Analyse B.C. Groundwater Level Data

context("gwl_zyp_test")

test_that("gwl_zyp_test returns data", {
  g_month <- get_gwl(wells = 309) %>%
    monthly_values() %>%
    make_well_ts()
  
  g_annual <- g_month %>%
    dplyr::select(-yearmonth) %>% 
    dplyr::group_by(EMS_ID, Well_Num, Year) %>%
    dplyr::summarize(nReadings = length(Well_Num),
                     mean_GWL = mean(med_GWL),
                     SD = sd(med_GWL),
                     med_GWL = median(med_GWL),
                     q95_GWL = quantile(med_GWL, 0.95))
  
  expect_silent(t <- gwl_zyp_test(g_annual, byID = "Well_Num", col = "mean_GWL"))
  expect_is(t, "data.frame")
  expect_gt(nrow(t), 0)
  expect_true(all(!is.na(t$intercept)))
  
  expect_silent(p <- gwl_area_plot(g_month, 
                                 trend = t$trend[1],
                                 intercept = t$intercept[1], sig = t$sig[1],
                                 trend_category = "Stable", mkperiod = "annual", showInterpolated = TRUE))
  expect_is(p, "ggplot")
})