/*======================================================================
 Duration Data Models: file = crime1.dat
 Unemployment spell in months
*/======================================================================

Read  ; Nobs = 2725 
	; Nvar = 16 
	; file = c:\econm\crime1.dat
	; Names = 2 $

/*
narr86    nfarr86   nparr86   pcnv      avgsen    tottime   ptime86   qemp86   
inc86     durat     black     hispan    born60    pcnvsq    pt86sq    inc86sq  
*/

Reject ; durat < 7 $   Now only 308 observations remain.

Create ; LogT = Log(durat) $

/*=================================================================
1. Nonparametric Life Table (no covariates)
*/=================================================================

SURVIVAL  ; lhs = durat 
	    ; plot $

/*=================================================================
2. (log-linear) Cox's Proportional Hazard Model (with Covariates)
*/=================================================================

SURVIVAL  ; lhs = durat 
          ; rhs = inc86, black, hispan ; plot $

/*=================================================================
3. Log-linear Parametric Hazard Models
*/=================================================================

SURVIVAL  ; lhs = LogT 
	    ; Model = Weibull
          ; rhs = One, inc86, black, hispan  $

SURVIVAL  ; lhs = LogT 
	    ; Model = Exponential 
          ; rhs = One, inc86, black, hispan  $

SURVIVAL  ; lhs = LogT 
	    ; Model = normal
          ; rhs = One, inc86, black, hispan  $

SURVIVAL  ; lhs = LogT 
	    ; Model = logistic
          ; rhs = One, inc86, black, hispan  $