亚洲免费在线-亚洲免费在线播放-亚洲免费在线观看-亚洲免费在线观看视频-亚洲免费在线看-亚洲免费在线视频

python 實現信息熵、條件熵、信息增益、基尼系數

系統 2327 0

注:該代碼為慕課網課程中老師講解

            
              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))
    

            
          

更多文章、技術交流、商務合作、聯系博主

微信掃碼或搜索:z360901061

微信掃一掃加我為好友

QQ號聯系: 360901061

您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描下面二維碼支持博主2元、5元、10元、20元等您想捐的金額吧,狠狠點擊下面給點支持吧,站長非常感激您!手機微信長按不能支付解決辦法:請將微信支付二維碼保存到相冊,切換到微信,然后點擊微信右上角掃一掃功能,選擇支付二維碼完成支付。

【本文對您有幫助就好】

您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描上面二維碼支持博主2元、5元、10元、自定義金額等您想捐的金額吧,站長會非常 感謝您的哦!!!

發表我的評論
最新評論 總共0條評論
主站蜘蛛池模板: 欧美a视频在线观看 | 不卡一级aaa全黄毛片 | 91精品国产综合久久久久 | 国产自产拍精品视频免费看 | 亚洲最大黄色网址 | 蕾丝视频在线看片国产 | 一区免费在线观看 | 欧美一级刺激毛片 | 久久视频精品 | 国产亚洲女在线精品 | 狠狠婷婷| 日本一级毛片片在线播放 | 草草影院1 | 成人日韩欧美 | 日日操日日爽 | 一级毛片免费一级直接观看 | 俄罗斯一级毛片aaaa | 在线观看视频一区二区 | 青草青草久热精品视频在线观看 | 亚洲九九色 | 这里只有精品免费视频 | 国产精品丝袜在线 | 免费观看黄色网 | 国产一区二区三区不卡免费观看 | 国产精品久久久久无码av | 精品欧美一区二区精品久久 | 青青青国产免费手机视频在线观看 | 看真人一级毛片 | 国产欧美一区二区精品仙草咪 | 欧美久久精品一级c片片 | 亚洲天堂区 | 色偷偷亚洲天堂 | 视频二区 中文字幕 欧美 | 欧美大黄大色一级毛片出红 | www国产精品com| 久久久久久久国产高清 | 中文字幕久久久久一区 | 性做久久久久久免费观看 | 国产精品高清全国免费观看 | 日本一级特黄aa毛片免费观看 | 国产情侣普通话刺激对白 |