-------------
       log:  E:\intdemdata\JOP Aug 01 replication.smcl
  log type:  smcl
 opened on:  21 Mar 2003, 15:28:34

. *This log file replicates the results reported in Dan Reiter, "Does Peace Nurture Democracy?"
. *Journal of Politics 63 (August) 2001, 935-948.  Analysis conducted using STATA 7.0.
. *Note that the results produced in Table 1 from p.943 are produced with a data set
. *different from the data set used to produce the results in Table 2 on p.945.
. set mem 70m
(71680k)

. *First, I reproduce the Table 1 results, which are the factors which hasten the transition of
. *democracies.
. use "E:\intdemdata\dtnew02.dta", clear

. stweib midtot loswar dregfrc6 lnrgdp growth islam if year>1959 & year<1993, nohr robust cluster(panid)

         failure _d:  dtrans
   analysis time _t:  dtctr
                 id:  panid


Fitting constant-only model:

Iteration 0:   log likelihood = -143.78171
Iteration 1:   log likelihood = -131.50681
Iteration 2:   log likelihood = -131.00303
Iteration 3:   log likelihood = -131.00139
Iteration 4:   log likelihood = -131.00139

Fitting full model:

Iteration 0:   log likelihood = -131.00139
Iteration 1:   log likelihood = -113.23526
Iteration 2:   log likelihood =  -109.3401
Iteration 3:   log likelihood = -109.29224
Iteration 4:   log likelihood = -109.29214
Iteration 5:   log likelihood = -109.29214

Weibull regression -- log relative-hazard form

No. of subjects =          167                     Number of obs   =      2564
No. of failures =           39
Time at risk    =         2564
                                                   Wald chi2(6)    =     48.52
Log likelihood  =   -109.29214                     Prob > chi2     =    0.0000

                            (standard errors adjusted for clustering on panid)
-------------+----------------------------------------------------------------
             |               Robust
          _t |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
      midtot |  -.2690823   .2049005    -1.31   0.189    -.6706799    .1325153
      loswar |   .3817912   1.548391     0.25   0.805    -2.652999    3.416581
    dregfrc6 |   4.812115   1.117968     4.30   0.000     2.620937    7.003293
      lnrgdp |   .4684231   .2147294     2.18   0.029     .0475612    .8892849
      growth |  -.0254716   .2172975    -0.12   0.907     -.451367    .4004237
       islam |  -.0004386   .0071002    -0.06   0.951    -.0143547    .0134775
       _cons |  -7.134767   1.542086    -4.63   0.000     -10.1572   -4.112334
-------------+----------------------------------------------------------------
       /ln_p |   -.557496   .1049813    -5.31   0.000    -.7632556   -.3517363
-------------+----------------------------------------------------------------
           p |   .5726412   .0601166                      .4661464    .7034656
         1/p |   1.746294   .1833283                      1.421534    2.145249
-------------+----------------------------------------------------------------

. *That replicates Model 1 of Table 1.
. stweib pbit loswar dregfrc6 lnrgdp growth islam if year>1959 & year<1989, nohr robust cluster(ccode)

         failure _d:  dtrans
   analysis time _t:  dtctr
                 id:  panid


Fitting constant-only model:

Iteration 0:   log likelihood = -97.899623
Iteration 1:   log likelihood = -88.703125
Iteration 2:   log likelihood = -88.317002
Iteration 3:   log likelihood = -88.314911
Iteration 4:   log likelihood = -88.314911

Fitting full model:

Iteration 0:   log likelihood = -88.314911
Iteration 1:   log likelihood =  -78.96665
Iteration 2:   log likelihood = -76.335099
Iteration 3:   log likelihood = -76.262878
Iteration 4:   log likelihood = -76.262731
Iteration 5:   log likelihood = -76.262731

Weibull regression -- log relative-hazard form

No. of subjects =          141                     Number of obs   =      2274
No. of failures =           23
Time at risk    =         2274
                                                   Wald chi2(6)    =     43.60
Log likelihood  =   -76.262731                     Prob > chi2     =    0.0000

                            (standard errors adjusted for clustering on ccode)
-------------+----------------------------------------------------------------
             |               Robust
          _t |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
        pbit |  -.0630291   .2075476    -0.30   0.761    -.4698149    .3437566
      loswar |   1.945241    1.18155     1.65   0.100    -.3705542    4.261037
    dregfrc6 |   5.744824   2.322617     2.47   0.013     1.192578    10.29707
      lnrgdp |   .6558053   .2402742     2.73   0.006     .1848765    1.126734
      growth |  -.0331497   .4990038    -0.07   0.947    -1.011179    .9448797
       islam |   .0056607   .0096413     0.59   0.557    -.0132359    .0245573
       _cons |  -8.612328   1.536652    -5.60   0.000    -11.62411   -5.600546
