opreg Y, exit(Exit) state(K) proxy(I) free(L) cvars(Year) vce(bootstrap, seed(1) rep(5))
Note: ignoring 2 observations with missing value for proxy variable

(running opreg on estimation sample)

Bootstrap replications (5)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5
.....


Olley-Pakes productivity estimator              Number of obs      =     30318
Group variable (i): Stkcd                       Number of groups   =      4914
Time variable (t): Year
                                                Obs per group: min =         1
                                                               avg =       6.2
                                                               max =        11

                               (Replications based on 4,914 clusters in Stkcd)
------------------------------------------------------------------------------
             |   Observed   Bootstrap                         Normal-based
             | coefficient  std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
Y            |
           K |   .0112469   .0262279     0.43   0.668    -.0401588    .0626526
           L |   .6383019   .0170709    37.39   0.000     .6048436    .6717602
        Year |   .0301264   .0026188    11.50   0.000     .0249936    .0352592
------------------------------------------------------------------------------
State:         K
Free:          L
Control:       Year
Proxy:         I

. gen OP1=Y-_b[L]*L-_b[K]*K
(8 missing values generated)

. sum lp1_1 OP1

    Variable |        Obs        Mean    Std. dev.       Min        Max
-------------+---------------------------------------------------------
       lp1_1 |     30,312    17.10338    1.073605   14.43997   21.44722
         OP1 |     30,312    16.45872    .9151992   13.80934   20.81245

问题不大 先算了再说 然后看看你算的结果 和一般期刊中是不是接近

换一种相似方法