--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  P:\Paper\FIR\final\replication\tc1fc1.txt
  log type:  text
 opened on:  21 Oct 2021, 13:29:41

. *
. * This is a replication program (stata 17.0) for:
. *
. * Kiyota, Kozo and Yoshinori Kurokawa (2021) "Factor Intensity Reversals Redux: Feenstra Is Right!"
. * Table C1 and Figure C1
. *
. * Please note that the different version of the software may generate slightly different results from the published version.
. *
. clear all

. set more off

. set matsize 5000
set matsize ignored.
    Matrix sizes are no longer limited by c(matsize) in modern Statas.  Matrix sizes are now limited by edition of Stata.  See limits for more details.

. use Kiyota_Kurokawa_2021_FIR, clear

. 
. *********************************************************************************************
. * Figure C1
. *********************************************************************************************
. 
. gen lnVA = ln(VA)
(45 missing values generated)

. gen logKL = log10(KL)
(30 missing values generated)

. gen logw = log10(wL/LH)
(30 missing values generated)

. keep if year == 1975 | year == 1985 | year == 1995 | year == 2005
(25,920 observations deleted)

. drop if prefecture == 0
(60 observations deleted)

. keep if ind >= 3 & ind <= 15
(376 observations deleted)

. twoway lfit logw logKL if ind == 3 || lfit logw logKL if ind == 4 || lfit logw logKL if ind == 5 || lfit logw logKL if ind == 6 || lfit logw logKL if ind == 7 || lfit logw lo
> gKL if ind == 8 || lfit logw logKL if ind == 9 || lfit logw logKL if ind == 10 || lfit logw logKL if ind == 11 || lfit logw logKL if ind == 12 || lfit logw logKL if ind == 13
>  || lfit logw logKL if ind == 14 || lfit logw logKL if ind == 15, by(year, note("") legend(off) yrescale xrescale) ytitle("Fitted value of log(w/r)") xtitle("log(K/L)")

. graph export fc1.pdf, replace
file fc1.pdf saved as PDF format

. *
. keep if year == 2005
(1,833 observations deleted)

. 
. *********************************************************************************************
. * Table C1
. *********************************************************************************************
. 
. preserve

. reg logw logKL, robust

Linear regression                               Number of obs     =        611
                                                F(1, 609)         =     312.96
                                                Prob > F          =     0.0000
                                                R-squared         =     0.3452
                                                Root MSE          =     .11284

------------------------------------------------------------------------------
             |               Robust
        logw | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
       logKL |   .1945736   .0109987    17.69   0.000     .1729737    .2161736
       _cons |   .2118034   .0127858    16.57   0.000     .1866938    .2369129
------------------------------------------------------------------------------

. outreg using tc1.csv, com bracket 3aster coefastr se replace

. nlcom (1/_b[logKL])

       _nl_1: 1/_b[logKL]

------------------------------------------------------------------------------
        logw | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |   5.139442   .2905175    17.69   0.000     4.570039    5.708846
------------------------------------------------------------------------------

. matrix M = r(b)

. matrix N = r(V)

. matrix m = vec(M)

. matrix n = vec(N)

. svmat m, names(sigma)

. svmat n, names(sigma_se)

. rename sigma1 sigma

. rename sigma_se1 sigma_se

. keep sigma sigma_se

. drop if sigma == .
(610 observations deleted)

. replace sigma_se = sigma_se^(1/2)
(1 real change made)

. restore

. *
. * by industry
. *
. forvalues i = 3/15 {
  2.    preserve
  3.    keep if ind == `i'
  4.    display "*** industry == `i' ***"
  5.    reg logw logKL, robust
  6.    outreg using tc1.csv, com bracket 3aster coefastr se append
  7.    nlcom (1/_b[logKL])
  8.    matrix M = r(b)
  9.    matrix N = r(V)
 10.    matrix m = vec(M)
 11.    matrix n = vec(N)
 12.    svmat m, names(sigma)
 13.    svmat n, names(sigma_se)
 14.    rename sigma1 sigma
 15.    rename sigma_se1 sigma_se
 16.    keep sigma sigma_se
 17.    drop if sigma == .
 18.    replace sigma_se = sigma_se^(1/2)
 19.    restore
 20.    }
(564 observations deleted)
*** industry == 3 ***

Linear regression                               Number of obs     =         47
                                                F(1, 45)          =      31.51
                                                Prob > F          =     0.0000
                                                R-squared         =     0.2402
                                                Root MSE          =     .07225

------------------------------------------------------------------------------
             |               Robust
        logw | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
       logKL |   .2751412   .0490151     5.61   0.000     .1764198    .3738627
       _cons |   .0091405   .0422528     0.22   0.830    -.0759611    .0942421
------------------------------------------------------------------------------

       _nl_1: 1/_b[logKL]

------------------------------------------------------------------------------
        logw | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |   3.634497    .647468     5.61   0.000     2.365483    4.903511
