Motivated by the work of Li and Meneveau [Phys. Rev. Lett. 95, 164502 (2005)], we propose and solve a model for the Lagrangian evolution of both longitudinal and transverse velocity and temperature increments for Boussinesq convection. From this model, the short-time evolution of an initially imposed Gaussian joint probability density function (PDF) of both velocity and temperature increments is computed analytically and the trend to non-Gaussian statistics shown in a quantitative way. Predictions for moments of the joint PDF are obtained and their behavior analyzed with respect to known experimental and numerical results. The obtained results do not depend on the model free parameters, a fact in favor of their robustness.