require(dplyr)
require(stringr)
require(purrr)
require(doParallel)
# define input data
data("chaco")
m <- gam(
ndvi ~
te(X,Y, by = tmax) +
s(X,Y, by = pr),
data = chaco,
method = "REML",
family = gaussian()
)
# examine the effect size
effect_size(m, 3)
Run the code above in your browser using DataLab