tests/testthat/test.paleocar_models.R

In bocinsky/paleocar: Paleoclimate Reconstruction from Tree Rings using Correlation Adjusted corRelation

library(raster)
testthat::context("paleocar tests")

testthat::test_that("Test that the paleocar function works", {
  
  # Force Raster to load large rasters into memory
  raster::rasterOptions(chunksize=2e+07,maxmemory=2e+08)
  
  ## Set the calibration period
  # Here, we use a 60 year period ending at 1983 
  # to maximize the number of dendro series.
  calibration.years <- 1924:1983
  
  ## Set the retrodiction years (AD)
  prediction.years <- 500:1400
  
  ######## END PARAMETERS
  
  ## Load spatial polygon for the boundary of Mesa Verde National Park in southwestern Colorado:
  data(mvnp)
  
  ## Get Tree-ring data from the ITRDB for 10-degree buffer around MVNP
  data(itrdb)
  
  ## Get 1/3 arc-second PRISM data for the VEPII north study area
  data(mvnp_prism)
  
  ## Run PaleoCAR over the MVNP brick
  mvnp_recon <- paleocar(predictands = mvnp_prism,
                         label = "mvnp_prism",
                         chronologies = itrdb,
                         calibration.years = calibration.years,
                         prediction.years = prediction.years,
                         out.dir = "./",
                         force.redo = T,
                         verbose = T)
  
  testthat::expect_is(mvnp_recon,"list")
  
  testthat::expect_named(mvnp_recon)
  
  unlink("mvnp_prism.models.Rds")
  unlink("mvnp_prism.prediction.Rds")
})