熵权法stata程序

• Edited by Linhao Cui. 公式部分参考：(https://wenku.baidu.com/view/a2600818bd64783e09122bf6.html)
• 仅供学习交流使用，如有疑问或者问题，欢迎质疑讨论，如有侵权 请联系撤回。

（1）数据标准化处理
$$\begin{gathered} 正向指标X_{ij}=\frac{X_{ij}-min(X_j)}{max(X_j)-min(X_j)} \\ 负指标X_{ij}=\frac{max(X_j)-X_{ij}}{max(X_j)-min(X_j)} \end{gathered}$$
(2)计算第i年第j项指标值的比重：
$$Y_{ij}=\frac{X_{ij}}{{\sum }_i^m{X}_{ij}}$$
(3)计算指标信息熵
$$e_j=-k\sum _i^m(Y_{ij}\times ln{Y}_{ij})$$
(4)计算信息熵冗余度：
$$d_j=1-e_j$$
(5) 计算指标权重：
$$W_i=d_j/\sum _j^n{d}_j$$
(6)单项指标评价得分
$$S_{ij}=W_i\times X_{ij}$$
其中，$X_{ij}$ 代表第i年第j项评价指标的数值，k= 1/lnm , m为评价年数，n为指标个数

**熵权法 
** step 1 指标的标准化处理
 by 地区,sort: egen m = count(年份)
global xlist " 高耗能工业利息 环保企业资产 高耗能行业资产 高耗能工业利息占比 环保企业资产占比 高耗能产业资产占比 绿色信贷发展水平 治污投资比 碳排放强度 绿色保险赔付率 绿色保险规模比"
foreach x of global xlist{
  bys 地区 年份, sort: egen minn`x' = min(`x') 
  bys 地区 年份, sort: egen maxx`x' = max(`x') 
  gen standard`x' = (`x' - minn`x')/ (maxx`x' - minn`x') if `x' > = 0 // 正向指标
  
}

** step 2 计算第i年第j项指标的比重
foreach x of global xlist {
  by 地区,sort: egen summm`x' = total(`x') 
  gen Y`x' = standard`x' / summm`x'
}

**step3 计算指标信息熵
foreach x of global xlist{
  bys 地区 年份, sort: egen he`x' = total(Y`x' * ln(Y`x'))
  gen e`x' = - log(m) * he`x'
}

**step4 计算信息熵冗余度
foreach x of global xlist{
  gen d`x' = 1 - e`x'
  by 地区,sort: egen uselesshe`x' = total(d`x')
  gen W`x' = d`x' / uselesshe`x'
}

**step5 计算指标综合得分
foreach x of global xlist {
  gen S`x' = W`x' * `x'
}