data jdp;
input lifetime @@;
life2k=(lifetime>1999);
cards;
1993 1993 1984 1997 1995 1991 1984 1996 1989 2001
;
run;
proc summary data=jdp;
 var life2k;
 output out=lab4c_out mean=p_hat;
run;
data lab4c_out;
 set lab4c_out;
 label p_ci_low='Lower Limit for p' p_ci_up='Upper Limit for p';
 n_tilde=10+1.96**2;
 p_tilde=(10*p_hat+.5*1.96**2)/n_tilde;
 p_ci_low1=p_tilde-1.96*sqrt(p_tilde*(1-p_tilde)/n_tilde);
 p_ci_low=max(p_ci_low1,0);
 p_ci_up1=p_tilde+1.96*sqrt(p_tilde*(1-p_tilde)/n_tilde);
 p_ci_up=min(p_ci_up1,1);
run;
proc print data=lab4c_out noobs label;
 var p_ci_low p_ci_up;
run;