data a; newp=.00005; do gen=1 to 10; seed=47639087; expon=2*ranuni(seed) + 4; n=round(10**expon,1); p=.newp; q=1-p; NumRR = round(n*p*p,1); NumSR = round(2*n*p*q,1); NumSS = round(n*q*q,1); newp=(NumRR + .5*NumSR) / (Numrr + numsr + numss); diffp=p-newp; newq=(NumSS + .5*NumSR)/ (Numrr + numsr + numss); diffq=q-newq; output; end;
... View more
