R codes for a flexible Monte Carlo sensitivity analysis approach for unmeasured confounding with multiple treatments and a binary outcome.