Hi folks,

My data set were analysed using

data mauricio; input x n y; datalines; 1735 04 2 1936 10 3 2000 05 1 . . . . . . . . .

proc genmod data=mauricio; model y/n = x / dist=bin type1 type3; output out=saida p=pre; run;

I want the sum below:

C = Sum [ (po/pre)**(y)*((1-po)/(1-pre))**(n-y) ]

in that 0.40 <= po <= 0.60.

and, I did:

data graf; set saida; do p=0.39 to 0.60 by 0.01; output; C=(pr/p)**(y)*((1-pr)/(1-p))**(n-y); end; run;

proc iml; reset print; use graf; read all var{x n y prop p pr C } ; list all; prp=pr'pr; quit;

but, I can't find the sum of C.

Any suggestions would be greatly appreciated

Silvano