Kinetic equations play a major rule in modeling large systems of interacting
particles. Uncertainties may be due to various reasons, like lack of knowledge
on the microscopic interaction details or incomplete informations at the
boundaries. These uncertainties, however, contribute to the curse of
dimensionality and the development of efficient numerical methods is a
challenge. In this paper we consider the construction of novel multi-scale
methods for such problems which, thanks to a control variate approach, are
capable to reduce the variance of standard Monte Carlo techniques