Sam had sent me this issue a few days ago, and this is exactly what we’re trying to decide now. See issue here.
Have you all landed on a set of outputs coming from a wrapper fitting function? I really liked the idea of passing back the stan_args, but also think we need to pass back the input data with the correct metadata mapped to it (e.g. the dates for example).
I would also be interested in making the summarized output in a format readily usable for scoringutils, we haven’t quite gotten to any evaluation modules in the package yet but intend to.