-------------+----------------------------------------------------------------
       /ln_p |  -.9418386   .0963077    -9.78   0.000    -1.130598    -.753079
-------------+----------------------------------------------------------------
           p |   .3899103   .0375514                      .3228401    .4709144
         1/p |   2.564693   .2469996                      2.123528    3.097509
-------------+----------------------------------------------------------------

. *That replicates Model 2 of Table 2.
. stweib pbitot loswar dregfrc6 lnrgdp growth islam if year>1959 & year<1989, nohr robust cluster(ccode)

         failure _d:  dtrans
   analysis time _t:  dtctr
                 id:  panid


Fitting constant-only model:

Iteration 0:   log likelihood = -97.899623
Iteration 1:   log likelihood = -88.703125
Iteration 2:   log likelihood = -88.317002
Iteration 3:   log likelihood = -88.314911
Iteration 4:   log likelihood = -88.314911

Fitting full model:

Iteration 0:   log likelihood = -88.314911
Iteration 1:   log likelihood = -78.972983
Iteration 2:   log likelihood = -76.359932
Iteration 3:   log likelihood = -76.288325
Iteration 4:   log likelihood = -76.288178
Iteration 5:   log likelihood = -76.288178

Weibull regression -- log relative-hazard form

No. of subjects =          141                     Number of obs   =      2274
No. of failures =           23
Time at risk    =         2274
                                                   Wald chi2(6)    =     42.31
Log likelihood  =   -76.288178                     Prob > chi2     =    0.0000

                            (standard errors adjusted for clustering on ccode)
-------------+----------------------------------------------------------------
             |               Robust
          _t |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
      pbitot |  -.0284425   .1483874    -0.19   0.848    -.3192766    .2623915
      loswar |   1.916897   1.309268     1.46   0.143    -.6492209    4.483015
    dregfrc6 |   5.733437   2.319435     2.47   0.013     1.187427    10.27945
      lnrgdp |   .6634493   .2457706     2.70   0.007     .1817477    1.145151
      growth |  -.0337533   .4957467    -0.07   0.946    -1.005399    .9378925
       islam |   .0055676   .0093721     0.59   0.552    -.0128014    .0239366
       _cons |  -8.679739   1.549911    -5.60   0.000    -11.71751   -5.641968
-------------+----------------------------------------------------------------
       /ln_p |  -.9358022   .0946136    -9.89   0.000    -1.121241    -.750363
-------------+----------------------------------------------------------------
           p |    .392271   .0371142                       .325875    .4721951
         1/p |   2.549258   .2411945                      2.117769    3.068662
-------------+----------------------------------------------------------------

. *That replicates Model 3 of Table 1.
. stweib iw dregfrc6 lnrgdp growth islam if year>1959 & year<1993, nohr robust cluster(ccode)

         failure _d:  dtrans
   analysis time _t:  dtctr
                 id:  panid


Fitting constant-only model:

Iteration 0:   log likelihood = -143.78171
Iteration 1:   log likelihood = -131.50681
Iteration 2:   log likelihood = -131.00303
Iteration 3:   log likelihood = -131.00139
Iteration 4:   log likelihood = -131.00139

Fitting full model:

Iteration 0:   log likelihood = -131.00139
Iteration 1:   log likelihood = -112.65672
Iteration 2:   log likelihood = -108.27519
Iteration 3:   log likelihood = -107.87288
Iteration 4:   log likelihood = -107.80652
Iteration 5:   log likelihood = -107.79056
Iteration 6:   log likelihood = -107.78696
Iteration 7:   log likelihood = -107.78618
Iteration 8:   log likelihood = -107.78604
Iteration 9:   log likelihood = -107.78603
Iteration 10:  log likelihood = -107.78603

Weibull regression -- log relative-hazard form

No. of subjects =          167                     Number of obs   =      2564
No. of failures =           39
Time at risk    =         2564
                                                   Wald chi2(5)    =   1299.80
Log likelihood  =   -107.78603                     Prob > chi2     =    0.0000

                            (standard errors adjusted for clustering on ccode)
-------------+----------------------------------------------------------------
             |               Robust
          _t |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
          iw |   -15.2468    .534836   -28.51   0.000    -16.29506   -14.19854
    dregfrc6 |   5.263675   1.084303     4.85   0.000      3.13848     7.38887
      lnrgdp |   .4528143   .1952429     2.32   0.020     .0701453    .8354833
      growth |  -.0111657   .2181786    -0.05   0.959    -.4387878    .4164565
       islam |   .0020841   .0073036     0.29   0.775    -.0122307     .016399
       _cons |  -7.205606   1.379117    -5.22   0.000    -9.908626   -4.502586
