熵权法STATA程序（修正）

鉴于很多童鞋和我探讨关于上一版的熵权法问题，在这里我打算补充修正一下熵权法程序，补充如下：
（1）这一版可用于面板数据处理，处理方法是将面板数据视为截面，同一城市不同年份视为不同个体,这样保证了权重不会随时间变化。而且这并不影响我们指标合成的初衷。
（2）这一版使用城市环境相关数据为例子，相应数据可在城市统计年鉴搜索得到。
（3）熵权法只是计算指标的客观权重方法，下文程序第五步计算得分时大家可以使用不同的指标标准化方法，如功效系数法、当量准值法等，这里我用的是极差标准化，大家也可以同时考虑纵向比较和横向比较，并不唯一。
当然可能仍然存在很多问题，欢迎大家和我讨论。
by Linhao Cui

** step 1 指标的标准化处理egen m = count(id)
set trace on
global zlist " 一般工业固体废物综合利用率_百分比_全市_补齐 年末耕地总资源 工业二氧化硫去除量_吨_全市_补齐 工业废水排放达标量_万吨_全市_补齐 污水处理厂集中处理率_百分比_全市_补齐 工业烟粉尘去除量_吨_全市_补齐 生活垃圾无害化处理率_百分比_全市_补齐 生活污水处理率 环境污染治理占GDP比重 可吸入细颗粒物年平均浓度 空气质量达到及好于二级的天数占比 建成区绿化覆盖率_百分比_市辖区 地区生产总值增长率 科技教育比重 专利授权数 在岗职工平均工资_元_全市_补齐 每万人在校大学生数 第三产业增加值占GDP比重 每万人拥有公共汽车_辆_市辖区 体育场馆数 公共图书馆 博物馆数"foreach x of global zlist{egen minn`x' = min(`x') egen maxx`x' = max(`x') gen standard`x' = (`x' - minn`x')/(maxx`x' - minn`x') //正向指标}foreach x of varlist 售水量_万吨_市辖区_补齐 能源消费量 电_补齐{egen minn`x' = min(`x') egen maxx`x' = max(`x') gen standard`x' = (maxx`x' - `x')/(maxx`x' - minn`x') //负向指标}global xlist "售水量_万吨_市辖区_补齐 能源消费量 一般工业固体废物综合利用率_百分比_全市_补齐 年末耕地总资源 电_补齐 工业二氧化硫去除量_吨_全市_补齐 工业废水排放达标量_万吨_全市_补齐 污水处理厂集中处理率_百分比_全市_补齐 工业烟粉尘去除量_吨_全市_补齐 生活垃圾无害化处理率_百分比_全市_补齐 生活污水处理率 可吸入细颗粒物年平均浓度 空气质量达到及好于二级的天数占比 建成区绿化覆盖率_百分比_市辖区 地区生产总值增长率 专利授权数 在岗职工平均工资_元_全市_补齐 每万人在校大学生数 第三产业增加值占GDP比重 每万人拥有公共汽车_辆_市辖区 体育场馆数 公共图书馆 博物馆数"** step 2 计算第i年第j项指标的比重
foreach x of global xlist {egen summm`x' = total(`x') gen Y`x' = standard`x'/summm`x'
}**step3 计算指标信息熵
foreach x of global xlist{egen he`x' = total(Y`x' * ln(Y`x'))gen e`x' = -(1/ln(m)) * he`x'
}**step4 计算信息熵冗余度
foreach x of global xlist{gen d`x' = 1 - e`x'  //stata貌似循环有问题 如遇到not found问题不能运行可尝试手动计算 egen uselesshe = rowtotal(d售水量_万吨_市辖区_补齐 d能源消费量 d一般工业固体废物综合利用率_百分比_全市_补齐 d年末耕地总资源 d电_补齐 d工业二氧化硫去除量_吨_全市_补齐 d工业废水排放达标量_万吨_全市_补齐 d污水处理厂集中处理率_百分比_全市_补齐 d工业烟粉尘去除量_吨_全市_补齐 d生活垃圾无害化处理率_百分比_全市_补齐 d生活污水处理率  d可吸入细颗粒物年平均浓度 d空气质量达到及好于二级的天数占比 d建成区绿化覆盖率_百分比_市辖区 d地区生产总值增长率  d专利授权数 d在岗职工平均工资_元_全市_补齐 d每万人在校大学生数 d第三产业增加值占GDP比重 d每万人拥有公共汽车_辆_市辖区 d体育场馆数 d公共图书馆 d博物馆数)gen W`x' = d`x'/uselesshe
}**step5 计算指标综合得分
foreach x of global xlist {gen S`x' = W`x' * `x'
}

熵权法STATA程序(第二版修正)相关推荐

基于熵权法对Topsis模型的修正
基于熵权法对Topsis模型的修正有n个要评价的对象,m个评价指标的标准化矩阵,可以使用层次分析法给这m个评价指标确定权重 ∑j=1mωj=1\sum_{j=1}^m{\omega_j}=1 j=1 ...
数学建模学习笔记（2）：TOPSIS方法（优劣解距离法）和熵权法修正
文章目录 TOPSIS方法概述 TOPSIS方法步骤(重点) 熵权法对TOPSIS模型的修正熵权法的步骤(重点) TOPSIS方法由C.L.Hwang和K.Yoon在1981年首次提出,在国内常简称 ...
综合评价分析法：熵权法matlab实现+层次分析法matlab实现
综合评价分析法:熵权法matlab实现+层次分析法matlab实现层次分析法matlab实现 clc;clear%层次分析法matlab实现A=[1 1/9 1/2 1/4 1/6 9 1 6 5 ...
State详细代码和介绍-Topsis熵权法评分（含视频教程）、Stata灰色关联度分析
一.topsis熵权法评分 1.数据来源:自主计算 2.时间跨度:无 3.区域范围:无 4.指标说明: 基于熵权法的topsis综合评分方法,结合了熵权法的客观求权重,以及topsis法的综合评分本教 ...
基于熵权法对TOPSIS法模型的修正
最近在学习数学建模,在B站发现一个特别不错的课程,讲的很全面,常考的算法都有涉及到:清风数学建模本文将结合熵权法介绍TOPSIS法,并将淡化原理的推导,更侧重于具体应用. TOPSIS法概述 TO ...
帮我用python写一个程序输入三个列表数据使用熵权法计算这三组数据的权重并输出...
好的,以下是用Python实现熵权法计算三组数据权重的代码: import numpy as np# 定义熵权法函数 def entropy_weight(x):x = np.array(x)# 计算 ...
《零基础数学建模》——TOPSIS+熵权法
目录前言一.TOPSIS法(优劣解距离法) 1.模型原理 2.建模步骤二.模型实现第一步:将原始矩阵正向化第二步:正向化矩阵标准化第三步:计算得分并归一化四.TOPSIS模型的总结与 ...
熵权法 —— python
目录一.熵权法介绍二.熵权法赋权步骤 1.指标正向化 mapminmax介绍 2.数据标准化 3.计算信息熵 4.计算权重以及得分四.实例分析 1.导入相关库 2.读取数据 2.指标正向化 2. ...
熵权法与Apriori算法对较多数据种类数据的处理
前言:前两天女朋友问我他们论文的时候看到他们用了基于熵权的TOPSIS模型做分析,由于本人没学过有关统计方面的知识,所以就去现学了一下.好巧不巧,在学校上课的时候老师提到了Apriori算法,突然联想 ...
数学建模笔记——评价类模型之熵权法
嗯,这次讲一讲熵权法,一种通过样本数据确定评价指标权重的方法. 熵权法之前我们提到了TOPSIS方法,用来处理有数据的评价类模型.TOPSIS方法还蛮简单的,大概就三步. 将原始数据矩阵正向化. 也 ...

熵权法STATA程序(第二版修正)

熵权法STATA程序（修正）

熵权法STATA程序(第二版修正)相关推荐

最新文章

热门文章

熵权法STATA程序(第二版修正)

熵权法STATA程序 （修正）

熵权法STATA程序(第二版修正)相关推荐

最新文章

热门文章

熵权法STATA程序（修正）