------------------------------------------------------------------------------
(46 observations deleted)
(1 real change made)
(564 observations deleted)
*** industry == 4 ***

Linear regression                               Number of obs     =         47
                                                F(1, 45)          =       1.85
                                                Prob > F          =     0.1809
                                                R-squared         =     0.0637
                                                Root MSE          =      .0802

------------------------------------------------------------------------------
             |               Robust
        logw | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
       logKL |   .0754355   .0555029     1.36   0.181    -.0363531    .1872241
       _cons |   .1295955   .0458888     2.82   0.007     .0371706    .2220204
------------------------------------------------------------------------------

       _nl_1: 1/_b[logKL]

------------------------------------------------------------------------------
        logw | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |   13.25636   9.753592     1.36   0.174    -5.860326    32.37305
------------------------------------------------------------------------------
(46 observations deleted)
(1 real change made)
(564 observations deleted)
*** industry == 5 ***

Linear regression                               Number of obs     =         47
                                                F(1, 45)          =       3.27
                                                Prob > F          =     0.0771
                                                R-squared         =     0.0818
                                                Root MSE          =     .07942

------------------------------------------------------------------------------
             |               Robust
        logw | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
       logKL |  -.0814612   .0450293    -1.81   0.077    -.1721549    .0092326
       _cons |   .4817293   .0510292     9.44   0.000     .3789512    .5845075
------------------------------------------------------------------------------

       _nl_1: 1/_b[logKL]

------------------------------------------------------------------------------
        logw | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |  -12.27579   6.785695    -1.81   0.070    -25.57551    1.023931
------------------------------------------------------------------------------
(46 observations deleted)
(1 real change made)
(564 observations deleted)
*** industry == 6 ***

Linear regression                               Number of obs     =         47
                                                F(1, 45)          =       0.94
                                                Prob > F          =     0.3377
                                                R-squared         =     0.0371
                                                Root MSE          =     .08134

------------------------------------------------------------------------------
             |               Robust
        logw | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
       logKL |   .0580607    .059915     0.97   0.338    -.0626143    .1787358
       _cons |   .4463262   .0880575     5.07   0.000     .2689692    .6236832
------------------------------------------------------------------------------

       _nl_1: 1/_b[logKL]

------------------------------------------------------------------------------
        logw | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |   17.22334   17.77341     0.97   0.333     -17.6119    52.05859
------------------------------------------------------------------------------
(46 observations deleted)
(1 real change made)
(564 observations deleted)
*** industry == 7 ***

Linear regression                               Number of obs     =         47
                                                F(1, 45)          =       0.00
                                                Prob > F          =     0.9848
                                                R-squared         =     0.0000
                                                Root MSE          =     .08289

------------------------------------------------------------------------------
             |               Robust
        logw | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
       logKL |  -.0006177   .0323463    -0.02   0.985    -.0657665    .0645311
       _cons |    .617593   .0625806     9.87   0.000     .4915491    .7436369
------------------------------------------------------------------------------

       _nl_1: 1/_b[logKL]

------------------------------------------------------------------------------
        logw | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |  -1618.931   84777.73    -0.02   0.985    -167780.2    164542.4
------------------------------------------------------------------------------
(46 observations deleted)
(1 real change made)
(564 observations deleted)
*** industry == 8 ***

Linear regression                               Number of obs     =         47
                                                F(1, 45)          =       0.88
                                                Prob > F          =     0.3532
                                                R-squared         =     0.0213
                                                Root MSE          =       .082

------------------------------------------------------------------------------
             |               Robust
        logw | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
       logKL |   .0657177    .070048     0.94   0.353    -.0753662    .2068015
       _cons |   .3248928   .0626877     5.18   0.000     .1986333    .4511523
------------------------------------------------------------------------------

       _nl_1: 1/_b[logKL]

------------------------------------------------------------------------------
        logw | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |    15.2166   16.21926     0.94   0.348    -16.57257    47.00578
------------------------------------------------------------------------------
(46 observations deleted)
(1 real change made)
(564 observations deleted)
*** industry == 9 ***

Linear regression                               Number of obs     =         47
                                                F(1, 45)          =       0.29
                                                Prob > F          =     0.5954
                                                R-squared         =     0.0078
                                                Root MSE          =     .08256

------------------------------------------------------------------------------
             |               Robust
        logw | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
       logKL |   .0270447   .0505721     0.53   0.595    -.0748128    .1289022
       _cons |   .4725345   .0680478     6.94   0.000     .3354792    .6095899
------------------------------------------------------------------------------

       _nl_1: 1/_b[logKL]

------------------------------------------------------------------------------
        logw | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |   36.97582   69.14277     0.53   0.593    -98.54152    172.4931
------------------------------------------------------------------------------
(46 observations deleted)
(1 real change made)
(564 observations deleted)
*** industry == 10 ***

