*change name of path to read in fulldata.dat . infile pubid entry exit event p1 p2 p3 p4 per1234 using C:\mikefolder\mikeFolder\survivalme > thods\programsforweb\Figures\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 . streg p1 p2 p3 p4, distribution(weibull) noconstant nohr time failure _d: event analysis time _t: exit enter on or after: time entry Fitting full model: Iteration 0: log likelihood = -2687.2908 (not concave) Iteration 1: log likelihood = -2338.5997 Iteration 2: log likelihood = -2274.849 Iteration 3: log likelihood = -2272.8382 Iteration 4: log likelihood = -2272.7926 Iteration 5: log likelihood = -2272.7926 Weibull regression -- accelerated failure-time form No. of subjects = 2314 Number of obs = 2314 No. of failures = 1053 Time at risk = 4478.967999 Wald chi2(4) = 1824.55 Log likelihood = -2272.7926 Prob > chi2 = 0.0000 ------------------------------------------------------------------------------ _t | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- p1 | .553122 .0458155 12.07 0.000 .4633252 .6429187 p2 | .7069708 .0427269 16.55 0.000 .6232276 .790714 p3 | 1.888003 .0897903 21.03 0.000 1.712018 2.063989 p4 | 2.862724 .0873612 32.77 0.000 2.691499 3.033949 -------------+---------------------------------------------------------------- /ln_p | .1048529 .0270345 3.88 0.000 .0518662 .1578396 -------------+---------------------------------------------------------------- p | 1.110547 .0300231 1.053235 1.170978 1/p | .900457 .0243434 .8539867 .9494559 ------------------------------------------------------------------------------ *the command will fit a conventional weibull model to the full data for all *four periods providing estimated betas for p1-p4, sigma (1/p), and log likelihood . exit, clear