proc phreg data=auswertung;
class patid event sma;
*model zeit*event(0,2)=sma_c/rl;
model zeit*event(0,2)=sma/rl;
array zz{*} zeit_01 zeit_02 zeit_11 zeit_12 zeit_21 zeit_22 zeit_31 zeit_32;
array tt{*} therapie_01 therapie_02 therapie_11 therapie_12 therapie_21 therapie_22 therapie_31 therapie_32;
sma_c=0;
do i=1 to (dim(zz)-1);
	if ((zeit-1)<zz[i+1] and zz[i]<zeit) then sma_c=sma_c+tt[i+1]*(min(zeit,zz[i])-max(zeit-1,zz[i+1]));
end;
if sma_c=0 then sma=0;else if 0<sma_c<0.5 then sma=1;else if 0.5<=sma_c<1 then sma=2;else if sma_c=1 then sma=3;
hazardratio sma /diff=pairwise;
id patid;
*format sma_c th_lag4k.;
title'fracture risk';
run;