---------------------------------------------------------------------------------------------------------- log: C:\WDDtext\log\7.7.Framingham.log log type: text opened on: 29 Aug 2001, 17:55:59 . * 7.7.Framingham.log . * . * Proportional hazards regression analysis of the effect of gender and . * baseline diastolic blood pressure (DBP) on coronary heart disease (CHD) . * adjusted for age, body mass index (BMI) and serum cholesterol (SCL). . * . set memory 2000 (2000k) . use C:\WDDtext\2.20.Framingham.dta, clear . set textsize 120 . * . * Univariate analysis of the effect of DBP on CHD . * . graph dbp, bin(50) freq xlabel(40,60 to 140) xtick(50,70 to 150) ylabel(0,100 to 600) ytick(50,100 to > 550) gap(4) . more . generate dbpgr = recode(dbp,60,70,80,90,100,110,111) . tabulate dbpgr chdfate | Coronary Heart | Disease dbpgr | Censored CHD | Total -----------+----------------------+---------- 60 | 132 18 | 150 70 | 592 182 | 774 80 | 1048 419 | 1467 90 | 863 404 | 1267 100 | 417 284 | 701 110 | 125 110 | 235 111 | 49 56 | 105 -----------+----------------------+---------- Total | 3226 1473 | 4699 . label define dbp 60 "DBP<= 60" 70 "60 " 111 "110< DBP" . label values dbpgr dbp . generate time= followup/365.25 . label variable time "Follow-up in Years" . stset time, failure(chdfate) failure event: chdfate ~= 0 & chdfate ~= . obs. time interval: (0, time] exit on or before: failure ------------------------------------------------------------------------------ 4699 total obs. 0 exclusions ------------------------------------------------------------------------------ 4699 obs. remaining, representing 1473 failures in single record/single failure data 103710.1 total analysis time at risk, at risk from t = 0 earliest observed entry t = 0 last observed exit t = 32 . sts graph , by(dbpgr) xlabel(0,5 to 30) ylabel(0,.2 to 1) ytick(.1,.3 to .9) l1title(Proportion With > out CHD) gap(6) noborder failure _d: chdfate analysis time _t: time . more . sts test dbpgr failure _d: chdfate analysis time _t: time Log-rank test for equality of survivor functions | Events Events dbpgr | observed expected ---------+------------------------- DBP<= 60 | 18 53.63 60chi2 = 0.0000 . sts test dbpgr if dbpgr == 60 | dbpgr == 70 failure _d: chdfate analysis time _t: time Log-rank test for equality of survivor functions | Events Events dbpgr | observed expected ---------+------------------------- DBP<= 60 | 18 32.58 60chi2 = 0.0052 . sts test dbpgr if dbpgr == 70 | dbpgr == 80 failure _d: chdfate analysis time _t: time Log-rank test for equality of survivor functions | Events Events dbpgr | observed expected ---------+------------------------- 60chi2 = 0.0028 . sts test dbpgr if dbpgr == 80 | dbpgr == 90 failure _d: chdfate analysis time _t: time Log-rank test for equality of survivor functions | Events Events dbpgr | observed expected ---------+------------------------- 70chi2 = 0.0090 . sts test dbpgr if dbpgr == 90 | dbpgr == 100 failure _d: chdfate analysis time _t: time Log-rank test for equality of survivor functions | Events Events dbpgr | observed expected ---------+------------------------- 80chi2 = 0.0000 . sts test dbpgr if dbpgr == 100 | dbpgr == 110 failure _d: chdfate analysis time _t: time Log-rank test for equality of survivor functions | Events Events dbpgr | observed expected ---------+------------------------- 90DBP100 | 284 306.92 100BP110 | 110 87.08 ---------+------------------------- Total | 394 394.00 chi2(1) = 7.77 Pr>chi2 = 0.0053 . sts test dbpgr if dbpgr == 110 | dbpgr == 111 failure _d: chdfate analysis time _t: time Log-rank test for equality of survivor functions | Events Events dbpgr | observed expected ---------+------------------------- 100BP110 | 110 122.95 110< DBP | 56 43.05 ---------+------------------------- Total | 166 166.00 chi2(1) = 5.29 Pr>chi2 = 0.0215 . xi: stcox i.dbpgr i.dbpgr _Idbpgr_1-7 (_Idbpgr_1 for dbpgr==60 omitted) failure _d: chdfate analysis time _t: time Iteration 0: log likelihood = -11834.856 Iteration 1: log likelihood = -11761.466 Iteration 2: log likelihood = -11726.108 Iteration 3: log likelihood = -11723.967 Iteration 4: log likelihood = -11723.942 Iteration 5: log likelihood = -11723.942 Refining estimates: Iteration 0: log likelihood = -11723.942 Cox regression -- Breslow method for ties No. of subjects = 4699 Number of obs = 4699 No. of failures = 1473 Time at risk = 103710.0917 LR chi2(6) = 221.83 Log likelihood = -11723.942 Prob > chi2 = 0.0000 ------------------------------------------------------------------------------ _t | _d | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _Idbpgr_2 | 1.968764 .486453 2.74 0.006 1.213037 3.195312 _Idbpgr_3 | 2.557839 .6157326 3.90 0.000 1.595764 4.099941 _Idbpgr_4 | 3.056073 .7362768 4.64 0.000 1.905856 4.900466 _Idbpgr_5 | 4.53703 1.103093 6.22 0.000 2.817203 7.306767 _Idbpgr_6 | 6.291702 1.600738 7.23 0.000 3.821246 10.35932 _Idbpgr_7 | 9.462228 2.566611 8.29 0.000 5.560408 16.10201 ------------------------------------------------------------------------------ . * . * Univariate analysis of the effect of gender on CHD . * . sts graph , by(sex) xlabel(0,5,10,15,20,25,30) ylabel(0,.1,.2,.3,.4,.5) failure l1title("Cumulative C > HD Morbidity") gap(3) noborder failure _d: chdfate analysis time _t: time . more . sts test sex failure _d: chdfate analysis time _t: time Log-rank test for equality of survivor functions | Events Events sex | observed expected ------+------------------------- Men | 823 589.47 Women | 650 883.53 ------+------------------------- Total | 1473 1473.00 chi2(1) = 154.57 Pr>chi2 = 0.0000 . generate male = sex==1 . tabulate male sex | Sex male | Men Women | Total -----------+----------------------+---------- 0 | 0 2650 | 2650 1 | 2049 0 | 2049 -----------+----------------------+---------- Total | 2049 2650 | 4699 . stcox male failure _d: chdfate analysis time _t: time Iteration 0: log likelihood = -11834.856 Iteration 1: log likelihood = -11759.624 Iteration 2: log likelihood = -11759.553 Refining estimates: Iteration 0: log likelihood = -11759.553 Cox regression -- Breslow method for ties No. of subjects = 4699 Number of obs = 4699 No. of failures = 1473 Time at risk = 103710.0917 LR chi2(1) = 150.61 Log likelihood = -11759.553 Prob > chi2 = 0.0000 ------------------------------------------------------------------------------ _t | _d | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- male | 1.900412 .0998308 12.22 0.000 1.714482 2.106504 ------------------------------------------------------------------------------ . * . * Fit multiplicative model of DBP and gender on risk of CHD . * . xi: stcox i.dbpgr male i.dbpgr _Idbpgr_1-7 (_Idbpgr_1 for dbpgr==60 omitted) failure _d: chdfate analysis time _t: time Iteration 0: log likelihood = -11834.856 Iteration 1: log likelihood = -11695.048 Iteration 2: log likelihood = -11659.586 Iteration 3: log likelihood = -11657.434 Iteration 4: log likelihood = -11657.409 Iteration 5: log likelihood = -11657.409 Refining estimates: Iteration 0: log likelihood = -11657.409 Cox regression -- Breslow method for ties No. of subjects = 4699 Number of obs = 4699 No. of failures = 1473 Time at risk = 103710.0917 LR chi2(7) = 354.89 Log likelihood = -11657.409 Prob > chi2 = 0.0000 ------------------------------------------------------------------------------ _t | _d | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _Idbpgr_2 | 1.911621 .4723633 2.62 0.009 1.177793 3.102662 _Idbpgr_3 | 2.429787 .585021 3.69 0.000 1.515737 3.895044 _Idbpgr_4 | 2.778377 .6697835 4.24 0.000 1.732176 4.456464 _Idbpgr_5 | 4.060083 .9879333 5.76 0.000 2.520075 6.541184 _Idbpgr_6 | 5.960225 1.516627 7.02 0.000 3.619658 9.814262 _Idbpgr_7 | 9.181868 2.490468 8.17 0.000 5.395767 15.6246 male | 1.833729 .0968002 11.49 0.000 1.653489 2.033616 ------------------------------------------------------------------------------ . lincom _Idbpgr_2 + male , hr ( 1) _Idbpgr_2 + male = 0.0 ------------------------------------------------------------------------------ _t | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- (1) | 3.505395 .8837535 4.98 0.000 2.138644 5.7456 ------------------------------------------------------------------------------ . lincom _Idbpgr_3 + male , hr ( 1) _Idbpgr_3 + male = 0.0 ------------------------------------------------------------------------------ _t | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- (1) | 4.455571 1.0938 6.09 0.000 2.753859 7.208835 ------------------------------------------------------------------------------ . lincom _Idbpgr_4 + male , hr ( 1) _Idbpgr_4 + male = 0.0 ------------------------------------------------------------------------------ _t | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- (1) | 5.094791 1.248052 6.65 0.000 3.152177 8.234591 ------------------------------------------------------------------------------ . lincom _Idbpgr_5 + male , hr ( 1) _Idbpgr_5 + male = 0.0 ------------------------------------------------------------------------------ _t | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- (1) | 7.445093 1.838249 8.13 0.000 4.588832 12.0792 ------------------------------------------------------------------------------ . lincom _Idbpgr_6 + male , hr ( 1) _Idbpgr_6 + male = 0.0 ------------------------------------------------------------------------------ _t | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- (1) | 10.92944 2.829359 9.24 0.000 6.580259 18.15318 ------------------------------------------------------------------------------ . lincom _Idbpgr_7 + male , hr ( 1) _Idbpgr_7 + male = 0.0 ------------------------------------------------------------------------------ _t | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- (1) | 16.83706 4.643535 10.24 0.000 9.806438 28.90822 ------------------------------------------------------------------------------ . display 2*(11723.942 -11657.409) 133.066 . display chi2tail(1,133.066) 8.746e-31 . * . * Fit model of DBP and gender on risk of CHD using interaction terms . * . xi: stcox i.dbpgr*i.male i.dbpgr _Idbpgr_1-7 (_Idbpgr_1 for dbpgr==60 omitted) i.male _Imale_0-1 (naturally coded; _Imale_0 omitted) i.dbpgr*i.male _IdbpXmal_#_# (coded as above) failure _d: chdfate analysis time _t: time Iteration 0: log likelihood = -11834.856 Iteration 1: log likelihood = -11689.462 Iteration 2: log likelihood = -11648.634 Iteration 3: log likelihood = -11646.81 Iteration 4: log likelihood = -11646.794 Iteration 5: log likelihood = -11646.794 Refining estimates: Iteration 0: log likelihood = -11646.794 Cox regression -- Breslow method for ties No. of subjects = 4699 Number of obs = 4699 No. of failures = 1473 Time at risk = 103710.0917 LR chi2(13) = 376.12 Log likelihood = -11646.794 Prob > chi2 = 0.0000 ------------------------------------------------------------------------------ _t | _d | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _Idbpgr_2 | 1.82731 .6428651 1.71 0.087 .9169625 3.64144 _Idbpgr_3 | 2.428115 .8298216 2.60 0.009 1.2427 4.744299 _Idbpgr_4 | 3.517929 1.201355 3.68 0.000 1.801384 6.870179 _Idbpgr_5 | 4.693559 1.628053 4.46 0.000 2.378188 9.263141 _Idbpgr_6 | 7.635131 2.736437 5.67 0.000 3.782205 15.41302 _Idbpgr_7 | 13.62563 5.067901 7.02 0.000 6.572973 28.24565 _Imale_1 | 2.372645 1.118489 1.83 0.067 .9418198 5.977199 _IdbpXma~2_1 | 1.058632 .5235583 0.12 0.908 .4015814 2.79072 _IdbpXma~3_1 | .9628061 .4637697 -0.08 0.937 .3745652 2.474858 _IdbpXma~4_1 | .6324678 .3047828 -0.95 0.342 .2459512 1.626402 _IdbpXma~5_1 | .7437487 .3621623 -0.61 0.543 .2863787 1.931576 _IdbpXma~6_1 | .6015939 .3059896 -1.00 0.318 .2220014 1.630239 _IdbpXma~7_1 | .401376 .2205419 -1.66 0.097 .1367245 1.178302 ------------------------------------------------------------------------------ . display 2*( 11657.409 -11646.794) 21.23 . display chi2tail(6, 21.23) .00166794 . lincom _Idbpgr_2 + _Imale_1 + _IdbpXmal_2_1, hr ( 1) _Idbpgr_2 + _Imale_1 + _IdbpXmal_2_1 = 0.0 ------------------------------------------------------------------------------ _t | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- (1) | 4.589761 1.595446 4.38 0.000 2.322223 9.071437 ------------------------------------------------------------------------------ . lincom _Idbpgr_3 + _Imale_1 + _IdbpXmal_3_1, hr ( 1) _Idbpgr_3 + _Imale_1 + _IdbpXmal_3_1 = 0.0 ------------------------------------------------------------------------------ _t | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- (1) | 5.546778 1.883042 5.05 0.000 2.851487 10.78972 ------------------------------------------------------------------------------ . lincom _Idbpgr_4 + _Imale_1 + _IdbpXmal_4_1, hr ( 1) _Idbpgr_4 + _Imale_1 + _IdbpXmal_4_1 = 0.0 ------------------------------------------------------------------------------ _t | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- (1) | 5.27908 1.795204 4.89 0.000 2.710807 10.28059 ------------------------------------------------------------------------------ . lincom _Idbpgr_5 + _Imale_1 + _IdbpXmal_5_1, hr ( 1) _Idbpgr_5 + _Imale_1 + _IdbpXmal_5_1 = 0.0 ------------------------------------------------------------------------------ _t | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- (1) | 8.282496 2.831657 6.18 0.000 4.23788 16.18728 ------------------------------------------------------------------------------ . lincom _Idbpgr_6 + _Imale_1 + _IdbpXmal_6_1, hr ( 1) _Idbpgr_6 + _Imale_1 + _IdbpXmal_6_1 = 0.0 ------------------------------------------------------------------------------ _t | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- (1) | 10.89815 3.935989 6.61 0.000 5.369473 22.11942 ------------------------------------------------------------------------------ . lincom _Idbpgr_7 + _Imale_1 + _IdbpXmal_7_1, hr ( 1) _Idbpgr_7 + _Imale_1 + _IdbpXmal_7_1 = 0.0 ------------------------------------------------------------------------------ _t | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- (1) | 12.976 5.253189 6.33 0.000 5.868692 28.69064 ------------------------------------------------------------------------------ . * . * Adjust model for age, BMI and SCL . * . xi: stcox i.dbpgr*i.male age i.dbpgr _Idbpgr_1-7 (_Idbpgr_1 for dbpgr==60 omitted) i.male _Imale_0-1 (naturally coded; _Imale_0 omitted) i.dbpgr*i.male _IdbpXmal_#_# (coded as above) failure _d: chdfate analysis time _t: time Iteration 0: log likelihood = -11834.856 Iteration 1: log likelihood = -11555.813 Iteration 2: log likelihood = -11518.758 Iteration 3: log likelihood = -11517.258 Iteration 4: log likelihood = -11517.247 Iteration 5: log likelihood = -11517.247 Refining estimates: Iteration 0: log likelihood = -11517.247 Cox regression -- Breslow method for ties No. of subjects = 4699 Number of obs = 4699 No. of failures = 1473 Time at risk = 103710.0917 LR chi2(14) = 635.22 Log likelihood = -11517.247 Prob > chi2 = 0.0000 ------------------------------------------------------------------------------ _t | _d | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _Idbpgr_2 | 1.708833 .6012239 1.52 0.128 .8574694 3.405496 _Idbpgr_3 | 2.036418 .6963299 2.08 0.038 1.041857 3.980389 _Idbpgr_4 | 2.54474 .8704832 2.73 0.006 1.301582 4.975254 _Idbpgr_5 | 2.886965 1.004923 3.05 0.002 1.459308 5.711315 _Idbpgr_6 | 4.676467 1.681721 4.29 0.000 2.311076 9.462841 _Idbpgr_7 | 8.149445 3.040448 5.62 0.000 3.922439 16.93167 _Imale_1 | 2.126092 1.002396 1.60 0.110 .8438459 5.356747 _IdbpXma~2_1 | 1.081597 .5349175 0.16 0.874 .4102917 2.851268 _IdbpXma~3_1 | 1.078892 .5197532 0.16 0.875 .4196762 2.773588 _IdbpXma~4_1 | .8176247 .3942289 -0.42 0.676 .3177866 2.103645 _IdbpXma~5_1 | 1.063911 .5185957 0.13 0.899 .4092542 2.765778 _IdbpXma~6_1 | .8651554 .4405249 -0.28 0.776 .3189145 2.347005 _IdbpXma~7_1 | .5022462 .2760492 -1.25 0.210 .1710296 1.474898 age | 1.055366 .003514 16.18 0.000 1.048501 1.062276 ------------------------------------------------------------------------------ . display 2*(11646.794 -11517.247) 259.094 . display chi2tail(1,259.094) 2.704e-58 . xi: stcox i.dbpgr*i.male age bmi i.dbpgr _Idbpgr_1-7 (_Idbpgr_1 for dbpgr==60 omitted) i.male _Imale_0-1 (naturally coded; _Imale_0 omitted) i.dbpgr*i.male _IdbpXmal_#_# (coded as above) failure _d: chdfate analysis time _t: time Iteration 0: log likelihood = -11825.305 Iteration 1: log likelihood = -11529.516 Iteration 2: log likelihood = -11492.307 Iteration 3: log likelihood = -11490.745 Iteration 4: log likelihood = -11490.733 Iteration 5: log likelihood = -11490.733 Refining estimates: Iteration 0: log likelihood = -11490.733 Cox regression -- Breslow method for ties No. of subjects = 4690 Number of obs = 4690 No. of failures = 1472 Time at risk = 103575.6441 LR chi2(15) = 669.14 Log likelihood = -11490.733 Prob > chi2 = 0.0000 ------------------------------------------------------------------------------ _t | _d | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _Idbpgr_2 | 1.610195 .5667799 1.35 0.176 .807719 3.20994 _Idbpgr_3 | 1.842678 .6309583 1.79 0.074 .9418602 3.605061 _Idbpgr_4 | 2.191021 .7518132 2.29 0.022 1.118332 4.292621 _Idbpgr_5 | 2.323355 .8136951 2.41 0.016 1.169511 4.615587 _Idbpgr_6 | 3.772948 1.365324 3.67 0.000 1.856328 7.668436 _Idbpgr_7 | 5.930829 2.23934 4.71 0.000 2.82958 12.43108 _Imale_1 | 2.031802 .9580764 1.50 0.133 .8063167 5.119849 _IdbpXma~2_1 | 1.114935 .5514427 0.22 0.826 .4229103 2.939346 _IdbpXma~3_1 | 1.113261 .5363989 0.22 0.824 .4329775 2.862389 _IdbpXma~4_1 | .8417087 .4059036 -0.36 0.721 .3270999 2.165924 _IdbpXma~5_1 | 1.13761 .5547716 0.26 0.791 .4374145 2.958652 _IdbpXma~6_1 | .9195821 .4684693 -0.16 0.869 .3388105 2.495883 _IdbpXma~7_1 | .5731364 .3153943 -1.01 0.312 .1949152 1.685273 age | 1.054176 .0035064 15.86 0.000 1.047326 1.061071 bmi | 1.040512 .0068282 6.05 0.000 1.027215 1.053982 ------------------------------------------------------------------------------ . display 2*(11517.247 -11490.733 ) 53.028 . display chi2tail(1,53.028) 3.288e-13 . xi: stcox i.dbpgr*i.male age bmi scl, mgale(mg) i.dbpgr _Idbpgr_1-7 (_Idbpgr_1 for dbpgr==60 omitted) i.male _Imale_0-1 (naturally coded; _Imale_0 omitted) i.dbpgr*i.male _IdbpXmal_#_# (coded as above) failure _d: chdfate analysis time _t: time Iteration 0: log likelihood = -11758.886 Iteration 1: log likelihood = -11419.565 Iteration 2: log likelihood = -11383.554 Iteration 3: log likelihood = -11382.142 Iteration 4: log likelihood = -11382.132 Iteration 5: log likelihood = -11382.132 Refining estimates: Iteration 0: log likelihood = -11382.132 Cox regression -- Breslow method for ties No. of subjects = 4658 Number of obs = 4658 No. of failures = 1465 Time at risk = 102895.0938 LR chi2(16) = 753.51 Log likelihood = -11382.132 Prob > chi2 = 0.0000 ------------------------------------------------------------------------------ _t | _d | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _Idbpgr_2 | 1.514961 .5334695 1.18 0.238 .7597392 3.020916 _Idbpgr_3 | 1.654264 .5669665 1.47 0.142 .8450299 3.238451 _Idbpgr_4 | 1.911763 .6566924 1.89 0.059 .9750921 3.748199 _Idbpgr_5 | 1.936029 .6796612 1.88 0.060 .9729479 3.852425 _Idbpgr_6 | 3.097614 1.123672 3.12 0.002 1.521425 6.306727 _Idbpgr_7 | 5.269096 1.988701 4.40 0.000 2.514603 11.04086 _Imale_1 | 1.984033 .9355668 1.45 0.146 .7873473 4.999554 _IdbpXma~2_1 | 1.173058 .5802796 0.32 0.747 .4448907 3.09304 _IdbpXma~3_1 | 1.18152 .5693995 0.35 0.729 .4594405 3.038457 _IdbpXma~4_1 | .8769476 .4230106 -0.27 0.785 .3407078 2.257175 _IdbpXma~5_1 | 1.265976 .6179759 0.48 0.629 .4863156 3.295585 _IdbpXma~6_1 | 1.023429 .5215766 0.05 0.964 .3769245 2.778823 _IdbpXma~7_1 | .6125694 .3371363 -0.89 0.373 .2082976 1.801467 age | 1.04863 .003559 13.99 0.000 1.041677 1.055628 bmi | 1.038651 .0070125 5.62 0.000 1.024998 1.052487 scl | 1.005788 .0005883 9.87 0.000 1.004635 1.006941 ------------------------------------------------------------------------------ . display 2*(11490.733 -11382.132) 217.202 . display chi2tail(1,217.202) 3.687e-49 . lincom _Idbpgr_2 + _Imale_1 + _IdbpXmal_2_1, hr ( 1) _Idbpgr_2 + _Imale_1 + _IdbpXmal_2_1 = 0.0 ------------------------------------------------------------------------------ _t | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- (1) | 3.525899 1.22701 3.62 0.000 1.782593 6.974087 ------------------------------------------------------------------------------ . lincom _Idbpgr_3 + _Imale_1 + _IdbpXmal_3_1, hr ( 1) _Idbpgr_3 + _Imale_1 + _IdbpXmal_3_1 = 0.0 ------------------------------------------------------------------------------ _t | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- (1) | 3.877883 1.319887 3.98 0.000 1.990112 7.556349 ------------------------------------------------------------------------------ . lincom _Idbpgr_4 + _Imale_1 + _IdbpXmal_4_1, hr ( 1) _Idbpgr_4 + _Imale_1 + _IdbpXmal_4_1 = 0.0 ------------------------------------------------------------------------------ _t | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- (1) | 3.326263 1.136427 3.52 0.000 1.702712 6.497885 ------------------------------------------------------------------------------ . lincom _Idbpgr_5 + _Imale_1 + _IdbpXmal_5_1, hr ( 1) _Idbpgr_5 + _Imale_1 + _IdbpXmal_5_1 = 0.0 ------------------------------------------------------------------------------ _t | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- (1) | 4.862796 1.673288 4.60 0.000 2.477352 9.545188 ------------------------------------------------------------------------------ . lincom _Idbpgr_6 + _Imale_1 + _IdbpXmal_6_1, hr ( 1) _Idbpgr_6 + _Imale_1 + _IdbpXmal_6_1 = 0.0 ------------------------------------------------------------------------------ _t | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- (1) | 6.289755 2.284736 5.06 0.000 3.086293 12.8183 ------------------------------------------------------------------------------ . lincom _Idbpgr_7 + _Imale_1 + _IdbpXmal_7_1, hr ( 1) _Idbpgr_7 + _Imale_1 + _IdbpXmal_7_1 = 0.0 ------------------------------------------------------------------------------ _t | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- (1) | 6.403836 2.608058 4.56 0.000 2.882542 14.22672 ------------------------------------------------------------------------------ . * . * Perform Cox-Snell generalized residual analysis . * . predict cs, csnell (41 missing values generated) . stset cs, failure(chdfate) failure event: chdfate ~= 0 & chdfate ~= . obs. time interval: (0, cs] exit on or before: failure ------------------------------------------------------------------------------ 4699 total obs. 41 event time missing (cs==.) PROBABLE ERROR ------------------------------------------------------------------------------ 4658 obs. remaining, representing 1465 failures in single record/single failure data 1465 total analysis time at risk, at risk from t = 0 earliest observed entry t = 0 last observed exit t = 2.833814 . sts generate km=s . generate es = exp(-cs) (41 missing values generated) . sort cs . graph km es cs, connect(ll) symbol(..) xlabel(0 .5 to 2.5) ylabel(0 .2 to 1.0) xtick(.25 .75 to 2.7 > 5) ytick(.1 .2 to 1) . log close log: C:\WDDtext\log\7.7.Framingham.log log type: text closed on: 29 Aug 2001, 17:57:09 -------------------------------------------------------------------------------------------------------