Analysis of stacking disorder in ice I using pair distribution functions