context("Test for lets.addvar")
data(PAM)
data(temp)
temp <- unwrap(temp)

tempstack <- c(temp, temp)

test_that("lets.addvar works fine", {
  
  
  PAM_temp_mean <- lets.addvar(PAM, temp)
  expect_true(is.matrix(PAM_temp_mean))
  expect_true(ncol(as.matrix(PAM_temp_mean)) == (ncol(PAM[[1]]) + 1))
  
})

test_that("lets.addvar works fine, different fun", {
  
  
  PAM_temp_mean <- lets.addvar(PAM, temp, fun = sd)
  expect_true(is.matrix(PAM_temp_mean))
  expect_true(ncol(as.matrix(PAM_temp_mean)) == (ncol(PAM[[1]]) + 1))

})

test_that("lets.addvar works fine, onlyvar = TRUE", {
  
  
  PAM_temp_mean <- lets.addvar(PAM, temp, onlyvar = TRUE)
  expect_true(is.matrix(PAM_temp_mean))
  expect_true(ncol(as.matrix(PAM_temp_mean)) == 1)

})

test_that("lets.addvar works fine, multiple rasters", {
  
  
  PAM_temp_mean <- lets.addvar(PAM, tempstack)
  expect_true(is.matrix(PAM_temp_mean))
  expect_true(ncol(as.matrix(PAM_temp_mean)) == (ncol(PAM[[1]]) + 2))
  
})