Sampling Plans and Initial Condition Problems For Continuous Time - - PowerPoint PPT Presentation

sampling plans and initial condition problems for
SMART_READER_LITE
LIVE PREVIEW

Sampling Plans and Initial Condition Problems For Continuous Time - - PowerPoint PPT Presentation

Sampling Plans and Initial Condition Problems For Continuous Time Duration Models James J. Heckman University of Chicago Econ 312, Spring 2019 Heckman Sampling Plans Sampling Plans and Initial Condition Problems For Continuous Time Duration


slide-1
SLIDE 1

Sampling Plans and Initial Condition Problems For Continuous Time Duration Models

James J. Heckman University of Chicago Econ 312, Spring 2019

Heckman Sampling Plans

slide-2
SLIDE 2

Sampling Plans and Initial Condition Problems For Continuous Time Duration Models

James J. Heckman University of Chicago Econ 312, Spring 2019

Heckman Sampling Plans

slide-3
SLIDE 3

Sampling plans and initial condition problems: Duration Models Consider a random sample of unemployment spells in progress. For sampled spells, one of the following duration times may be

  • bserved:
  • time in state up to sampling date (Tb) (recall of time spent)
  • time in state after sampling date (Ta) (prospective sampling

forward)

  • total time in completed spell observed at origin of sample

(Tc = Ta + Tb)

Heckman Sampling Plans

slide-4
SLIDE 4

Duration of spells beginning after the origin date of the sample, denoted Td, are not subject to initial condition problems. The intake rate at time −tb (assuming sample occurs at time 0: the proportion of the population entering a spell at −tb. Assume:

  • A time homogenous environment, i.e. constant intake rate,

k(−tb) = k, ∀b

  • A model without observed or unobserved explanatory variables.
  • No right censoring, so Tc = Ta + Tb
  • Underlying distribution f (x) is nondefective
  • m =

0 (x)dx < ∞

Heckman Sampling Plans

slide-5
SLIDE 5

The proportion of the population experiencing a spell at t = 0, the

  • rigin date of the sample, is

P0 = ∞ k(−tb)(1 − F(tb))dtb = k ∞ (1 − F(tb))dtb = k

  • tb(1 − F(tb))|∞

0 −

∞ tbd(1 − F(tb))

  • =

k ∞ tbf (tb)dtb = km where 1 − F(tb) is the probability the spell lasts from −tb to 0 (or equivalently, from 0 to −tb).

Heckman Sampling Plans

slide-6
SLIDE 6

So the density of a spell of length tb interrupted at the beginning of the sample (t = 0) is g(tb) = proportion surviving til t = 0 from batch tb total surviving til t = 0 = k(−tb)(1 − F(tb)) P0 = 1 − F(tb) m = f (tb) Notice: g is the distribution of Tb in the population constructed by sampling rule of source population. Distinguish from F : cdf of the true population. G : cdf of the sampled spells.

Heckman Sampling Plans

slide-7
SLIDE 7

The probability that a spell lasts until tc given that it has lasted from −tb to 0, is the conditioned density of tc given 0 < tb < tc. f (tc|t > tb > 0) = f (tc) 1 − F(tb); tc ≥ tb ≥ 0 So the density of a spell in the sampled population that lasts, tc is g(tc) = tc f (tc|t ≥ tb)f (t ≥ tb)dtb = tc f (tc) m dtb = f (tc)tc m

Heckman Sampling Plans

slide-8
SLIDE 8

Likewise, the density of a sampled spell that lasts until ta is g(ta) = ∞ f (ta + tb|tb)Pr(t ≥ ta ≥ 0))dtb = ∞ f (ta + tb) m dtb = 1 m ∞

ta

f (tb)dtb = 1 − F(ta) m (Stationarity, mirror images have some densities). So the functional form of f (tb) = f (ta): Consequences of stationarity.

Heckman Sampling Plans

slide-9
SLIDE 9

Some useful results that follow from this model:

1 If f (t) = θe−tθ, then g(tb) = θe−tbθ and g(ta) = θe−taθ.

Proof: f (t) = θe−tθ → m = 1 θ, F(t) = 1 − e−tθ → g(ta) = 1 − F(t) m = θe−tθ

Heckman Sampling Plans

slide-10
SLIDE 10

1 E(Ta) = m

2 (1 + σ2 m2). Proof: E(Ta) =

  • taf (ta)dta =
  • ta

1 − G(ta) m dta = 1 m 1 2t2

a(1 − F(ta))|∞ 0 −

1 2t2

ad(1 − F(ta))

  • =

1 m 1 2t2

aF(ta)dta = 1

2m[var(ta) + E 2(ta)] = 1 2m[σ2 + m2]

Heckman Sampling Plans

slide-11
SLIDE 11

1 E(Tb) = m

2 (1 + σ2 m2). Proof: See proof of Proposition 2.

2 E(Tc) = m(1 + σ2

m2). Proof: E(Tc) = t2

c F(tc)

m dtc = 1 m(var(tc) + E 2(tc)) → E(Tc) = 2E(Ta) = 2E(Tb), E(Tc) > m unless σ2 = 0

Heckman Sampling Plans

