/* Historical Control Program for Binomial Responses Input Data rc=Number of Successes (Responses) in the Historical Control Sample nc=Number of Patients in the Historical Control Sample pexp=Hypothesized Response Rate for the Experimental Therapy (>rc/nc) alpha=Type I error Power=Power Output: ne=Sample size for experimental group re=Minimum number of responses to achieve P-val1 then pexp=.01*pexp; if alpha>1 then alpha=.01*alpha;if power>1 then power=.01*power; za=-probit(alpha); aa: ne=ne+1; re=0; bb:re=re+1;pe=re/ne; pstar=(re+rc)/(ne+nc); * if re999.6 then put @10 rc @16 nc @25 pexp @36 alpha @44 power @52 'No soln or N>999'; run;