Where is code for the model optimization?

[radfordneal]

The paper and supplemental materials mention that parameter values for the ODE model were found by a combination of latin hypercube sampling (in R) and Nelder-Mead optimization (in Mathematica). But I can't find the code used to do this.

I'm aware of Mathematica code at https://github.com/c2-d2/postpandemic, but that code seems to just simulate future scenarios with specified parameter values.

Thanks.