******************************************************;
title 'example of epidemiology 2003 APR';
**********************************;
data adata;
infile '/?????/???/chd.txt';
  input (chd cat age ecg) (3.);
******* atarashii hensu no keisan rei *********************************;
*tal=0;
*tal=ia1+ia2+ia3+ia4+ia5+ia6+ia7+ia8+ia9+ia10+ia11+ia12+ia13;
****************************************;
*run;
*********************;
proc logistic data = adata simple descending;
 model chd = cat age ecg
     /details plrl;
run;
****************************************;
data adata;
infile '/?????/???/6mp.txt';
  input (id sex age trt wek obs) (3.);
proc lifetest plots=(s,ls,lls);
   time wek*obs(1);
    strata trt;
*   test sex age;
run;
proc lifetest plots=(s,ls,lls) method=act;
   time wek*obs(1);
    strata age(25 45);
*   test trt sex;
run;
***************************************;
data cox;
infile '/?????/???/6mp.txt';
  input (id sex age trt wek obs) (3.);
if age<25 then cage=0;
if age>=25 then cage=1;
if age>=45 then cage=2;
run;
proc phreg simple;
   model wek*obs(1)=sex cage trt
                    /selection = stepwise slentry=.25 slstay=.15 details;
run;
***************************************;