*change name of path to read in fulldata.dat
> infile pubid entry exit event p1 p2 p3 p4 per1234 using E:\GeneralGamma\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(gamma) noconstant                                          
                                                                                             
         failure _d:  event                                                                  
   analysis time _t:  exit                                                                   
  enter on or after:  time entry                                                             
Fitting full model:                                                                          
                                                                                             
Iteration 0:   log likelihood = -3146.5895  (not concave)                                    
Iteration 1:   log likelihood =  -2464.221  (not concave)                                    
Iteration 2:   log likelihood = -2355.2306                                                   
Iteration 3:   log likelihood = -2277.4682                                                   
Iteration 4:   log likelihood = -2269.6253                                                   
Iteration 5:   log likelihood = -2269.3496                                                   
Iteration 6:   log likelihood = -2269.3491                                                   
Iteration 7:   log likelihood = -2269.3491                                                   
                                                                                             
Gamma 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)    =   1356.69               
Log likelihood  =   -2269.3491                     Prob > chi2     =    0.0000               
                                                                                             
------------------------------------------------------------------------------               
          _t |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]               
-------------+----------------------------------------------------------------               
          p1 |    .642839   .0552769    11.63   0.000     .5344982    .7511797               
          p2 |   .7937214   .0528871    15.01   0.000     .6900646    .8973783               
          p3 |   1.961447    .096863    20.25   0.000     1.771599    2.151295               
          p4 |   2.914392    .090706    32.13   0.000     2.736612    3.092173               
-------------+----------------------------------------------------------------               
     /ln_sig |   -.217642   .0524492    -4.15   0.000    -.3204405   -.1148435               
      /kappa |   1.270146    .109724    11.58   0.000     1.055091    1.485201               
-------------+----------------------------------------------------------------               
       sigma |   .8044134   .0421908                      .7258292    .8915056               
------------------------------------------------------------------------------               
                  
*the command will fit a six-parameter conventional general gamma model to the full
*data for all four periods producing estimated betas for p1-p4, sigma (1/p), lambda (kappa)
*and log likelihood                                                                           

. exit, clear