slide-12
SLIDE 12

Some Additional Results: h(t) = hazard : h(t) = f (t) 1 − F(t).

1 h′(t) > 0 → E(Ta) = E(Tb) < m. Proof: See Barlow and

Proschan.

2 h′(t) < 0 → E(Ta) = E(Tb) > m. Proof: See Barlow and

Proschan.

Heckman Sampling Plans

slide-13
SLIDE 13

Examples

Heckman Sampling Plans

slide-14
SLIDE 14

Specification of the Distribution

Weibull Distribution

  • Parameters: λ > 0, k > 0
  • Probability Density Function (PDF):

λ k t λ k−1 exp

t k k

  • Cumulative Density Function:

1 − exp

t k k

  • Set of Parameters: 

   λ1, k1 = 0.5 λ2, k1 = 1.0 λ3, k1 = 2.0 λ3, k1 = 3.0     , respectively

Heckman Sampling Plans

slide-15
SLIDE 15

Basic Distribution Graphs

& &

0.5 1 1.5 2 0.5 1 1.5 2 2.5 3 t P D F

  • f

t h e S p e l l s : W e i b u l l D i s t r i b u t i

  • n

s Weibull Distribution λ = 0.1, k = 0.5 Weibull Distribution λ = 0.5, k = 1.0 Weibull Distribution λ = 0.5, k = 2.0 Weibull Distribution λ = 1.0, k = 3.0 0.5 1 1.5 2 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 t C D F

  • f

t h e D i s t r i b u t i

  • n

: W e i b u l l Weibull Distribution λ = 0.1, k = 0.5 Weibull Distribution λ = 0.5, k = 1.0 Weibull Distribution λ = 0.5, k = 2.0 Weibull Distribution λ = 1.0, k = 3.0

Heckman Sampling Plans

slide-16
SLIDE 16

Basic Duration Graphs

& !" &

0.5 1 1.5 2 1 2 3 4 5 6 7 8 9 10 t Hazard Function of the Distribution: Weibull Weibull Distribution λ = 0.1, k = 0.5 Weibull Distribution λ = 0.5, k = 1.0 Weibull Distribution λ = 0.5, k = 2.0 Weibull Distribution λ = 1.0, k = 3.0 0.5 1 1.5 2 1 2 3 4 5 6 7 8 9 10 t Integrated Hazard Function of the Distribution: Weibull Weibull Distribution λ = 0.1, k = 0.5 Weibull Distribution λ = 0.5, k = 1.0 Weibull Distribution λ = 0.5, k = 2.0 Weibull Distribution λ = 1.0, k = 3.0

Heckman Sampling Plans

slide-17
SLIDE 17

Observed and Original Distribution for Tb (Example 1)

0.5 1 1.5 2 0.5 1 1.5 2 2.5 3 t Observed (T b) and Original PDFs of the Spells The Observed PDF of Spells (T

b)

The Original PDF (Weibull Distribution λ = 0.1, k = 0.5)

Heckman Sampling Plans

slide-18
SLIDE 18

Observed and Original Distribution for Tb (Example 3)

0.5 1 1.5 0.5 1 1.5 2 2.5 t Observed (T b) and Original PDFs of the Spells The Observed PDF of Spells (T

b)

The Original PDF (Weibull Distribution λ = 0.5, k = 2.0)

Heckman Sampling Plans

slide-19
SLIDE 19

Observed and Original Distribution for Tb (Example 4)

0.5 1 1.5 2 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 t Observed (T b) and Original PDFs of the Spells The Observed PDF of Spells (T

b)

The Original PDF (Weibull Distribution λ = 1.0, k = 3.0)

Heckman Sampling Plans

slide-20
SLIDE 20

Observed and Original Distribution for Tc (Example 1)

0.5 1 1.5 2 0.5 1 1.5 2 2.5 3 t O b s e r v e d ( T c ) a n d O r i g i n a l P D F s

  • f

t h e S p e l l s The Observed PDF of Spells (T

c)

The Original PDF (Weibull Distribution λ = 0.1, k = 0.5)

Heckman Sampling Plans

slide-21
SLIDE 21

Observed and Original Distribution for Tc (Example 2)

0.5 1 1.5 2 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 t O b s e r v e d ( T c ) a n d O r i g i n a l P D F s

  • f

t h e S p e l l s The Observed PDF of Spells (T

c)

The Original PDF (Weibull Distribution λ = 0.5, k = 1.0)

  • Heckman

Sampling Plans

slide-22
SLIDE 22

Observed and Original Distribution for Tc (Example 3)

0.5 1 1.5 0.5 1 1.5 2 2.5 t Observed (T c) and Original PDFs of the Spells The Observed PDF of Spells (T

c)

The Original PDF (Weibull Distribution λ = 0.5, k = 2.0)

  • Heckman

Sampling Plans

slide-23
SLIDE 23

Observed and Original Distribution for Tc (Example 4)

0.5 1 1.5 2 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 t Observed (T c) and Original PDFs of the Spells The Observed PDF of Spells (T

c)

The Original PDF (Weibull Distribution λ = 1.0, k = 3.0)

Heckman Sampling Plans