Linear regression                               Number of obs     =         47
                                                F(1, 45)          =       2.35
                                                Prob > F          =     0.1320
                                                R-squared         =     0.0456
                                                Root MSE          =     .08097

------------------------------------------------------------------------------
             |               Robust
        logw | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
       logKL |     .12447   .0811377     1.53   0.132    -.0389498    .2878898
       _cons |   .2741357   .0483434     5.67   0.000     .1767672    .3715042
------------------------------------------------------------------------------

       _nl_1: 1/_b[logKL]

------------------------------------------------------------------------------
        logw | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |   8.034066   5.237134     1.53   0.125    -2.230528    18.29866
------------------------------------------------------------------------------
(46 observations deleted)
(1 real change made)
(564 observations deleted)
*** industry == 11 ***

Linear regression                               Number of obs     =         47
                                                F(1, 45)          =       6.84
                                                Prob > F          =     0.0121
                                                R-squared         =     0.1602
                                                Root MSE          =     .07596

------------------------------------------------------------------------------
             |               Robust
        logw | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
       logKL |   .2255237   .0862347     2.62   0.012     .0518381    .3992094
       _cons |   .2288853    .080208     2.85   0.007     .0673382    .3904325
------------------------------------------------------------------------------

       _nl_1: 1/_b[logKL]

------------------------------------------------------------------------------
        logw | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |   4.434123   1.695499     2.62   0.009     1.111005    7.757241
------------------------------------------------------------------------------
(46 observations deleted)
(1 real change made)
(564 observations deleted)
*** industry == 12 ***

Linear regression                               Number of obs     =         47
                                                F(1, 45)          =       1.75
                                                Prob > F          =     0.1920
                                                R-squared         =     0.0606
                                                Root MSE          =     .08033

------------------------------------------------------------------------------
             |               Robust
        logw | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
       logKL |   .0882832   .0666442     1.32   0.192    -.0459451    .2225114
       _cons |   .3356892   .0792015     4.24   0.000     .1761692    .4952092
------------------------------------------------------------------------------

       _nl_1: 1/_b[logKL]

------------------------------------------------------------------------------
        logw | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |   11.32719   8.550792     1.32   0.185    -5.432056    28.08643
------------------------------------------------------------------------------
(46 observations deleted)
(1 real change made)
(564 observations deleted)
*** industry == 13 ***

Linear regression                               Number of obs     =         47
                                                F(1, 45)          =      19.94
                                                Prob > F          =     0.0001
                                                R-squared         =     0.2751
                                                Root MSE          =     .07057

------------------------------------------------------------------------------
             |               Robust
        logw | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
       logKL |   .1781138   .0398911     4.47   0.000      .097769    .2584585
       _cons |   .3076581    .035781     8.60   0.000     .2355915    .3797248
------------------------------------------------------------------------------

       _nl_1: 1/_b[logKL]

------------------------------------------------------------------------------
        logw | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |   5.614389   1.257422     4.47   0.000     3.149888     8.07889
------------------------------------------------------------------------------
(46 observations deleted)
(1 real change made)
(564 observations deleted)
*** industry == 14 ***

Linear regression                               Number of obs     =         47
                                                F(1, 45)          =       3.49
                                                Prob > F          =     0.0682
                                                R-squared         =     0.0817
                                                Root MSE          =     .07943

------------------------------------------------------------------------------
             |               Robust
        logw | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
       logKL |   .0554996   .0297015     1.87   0.068    -.0043223    .1153215
       _cons |   .3856056   .0284437    13.56   0.000      .328317    .4428942
------------------------------------------------------------------------------

       _nl_1: 1/_b[logKL]

------------------------------------------------------------------------------
        logw | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |   18.01815    9.64271     1.87   0.062    -.8812143    36.91751
------------------------------------------------------------------------------
(46 observations deleted)
(1 real change made)
(564 observations deleted)
*** industry == 15 ***

Linear regression                               Number of obs     =         47
                                                F(1, 45)          =      14.85
                                                Prob > F          =     0.0004
                                                R-squared         =     0.2627
                                                Root MSE          =     .07117

------------------------------------------------------------------------------
             |               Robust
        logw | Coefficient  std. err.      t    P>|t|     [95% conf. interval]
-------------+----------------------------------------------------------------
       logKL |    .265438   .0688699     3.85   0.000     .1267269     .404149
       _cons |    .131817   .0548166     2.40   0.020     .0214108    .2422233
------------------------------------------------------------------------------

       _nl_1: 1/_b[logKL]

------------------------------------------------------------------------------
        logw | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
       _nl_1 |   3.767359   .9774694     3.85   0.000     1.851554    5.683164
------------------------------------------------------------------------------
(46 observations deleted)
(1 real change made)

. 
. log close
      name:  <unnamed>
       log:  P:\Paper\FIR\final\replication\tc1fc1.txt
  log type:  text
 closed on:  21 Oct 2021, 13:29:47
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
