  We constrain deviations from general relativity (GR) including both redshift
and scale dependencies in the modified gravity (MG) parameters. In particular,
we employ the under-used binning approach and compare the results to functional
forms. We use available datasets such as Cosmic Microwave Background (CMB) from
Planck 2018, Baryonic Acoustic Oscillations (BAO) and Redshift Space
Distortions (BAO/RSD) from the BOSS DR12, the 6DF Galaxy Survey, the SDSS DR7
Main Galaxy Sample, the correlation of Lyman-$\alpha$ forest absorption and
quasars from SDSS-DR14, Supernova Type Ia (SNe) from the Pantheon compilation,
and DES Y1 data. Moreover, in order to maximize the constraining power from
available datasets, we analyze MG models where we alternatively set some of the
MG parameters to their GR values and vary the others. Using functional forms,
we find an up to 3.5-$\sigma$ tension with GR in $\Sigma$ (while $\mu$ is
fixed) when using Planck+SNe+BAO+RSD; this goes away when lensing data is
included, i.e. CMB lensing and DES (CMBL+DES). Using different binning methods,
we find that a tension with GR above 2-$\sigma$ in the (high-z, high-k) bin is
persistent even when including CMBL+DES to Planck+SNe+BAO+RSD. Also, we find
another tension above 2-$\sigma$ in the (low-z, high-k) bin, but that can be
reduced with the addition of lensing data. Furthermore, we perform a model
comparison using the Deviance Information Criterion statistical tool and find
that the MG model ($\mu=1$, $\Sigma$) is weakly favored by the data compared to
$\Lambda$CDM, except when DES data is included. Another noteworthy result is
that we find that the binning methods do not agree with the widely-used
functional parameterization where the MG parameters are proportional to
$\Omega_{\text{DE}}(a)$, and this is clearly apparent in the high-z and high-k
regime where this parameterization underestimates the deviations from GR.
