*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\fulldata.dat
stset exit, failure(event) enter(entry)
streg p1 p2 p3 p4, distribution(gamma) noconstant
*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