normal probability plot

[shahd]

Spoiler

Hi, I am trying to plot the normal probability plot of deviance could anyone help with that please 
data rate;
input x1 x2 x3 x4 rate;
datalines;
0 0 0 0 12.60
1 0 0 0 18.70
0 1 0 0 3.28
1 1 0 0 3.44
0 0 1 0 4.98
1 0 1 0 5.70
0 1 1 0 9.97
1 1 1 0 9.07
0 0 0 1 2.07
1 0 0 1 2.44
0 1 0 1 4.09
1 1 0 1 4.53
0 0 1 1 7.77
1 0 1 1 9.43
0 1 1 1 11.75
1 1 1 1 16.30

proc genmod;
model rate = x1 x2 x3 x4/ dist =gamma link = log type3 ;
run;

ods graphics on;

proc genmod data=rate plots(unpack)=all;

model rate=x1 x2 x3 x4/dist=gamma link=log type1 typ3;

run;

ods graphics on;

proc genmod data=rate plots=resdev;

model rate=x1 x2 x3 x4/dist=gamma link=log type1 typ3;

run;
output out=resdev=deviance;
proc univariate data=mydata;
var deviance;
probplot;
run;

[PeterClemmensen]

What do you mean by "normal probability plot of deviance"?

[Ksharp]

@Rick_SAS have written many blog about it. Check PDF() function.

[Ksharp]

data have;
do q=-3 to 3 by 0.1;
 d=pdf('normal',q);output;
end;
run;
proc sgplot data=have;
series x=q y=d;
run;