This model assumes that the outcome Y is normally distributed conditional on S.1 and Z, with mean determined by the model formula. It also assumes that larger values of Y are more indicative of poor outcomes, e.g., blood pressure.
risk_continuous(model = Y ~ S.1 * Z, D = 5000)
Formula specifying the risk model for the mean
number of samples for the simulated annealing integration