libname a xport '<jtue.sarkooma>analyysi020712.tpt';
options ls=78 formdlim='-';
proc copy in=a out=work memtype=data;
run;

data analyysi;
   set analyysi;
   if stratase ne 9 then delete;

if vastattu = . then vastattu = 0;

if whoteq ne . then teqk1 = 0;
if whoteq ne . then teqk2 = 0;
if whoteq ne . then teqk3 = 0;
if whoteq ne . then teqk4 = 0;
if whoteq ne . then teqk5 = 0;

if                 whoteq lt 16.682439482 then teqk1  = 1;
if 16.682439482 le whoteq lt 24.611149088 then teqk2  = 1;
if 24.611149088 le whoteq lt 32.55454537  then teqk3  = 1;
if 32.55454537  le whoteq lt 48.255828587 then teqk4  = 1;
if 48.255828587 le whoteq le 145.45739916 then teqk5  = 1;

if teqk1=1 then teqk=1;
if teqk2=1 then teqk=2;
if teqk3=1 then teqk=3;
if teqk4=1 then teqk=4;
if teqk5=1 then teqk=5;

if asukkait ne . then asukkai1 = 0;
if asukkait ne . then asukkai2 = 0;
if asukkait ne . then asukkai3 = 0;
if asukkait ne . then asukkai4 = 0;

if           asukkait < 13750 then asukkai1 = 1;
if 13750  le asukkait < 56664 then asukkai2 = 1;
if 56664  le asukkait < 191247 then asukkai3 = 1;
if 191247 le asukkait le 525031 then asukkai4 = 1;

bmi = paino / pituus / pituus * 10000;

mkala = 1.180840432 +
      PETOKALA * 0.285437176 +
      MUIKKUA  * 0.120796063 +
      SIS_VESI * 0.081192631 +
      KIRJOLOH * 0.432524410 +
      SILAKKAA * 0.201481292 +
      IT_MEREN * 0.186820358 +
      MUUTA_IT *-0.096115539 +
      PAKASTEK * 0.105839909 +
      KALAS_IL *-0.059836406 +
      VALTAMER * 0.095303307 +
      _YRI_ISI * 0.033074683;

mkaladio = 1.180840432 +
      SILAKKAA * 0.201481292 +
      IT_MEREN * 0.186820358;

mkalamuu = 1.180840432 +
      PETOKALA * 0.285437176 +
      MUIKKUA  * 0.120796063 +
      SIS_VESI * 0.081192631 +
      KIRJOLOH * 0.432524410 +
      MUUTA_IT *-0.096115539 +
      PAKASTEK * 0.105839909 +
      KALAS_IL *-0.059836406 +
      VALTAMER * 0.095303307 +
      _YRI_ISI * 0.033074683;

if silakkaa = 0 then jsilakka = 0;
if silakkaa = 1 then jsilakka = 0.5;
if silakkaa = 2 then jsilakka = 1.5;
if silakkaa = 3 then jsilakka = 4;
if silakkaa = 4 then jsilakka = 8;
if silakkaa = 5 then jsilakka = 20;

if it_meren = 0 then jit_mere = 0;
if it_meren = 1 then jit_mere = 0.5;
if it_meren = 2 then jit_mere = 1.5;
if it_meren = 3 then jit_mere = 4;
if it_meren = 4 then jit_mere = 8;
if it_meren = 5 then jit_mere = 20;

if petokala = 0 then jpetokal = 0;
if petokala = 1 then jpetokal = 0.5;
if petokala = 2 then jpetokal = 1.5;
if petokala = 3 then jpetokal = 4;
if petokala = 4 then jpetokal = 8;
if petokala = 5 then jpetokal = 20;

if kirjoloh = 0 then jkirjolo = 0;
if kirjoloh = 1 then jkirjolo = 0.5;
if kirjoloh = 2 then jkirjolo = 1.5;
if kirjoloh = 3 then jkirjolo = 4;
if kirjoloh = 4 then jkirjolo = 8;
if kirjoloh = 5 then jkirjolo = 20;

jkaladio = jsilakka + jit_mere;
jkalamuu = jpetokal + jkirjolo;

if jkaladio ne . then jkaladi1 = 0;
if jkaladio ne . then jkaladi2 = 0;
if jkaladio ne . then jkaladi3 = 0;
if jkaladio ne . then jkaladi4 = 0;

