Stat Computing > Stata > Textbook Examples > Applied Longitudinal Data Analysis

Stata Textbook Examples
Applied Longitudinal Data Analysis: Modeling Change and Event Occurrence
by Judith D. Singer and John B. Willett
Chapter 12:  Extending the Discrete-time Hazard Model

Table 12.2, page 413.

use http://www.ats.ucla.edu/stat/stata/examples/alda/data/tenure_pp, clear

/* compute powers of period for polynomial models */
generate p2 = period^2
generate p3 = period^3
generate p4 = period^4
generate p5 = period^5
 
quietly tabulate id
global nid = r(r)
  
/* constant model */ 
quietly logit event
 
/* display deviance aic bic */
display -2*e(ll) "  " -2*(e(ll)-(e(df_m)+1)) "   " -2*(e(ll)-(log($nid)/2*(e(df_m)+1)))
1037.5652  1039.5652   1043.1259
 
/* linear model */
quietly logit event period
display -2*e(ll) "  " -2*(e(ll)-(e(df_m)+1)) "   " -2*(e(ll)-(log($nid)/2*(e(df_m)+1)))
867.46187  871.46187   875.02255
 
/* quadratic model */
quietly logit event period p2
display -2*e(ll) "  " -2*(e(ll)-(e(df_m)+1)) "   " -2*(e(ll)-(log($nid)/2*(e(df_m)+1)))
836.30407  842.30407   852.98612
 
/* cubic model */
quietly logit event period p2 p3
display -2*e(ll) "  " -2*(e(ll)-(e(df_m)+1)) "   " -2*(e(ll)-(log($nid)/2*(e(df_m)+1)))
833.17249  841.17249   855.41521
 
/* fourth order model */
quietly logit event period p2 p3 p4
display -2*e(ll) "  " -2*(e(ll)-(e(df_m)+1)) "   " -2*(e(ll)-(log($nid)/2*(e(df_m)+1)))
832.74267  842.74267   860.54608
 
quietly logit event period p2 p3 p4 p5
display -2*e(ll) "  " -2*(e(ll)-(e(df_m)+1)) "   " -2*(e(ll)-(log($nid)/2*(e(df_m)+1)))
832.73216  844.73216   866.09625
 
/* general model */
quietly logit event period d1-d9
display -2*e(ll) "  " -2*(e(ll)-(e(df_m)+1)) "   " -2*(e(ll)-(log($nid)/2*(e(df_m)+1)))
831.2038  849.2038   881.24993

UCLA Researchers are invited to our Statistical Consulting Services
We recommend others to our list of Other Resources for Statistical Computing Help
These pages are Copyrighted (c) by UCLA Academic Technology Services