We achieve a similar but improved result compared to MS14 using a more involved comparison.
We use a larger \(k\) -range, fitting out to \(k=1\ h/\rm {Mpc}\)  using the matter power spectrum from simulations.
We also employ a more diverse range of baryonic physics models by multiplying the dark-matter only CrowCanyon matter power at \(z=0\)  by the ratio \(P_{\rm {baryon}}/P_{\rm {DMO}}\)  for 13 different models from [1]}.
Figure REF  shows fits to dark-matter-only power spectra and power spectra including the effects of baryons for several models with \(A_{0}\)  and \(R\)  fixed to their DMO values.
Clearly the \(R_{n,nh}\)  parameters are flexible enough to accurately account for feedback, and marginalizing over them should remove biases in cosmological parameter constraints.
We note that while all models are fit to the \(1\%\)  level down to \(k=1\ h/\rm {Mpc}\) , the simulations with the largest deviations in the \(R_{nh}\)  parameters from the DMO case are Illustris and Horizon-AGN.
This is the case for Illustris since it has been shown to have an unrealistically strong feedback model (compared to other hydrodynamic simulations) in terms of its effects on the power spectrum for \(k\le 1\)  due to low baryon fraction (compared with the observed value in galaxy groups) in high mass halos [2]}.
For Horizon-AGN, there is a large-scale \(1\%\)  excess of the \(P_{\rm {baryon}}/P_{\rm {DMO}}\)  ratio above unity that causes a relatively large change in the HZPT \(R_{nh}\)  parameters (the source of this deviation is discussed in Appendix A of [3]}, and may not be physical).
<FIGURE>