面板数据熵权法的Stata程序

之前，推荐过①熵(shang)值法, 为多指标综合评价提供依据的工具，②模糊综合评价法是什么？如何通过软件实现, 怎么去解读他的结果呢？③层次分析法AHP是什么? 决策评价大法，④全局主成分分析GPCA是什么？为什么比传统PCA更好, 具体如何使用呢？⑤熵权法代码, entropy weight method了解一下！

*自定义函数

cap program drop EWM

program define EWM

args var statue rn

var：待处理的变量，statue=±1分别表示正向指标和负向指标，rn:r年与n个观测值的乘积

quietly{

*===第一步：标准化 var'_sd

sumvar'

scalar min=r(min)

scalar max=r(max)

g var'_sd=(var'-min)/(max-min)

if statue'==-1{

noisily dis as error "负向指标"

replacevar'sd=1-`var'sd

}

*===第二步：指标归一化 var'_p

 genvar'sds= `var'sd+0.00000001  //添加偏移量

egen var'_sds_sum=sum(var'sds)gen var'_sdw=var'sds / `var'sds_sum

*===第三步：计算信息熵 var'_e

genvar'sij=-1/ln(rn')*var'sdw*ln(var'_sdw)

egenvar's=sum(`var'sij)

*===第四步：计算信息熵的冗余度 var'_d

genvar'g=1-`var's

*===第五步：清除多余变量 只保留var'_sdvar'ddrop var'_sdsvar'sdw `var'sds_sum var'_sijvar'_s }

end

*主函数

*【1】求`var'_d

*调用EWM程序，依次输入三个参数var statue rn, 注：负向指标会提示文字*放入指标（除第5个为负向指标外，其他都是正向指标

EWM var1 1 270

EWM var2 1 270

EWM var3 1 270

EWM var4 1 270

EWM var5 -1 270

EWM var6 1 270

EWM var7 1 270

EWM var8 1 270

EWM var9 1 270

EWM var10 1 270

【2】差异系数归一化，求权重w

gen sum_d=var1_d+var2_d+var3_d+var4_d+var5_d+var6_d+var7_d+var8_d+var9_d+var10

forvalues i=1/10{

gen wi'=vari'_d/sum_d

}drop sum_d

*展示权重

*list w* in 1

*【3】最终变量测度

gen Y = var1_sdw1+var2_sdw2+var3_sdw3+var4_sdw4+var5_sdw5+var6_sdw6+var7_sdw7+var8_sdw8+var9_sdw9+var10_sdw10

drop _

推荐内容：1.主成分分析/因子分析/判别分析，见过的最通俗易懂的解释，2.面板数据聚类, 因子分析和主成分分析咋做? 3.主成分分析和因子分析 PC, FA 及运行结果，4.主成分分析PCA理论讲解与程序实现指南，5.因子分析法FA理论讲解与程序实现指南，6.GVAR, 全局VAR模型是什么？该如何用软件实现, 有哪些研究文献和最新进展！

下面短链接文章属于合集，可以收藏起来阅读，不然以后都找不到了。