Based on the very general assumptions on evolution, we have modeled the statistics of SNRs evolving at various initial and environmental conditions. The method is applied to M33, for which the value of the filling factor of the hot phase of the ISM is estimated to be ~ 10%.