The model: Person 1 (birth occurs in the first interval): Which leads to:

1. Single and Multiple Spell Discrete Time Hazards Models with Parametric and Non-Parametric Corrections for Unobserved HeterogeneityDavid K. Guilkey

2. Demographic Applications:Single Spell1. Time until death2. Time until retirement3. Time until first marriage4. Time until first birthMultiple Spell1. Time until birth of each child2. Duration of each spell of employmentWe will use time until first birth and the timing of subsequent births as an example throughout the presentation.

3. The variable of interest is: P(t ≤ T < t+n | T > t)This is the conditional probability that an individual experiences the event between t and t+n given that she has not experienced the event until that time.Example:The dependent variable is the timing of a first birth. Suppose the discrete time interval is a year and we observe each woman from the beginning of her child bearing years:0…..1…..2…..3Consider three cases:Person 1: Has a birth in year 1 (time 0 may be age 12)Person 2: Has a birth in year 2Person 3: Still has not had a birth at the end of the observation period

4. Some important notes:1. Since we are following the woman from the beginning of her child bearing years, we have eliminated the possibility of left censoring (the event occurs before the observation period).2. Left censoring combined with unobserved heterogeneity introduces bias into the estimation results. The correction requires the estimation of an “initial conditions” equation similar to Heckman selection equation which are well known to yield unstable parameter estimates.3. The third person is right censored. However, right censoring is easily handled as part of the estimation process.4. As will be seen below, the dependent variable in a discrete time hazard model is dichotomous. Can use probit, logit or complementary log log (cloglog) models. Logit and cloglog are most often used. I use logit since one of the software packages needs logit – results were nearly the same for cloglog in models where software allowed for both (STATA).

5. The model:Person 1 (birth occurs in the first interval):Which leads to:

6. Person 2 (No birth in the first year and a birth in the second year):Joint probability is:

7. Person 3: (No births in the observation period)Estimation:Time 1: 3 observationsTime 2: 2 observationsTime 3: 1 observationThe three sets of coefficients could be estimated in three separate logits for the set of individuals at risk. This is true since there is no unobserved heterogeneity that links the three time periods together.

8. Duration dependenceThis is a concept similar to state dependence in a standard panel data model.Duration dependence occurs when the value of the hazard at any point in time depends on the amount of time that has already elapsed.Relates to the propensity of a state towards self-perpetuationExamples:Mortality – hazard increases with time regardless of the values of the other covariatesUnemployment duration – hazard of finding employment may decrease as the length of the unemployment spell increases

9. Modeling Duration DependenceIn our current model, duration dependence is captured by the intercept terms in the equation since they are allowed to differ at each point in time.To see more clearly, assume that the effects of the covariates is the same at each point in time (the β’s are the same in the previous equations).Now define T1ti=1 if if t=1 and 0 otherwise – with T2ti and T3ti defined similarlyThen we can write (no constant in the model):Which allows for a very flexible pattern of duration dependence – can be non-linear for example

10. A less flexible pattern that requires the estimation of fewer parameters is:In our example, we will be examining the birth hazard starting all women at age 10 and so age and duration dependence are not separately identified.A parametric model which allows for non-linear duration dependence is:

11. Empirical ExampleData from Indonesia Family Life Survey.We first examine timing of first birth – women followed from age 10 until first birth.Data set up:

12. Simple Models (linear and non-linear duration dependence):

13. Non-parametric Duration Dependence (using duration or age dummies):

14. Duration Dependence and Unobserved HeterogeneityReview of dynamic panel data model:where we have a time varying error and a persistent error (sometimes referred to as time invariant unobserved heterogeneity)Define:Then:

15. Alternative model (state dependence):where |α|<1Now:It is very difficult to distinguish between the models – so we use the hybrid model:A problem is that this model is more difficult to estimate – neither ordinary least squares nor fixed effects methods yield consistent estimators – use maximum likelihood (with initial conditions problem) or instrumental variables.

16. Return to first example and unobserved heterogeneity (using person 2 as the example):Person 2 (No birth in the first year and a birth in the second year):We can no longer estimate parameters time period by time period – due to selection on unobservables (just as in standard Heckman selectivity model)Joint probability is now:

17. The unconditional joint probability is:Most commonly used distributional assumption for the unobserved heterogeneity is the normal distribution. The integral is approximated using Hermite point and weights (simply looked up in a table for the normal distribution):K is the number of interpolation points – more accurate to add more but slower (STATA default is 12 – frequently not enough for rare events)Heckman-Singer approach: Do not assume a distribution – directly estimate the points and weights as part of the maximum likelihood estimation process – referred to as the discrete factor approximation.

18. Identification “it is somewhat heroic to think that we can distinguish between duration dependence and unobserved heterogeneity when we only observe a single cycle for each agent” (Wooldridge – page 705)Example:Model with no censoring estimated by OLS.Can identify both using functional form – but the model parameter estimates are frequently unstable.

19. ExamplesAssume normality:

20. Cannot directly compare the coefficients with and without heterogeneity correction because of possible scale differences for discrete dependent variable models. However, scale effects can be removed if you compare ratios of coefficients:Without unobserved heterogeneity:With Unobserved heterogeneity:

21. Use Discrete Factor Method

22. Multiple Spell Discrete Time Hazards ModelsModel with no unobserved heterogeneity:Allow for M births:With no heterogeneity, estimate M+1 single spell hazards models (or fully interacted model).Results for fully interacted model (m=0,1,2,3,4):

24. Simple Model with coefficients restricted to be the same (using all available births for all women):

25. Add unobserved heterogeneity to the model:In order to use STATA, must assume a restrictive form of unobserved heterogeneity for both parametric and non-parametric forms.Parametric:More flexible specification would be:where Σ is m x m

26. Estimate assuming normally distributed unobserved heterogeneity (restrict coefficients across births):

27. Use the discrete factor model (restrict coefficients across births):

28. Normally distributed unobserved heterogeneity (unrestricted coefficients):

29. Discrete factor model with two points of support (unrestricted coefficients):

30. Add non-parametric unobserved heterogeneity with three points of support – unrestricted across equations using fortran:

31. Continued:

32. Continued:

33. Can use likelihood ratio test to compare model without heterogeneity to:1. Discrete factor model with two points of support using STATA where we have a restricted form of heterogeneity2. Discrete factor model with three points of support and unrestricted heterogeneity using fortran.Tests sequentially reject the simpler models with p levels close to zero.