if        jkaladio < 0.5 then jkaladi1 = 1;
if 0.5 le jkaladio < 1   then jkaladi2 = 1;
if 1   le jkaladio < 2   then jkaladi3 = 1;
if 2   le jkaladio le 12 then jkaladi4 = 1;

if jkalamuu ne . then jkalamu1 = 0;
if jkalamuu ne . then jkalamu2 = 0;
if jkalamuu ne . then jkalamu3 = 0;
if jkalamuu ne . then jkalamu4 = 0;

if        jkalamuu < 1.5 then jkalamu1 = 1;
if 1.5 le jkalamuu < 2   then jkalamu2 = 1;
if 2   le jkalamuu < 3.5 then jkalamu3 = 1;
if 3.5 le jkalamuu le 24 then jkalamu4 = 1;

if dgluokka = 4 or dgluokka = 5 or dgluokka > 6 then dgluokka = 13;

if kouluv ne . then kouluv1 = 0;
if kouluv ne . then kouluv2 = 0;
if kouluv ne . then kouluv3 = 0;
if kouluv ne . then kouluv4 = 0;

if         kouluv <  9    then kouluv1 = 1;
if 9    le kouluv <  11.5 then kouluv2 = 1;
if 11.5 le kouluv <  15   then kouluv3 = 1;
if 15   le kouluv le 22   then kouluv4 = 1;

if kouluv ne . then koulu1 = 0;
if kouluv ne . then koulu2 = 0;
if kouluv ne . then koulu3 = 0;

if       kouluv   <10 then koulu1 = 1;
if 10 le kouluv le 12 then koulu2 = 1;
if 12 <  kouluv       then koulu3 = 1;

if pituus ne . then pituus1 = 0;
if pituus ne . then pituus2 = 0;
if pituus ne . then pituus3 = 0;
if pituus ne . then pituus4 = 0;

if        pituus < 163     then pituus1 = 1;
if 163 le pituus < 170     then pituus2 = 1;
if 170 le pituus < 178     then pituus3 = 1;
if 178 le pituus le 197    then pituus4 = 1;

if paino ne . then paino1 = 0;
if paino ne . then paino2 = 0;
if paino ne . then paino3 = 0;
if paino ne . then paino4 = 0;

if        paino < 65 then paino1 = 1;
if 65 le  paino < 74 then paino2 = 1;
if 74 le  paino < 84 then paino3 = 1;
if 84 le  paino le 125 then paino4 = 1;

if bmi ne . then bmi1 = 0;
if bmi ne . then bmi2 = 0;
if bmi ne . then bmi3 = 0;
if bmi ne . then bmi4 = 0;

if                 bmi < 23.163373273 then bmi1 = 1;
if 23.163373273 le bmi < 24.977043159 then bmi2 = 1;
if 24.977043159 le bmi < 28.256989887 then bmi3 = 1;
if 28.256989887 le bmi le 39.442073936 then bmi4 = 1;

run;

data analyysi;
set analyysi;

if kalaa ne . then kalaa1 = 0;
if kalaa ne . then kalaa2 = 0;
if kalaa ne . then kalaa3 = 0;
if kalaa ne . then kalaa4 = 0;

if kalaa = 1 then kalaa1=1;
if kalaa = 2 then kalaa2=1;
if kalaa = 3 then kalaa3=1;
if 3 lt kalaa le 6 then kalaa4 = 1;

if mkala ne . then mkala1 = 0;
if mkala ne . then mkala2 = 0;
if mkala ne . then mkala3 = 0;
if mkala ne . then mkala4 = 0;

if                mkala < 2.524308051 then mkala1 = 1;
if 2.524308051 le mkala < 3.061359395 then mkala2 = 1;
if 3.061359395 le mkala < 3.463804967 then mkala3 = 1;
if 3.463804967 le mkala le 4.723833506 then mkala4 = 1;

if mkaladio ne . then mkaladi1 = 0;
if mkaladio ne . then mkaladi2 = 0;
if mkaladio ne . then mkaladi3 = 0;
if mkaladio ne . then mkaladi4 = 0;

