Not yet reviewed by Pith; the record is open.

Nongaussian statistics are a powerful discriminant between inflationary models, particularly those with noncanonical kinetic terms. Focusing on theories where the Lagrangian is an arbitrary Lorentz-invariant function of a scalar field and its first derivatives, we review and extend the calculation of the observable three-point function. We compute the "next-order" slow-roll corrections to the bispectrum in closed form, and obtain quantitative estimates of their magnitude in DBI and power-law k-inflation. In the DBI case our results enable us to estimate corrections from the shape of the potential and the warp factor: these can be of order several tens of percent. We track the possible sources of large logarithms which can spoil ordinary perturbation theory, and use them to obtain a general formula for the scale dependence of the bispectrum. Our result satisfies the next-order version of Maldacena's consistency condition and an equivalent consistency condition for the scale dependence. We identify a new bispectrum shape available at next-order, which is similar to a shape encountered in Galileon models. If fNL is sufficiently large this shape may be independently detectable.