Peter Boorman

Thanks for getting back to me so quickly! I tried adding the `INSTRUME="XIFU"` and `TELESCOP="ATHENA"` keywords to the primary header (as they were only present in the spectrum header previously),...

Thanks! I replaced the `datasets` line and also converted it to a list to match the form of `datasets` from the previous line. Now I get a new error due...

Thanks! I had a go at tweaking the store_chain() function. A version I wrote that appears to work for models with one data group is: ``` def store_chain_universal(chainfilename, transformations, posterior,...

Thanks! After plotting it, the best-fit for all datasets seems to be below the data systematically.

Thanks a lot for your help, Johannes. Increasing the custom log-Gaussian prior standard deviation allowed sufficient dynamic range for all parameters in the fit to generate a corner plot.