data <- data.frame(tmin = rnorm(10), tmin.l1 = rnorm(10),
pdsi = rnorm(10), urban = rnorm(10))
# Extract covariate names
covs <- extract_names(data, pattern = c("tmin", "pdsi", "urban"))
# Combine "tmin" and "pdsi" into all possible pairings
cov_multi(covariates = covs, pattern = c("tmin", "pdsi"))
# Combine "tmin" and "urban", treating "urban" as an exact match
cov_multi(covariates = covs, pattern = "tmin", name = "urban")
# Use output as input to write_inla_formulas()
combined_covs <- cov_multi(covariates = covs, pattern = c("tmin", "pdsi"))
formulas <- write_inla_formulas(outcome = "cases", covariates = combined_covs)
Run the code above in your browser using DataLab