Anglada-Escude and Tuomi question the statistical rigor of our analysis while ignoring the stellar activity aspects that we present. Although we agree that improvements in multiparametric radial velocity (RV) modeling are necessary for the detection of Earth-mass planets, the key physical points we raised were not challenged. We maintain that activity on Gliese 581 induces RV shifts that were interpreted as exoplanets.