/* Elinaikamallien SAS-demoja */ /* Kopioi data omaan hakemistoon. Alla C:\MyTemp\ea_sas */ /* Annetaan oletushakemiston polku, jossa data sijaitsee. */ libname d 'C:\MyTemp\ea_sas'; /* Luetaan tekstitiedosto enstul.dat SASiin, nimetään muuttujat ja annetaan aineistolle nimeksi lapset. Aineisto tallentuu omaan hakemistoon. */ data d.lapset; infile 'C:\MyTemp\ea_sas\enstul.dat'; input num time delta sukup tup lapsia; run; /* Estimoidaan välttöfkt KM-menetelmällä; piirretään estimaatti sekä (t,-log S) ja (log t, log (-log S))-käyrät */ proc lifetest data=d.lapset method=km plots=(s,ls,lls); time time*delta(0); run; /* Lasketaan luottamusvälit KM-estimaatille ja plotataan ne. Vaihtoehtoisia tapoja: conftype=log tai conftype=loglog */ proc lifetest data=d.lapset; time time*delta(0); survival out=d.Out1 conftype=linear; run; goptions reset=all; symbol1 c=black i=stepjll; symbol2 c=red i=stepjl line=21; symbol3 c=red i=stepjl line=21; proc gplot data = d.Out1; plot (survival sdf_lcl sdf_ucl)*time / overlay; run; /* Ryhmitellään aineisto ja plotataan edellisten lisäksi vaara- ja tiheysfunktiot */ proc lifetest data=d.lapset method=lt intervals=(0 to 25 by 1) plots=(s,h,p,ls,lls); time time*delta(0); survival out=d.Out1 conftype=linear; run; /* Esim. vaaran plottaus porrasfunktiona */ proc gplot data = d.Out1; plot (hazard)*time; run; /* Välttökäyräestimaatit LAPSIA-mjan eri luokissa, ja erojen testaus */ goptions reset=all; proc lifetest data=d.lapset plots=(s); time time*delta(0); strata lapsia; run; /* Välttökäyräestimaatit sukupuolittain */ proc lifetest data=d.lapset plots=(s); time time*delta(0); strata sukup; run; /* Erojen testus Peto-Peton -testillä */ proc lifetest data=d.lapset plots=(s); time time*delta(0); strata sukup / test=(peto); run; /* Testataan tupakoinnin vaikutusta, eliminoimalla sukupuolen mahdollinen vaikutus */ proc lifetest data=d.lapset plots=(s); time time*delta(0); strata sukup / group=tup; run;