------------------------------------------------------------------------------------------------- log: C:\MyDocs\MPH\Text\SecondEdition\WebDoFiles\7.6.Framingham.log log type: text opened on: 3 Jan 2008, 17:22:59 . * 7.6.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 Current memory allocation current memory usage settable value description (1M = 1024k) -------------------------------------------------------------------- set maxvar 5000 max. variables allowed 1.909M set memory 2000k max. data space 1.953M set matsize 400 max. RHS vars in models 1.254M ----------- 5.116M . use C:\WDDtext\2.20.Framingham.dta, clear . * . * Univariate analysis of the effect of DBP on CHD . * . histogram dbp, start(39) width(2) frequency ylabel(0(100)500) /// > ymtick(0(25)525) xlabel(40(20)140) xmtick(40(5)140) /// > ytitle(Number of Study Subjects) (bin=55, start=39, width=2) . 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 | 1,048 419 | 1,467 90 | 863 404 | 1,267 100 | 417 284 | 701 110 | 125 110 | 235 111 | 49 56 | 105 -----------+----------------------+---------- Total | 3,226 1,473 | 4,699 . label define dbp 60 " DBP <= 60" 70 "60 < DBP <= 70" /// > 80 "70 < DBP <= 80" 90 "80 < DBP <= 90" /// > 100 "90 < DBP <= 100" 110 "100 < DBP <= 110" /// > 111 "110 < DBP" . label variable dbpgr "DBP level" . 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) ytitle(Proportion Without CHD) /// > ylabel(0(.1)1, angle(0)) xlabel(0(5)30) xmtick(1(1)32) /// > xtitle(Years of Follow-up) plot2opts(color(gray)) /// > plot3opts(lpattern(dash)) /// > plot4opts(lpattern(dash)) /// > plot5opts(lpattern(longdash)) /// > plot6opts(lpattern(longdash)) /// > plot7opts(lpattern("_-")) /// > title(" ",size(zero)) legend(ring(0) position(7) col(1) /// > order(1 "DBP <= 60" 2 "60 < DBP <= 70" /// > 3 "70 < DBP <= 80" 4 "80 < DBP <= 90" /// > 5 "90 < DBP <= 100" 6 "100 < DBP <= 110" 7 "110 < DBP")) 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 60 < DBP <= 70 | 182 275.72 70 < DBP <= 80 | 419 489.41 80 < DBP <= 90 | 404 395.62 90 < DBP <= 100 | 284 187.97 100 < DBP <= 110 | 110 52.73 110 < DBP | 56 17.94 -----------------+------------------------- Total | 1473 1473.00 chi2(6) = 259.71 Pr>chi2 = 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 60 < DBP <= 70 | 182 167.42 ----------------+------------------------- Total | 200 200.00 chi2(1) = 7.80 Pr>chi2 = 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 ----------------+------------------------- 60 < DBP <= 70 | 182 217.24 70 < DBP <= 80 | 419 383.76 ----------------+------------------------- Total | 601 601.00 chi2(1) = 8.96 Pr>chi2 = 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 ----------------+------------------------- 70 < DBP <= 80 | 419 456.22 80 < DBP <= 90 | 404 366.78 ----------------+------------------------- Total | 823 823.00 chi2(1) = 6.82 Pr>chi2 = 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 -----------------+------------------------- 80 < DBP <= 90 | 404 466.63 90 < DBP <= 100 | 284 221.37 -----------------+------------------------- Total | 688 688.00 chi2(1) = 26.18 Pr>chi2 = 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 -----------------+------------------------- 90 < DBP <= 100 | 284 306.92 100 < DBP <= 110 | 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 -----------------+------------------------- 100 < DBP <= 110 | 110 122.95 110 < DBP | 56 43.05 -----------------+------------------------- Total | 166 166.00 chi2(1) = 5.29 Pr>chi2 = 0.0215 . * . * Hazard regression analysis with dbp - 60 as a continuous . * covariate. . * . generate dbp60 = dbp - 60 . stcox dbp60 failure _d: chdfate analysis time _t: time Iteration 0: log likelihood = -11834.856 Iteration 1: log likelihood = -11716.395 Iteration 2: log likelihood = -11713.818 Iteration 3: log likelihood = -11713.816 Refining estimates: Iteration 0: log likelihood = -11713.816 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) = 242.08 Log likelihood = -11713.816 Prob > chi2 = 0.0000 ------------------------------------------------------------------------------ _t | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- dbp60 | 1.032064 .0019926 16.35 0.000 1.028166 1.035977 ------------------------------------------------------------------------------ . estimates store simple . predict relhaz1, hr . predict loghaz, xb . predict se, stdp . generate logcil = loghaz - 1.96*se . generate logciu = loghaz + 1.96*se . generate cil = exp(logcil) . generate ciu = exp(logciu) . sort dbp . * . * Plot risk of CHD for patients with a given baseline DBP . * relative to patients with a basline DBP of 60. Show . * the 95% confidence band for this plot. . * . twoway rarea cil ciu dbp, color(yellow) /// > || line relhaz1 dbp /// > , legend(off) ytitle(Relative risk of CHD) /// > ymtick(1(1)16) ylabel(0(3)15) xlabel(40(10)140) /// > xmtick(45(5)145) . more . twoway rarea cil ciu dbp, color(yellow) lwidth(none) /// > || line relhaz1 dbp /// > , legend(off) ymtick(1(1)16) ylabel(0(3)15) /// > xlabel(40(10)140) xmtick(45(5)145) ytitle(Relative risk of CHD) . more . * . * Restricted cubic spline model of the effect of DBP on CHD . * risk. Use a 5 knot model with default knot values. . * . mkspline _Sdbp60 = dbp60, cubic displayknots | knot1 knot2 knot3 knot4 knot5 -------------+------------------------------------------------------- dbp60 | 4 14 20 29.5 45 . stcox _S*, nohr failure _d: chdfate analysis time _t: time Iteration 0: log likelihood = -11834.856 Iteration 1: log likelihood = -11756.455 Iteration 2: log likelihood = -11715.422 Iteration 3: log likelihood = -11711.467 Iteration 4: log likelihood = -11711.393 Iteration 5: log likelihood = -11711.393 Refining estimates: Iteration 0: log likelihood = -11711.393 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(4) = 246.93 Log likelihood = -11711.393 Prob > chi2 = 0.0000 ------------------------------------------------------------------------------ _t | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _Sdbp601 | .0618603 .016815 3.68 0.000 .0289035 .094817 _Sdbp602 | -.2268319 .1120642 -2.02 0.043 -.4464737 -.0071902 _Sdbp603 | .93755 .4547913 2.06 0.039 .0461754 1.828925 _Sdbp604 | -.982937 .4821521 -2.04 0.041 -1.927938 -.0379362 ------------------------------------------------------------------------------ . lrtest simple . Likelihood-ratio test LR chi2(3) = 4.85 (Assumption: simple nested in .) Prob > chi2 = 0.1835 . * . * Let's try a three knot model . * . drop _S* . mkspline _Sdbp60 = dbp60, cubic nknots(3) displayknots | knot1 knot2 knot3 -------------+--------------------------------- dbp60 | 8 20 40 . stcox _S*, nohr failure _d: chdfate analysis time _t: time Iteration 0: log likelihood = -11834.856 Iteration 1: log likelihood = -11746.134 Iteration 2: log likelihood = -11715.117 Iteration 3: log likelihood = -11713.647 Iteration 4: log likelihood = -11713.643 Iteration 5: log likelihood = -11713.643 Refining estimates: Iteration 0: log likelihood = -11713.643 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(2) = 242.43 Log likelihood = -11713.643 Prob > chi2 = 0.0000 ------------------------------------------------------------------------------ _t | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _Sdbp601 | .0347213 .0057337 6.06 0.000 .0234835 .0459592 _Sdbp602 | -.0041479 .0070762 -0.59 0.558 -.0180169 .0097212 ------------------------------------------------------------------------------ . lrtest simple . Likelihood-ratio test LR chi2(1) = 0.35 (Assumption: simple nested in .) Prob > chi2 = 0.5560 . more . * . * The range of DBPs above the last knot is 40 mm Hg. . * Let's add another knot at DBP = 120 to try to avoid . * excessive model extrapolation for high DBPs . * . drop _S* . mkspline _Sdbp60 = dbp60, cubic knots(8 20 40 60) . stcox _S*, nohr failure _d: chdfate analysis time _t: time Iteration 0: log likelihood = -11834.856 Iteration 1: log likelihood = -11755.789 Iteration 2: log likelihood = -11717.419 Iteration 3: log likelihood = -11713.066 Iteration 4: log likelihood = -11712.891 Iteration 5: log likelihood = -11712.89 Refining estimates: Iteration 0: log likelihood = -11712.89 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(3) = 243.93 Log likelihood = -11712.89 Prob > chi2 = 0.0000 ------------------------------------------------------------------------------ _t | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _Sdbp601 | .029958 .0068727 4.36 0.000 .0164878 .0434282 _Sdbp602 | .0229262 .0332729 0.69 0.491 -.0422876 .0881399 _Sdbp603 | -.0569008 .0676407 -0.84 0.400 -.189474 .0756725 ------------------------------------------------------------------------------ . lrtest simple . Likelihood-ratio test LR chi2(2) = 1.85 (Assumption: simple nested in .) Prob > chi2 = 0.3963 . more . drop loghaz se logcil logciu cil ciu . predict loghaz, xb . predict relhaz4, hr . predict se, stdp . generate logcil = loghaz - 1.96*se . generate logciu = loghaz + 1.96*se . generate cil = exp(logcil) . generate ciu = exp(logciu) . * . * Plot relative risks of CHD with 95% confidence bands . * for this model. . * . twoway rarea cil ciu dbp, color(yellow) /// > || line relhaz1 dbp, color(white) /// > || line relhaz4 dbp, lpattern(solid) /// > , legend(off) ytitle(Relative risk of CHD) /// > ymtick(1(1)11) ylabel(0(2)10) xlabel(40(10)140) /// > xmtick(45(5)145) . more . replace ciu = min(ciu, 11.41) (8 real changes made) . replace relhaz1 = min(relhaz1,12.5) (3 real changes made) . twoway rarea cil ciu dbp, color(yellow) lwidth(none) /// > || line relhaz1 dbp, color(white) /// > || line relhaz4 dbp, lpattern(solid) /// > , legend(off) ytitle(Relative risk of CHD) /// > ymtick(1(1)11) ylabel(0(2)10) /// > xlabel(40(10)140) xmtick(45(5)145) . more . * . * Model effects of DBP on CHD risk using a categroical . * variable for DBP . * . xi: stcox i.dbpgr i.dbpgr _Idbpgr_60-111 (naturally coded; _Idbpgr_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 | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _Idbpgr_70 | 1.968764 .486453 2.74 0.006 1.213037 3.195312 _Idbpgr_80 | 2.557839 .6157326 3.90 0.000 1.595764 4.099941 _Idbpgr_90 | 3.056073 .7362768 4.64 0.000 1.905856 4.900466 _Idbpgr_100 | 4.53703 1.103093 6.22 0.000 2.817203 7.306767 _Idbpgr_110 | 6.291702 1.600738 7.23 0.000 3.821246 10.35932 _Idbpgr_111 | 9.462228 2.566611 8.29 0.000 5.560408 16.10201 ------------------------------------------------------------------------------ . * . * Store estimates from this model for future likelihood ratio . * tests (tests of change in model deviance). . * . estimates store _dbpgr . * . * Univariate analysis of the effect of gender on CHD . * . sts graph, by(sex) failure ylabel(0(.1).5) /// > ytitle(Cumulative CHD Morbidity) xlabel(0(5)30) /// > xtitle(Years of Follow-up) title(" ", ring(0) position(11)) /// > legend(off) 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 if !missing(sex) . tabulate male sex | Sex male | Men Women | Total -----------+----------------------+---------- 0 | 0 2,650 | 2,650 1 | 2,049 0 | 2,049 -----------+----------------------+---------- Total | 2,049 2,650 | 4,699 . 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 | 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_60-111 (naturally coded; _Idbpgr_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 | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _Idbpgr_70 | 1.911621 .4723633 2.62 0.009 1.177793 3.102662 _Idbpgr_80 | 2.429787 .585021 3.69 0.000 1.515737 3.895044 _Idbpgr_90 | 2.778377 .6697835 4.24 0.000 1.732176 4.456464 _Idbpgr_100 | 4.060083 .9879333 5.76 0.000 2.520075 6.541184 _Idbpgr_110 | 5.960225 1.516627 7.02 0.000 3.619658 9.814262 _Idbpgr_111 | 9.181868 2.490468 8.17 0.000 5.395767 15.6246 male | 1.833729 .0968002 11.49 0.000 1.653489 2.033616 ------------------------------------------------------------------------------ . lrtest _dbpgr . Likelihood-ratio test LR chi2(1) = 133.07 (Assumption: _dbpgr nested in .) Prob > chi2 = 0.0000 . estimates store dbp_male . lincom _Idbpgr_70 + male , hr ( 1) _Idbpgr_70 + male = 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_80 + male , hr ( 1) _Idbpgr_80 + male = 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_90 + male , hr ( 1) _Idbpgr_90 + male = 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_100 + male , hr ( 1) _Idbpgr_100 + male = 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_110 + male , hr ( 1) _Idbpgr_110 + male = 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_111 + male , hr ( 1) _Idbpgr_111 + male = 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 ------------------------------------------------------------------------------ . * . * Fit model of DBP and gender on risk of CHD using interaction terms . * . xi: stcox i.dbpgr*male i.dbpgr _Idbpgr_60-111 (naturally coded; _Idbpgr_60 omitted) i.dbpgr*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 | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _Idbpgr_70 | 1.82731 .6428651 1.71 0.087 .9169625 3.64144 _Idbpgr_80 | 2.428115 .8298216 2.60 0.009 1.2427 4.744299 _Idbpgr_90 | 3.517929 1.201355 3.68 0.000 1.801384 6.870179 _Idbpgr_100 | 4.693559 1.628053 4.46 0.000 2.378188 9.263141 _Idbpgr_110 | 7.635131 2.736437 5.67 0.000 3.782205 15.41302 _Idbpgr_111 | 13.62563 5.067901 7.02 0.000 6.572973 28.24565 male | 2.372645 1.118489 1.83 0.067 .9418198 5.977199 _IdbpXmal_70 | 1.058632 .5235583 0.12 0.908 .4015814 2.79072 _IdbpXmal_80 | .9628061 .4637697 -0.08 0.937 .3745652 2.474858 _IdbpXmal_90 | .6324678 .3047828 -0.95 0.342 .2459512 1.626402 _IdbpXma~100 | .7437487 .3621623 -0.61 0.543 .2863787 1.931576 _IdbpXma~110 | .6015939 .3059896 -1.00 0.318 .2220014 1.630239 _IdbpXma~111 | .401376 .2205419 -1.66 0.097 .1367245 1.178302 ------------------------------------------------------------------------------ . lrtest dbp_male . Likelihood-ratio test LR chi2(6) = 21.23 (Assumption: dbp_male nested in .) Prob > chi2 = 0.0017 . estimates store dbp_maleInteract . lincom _Idbpgr_70 + male + _IdbpXmal_70, hr ( 1) _Idbpgr_70 + male + _IdbpXmal_70 = 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_80 + male + _IdbpXmal_80, hr ( 1) _Idbpgr_80 + male + _IdbpXmal_80 = 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_90 + male + _IdbpXmal_90, hr ( 1) _Idbpgr_90 + male + _IdbpXmal_90 = 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_100 + male + _IdbpXmal_100, hr ( 1) _Idbpgr_100 + male + _IdbpXmal_100 = 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_110 + male + _IdbpXmal_110, hr ( 1) _Idbpgr_110 + male + _IdbpXmal_110 = 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_111 + male + _IdbpXmal_111, hr ( 1) _Idbpgr_111 + male + _IdbpXmal_111 = 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*male age i.dbpgr _Idbpgr_60-111 (naturally coded; _Idbpgr_60 omitted) i.dbpgr*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 | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _Idbpgr_70 | 1.708833 .6012239 1.52 0.128 .8574694 3.405496 _Idbpgr_80 | 2.036418 .6963299 2.08 0.038 1.041857 3.980389 _Idbpgr_90 | 2.54474 .8704832 2.73 0.006 1.301582 4.975254 _Idbpgr_100 | 2.886965 1.004923 3.05 0.002 1.459308 5.711315 _Idbpgr_110 | 4.676467 1.681721 4.29 0.000 2.311076 9.462841 _Idbpgr_111 | 8.149445 3.040448 5.62 0.000 3.922439 16.93167 male | 2.126092 1.002396 1.60 0.110 .8438459 5.356747 _IdbpXmal_70 | 1.081597 .5349175 0.16 0.874 .4102917 2.851268 _IdbpXmal_80 | 1.078892 .5197532 0.16 0.875 .4196762 2.773588 _IdbpXmal_90 | .8176247 .3942289 -0.42 0.676 .3177866 2.103645 _IdbpXma~100 | 1.063911 .5185957 0.13 0.899 .4092542 2.765778 _IdbpXma~110 | .8651554 .4405249 -0.28 0.776 .3189145 2.347005 _IdbpXma~111 | .5022462 .2760492 -1.25 0.210 .1710296 1.474898 age | 1.055366 .003514 16.18 0.000 1.048501 1.062276 ------------------------------------------------------------------------------ . lrtest dbp_maleInteract . Likelihood-ratio test LR chi2(1) = 259.09 (Assumption: dbp_maleInte~t nested in .) Prob > chi2 = 0.0000 . xi: stcox i.dbpgr*male age if !missing(bmi) & !missing(scl) i.dbpgr _Idbpgr_60-111 (naturally coded; _Idbpgr_60 omitted) i.dbpgr*male _IdbpXmal_# (coded as above) failure _d: chdfate analysis time _t: time Iteration 0: log likelihood = -11758.886 Iteration 1: log likelihood = -11482.835 Iteration 2: log likelihood = -11446.25 Iteration 3: log likelihood = -11444.769 Iteration 4: log likelihood = -11444.759 Iteration 5: log likelihood = -11444.759 Refining estimates: Iteration 0: log likelihood = -11444.759 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(14) = 628.25 Log likelihood = -11444.759 Prob > chi2 = 0.0000 ------------------------------------------------------------------------------ _t | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _Idbpgr_70 | 1.734122 .6101188 1.56 0.118 .8701616 3.455884 _Idbpgr_80 | 2.064992 .7060844 2.12 0.034 1.056493 4.036179 _Idbpgr_90 | 2.551707 .8728997 2.74 0.006 1.305112 4.989003 _Idbpgr_100 | 2.799698 .9759998 2.95 0.003 1.413757 5.544312 _Idbpgr_110 | 4.745228 1.706296 4.33 0.000 2.345205 9.601374 _Idbpgr_111 | 8.121299 3.03002 5.61 0.000 3.908823 16.87349 male | 2.145405 1.011508 1.62 0.105 .8515059 5.405436 _IdbpXmal_70 | 1.061787 .5251211 0.12 0.904 .4027762 2.79905 _IdbpXmal_80 | 1.056844 .5091531 0.11 0.909 .4110831 2.717016 _IdbpXmal_90 | .8005036 .3860191 -0.46 0.644 .3110975 2.059823 _IdbpXma~100 | 1.079666 .5267195 0.16 0.875 .4149804 2.808999 _IdbpXma~110 | .8434902 .4294899 -0.33 0.738 .3109307 2.288213 _IdbpXma~111 | .4980528 .2737471 -1.27 0.205 .1695999 1.462599 age | 1.055745 .0035312 16.22 0.000 1.048847 1.062689 ------------------------------------------------------------------------------ . estimates store dbp_maleInteract_age . xi: stcox i.dbpgr*male age bmi scl i.dbpgr _Idbpgr_60-111 (naturally coded; _Idbpgr_60 omitted) i.dbpgr*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 | Haz. Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _Idbpgr_70 | 1.514961 .5334695 1.18 0.238 .7597392 3.020916 _Idbpgr_80 | 1.654264 .5669665 1.47 0.142 .8450299 3.238451 _Idbpgr_90 | 1.911763 .6566924 1.89 0.059 .9750921 3.748199 _Idbpgr_100 | 1.936029 .6796612 1.88 0.060 .9729479 3.852425 _Idbpgr_110 | 3.097614 1.123672 3.12 0.002 1.521425 6.306727 _Idbpgr_111 | 5.269096 1.988701 4.40 0.000 2.514603 11.04086 male | 1.984033 .9355668 1.45 0.146 .7873473 4.999554 _IdbpXmal_70 | 1.173058 .5802796 0.32 0.747 .4448907 3.09304 _IdbpXmal_80 | 1.18152 .5693995 0.35 0.729 .4594405 3.038457 _IdbpXmal_90 | .8769476 .4230106 -0.27 0.785 .3407078 2.257175 _IdbpXma~100 | 1.265976 .6179759 0.48 0.629 .4863156 3.295585 _IdbpXma~110 | 1.023429 .5215766 0.05 0.964 .3769245 2.778823 _IdbpXma~111 | .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 ------------------------------------------------------------------------------ . lrtest dbp_maleInteract_age . Likelihood-ratio test LR chi2(2) = 125.25 (Assumption: dbp_maleInte~e nested in .) Prob > chi2 = 0.0000 . lincom _Idbpgr_70 + male + _IdbpXmal_70, hr ( 1) _Idbpgr_70 + male + _IdbpXmal_70 = 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_80 + male + _IdbpXmal_80, hr ( 1) _Idbpgr_80 + male + _IdbpXmal_80 = 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_90 + male + _IdbpXmal_90, hr ( 1) _Idbpgr_90 + male + _IdbpXmal_90 = 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_100 + male + _IdbpXmal_100, hr ( 1) _Idbpgr_100 + male + _IdbpXmal_100 = 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_110 + male + _IdbpXmal_110, hr ( 1) _Idbpgr_110 + male + _IdbpXmal_110 = 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_111 + male + _IdbpXmal_111, hr ( 1) _Idbpgr_111 + male + _IdbpXmal_111 = 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 ------------------------------------------------------------------------------ . . log close log: C:\MyDocs\MPH\Text\SecondEdition\WebDoFiles\7.6.Framingham.log log type: text closed on: 3 Jan 2008, 17:24:05 -----------------------------------------------------------------------------------------------