tests/testthat/test_so_central_or.R

In alastairrushworth/smnet: Smoothing for Stream Network Data

context('Testing SO Central OR')

mf_location <- paste0(normalizePath('testdata'), "/SoCentralOR.ssn")
mf          <- suppressWarnings(SSN::importSSN(mf_location, predpts = "preds", o.write = TRUE))
adjacency   <- get_adjacency(mf_location, netID = 1)


# test that maxit works properly
test_that("Linear model with network", {
  # simple gam
  mod_smn     <- 
    smnet(
      formula = ELEV ~ CUMDRAINAG + network(adjacency = adjacency, weight = "areaPI"),
      data.object = mf, netID = 1, control = list(verbose = FALSE))
  expect_equal(class(mod_smn), 'smnet')
  expect_equal(names(summary(mod_smn, verbose = FALSE)), c("linear.terms", "smooth.terms"))
  plot(mod_smn, type = "covariates")
  plot(mod_smn, type = "segments", weight = 0.2, shadow = 2, 
       xlim = c(1060000, 1130000), ylim = c(1380000, 1480000))
  plot(mod_smn, type = "full", weight = 0.2, shadow = 2, 
       xlim = c(1060000, 1130000), ylim = c(1380000, 1480000))
  expect_equal(
    names(predict(mod_smn, newdata = getSSNdata.frame(mf, "preds"))),
    c('predictions', 'predictions.se')
  )
  expect_equal(
    names(predict(mod_smn, newdata = getSSNdata.frame(mf, "Obs"))),
    c('predictions', 'predictions.se')
  )
})