数据库：工业企业数据库，用LP法在计算TFP时，我的命令如下：
按产业分别设置生产函数，所有按产业类别计算TFP
共562个产业类别

forvalues i=1/562 {
levpet lny if industryid ==`i', free(lnL) proxy(lnm) capital(lnK) valueadded  reps(250)
predict TFP_LP if industryid ==`i', omega
}

计算结果如下：
Levinsohn-Petrin productivity estimator
Dependent variable represents value added.      Number of obs      =     25635
Group variable (i): id                          Number of groups   =    516347
Time variable (t): year
                                                Obs per group: min =         1
                                                               avg =       3.6
                                                               max =        10
------------------------------------------------------------------------------
         lny |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         lnL |   .1120306     .00625    17.92   0.000     .0997809    .1242804
         lnK |   .1467407   .0101792    14.42   0.000     .1267899    .1666915
------------------------------------------------------------------------------
Wald test of constant returns to scale: Chi2 = 4326.34 (p = 0.0000).

分组数与样本数明显与我想要的不符；

另外：在计算TFP时，没法按产业来取系数，所有的样本都用一个系数估算，与我要求的不符。