*-几何加权乘积

*-Problem: 几何加权实际汇率的计算(以出口额为权重)

  *---------------------------------------------------------
  *         i=N                             i=N
  *  EER = PRODUCT(exchange_i^weight_i),    SUM(weight_i)=1
  *         i=1                             i=1   
  
  *      = (x1^w1)*(x2^w2)*...*(xn^wn)
  *
  *  where, x=exchange, w=weight
  *---------------------------------------------------------
  
*-数据概况:
   
  *-firm1 的产品出口到三个国家(country=1,2,3), export 是对应的出口额
  * 则 firm1 对 country1 的出口权重为 weight11=60/(60+50+80)

    clear
    input firm country exch  export
            1   1       8      60
            1   2       0.03   50
            1   3      10      80
            2   2       0.03   68
            2   4       3      80
            3   2       0.03  200
            3   3      10      90
            3   4       3     250
            3   5       1.1   120
    end

    bysort firm: egen weight = pc(export), prop  // 权重: weight
   
   
  *-----------
  *-Method 1:  通过对数变换转变成求和问题，进而exp变换还原成连乘积
  
    *--------------------------------------------------
    *-Note: 原理
    *
    *    ln(EER) = ln[(x1^w1)*(x2^w2)*(x3^w3)]
    *            = w1*ln(x1) + w2*ln(x2) + w3*ln(x3)
    *            = Y
    *  so
    *       EER  = exp(Y)
    *--------------------------------------------------
    gen lnexch = ln(exch)
    gen lnexch_x_weight = lnexch*weight
    bysort firm: egen sum_exch = sum(lnexch_x_weight)
    gen EER1 = exp(sum_exch)

  *-----------
  *-Method 2:  直接计算连乘积
  
    net install dm71.pkg, replace    // 安装外部命令 prod
    gen exch_weight = exch^weight
    egen EER2 = prod(exch_weight), by(firm)

  *-查验结果
    br firm-weight EER*

firm	country	exch	export	weight	EER1	EER2
1	1	8	60	.31578947	2.0206048	2.0206048
1	2	.03	50	.26315789	2.0206048	2.0206048
1	3	10	80	.42105263	2.0206048	2.0206048
2	2	.03	68	.45945946	.36157828	.36157828
2	4	3	80	.54054054	.36157828	.36157828
3	2	.03	200	.3030303	.72969088	.72969088
3	3	10	90	.13636364	.72969088	.72969088
3	4	3	250	.37878788	.72969088	.72969088
3	5	1.1	120	.18181818	.72969088	.72969088

连老师，为什么在执行  bysort firm: egen weight = pc(export), prop  时老是提示option prop not allowed？我已经检查过了 语法没有任何问题 egen的帮助文件也是如你所说的那样来操作的，去掉prop是可以操作的，这是是什么情况？

我测试过，没有问题。

破解的软件不可靠，用stata12不行，用stata11又可以了。谢谢连老师！