  We present constraints on canonical single-field inflation derived from WMAP
five year, ACBAR, QUAD, BICEP data combined with the halo power spectrum from
SDSS LRG7. Models with a non-scale-invariant spectrum and a red tilt n_s < 1
are now preferred over the Harrison-Zel'dovich model (n_s = 1, tensor-to-scalar
ratio r = 0) at high significance. Assuming no running of the spectral indices,
we derive constraints on the parameters (n_s, r) and compare our results with
the predictions of simple inflationary models. The marginalised credible
intervals read n_s = 0.962^{+0.028}_{-0.026} and r < 0.17 (at 95% confidence
level). Interestingly, the 68% c.l. contours favour mainly models with a convex
potential in the observable region, but the quadratic potential model remains
inside the 95% c.l. contours. We demonstrate that these results are robust to
changes in the datasets considered and in the theoretical assumptions made. We
then consider a non-vanishing running of the spectral indices by employing
different methods, non-parametric but approximate, or parametric but exact.
With our combination of CMB and LSS data, running models are preferred over
power-law models only by a Delta chi^2 ~ 5.8, allowing inflationary stages
producing a sizable negative running -0.063^{+0.061}_{-0.049} and larger
tensor-scalar ratio r < 0.33 at the 95% c.l. This requires large values of the
third derivative of the inflaton potential within the observable range. We
derive bounds on this derivative under the assumption that the inflaton
potential can be approximated as a third order polynomial within the observable
range.