-------------+----------------------------------------------------------------
       /ln_p |  -.5761715   .1022734    -5.63   0.000    -.7766237   -.3757193
-------------+----------------------------------------------------------------
           p |   .5620461   .0574824                      .4599563    .6867951
         1/p |   1.779214   .1819663                      1.456038    2.174119
-------------+----------------------------------------------------------------

. *That replicates Model 4 of Table 1.
. *Now, on to Table 2, which uses a different data set to analyze the factors which
. *hasten democratic breakdown.
. clear

. use "E:\intdemdata\dcnew01.dta", clear

. stweib  midtot  dregfrc6 lnrgdp growth islam if year>1959 & year<1993, robust cluster(ccode) nohr

         failure _d:  autrev
   analysis time _t:  dsctr
                 id:  panelid


Fitting constant-only model:

Iteration 0:   log likelihood = -54.304906
Iteration 1:   log likelihood = -51.384946
Iteration 2:   log likelihood = -51.339033
Iteration 3:   log likelihood = -51.339021
Iteration 4:   log likelihood = -51.339021

Fitting full model:

Iteration 0:   log likelihood = -51.339021
Iteration 1:   log likelihood = -36.201247
Iteration 2:   log likelihood = -31.677484
Iteration 3:   log likelihood = -31.332918
Iteration 4:   log likelihood = -31.332092
Iteration 5:   log likelihood = -31.332092

Weibull regression -- log relative-hazard form

No. of subjects =           82                     Number of obs   =      1277
No. of failures =           16
Time at risk    =         1277
                                                   Wald chi2(5)    =     48.38
Log likelihood  =   -31.332092                     Prob > chi2     =    0.0000

                            (standard errors adjusted for clustering on ccode)
-------------+----------------------------------------------------------------
             |               Robust
          _t |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
      midtot |  -.2490467   .4368666    -0.57   0.569     -1.10529    .6071961
    dregfrc6 |   4.073772   3.488994     1.17   0.243    -2.764529    10.91207
      lnrgdp |  -1.860991   .3713526    -5.01   0.000    -2.588828   -1.133153
      growth |  -.3851388   1.277004    -0.30   0.763    -2.888021    2.117743
       islam |   .0184209   .0072628     2.54   0.011      .004186    .0326557
       _cons |    8.09606   1.956796     4.14   0.000      4.26081    11.93131
-------------+----------------------------------------------------------------
       /ln_p |   .1903058   .1291164     1.47   0.141    -.0627577    .4433692
-------------+----------------------------------------------------------------
           p |   1.209619   .1561817                       .939171    1.557947
         1/p |   .8267063   .1067413                      .6418702    1.064769
-------------+----------------------------------------------------------------

. *That replicates Model 5 of Table 2.
. stweib pbit dregfrc6 lnrgdp growth islam if year>1959 & year<1989, robust cluster(ccode) nohr

         failure _d:  autrev
   analysis time _t:  dsctr
                 id:  panelid


Fitting constant-only model:

Iteration 0:   log likelihood = -43.063836
Iteration 1:   log likelihood = -40.678014
Iteration 2:   log likelihood = -40.643988
Iteration 3:   log likelihood = -40.643976
Iteration 4:   log likelihood = -40.643976

Fitting full model:

Iteration 0:   log likelihood = -40.643976
Iteration 1:   log likelihood = -27.741178
Iteration 2:   log likelihood = -23.919568
Iteration 3:   log likelihood = -23.586742
Iteration 4:   log likelihood = -23.585435
Iteration 5:   log likelihood = -23.585435

Weibull regression -- log relative-hazard form

No. of subjects =           63                     Number of obs   =      1058
No. of failures =           13
Time at risk    =         1058
                                                   Wald chi2(5)    =     60.25
Log likelihood  =   -23.585435                     Prob > chi2     =    0.0000

                            (standard errors adjusted for clustering on ccode)
-------------+----------------------------------------------------------------
             |               Robust
          _t |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
        pbit |  -.7582944   .6604546    -1.15   0.251    -2.052762    .5361728
    dregfrc6 |   1.731073   3.523074     0.49   0.623    -5.174026    8.636172
      lnrgdp |  -1.808032   .2632866    -6.87   0.000    -2.324064      -1.292
      growth |   .0235162   1.073567     0.02   0.983    -2.080636    2.127668
       islam |   .0191514   .0084254     2.27   0.023     .0026379    .0356649
       _cons |   7.760785   1.458071     5.32   0.000     4.903018    10.61855
-------------+----------------------------------------------------------------
       /ln_p |    .306264   .1604111     1.91   0.056     -.008136    .6206641
-------------+----------------------------------------------------------------
           p |   1.358341    .217893                       .991897    1.860163
         1/p |   .7361922   .1180934                      .5375873    1.008169