if                mkaladio <  1.382321724 then mkaladi1 = 1;
if 1.382321724 le mkaladio <  1.569142082 then mkaladi2 = 1;
if 1.569142082 le mkaladio <  1.770623374 then mkaladi3 = 1;
if 1.770623374 le mkaladio le 2.547226674 then mkaladi4 = 1;

if mkalamuu ne . then mkalamu1 = 0;
if mkalamuu ne . then mkalamu2 = 0;
if mkalamuu ne . then mkalamu3 = 0;
if mkalamuu ne . then mkalamu4 = 0;

if                mkalamuu <  2.231277899 then mkalamu1 = 1;
if 2.231277899 le mkalamuu <  2.633878211 then mkalamu2 = 1;
if 2.633878211 le mkalamuu <  3.01143575 then mkalamu3 = 1;
if 3.01143575  le mkalamuu le 4.510720146 then mkalamu4 = 1;

if tstatus ne . then tstatus1 = 0;
if tstatus ne . then tstatus2 = 0;
if tstatus ne . then tstatus3 = 0;

if tstatus = 1 then tstatus1 = 1;
if tstatus = 2 then tstatus2 = 1;
if tstatus = 3 then tstatus3 = 1;

if tstatus ne . then current = 0;
if tstatus ne . then former  = 0;
if tstatus ne . then never   = 0;
if tstatus = 1 then current = 1;
if tstatus = 2 then former  = 1;
if tstatus = 3 then never   = 1;

if packyear ne . then packyea1 = 0;
if packyear ne . then packyea2 = 0;
if packyear ne . then packyea3 = 0;
if packyear ne . then packyea4 = 0;

if        packyear = 0      then packyea1 = 1;
if 0    < packyear < 0.1    then packyea2 = 1;
if 0.1 le packyear < 10     then packyea3 = 1;
if 10  le packyear le 112.5 then packyea4 = 1;

if alko ne . then alko1 = 0;
if alko ne . then alko2 = 0;
if alko ne . then alko3 = 0;
if alko ne . then alko4 = 0;

if      alko < 3 then alko1 = 1;
if 3 le alko < 4 then alko2 = 1;
if 4 le alko < 6 then alko3 = 1;
if 6 le alko le 10 then alko4 = 1;

if alkous ne . then alkous1 = 0;
if alkous ne . then alkous2 = 0;
if alkous ne . then alkous3 = 0;
if alkous ne . then alkous4 = 0;

if        alkous < 0.1 then alkous1 = 1;
if 0.1 le alkous < 0.5 then alkous2 = 1;
if 0.5 le alkous <   1 then alkous3 = 1;
if   1 le alkous le  7 then alkous4 = 1;

if alkokulu ne . then alkokul1 = 0;
if alkokulu ne . then alkokul2 = 0;
if alkokulu ne . then alkokul3 = 0;
if alkokulu ne . then alkokul4 = 0;

if        alkokulu < 0.1 then alkokul1 = 1;
if 0.1 le alkokulu <   1 then alkokul2 = 1;
if   1 le alkokulu <   3 then alkokul3 = 1;
if   3 le alkokulu <  84 then alkokul4 = 1;

if LIUOTTI + MAALIT + FORMALDE + PESUAINE + METALLIT + MUU1 + MUU2 = 7 then kemikmuu = 0;
if LIUOTTI  = 2 then kemikmuu = 1;
if MAALIT   = 2 then kemikmuu = 1;
if FORMALDE = 2 then kemikmuu = 1;
if PESUAINE = 2 then kemikmuu = 1;
if METALLIT = 2 then kemikmuu = 1;
if MUU1     = 2 then kemikmuu = 1;
if MUU2     = 2 then kemikmuu = 1;

if leikkaus ne . then leikkau1 = 0;
if leikkaus ne . then leikkau2 = 0;
if leikkaus ne . then leikkau3 = 0;
if leikkaus ne . then leikkau4 = 0;

if               leikkaus < '08nov97'd then leikkau1 = 1;
if '08nov97'd le leikkaus < '20mar98'd then leikkau2 = 1;
if '20mar98'd le leikkaus < '17oct98'd then leikkau3 = 1;
if '17oct98'd le leikkaus le '13dec99'd then leikkau4 = 1;

survival = 2 - tapver;

run;

data analyysi;
set analyysi;

