test_that("plot method for an object of class PredictLoad resturns expected plot structure.", {
burd <- ReadInDataFromR(x.C = burdRC, x.Q = burdRQ)
loaddata <- CreateData(Q = burd$Q, Conc = burd$Conc,
date.range = list(model = c("1973-12-02", "2015-06-30"),
pred = c("1973-12-02", "2015-06-30"),
hour = FALSE),
samp.unit = "day", Ytype = "WY", Qflush = 0.9,
Reg = list(type = "none", rainfall = NULL, date = NULL))
mod1 <- FitModel(x = loaddata$CQ, parms = list(flow = "quadratic", seasonal = TRUE,
RFlimb = FALSE,
MA = c(MA1day = FALSE, MA2days = FALSE,
MAweek = TRUE, MAmonth = TRUE,
MA6months = TRUE, MA12months = TRUE),
trend = FALSE, correlation = FALSE))
predLoad <- predictL(object = mod1, objfix = mod1$gam, x = loaddata,
flow.error = list(me = 0, ce = 0),
samp.unit = "day", pvalue = 0.2)
mod1pL <- plot(predLoad$loadest, type = "annual", Conc = "TSS", scale = "Mt")
# Expectations
expect_type(mod1pL, "list")
expect_named(mod1pL, c("pL1", "pFWC1"))
expect_s3_class(mod1pL$pL1, "arrangelist")
expect_s3_class(mod1pL$pFWC1, "gg")
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.