stata熵值法编程，程序有点问题，有偿请教大神 - Stata专版

0关注
2粉丝

已卖：29份资源

博士生

2%

还不是VIP/贵宾

-

0%

威望: 0 级
论坛币: 19 个
通用积分: 5.0615
学术水平: 0 点
热心指数: 0 点
信用等级: 0 点
经验: 406 点
帖子: 29
精华: 0
在线时间: 412 小时
注册时间: 2009-3-4
最后登录: 2025-11-7

楼主

hanxinyue1983 发表于 2019-2-17 20:14:50 |AI写论文

是否 +2 论坛币

k人参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我修改了一个程序包计算熵值，但是运行后显示找不到x_1，不知道问题出在哪里，请大神帮忙看看
程序包是为了实现统一各指标权重的基础上，区分正负指标计算熵值
问题可能出现在这部分程序：

哪位大神知道问题出在哪里，可以有偿帮助修改
apt prog drop szf1
prog szf1
set more off
preserve
quietly {
local m = `1'
local n = `2'
local n2 = `n'+2
if `m'>`n' {
mat A = J(`m',3,.)
}
else {
mat A = J(`n',3,.)
}
//所有指标
global all_var $positive_var $negative_var
//年份
qui sum year
global min_year=r(min)
global max_year=r(max)

forvalues year=$min_year / $max_year{
use data.dta, clear
keep if year==`year'

//标准化数据正向指标
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
}
//
forvalues i = 3/`n2' {
local j = `i'-2
capt drop bss`j'
ge bss`j' = `x_`i''
su bss`j'
ge bs`j' = bss`j'/r(sum)
}