基於鄰域粗糙集演算法python實現
還在跑程式,我這渣筆記本CPU估計要跑很久,不如來貼個程式碼吧。
這個是之前一篇粗糙集相對屬性約簡python實現的一個小小改進
當中距離函式使用的是歐幾里得距離公式,入坑沒多久,程式碼寫的不是很好,請見諒。
# _*_coding:utf-8 _*_
#@Time :2018/12/22 上午10:45
#@Author :we2swing
#@FileName: test.py
import pandas as pd
import numpy as np
import math
import matplotlib.pyplot as plt
def basic_set(df):
basic = {}
for i in df.drop_duplicates().values.tolist():
basic[str(i)] = []
for j, k in enumerate(df.values.tolist()):
if k == i:
basic[str(i)].append(j)
return basic
def Euclidean(x_data):
n = x_data.shape[1]#代表列數
m = x_data.shape[0]#代表行數
num = [0] * n
num2 = np.zeros((m,m))
#轉為list
arr = np.array(x_data)
counti,countj,countk,countz = 0,0,0,0
#遍歷每一行的元素
for j in arr:
countk = 0
for k in arr:
if countj != countk:
countz = 0
# 遍歷每一行的每一個元素
for z in k:
temp = np.square(z - j[countz])
if countz == 0:
num[countz] = temp
else:
#把每一行中每個元素的差值平方相加
num[countz] = num[countz-1] + temp
countz = countz + 1
else:
num[countz-1] = 0
num2[countj][countk] = round(math.sqrt(num[countz-1]),1)
countk = countk + 1
countj = countj + 1
return num2
def key_basic(num,k):
# k = 4 # 設δ=k
# k = [4.0,4.1,4.2,4.3,4.4,4.5,4.6,4.7,4.8,4.9,5.0]
num1 = {}
counti, countj = 0, 0
for i in num:
num1[counti] = []
countj = 0
for j in i:
if j <= k:
num1[counti].append(countj)
countj = countj + 1
counti = counti + 1
return num1
def compare(data):
# k = [4.0, 4.1, 4.2, 4.3, 4.4, 4.5, 4.6, 4.7, 4.8, 4.9, 5.0]
global t1
t1 = 4.0
data = data.dropna(axis=0, how='any')
x_data = data.drop(['judge'], axis=1)
y_data = data.loc[:, 'judge']
n = x_data.shape[1] # 代表列數
m = x_data.shape[0] # 代表行數
# 決策屬性基本集
y_basic_set = sorted([v for k, v in basic_set(y_data).items()])
num = Euclidean(x_data)
# print(num)
x_basic_set = [v for k, v in key_basic(num,t1).items()]
# γC(D)
pos = []
for i in x_basic_set:
for j in y_basic_set:
if set(i).issubset(j):
pos.append(i)
# pos.sort()#排序
r_x_y = len(pos) / len(data) # 也許可以寫一個card函式
print('依賴度r_x_(y):', r_x_y)
# 計算每一個γ(C-a)(D)
# 獲取條件屬性個數為總下標存入集合
# 收集屬性重要性
imp_attr = []
columns_num = list(range(len(x_data.columns)))
for i in columns_num:
c = columns_num.copy()
c.remove(i)
u = data.iloc[:, c]
num = Euclidean(u)
u = sorted([v for k, v in key_basic(num,t1).items()])
# γC(D)
pos_va = []
for k in u:
for j in y_basic_set:
if set(k).issubset(j):
pos_va.append(k)
r_x_y_2 = len(pos_va) / len(data)
r_diff = round((r_x_y - r_x_y_2), 4)
imp_attr.append(r_diff)
dict_imp = {}
for o, p in enumerate(imp_attr):
dict_imp[data.columns[o]] = p
result = dict_imp
# print(imp_attr)
return result
def sepdata(data):
# 獲取資料長度
len = data.iloc[:, 0].size
# 將資料劃分
if len % 100 != 0:
if len > 100:
num = len // 100 + 1
else:
num = 1
else:
if len > 100:
num = int(len / 100)
else:
num = 1
arr = [[]] * num
count = 0
for i in arr:
# 如果數少於100或者最後一部分數少於100,則放入一個由數長決定的陣列
if num == 1:
arr[count] = data.iloc[0:len] # 取100開始,取
elif count == num - 1:
arr[count] = data.iloc[100 * count:len]
else:
arr[count] = data.iloc[100 * count:(count + 1) * 100]
count = count + 1
sorted_dict_imp = [[]] * num
total = [0] * 27
title = ['C1', 'C2', 'C3', 'C4', 'C5', 'C6', 'C7', 'C8', 'C9', 'C10', 'C11', 'C12', 'C13', 'C14', 'C15', 'C16',
'C17', 'C18', 'C19', 'C20', 'C21', 'C22', 'C23', 'C24', 'C25', 'C26', 'C27']
# total = [0] * 16
# title = ['C1', 'C2', 'C3', 'C4', 'C5', 'C6', 'C7', 'C8', 'C9', 'C10', 'C11', 'C12', 'C13', 'C14', 'C15', 'C16']
print(sorted_dict_imp)
count = 0
for i in arr:
print('-------------------------------------第%d個資料集資料-----------------------------------------' % (count + 1))
sorted_dict_imp[count] = compare(i)
count = count + 1
count1 = 0
# 將dict的key為C1-Cn的value存入total中儲存,並且相加
for i in sorted_dict_imp:
count = 0
if count1 == 0:
for j in title:
total[count] = i.get(j)
count = count + 1
else:
for z in title:
total[count] = i.get(z) + total[count]
count = count + 1
count1 = count1 + 1
# 輸出最終C1-Cn的結果
count = 0
for i in title:
print(i, ':', round(total[count], 4))
count = count + 1
figure(total)
def figure(opt):
count = 0
plt.rcParams['font.sans-serif'] = ['SimHei'] # 用來正常顯示中文標籤
title = ['C1', 'C2', 'C3', 'C4', 'C5', 'C6', 'C7', 'C8', 'C9', 'C10', 'C11', 'C12', 'C13', 'C14', 'C15', 'C16',
'C17', 'C18', 'C19', 'C20', 'C21', 'C22', 'C23', 'C24', 'C25', 'C26', 'C27']
plt.scatter(title, opt, c='b', s=7)
plt.plot(title, opt, linewidth=3, alpha=0.7)
plt.show()
def main():
data = pd.read_csv(filepath_or_buffer='/Users/we2swing/Desktop/111.csv')
sepdata(data)
if __name__ == '__main__':
main()
裡面很多程式碼延續的是之前上面傳送門裡文章的程式碼。
相關文章
- 推薦系統實踐 0x06 基於鄰域的演算法(1)演算法
- 推薦系統實踐 0x07 基於鄰域的演算法(2)演算法
- k近鄰演算法python實現 -- 《機器學習實戰》演算法Python機器學習
- 用Python實現K-近鄰演算法Python演算法
- Python 圖_系列之基於鄰接炬陣實現廣度、深度優先路徑搜尋演算法Python演算法
- 基於Python實現MapReducePython
- 《機器學習實戰》k最近鄰演算法(K-Nearest Neighbor,Python實現)機器學習演算法RESTPython
- k近鄰演算法的實現:kd樹演算法
- 基於BKDRhash實現Hash演算法演算法
- 基於Python的Akka實現Python
- Python實現一條基於POS演算法的區塊鏈Python演算法區塊鏈
- Python基於TCP實現聊天功能PythonTCP
- Python:K-近鄰演算法Python演算法
- 機器學習實戰----k值近鄰演算法(Python語言)機器學習演算法Python
- 排序演算法總結(基於Java實現)排序演算法Java
- 基於 Redis 實現基本搶紅包演算法Redis演算法
- 基於Python的Xgboost模型實現Python模型
- python基於opencv 實現影像時鐘PythonOpenCV
- K-近鄰演算法介紹與程式碼實現演算法
- 基於Vue實現圖片在指定區域內移動Vue
- 基於Gossip流言演算法實現alertmanager高可用Go演算法
- Python圖片驗證碼降噪 — 8鄰域降噪Python
- 【機器學習演算法實現】logistic迴歸__基於Python和Numpy函式庫機器學習演算法Python函式
- 教你如何基於python實現指令碼加密Python指令碼加密
- python實現基於smtp傳送郵件Python
- 基於Python實現的口罩佩戴檢測Python
- 基於canvas剪輯區域功能實現橡皮擦效果Canvas
- 機器學習經典分類演算法 —— k-近鄰演算法(附python實現程式碼及資料集)機器學習演算法Python
- K均值演算法基於CUDA環境的實現演算法
- 基於令牌桶演算法實現一個限流器演算法
- 基於Python的tf-idf演算法實現:以《笑傲江湖》為例Python演算法
- 基於面繪製的MC演算法以及基於體繪製的 Ray-casting 實現Dicom影像的三維重建(python實現)演算法ASTPython
- 基於PHP + TRIE樹實現敏感詞過濾演算法PHP演算法
- 如何基於MindSpore實現萬億級引數模型演算法?模型演算法
- 基於MeanShift的目標跟蹤演算法、實現演算法
- 基於C++和OpenCv的SIFT_影象區域性特徵檢測演算法程式碼的實現C++OpenCV特徵演算法
- 機器學習實戰筆記-k近鄰演算法機器學習筆記演算法
- k-鄰近演算法實現約會網站的配對效果演算法網站