observations and avoid integrating over the redshift $z$. The corresponding $\chi^2$ for the OHD is defined as $\chi^2_{OHD}=\sum^{36}_{i=1}[\frac{H_0E(z_i)-H_{obs}(z_i)}{\sigma_i}]^2$, where $H_{obs}(z_i)$ is the observed value of the Hubble parameter at a given redshift $z_i$. Another important and powerful constraint comes from the CMB anisotropy observations, and we use the Planck-2015 data sets to constrain the DGP braneworld scenario. For simplicity, we take the shift parameter $\mathcal{R}$ instead of the full data from the CMB anisotropy observations, because taking the entire CMB data to perform a global fitting will expend a great deal of power and computation time. The shift parameter $\mathcal{R}$ is defined as $\mathcal{R}=\sqrt{\Omega_{m0}}\int^{z_C}_0\frac{dz'}{E(z')}$, where $z_C$ is the redshift of recombination and $\Omega_{m0}$ is the present-day value of the matter density ratio parameter. As shown in a recent paper [@3], the shift parameter for the Planck-2015 data sets is $\mathcal{R}=1.7482\pm0.0048$ and $z_C=1089.90$. The corresponding $\chi^2$ for the CMB anisotropy observations can be expressed as $\chi^2_{CMB}(\theta)=[\frac{\mathcal{R}_{obs}-\mathcal{R}(\theta)}{\sigma_{\mathcal{R}}}]^2$, where $\mathcal{R}_{obs}$ and $\sigma_{\mathcal{R}}$ denote the values of the shift parameter and the corresponding $1\sigma$ statistical error, respectively. In addition, as in our previous works [@57; @58; @59], we still take a single data point from the single gravitational-wave (GW) event GW150914 [@61]. Note that, 