/* This program computes a Weibull duration model with covariate on the
   duration of strikes.  The covariate is unanticipated industrial production.
   See William H. Greene, Third ed., Econometric Analysis, example 20.17,
   p. 990.  */


options pagesize=60 linesize=74 nodate;

data strike;
  input year x dur;
cards;
1968 0.01138 7
1968 0.01138 9
1968 0.01138 13
1968 0.01138 14
1968 0.01138 26
1968 0.01138 29
1968 0.01138 52
1968 0.01138 130
1969 0.02299 9
1969 0.02299 37
1969 0.02299 41
1969 0.02299 49
1969 0.02299 52
1969 0.02299 119
1970 -0.03957 3
1970 -0.03957 17
1970 -0.03957 19
1970 -0.03957 28
1970 -0.03957 72
1970 -0.03957 99
1970 -0.03957 104
1970 -0.03957 114
1970 -0.03957 152
1970 -0.03957 153
1970 -0.03957 216
1971 -0.05467 15
1971 -0.05467 61
1971 -0.05467 98
1972 0.00535 2
1972 0.00535 25
1972 0.00535 85
1973 0.07427 3
1973 0.07427 10
1974 0.06450 1
1974 0.06450 2
1974 0.06450 3
1974 0.06450 3
1974 0.06450 3
1974 0.06450 4
1974 0.06450 8
1974 0.06450 11
1974 0.06450 22
1974 0.06450 23
1974 0.06450 27
1974 0.06450 32
1974 0.06450 33
1974 0.06450 35
1974 0.06450 43
1974 0.06450 43
1974 0.06450 44
1974 0.06450 100
1975 -0.10443 5
1975 -0.10443 49
1976 -0.00700 2
1976 -0.00700 12
1976 -0.00700 12
1976 -0.00700 21
1976 -0.00700 21
1976 -0.00700 27
1976 -0.00700 38
1976 -0.00700 42
1976 -0.00700 117
;
/*
proc lifereg data=strike;
  model dur = / dist = weibull;
  model dur = / dist = exponential;
  model dur = / dist = lnormal;
  model dur = / dist = llogistic;
*/
proc lifereg data=strike;
  model dur =  x / dist = weibull;

run;