Comapred to the results from my one-cell analytical model:

From an MCMC sampling over the scale_He parameter the posterior distribution is shown below, vertical lines show the mean and the 2-sigma interval:

If I sample T0 from the the scale_He distribution obtained and compare to the data:

Finally, for the fit I only compared against the Gaikwad T0 data since the effect of changing scale_He in the optical depth is not very significant but I’ll extend to fit against \(\tau_{eff}\) when doing the multi-parameter fit: