[编程问题求助] stata熵值法编程，程序有点问题，有偿请教大神 [推广有奖]

11楼

米高兄弟

发表于 2019-6-29 23:41:23

***********未进行分年进行整体熵值法********888
//正向指标
global positive_var x1 x2 x3
//负向指标
global negative_var x4 x5

//标准化数据正向指标
foreach i in $positive_var {
qui sum `i'
gen x_`i'=(`i'-r(min))/(r(max)-r(min))
replace x_`i'=0.00001 if x_`i'==0
}

//标准化数据负向指标
foreach i in $negative_var {
qui sum `i'
gen x_`i'=(r(max)-`i')/(r(max)-r(min))
replace x_`i'=0.00001 if x_`i'==0
}

//计算各指标的比重
global positive_var x_x1 x_x2 x_x3
global negative_var x_x4  x_x5
global all_var $positive_var $negative_var  //所有指标
global all_var x_x1 x_x2 x_x3 x_x4  x_x5

foreach i in $all_var{
      qui sum `i'
egen `i'_sum=sum(`i')
gen y_`i'=`i'/`i'_sum
}

//根据比重计算各分量的信息熵
gen n=_N

foreach i in $all_var {
gen y_lny_`i'=y_`i'*ln(y_`i')
}

//求和
global positive_var y_lny_x_x1 y_lny_x_x2 y_lny_x_x3
global negative_var y_lny_x_x4  y_lny_x_x5
global all_var $positive_var $negative_var  //所有指标
foreach i in $all_var {
egen y_lny_`i'_sum=sum(`i')
}
//计算各指标的贡献总量

foreach i in $all_var {
gen E_`i'= -1/ln(n)*y_lny_`i'_sum
}
//计算各指标权重
foreach i in $all_var {
gen d_`i'= 1-E_`i'
}
egen d_sum = rowtotal(d_*)
foreach i in $all_var {
gen W_`i'= d_`i'/d_sum
}
//计算综合得分
foreach i in $all_var {
gen Score_`i'= `i'*W_`i'
}

egen Score=rowtotal(Score_*)

keep id year $all_var Score
save data_`year', replace
}

已有 1 人评分	论坛币	学术水平	热心指数	信用等级	收起理由
admin_kefu	+ 20	+ 2	+ 2	+ 2	热心帮助其他会员

总评分: 论坛币 + 20 学术水平 + 2 热心指数 + 2 信用等级 + 2 查看全部评分

12楼

米高兄弟

发表于 2019-6-29 23:45:44

有位大咖卖了个命令，但她是分类回归的，花钱很多初学者都还是懵的，我初步研究了一下按照不分时间段，即同一年下如何按照步骤做的方法，供大家参考，都是穷学生，没钱，能帮就帮吧，看到这段命令就节约看不到就算了！楼主的命令是她部分的内容和另外一个老师的结合，都不合适，即没理解透彻，但也不怪，需要灰色关联命令的也可以联系我，可以模仿进行写循环命令，我因为循环不太熟悉用的最笨方法操作的！