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")
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.