I repeat the mcmc sampling using the published covariace matrices for the P(k) observations, the results are below

Posterior distribution

Flux Power Spectrum

Thermal Evolution

HI Effective Optical Depth

HeII Effective Optical Depth

HI Fraction

Thompson Optical Depth

Temperature evolution for reduced photoheating