感谢您的关注和支持，谢谢

感谢回复！！

两个软件都可以，最重要的是数据处理的逻辑和步骤要清晰。据我的经验来看，stata软件会更“智能”些，你可以保存步骤，下次直接用或修改后再用。

是呢，谢谢，
朋友，我还有个问题，就是第一步去量纲化以后，会出现0但是后面计算熵值又需要取对数，这个你是怎么处理的？我之前看论文写的是 所有的去量纲化的数值+0.0001   ，但是这几天看代码发现有的处理方法是 直接把0替换成0.00001，你觉得那种会合适一些？

直接用stata代码计算吧，会非常快的

谢谢回复，stata代码调好，确实很快

stata代码

从个人使用经验来看，stata使用起来更简单，并且在do文档中能快速便捷找到所存在问题。
具体来说操作步骤如下：
xtset id year
****标准化-正向指标********
global all_var        x2        x3        x4        x5        x6        x7        x9        x10        x14        x15        x16        x19
foreach i in $all_var{
egen min_`i'=min(`i')
egen max_`i'=max(`i')
gen s`i'=(`i'-min_`i')/(max_`i'-min_`i')  //正向指标
replace s`i'=0.0001 if s`i'==0
}

*****标准化-负向指标****
global neg_var x1 x8 x11        x12        x13 x17        x18       
foreach i in $neg_var{
egen min_`i'=min(`i')
egen max_`i'=max(`i')
gen s`i'=(max_`i'-`i')/(max_`i'-min_`i')  //正向指标
replace s`i'=0.0001 if s`i'==0
}

order id year x* min* max* s*

****p值********
forvalue i=1(1)19{
egen sums_`i'=sum(sx`i')
gen p`i'=sx`i'/sums_`i'
}
order id year x* min* max* s* sums* p*


****e/d********
forvalue i=1(1)19{
egen l`i'=sum(p`i'*ln(p`i'))
gen e`i'=-l`i'/ln(3)  //obs的个数，数据拉到底就知道了
gen d`i'=1-e`i'
}

order id year x* min* max* s* sums* p* l* e* d*

//求权重w-----------------------------------------------------------------------
forvalue i=1(1)19{
egen f`i'=rowtotal(d*)
gen w`i'=d`i'/f`i'
} //w就是权重

我喜欢用excel，简单好用