*change name of path to read in fulldata.dat . infile pubid entry exit event p1 p2 p3 p4 per1234 using C:\mikefolder\mikeFolder\survivalmethods\programsforweb\Figures_revise > dpaper\fulldata.dat (2314 observations read) . stset exit, failure(event) enter(entry) failure event: event != 0 & event < . obs. time interval: (0, exit] enter on or after: time entry exit on or before: failure ------------------------------------------------------------------------------ 2314 total obs. 0 exclusions ------------------------------------------------------------------------------ 2314 obs. remaining, representing 1053 failures in single record/single failure data 4478.968 total analysis time at risk, at risk from t = 0 earliest observed entry t = 0 last observed exit t = 8.953 *The command below "streg p2 p3 p4, d(gamma) anc(p2 p3 p4) anc2(p2 p3 p4)" will fit a saturated GG model for all periods *with one command as opposed to doing it separately for each period. . streg p2 p3 p4, d(gamma) anc(p2 p3 p4) anc2(p2 p3 p4) failure _d: event analysis time _t: exit enter on or after: time entry Fitting constant-only model: Iteration 0: log likelihood = -2740.1508 (not concave) Iteration 1: log likelihood = -2645.1735 (not concave) Iteration 2: log likelihood = -2400.6906 (not concave) Iteration 3: log likelihood = -2302.198 (not concave) Iteration 4: log likelihood = -2276.2273 (not concave) Iteration 5: log likelihood = -2270.3933 Iteration 6: log likelihood = -2262.6202 Iteration 7: log likelihood = -2261.6685 Iteration 8: log likelihood = -2261.6596 Iteration 9: log likelihood = -2261.6596 Fitting full model: Iteration 0: log likelihood = -2261.6596 Iteration 1: log likelihood = -2259.0287 (not concave) Iteration 2: log likelihood = -2255.2842 (not concave) Iteration 3: log likelihood = -2254.4294 Iteration 4: log likelihood = -2252.4243 Iteration 5: log likelihood = -2251.2608 Iteration 6: log likelihood = -2250.9031 Iteration 7: log likelihood = -2250.8108 Iteration 8: log likelihood = -2250.7673 Iteration 9: log likelihood = -2250.7658 Iteration 10: log likelihood = -2250.7658 Gamma regression -- accelerated failure-time form No. of subjects = 2314 Number of obs = 2314 No. of failures = 1053 Time at risk = 4478.967999 LR chi2(3) = 21.79 Log likelihood = -2250.7658 Prob > chi2 = 0.0001 ------------------------------------------------------------------------------ | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- _t | p2 | -.0244036 .0936539 -0.26 0.794 -.2079619 .1591546 p3 | 1.771969 .2573728 6.88 0.000 1.267528 2.27641 p4 | 2.397351 .4960183 4.83 0.000 1.425173 3.369529 _cons | .6736713 .0716252 9.41 0.000 .5332885 .814054 -------------+---------------------------------------------------------------- ln_sig | p2 | .0917806 .1048627 0.88 0.381 -.1137466 .2973078 p3 | .9545557 .4117426 2.32 0.020 .147555 1.761556 p4 | .5767601 .7730567 0.75 0.456 -.9384032 2.091923 _cons | -.2572982 .0878631 -2.93 0.003 -.4295067 -.0850896 -------------+---------------------------------------------------------------- kappa | p2 | -.5178993 .2276772 -2.27 0.023 -.9641384 -.0716602 p3 | -.8198126 .7261392 -1.13 0.259 -2.243019 .6033941 p4 | -.5240963 1.926999 -0.27 0.786 -4.300945 3.252752 _cons | 1.368894 .1890939 7.24 0.000 .998277 1.739511 ------------------------------------------------------------------------------ . exit, clear