  Anglada-Escud\'e 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.
