thats looking good
You should try to use either the corner or getdist packages to plot the 'corner' plots commonly used when showing correlated posteriors (cf. https://emcee.readthedocs.io/en/v2.2.1/user/line/) so that we can see what's up with the other uncertainties
This week I have been working towards a Markov-chain Monte Carlo (MCMC) approach for fitting Megastat cooldown data and obtaining estimates on various emissivity values. I started with a simple model, only simulating the radiative cooling from the inner shield to the test mass. I supply the test mass and inner shield temperature data, and the emissivity of Aquadag (coating both the TM and inner surface of IS) as the only fit parameter. I am using Stan for the modeling, and Attachment 1 is a copy of my stan model.
The results of the simple model are in Attachment 2. The emissivity of Aquadag is estimated as a gaussian centered around 0.7. I am still determining whether this result is "real", or if the simulation is simply returning my prior. I will look into this more before adding complexity to the model.