The stacking disorder in ice I samples is analysed in a quantitative fashion using pair distribution functions derived from neutron total scattering. This approach is found to have several advantages compared with the traditional analysis of powder diffraction data.