model_lr <- glm_bw(formula = Radiation ~ Pain + factor(Satisfaction) +
rcs(Tampascale,3) + Age + Duration + JobControl + JobDemands + SocialSupport,
data=lbpmilr_dev, p.crit = 0.05)
if (FALSE) {
stab_res <- stab_single(model_lr, start_model = TRUE, nboot=20, p.crit=0.05)
stab_res$bif
stab_res$bif_perc
stab_res$model_stab
}
Run the code above in your browser using DataLab