Multi-scale control variate methods for uncertainty quantification in kinetic equations