Date:         Fri, 12 Apr 2002 08:36:07 -0300
Reply-To:     Silvano Cesar da Costa <silvano@UEL.BR>
Sender:       "SAS(r) Discussion" <SAS-L@LISTSERV.UGA.EDU>
From:         Silvano Cesar da Costa <silvano@UEL.BR>
Subject:      Genmod - IML
Content-Type: text/plain

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