https://github.com/StatisticsHealthEconomics/multimcm/blob/8e7a88dd08053ae64bd8c2673350048eecbbcb9f/R/stan_code_blocks.R#L418 The questions is how to do this quickly and easily. Think keep the loop as-is and change to `S(time[i], ...` I will have to make some small changes to the plotting function to update this.