After fitting the three fake datasets, the three different reciprocal_dispersion values are basically the same at: 10.7, 10.8, and 10.8. I’m a bit surprised, because I thought the phi values would be recovered by the fits, but obviously not. Is there something in the fitting that is supposed to be able to recover the phi values? I’d like to follow the process of recovering the known inputs to the data generating process, but I don’t know how to do that with the available code. Thanks.