if  kouluv1=1 or  kouluv2=2 then  kouluvb=0; else if  kouluv3=1 or  kouluv4=1 then  kouluvb=1;
if asukkai1=1 or asukkai2=1 then asukkaib=0; else if asukkai3=1 or asukkai4=1 then asukkaib=1;
if     bmi1=1 or     bmi2=1 then     bmib=0; else if     bmi3=1 or     bmi4=1 then     bmib=1;
if   kalaa1=1 or   kalaa2=1 then   kalaab=0; else if   kalaa3=1 or   kalaa4=1 then   kalaab=1;
if mkaladi1=1 or mkaladi2=1 then mkaladib=0; else if mkaladi3=1 or mkaladi4=1 then mkaladib=1;
if mkalamu1=1 or mkalamu2=1 then mkalamub=0; else if mkalamu3=1 or mkalamu4=1 then mkalamub=1;
if tstatus =3 then tstatusb=0; else if tstatus = 1 or tstatus=2 then tstatusb=1;
if packyea1=1 or packyea2=1 then packyeab=0; else if packyea3=1 or packyea4=1 then packyeab=1;
if  alkous1=1 or  alkous2=1 then  alkousb=0; else if  alkous3=1 or  alkous4=1 then  alkousb=1;
if alkokul1=1 or alkokul2=1 then alkokulb=0; else if alkokul3=1 or alkokul4=1 then alkokulb=1;
if leikkau1=1 or leikkau2=1 then leikkaub=0; else if leikkau3=1 or leikkau4=1 then leikkaub=1;
if jkaladi1=1 or jkaladi2=1 then jkaladib=0; else if jkaladi3=1 or jkaladi4=1 then jkaladib=1;
if jkalamu1=1 or jkalamu2=1 then jkalamub=0; else if jkalamu3=1 or jkalamu4=1 then jkalamub=1;

if               leikkaus < '01jan97'd then leik=1;
if '01jan97'd le leikkaus < '01apr97'd then leik=2;
if '01apr97'd le leikkaus < '01jul97'd then leik=3;
if '01jul97'd le leikkaus < '01oct97'd then leik=4;
if '01oct97'd le leikkaus < '01jan98'd then leik=5;
if '01jan98'd le leikkaus < '01apr98'd then leik=6;
if '01apr98'd le leikkaus < '01jul98'd then leik=7;
if '01jul98'd le leikkaus < '01oct98'd then leik=8;
if '01oct98'd le leikkaus < '01jan99'd then leik=9;
if '01jan99'd le leikkaus < '01apr99'd then leik=10;
if '01apr99'd le leikkaus < '01jul99'd then leik=11;
if '01jul99'd le leikkaus < '01oct99'd then leik=12;
if '01oct99'd le leikkaus < '01jan00'd then leik=13;
if '01jan00'd le leikkaus              then leik=14;


run;

data vastattu;
   set analyysi;
   if vastattu = 0 then delete;
run;

data ruuhka;
   set analyysi;
   if leikkaus <'01jul97'd then delete;
   if leikkaus >'31dec98'd then delete;
run;

data ruuh_a;
   set analyysi;
   if leikkaus ge'01jul97'd and tapver = 1 then delete;
run;

data ruuh_b;
   set analyysi;
   if leikkaus <'01jul97'd and tapver = 1 then delete;
   if leikkaus >'31dec98'd and tapver = 1 then delete;
run;

data ruuh_c;
   set analyysi;
   if leikkaus le '31dec98'd and tapver = 1 then delete;
run;

data ruuh_ab;
   set analyysi;
   if leikkaus >'31dec98'd and tapver = 1 then delete;
run;

data ruuh_bc;
   set analyysi;
   if leikkaus <'01jul97'd and tapver = 1 then delete;
run;

data sadepois;
   set analyysi;
   if sadeaik=1 then delete;
run;

data ruuh_b2;
   set analyysi;
   if leikkaus <'01jul97'd and tapver = 1 then delete;
   if leikkaus >'28feb99'd and tapver = 1 then delete;
run;

data analyysi;
   set analyysi;
run;

title 'STRATASET 9';

