An extension of the classical linear Model Output Statistics (MOS) technique is proposed allowing for the post-processing of ensemble forecasts. In this new approach, the cost function on which the least square parameter estimation is based takes into account the presence of errors in both observations and model observables (referred to as Error-in-Variables MOS, EVMOS), unlike the classical linear MOS cost function whose implicit assumption is the absence of errors in the model observables. It allows for the maintenance of an appropriate variability for the corrected forecasts, even for long lead times and for providing a framework in which both deterministic and probabilistic forecasts can be corrected. The scheme is successfully tested for ensemble correction in the context of an idealized low-order chaotic system, the Lorenz atmospheric model, in the presence of model errors, and compared with a classical technique, known as the non-homogeneous Gaussian regression (NGR) method. The potential use of this approach is also briefly discussed.