-------------+----------------------------------------------------------------

. *That replicates Model 6 of Table 2.
. stweib pbitot dregfrc6 lnrgdp growth islam if year>1959 & year<1989, robust cluster(ccode) nohr

         failure _d:  autrev
   analysis time _t:  dsctr
                 id:  panelid


Fitting constant-only model:

Iteration 0:   log likelihood = -43.063836
Iteration 1:   log likelihood = -40.678014
Iteration 2:   log likelihood = -40.643988
Iteration 3:   log likelihood = -40.643976
Iteration 4:   log likelihood = -40.643976

Fitting full model:

Iteration 0:   log likelihood = -40.643976
Iteration 1:   log likelihood = -28.423601
Iteration 2:   log likelihood = -25.012677
Iteration 3:   log likelihood = -24.746118
Iteration 4:   log likelihood = -24.745211
Iteration 5:   log likelihood = -24.745211

Weibull regression -- log relative-hazard form

No. of subjects =           63                     Number of obs   =      1058
No. of failures =           13
Time at risk    =         1058
                                                   Wald chi2(5)    =     37.22
Log likelihood  =   -24.745211                     Prob > chi2     =    0.0000

                            (standard errors adjusted for clustering on ccode)
-------------+----------------------------------------------------------------
             |               Robust
          _t |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
      pbitot |  -.0836085   .2557773    -0.33   0.744    -.5849227    .4177057
    dregfrc6 |    1.67088   3.344377     0.50   0.617    -4.883978    8.225738
      lnrgdp |  -1.692386   .3346012    -5.06   0.000    -2.348193    -1.03658
      growth |  -.0746846   1.151656    -0.06   0.948    -2.331888    2.182519
       islam |   .0185068   .0087361     2.12   0.034     .0013843    .0356292
       _cons |   7.060131   2.035354     3.47   0.001     3.070911    11.04935
-------------+----------------------------------------------------------------
       /ln_p |   .2302917   .1641619     1.40   0.161    -.0914596    .5520431
-------------+----------------------------------------------------------------
           p |   1.258967   .2066744                      .9125982    1.736798
         1/p |   .7943018   .1303941                      .5757723    1.095773
-------------+----------------------------------------------------------------

. *That replicates Model 7 of Table 2.
. stweib iw dregfrc6 lnrgdp growth islam if year>1959 & year<1993, robust cluster(ccode) nohr

         failure _d:  autrev
   analysis time _t:  dsctr
                 id:  panelid


Fitting constant-only model:

Iteration 0:   log likelihood = -54.304906
Iteration 1:   log likelihood = -51.384946
Iteration 2:   log likelihood = -51.339033
Iteration 3:   log likelihood = -51.339021
Iteration 4:   log likelihood = -51.339021

Fitting full model:

Iteration 0:   log likelihood = -51.339021
Iteration 1:   log likelihood = -36.383223
Iteration 2:   log likelihood = -31.890671
Iteration 3:   log likelihood = -31.485739
Iteration 4:   log likelihood = -31.460983
Iteration 5:   log likelihood =  -31.45641
Iteration 6:   log likelihood = -31.455838
Iteration 7:   log likelihood = -31.455811
Iteration 8:   log likelihood = -31.455809

Weibull regression -- log relative-hazard form

No. of subjects =           82                     Number of obs   =      1277
No. of failures =           16
Time at risk    =         1277
                                                   Wald chi2(5)    =   1897.26
Log likelihood  =   -31.455809                     Prob > chi2     =    0.0000

                            (standard errors adjusted for clustering on ccode)
-------------+----------------------------------------------------------------
             |               Robust
          _t |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
          iw |  -20.74407   .6538103   -31.73   0.000    -22.02552   -19.46263
    dregfrc6 |   4.388851   3.335301     1.32   0.188    -2.148219    10.92592
      lnrgdp |   -1.89001   .3629545    -5.21   0.000    -2.601387   -1.178632
      growth |  -.4707213   1.204923    -0.39   0.696    -2.832327    1.890885
       islam |   .0190177   .0076708     2.48   0.013     .0039832    .0340522
       _cons |   8.125682   2.006607     4.05   0.000     4.192804    12.05856
-------------+----------------------------------------------------------------
       /ln_p |   .1919317   .1332162     1.44   0.150    -.0691673    .4530307
-------------+----------------------------------------------------------------
           p |   1.211588   .1614032                      .9331705    1.573073
         1/p |   .8253633   .1099518                      .6356986    1.071616
-------------+----------------------------------------------------------------

. *That replicates Model 8 from Table 2.
. log close
       log:  E:\intdemdata\JOP Aug 01 replication.smcl
  log type:  smcl
 closed on:  21 Mar 2003, 15:55:18
-------------