Standards for sensitivity analyses in DLNM studies. Use when varying degrees of freedom, maximum lags, confounding adjustments, or performing subgroup analyses.
| Parameter | Default | Sensitivity Range |
|---|---|---|
| df for exposure-response | 3-4 | 2-6 |
| df for lag-response | 3-4 | 2-6 |
| Maximum lag | 21 days | 14, 28, 30 |
| df/year for time trend | 7 | 4, 5, 6, 8, 9, 10 |
| Temperature adjustment | ns(temp, df=4) | df = 3, 5, 6; cross-basis |
# Systematic df sensitivity
results <- lapply(2:6, function(dfvar) {
cb_sens <- crossbasis(exposure, lag = maxlag,
argvar = list(fun = "ns", df = dfvar),
arglag = list(fun = "ns", df = 3))
model_sens <- glm(death ~ cb_sens + covariates, family = quasipoisson(), data = df)
pred_sens <- crosspred(cb_sens, model_sens, cen = cen_value, at = at_vals)
data.frame(df_var = dfvar, rr = pred_sens$allRRfit["30"])
})