proc phreg data=work.vastattu;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp kalaa2 kalaa3 kalaa4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp mkala2 mkala3 mkala4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp mkaladi2 mkaladi3 mkaladi4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp mkalamu2 mkalamu3 mkalamu4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp kouluv2 kouluv3 kouluv4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp pituus2 pituus3 pituus4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp paino2 paino3 paino4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp bmi2 bmi3 bmi4 
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp tstatus2 tstatus3
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp packyea2 packyea3 packyea4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp alko2 alko3 alko4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp alkous2 alkous3 alkous4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp alkokul2 alkokul3 alkokul4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)=
   sp kalaa2 kalaa3 kalaa4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)=
   sp mkala2 mkala3 mkala4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)=
   sp mkaladi2 mkaladi3 mkaladi4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)=
   sp mkalamu2 mkalamu3 mkalamu4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)=
   sp kouluv2 kouluv3 kouluv4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)=
   sp pituus2 pituus3 pituus4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)=
   sp paino2 paino3 paino4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)=
   sp bmi2 bmi3 bmi4 
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)=
   sp tstatus2 tstatus3
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)=
   sp packyea2 packyea3 packyea4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)=
   sp alko2 alko3 alko4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)=
   sp alkous2 alkous3 alkous4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)=
   sp alkokul2 alkokul3 alkokul4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.analyysi;
   model survival*tapver(0)= whoteq
   sp
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.analyysi;
   model survival*tapver(0)= whoteq

   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp kouluv2 kouluv3 kouluv4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.analyysi;
   model survival*tapver(0)= whoteq
   sp asukkai2 asukkai3 asukkai4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp bmi2 bmi3 bmi4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp kalaa2 kalaa3 kalaa4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp mkaladi2 mkaladi3 mkaladi4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp mkalamu2 mkalamu3 mkalamu4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp jkaladi2 jkaladi3 jkaladi4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp jkalamu2 jkalamu3 jkalamu4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.analyysi;
   model survival*tapver(0)= whoteq
   sp sadeaik
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp former current
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp packyea2 packyea3 packyea4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp alkous2 alkous3 alkous4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp alkokul2 alkokul3 alkokul4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp hyonteis
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp kasvinsu
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp kyllaste
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp kemikmuu
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.analyysi;
   model survival*tapver(0)= whoteq
   sp leikkau2 leikkau3 leikkau4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp kouluvb
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.analyysi;
   model survival*tapver(0)= whoteq
   sp asukkaib
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp bmib
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp kalaab
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp mkaladib
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp mkalamub
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp jkaladib
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp jkalamub
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp tstatusb
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp packyeab
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp alkousb
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= whoteq
   sp alkokulb
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.analyysi;
   model survival*tapver(0)= whoteq
   sp leikkaub
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.ruuhka;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.ruuhka;
   model survival*tapver(0)= whoteq
   sp
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp kouluv2 kouluv3 kouluv4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.ruuh_a;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.ruuh_b;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.ruuh_c;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.ruuh_ab;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.ruuh_bc;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.ruuh_a;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
 
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.ruuh_b;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5

   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.ruuh_c;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5

   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.ruuh_ab;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5

   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.ruuh_bc;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5

   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp koulu2 koulu3
   /ties=discrete rl;
   strata stratum;
run;

proc sort;
   by tapver leik;
run;

proc means;
   by tapver leik;
   var whoteq;
run;

data ruuh_a;
set ruuh_a;
run;

proc freq;
   tables teqk*tapver;
run;

data ruuh_b;
set ruuh_b;
run;

proc freq;
   tables teqk*tapver;
run;

data ruuh_c;
set ruuh_c;
run;

proc freq;
   tables teqk*tapver;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp jkaladi2 jkaladi3 jkaladi4
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.vastattu;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp jkalamu2 jkalamu3 jkalamu4
   /ties=discrete rl;
   strata stratum;
run;

data vastattu;
set vastattu;
run;

proc sort;
by tapver;
run;

proc means n mean p50 std data=vastattu;
   by tapver;
   var jkaladio jkalamuu;
run;

proc phreg data=work.analyysi;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp sadeaik
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.sadepois;
   model survival*tapver(0)= teqk2 teqk3 teqk4 teqk5
   sp
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.ruuh_b;
   model survival*tapver(0)= whoteq
   sp
   /ties=discrete rl;
   strata stratum;
run;

proc phreg data=work.ruuh_b2;
   model survival*tapver(0)= whoteq
   sp
   /ties=discrete rl;
   strata stratum;
run;