This function carries out a Wilson plot calculation based on either .ins + .hkl or .cif + .fcf. Various E-statistics are calculated.