proc sort data=xstd;
      by xgroup;
   run;

   ods graphics on;
   proc freq data=xstd order=data;
      tables Z / chisq testp=(0.0258 0.09685 0.13578 0.14344 0.20106 0.27451 0.12268)
                    plots(only)=deviationplot(type=dot);
      weight percent;
      by xgroup;
      title 'Chi square of theoretical distribution';
   run;
   ods graphics off;