Comparison of $f(R,T)$ Gravity with Multiple Datasets

We examine $f(R,T)=R+\lambda T^\epsilon$ gravity models by finding the best fit parameters for various values of $\epsilon$. This is accomplished by analyzing data from the cosmic microwave background, baryon acoustic oscillation observations, cosmic chronometer, and Type Ia supernovae introducing correlations and radiation effects to our previous work. We find the probability distribution for the best fit model around the value of $\epsilon=0.010^{+0.013}_{-0.021}$, consistent with the standard cosmological model value of $\epsilon = 0$.