Python计算信息熵实例
作者:Mr.Jcak 时间:2022-08-17 16:13:39
计算信息熵的公式:n是类别数,p(xi)是第i类的概率
假设数据集有m行,即m个样本,每一行最后一列为该样本的标签,计算数据集信息熵的代码如下:
from math import log
def calcShannonEnt(dataSet):
numEntries = len(dataSet) # 样本数
labelCounts = {} # 该数据集每个类别的频数
for featVec in dataSet: # 对每一行样本
currentLabel = featVec[-1] # 该样本的标签
if currentLabel not in labelCounts.keys(): labelCounts[currentLabel] = 0
labelCounts[currentLabel] += 1
shannonEnt = 0.0
for key in labelCounts:
prob = float(labelCounts[key])/numEntries # 计算p(xi)
shannonEnt -= prob * log(prob, 2) # log base 2
return shannonEnt
补充知识:python 实现信息熵、条件熵、信息增益、基尼系数
我就废话不多说了,大家还是直接看代码吧~
import pandas as pd
import numpy as np
import math
## 计算信息熵
def getEntropy(s):
# 找到各个不同取值出现的次数
if not isinstance(s, pd.core.series.Series):
s = pd.Series(s)
prt_ary = pd.groupby(s , by = s).count().values / float(len(s))
return -(np.log2(prt_ary) * prt_ary).sum()
## 计算条件熵: 条件s1下s2的条件熵
def getCondEntropy(s1 , s2):
d = dict()
for i in list(range(len(s1))):
d[s1[i]] = d.get(s1[i] , []) + [s2[i]]
return sum([getEntropy(d[k]) * len(d[k]) / float(len(s1)) for k in d])
## 计算信息增益
def getEntropyGain(s1, s2):
return getEntropy(s2) - getCondEntropy(s1, s2)
## 计算增益率
def getEntropyGainRadio(s1, s2):
return getEntropyGain(s1, s2) / getEntropy(s2)
## 衡量离散值的相关性
import math
def getDiscreteCorr(s1, s2):
return getEntropyGain(s1,s2) / math.sqrt(getEntropy(s1) * getEntropy(s2))
# ######## 计算概率平方和
def getProbSS(s):
if not isinstance(s, pd.core.series.Series):
s = pd.Series(s)
prt_ary = pd.groupby(s, by = s).count().values / float(len(s))
return sum(prt_ary ** 2)
######## 计算基尼系数
def getGini(s1, s2):
d = dict()
for i in list(range(len(s1))):
d[s1[i]] = d.get(s1[i] , []) + [s2[i]]
return 1-sum([getProbSS(d[k]) * len(d[k]) / float(len(s1)) for k in d])
## 对离散型变量计算相关系数,并画出热力图, 返回相关性矩阵
def DiscreteCorr(C_data):
## 对离散型变量(C_data)进行相关系数的计算
C_data_column_names = C_data.columns.tolist()
## 存储C_data相关系数的矩阵
import numpy as np
dp_corr_mat = np.zeros([len(C_data_column_names) , len(C_data_column_names)])
for i in range(len(C_data_column_names)):
for j in range(len(C_data_column_names)):
# 计算两个属性之间的相关系数
temp_corr = getDiscreteCorr(C_data.iloc[:,i] , C_data.iloc[:,j])
dp_corr_mat[i][j] = temp_corr
# 画出相关系数图
fig = plt.figure()
fig.add_subplot(2,2,1)
sns.heatmap(dp_corr_mat ,vmin= - 1, vmax= 1, cmap= sns.color_palette('RdBu' , n_colors= 128) , xticklabels= C_data_column_names , yticklabels= C_data_column_names)
return pd.DataFrame(dp_corr_mat)
if __name__ == "__main__":
s1 = pd.Series(['X1' , 'X1' , 'X2' , 'X2' , 'X2' , 'X2'])
s2 = pd.Series(['Y1' , 'Y1' , 'Y1' , 'Y2' , 'Y2' , 'Y2'])
print('CondEntropy:',getCondEntropy(s1, s2))
print('EntropyGain:' , getEntropyGain(s1, s2))
print('EntropyGainRadio' , getEntropyGainRadio(s1 , s2))
print('DiscreteCorr:' , getDiscreteCorr(s1, s1))
print('Gini' , getGini(s1, s2))
来源:https://blog.csdn.net/weixin_38314865/article/details/84073484
标签:Python,信息熵
![](/images/zang.png)
![](/images/jiucuo.png)
猜你喜欢
scrapy-redis源码分析之发送POST请求详解
2021-05-19 05:24:03
简化翻页提升可用性
2009-03-28 11:52:00
![](https://img.aspxhome.com/file/UploadPic/20093/28/wp_pagenavi_1_en-28s.png)
Python创建对称矩阵的方法示例【基于numpy模块】
2021-01-03 07:35:59
![](https://img.aspxhome.com/file/2023/6/96956_0s.jpg)
FFrpc python客户端lib使用解析
2023-09-06 00:12:20
pandas中merge()函数的用法解读
2023-10-02 08:49:46
扩展性很好的一个分页存储过程分享
2011-11-03 17:04:16
python皮尔逊相关性数据分析分析及实例代码
2021-03-12 13:23:34
![](https://img.aspxhome.com/file/2023/7/92197_0s.jpg)
如何“看看”服务器上的那个文件还在不在?
2009-11-02 20:17:00
浅谈python中对于json写入txt文件的编码问题
2022-01-28 05:08:58
Python中 map()函数的用法详解
2021-03-25 08:16:44
![](https://img.aspxhome.com/file/2023/0/73000_0s.png)
Python自然语言处理 NLTK 库用法入门教程【经典】
2022-10-06 04:36:54
python字典的遍历3种方法详解
2022-05-01 06:00:44
![](https://img.aspxhome.com/file/2023/6/79646_0s.png)
Python可视化神器pyecharts绘制水球图
2023-10-07 14:19:41
![](https://img.aspxhome.com/file/2023/4/102954_0s.png)
用javascript 转换外部链接样式
2007-09-17 10:48:00
jquery AJAX 三个发送状态 posting, error, success
2010-07-31 18:59:00
Python类和对象基础入门介绍
2022-04-03 12:14:48
python构造IP报文实例
2023-07-10 20:40:51
ASC码对照表
2008-08-07 13:07:00
快速解决jupyter notebook启动需要密码的问题
2023-09-27 23:56:49
![](https://img.aspxhome.com/file/2023/9/97039_0s.jpg)
如何导出python安装的所有模块名称和版本号到文件中
2022-07-26 15:38:59
![](https://img.aspxhome.com/file/2023/9/88839_0s.png)