數據分析之PageRank:希拉里郵件門事件,如何計算每個人物在郵件關系網絡中的影響力?

我們用pagerank來查看,各個點的影響力,以及用nexworkx中的graph來展示其關系。

如圖ABCD代表著4個網頁,其中箭頭指向表示,從某網頁跳轉到另一個網頁,那么我們來計算一下,這4個網頁的影響力是多少吧!

圖片解釋

pagerank算法原理就是:一個網頁的影響力 = 所有入鏈集合的頁面的加權影響力之和。

如何計算,公式自己去網絡上查找,我們用代碼演示一遍。

import matplotlib.pyplot as plt
import networkx as nx
G = nx.DiGraph()
# 有向圖之間邊的關系
edges = [("A", "B"), ("A", "C"), ("A", "D"), ("B", "A"), ("B", "D"), ("C", "A"), ("D", "B"), ("D", "C")]
for edge in edges:
    G.add_edge(edge[0], edge[1])
pagerank_list = nx.pagerank(G, alpha=1)#alpha為阻尼因子,alpha=0.85表示跳轉率為15%
print("pagerank 值是:", pagerank_list)
nx.set_node_attributes(G,values=pagerank_list,name='pagerank_list')

# 繪制圓環圖像
positions=nx.circular_layout(G)
#繪制反射圖像
# positions=nx.spring_layout(G)
# 繪制節點
nx.draw_networkx_nodes(G, positions,  alpha=0.4)
# 繪制邊
nx.draw_networkx_edges(G, positions,  alpha=0.2)
# 繪制節點的 label
nx.draw_networkx_labels(G, positions, font_size=10)
plt.show()

得到的結果:


Figure_1.png

點越是大,影響力越大,當數據夠大時會顯現出來。

pagerank 值是: {'A': 0.33333396911621094, 'B': 0.22222201029459634, 'C': 0.22222201029459634, 'D': 0.22222201029459634}

接下來我們看下如何利用pagerank和graph來展示希拉里的郵件。
一張圖顯示步驟:

307055050e005ba5092028a074a5c142.png
1、提取收件名和發件名

首先我們需要將所有的文件提取出來,轉化為字典模式,id對應姓名和真實身份,最后提取的emails.MetadataTo和emails.MetadataFrom要對應其真實身份。

import pandas as pd
import numpy as np
from collections import defaultdict
emails = pd.read_csv("E:/數據學習網站/PageRank-master/input/Emails.csv")
file = pd.read_csv("E:/數據學習網站/PageRank-master/input/Aliases.csv")
aliases={}
for index,row in file.iterrows():
    aliases[row['Alias']]=row['PersonId']
# print(aliases)

file = pd.read_csv("E:/數據學習網站/PageRank-master/input/Persons.csv")
persons={}
for index,row in file.iterrows():
    persons[row['Id']]=row['Name']

# 單詞間的轉換
def unify_name(name):
#大小寫,刪除標點和后綴
    name =str(name).lower()
    name=name.replace(",","").split("@")[0]
#通過key進行名字轉化,提取aliases名字對應的唯一身份,中間通過相同keys連接。
#如果名字再aliases的keys中所對應的的name一樣,那么返回persons中aliases,name的keys,所對應的name。返回name.
    if name in aliases.keys():
        return persons[aliases[name]]
    return name

#提取收件和發件名,apply()代表執行函數,用前面的文字代入參數
emails.MetadataTo=emails.MetadataTo.apply(unify_name)
emails.MetadataFrom=emails.MetadataFrom.apply(unify_name)

2、圖像的準備函數

利用nexworkx畫圖:邊,點和標簽的大小,圖像可是化

#spring_layout為放射性圖像
#circular_layout為環形圖像
def show_graph(graph,layout='spring_layout'):
    if layout=="circular_layout":
        position=nx.circular_layout(graph)
    else:
        position=nx.spring_layout(graph)
#node點,edge為邊,lable為標簽
#nodesize點大小設置,edgesize大小設置
    nodesize=[x['pagerank']*20000 for v,x in graph.node(data=True)]
    edgesize=[np.sqrt(e[2]['weight']) for e in graph.edges(data=True)]
#開始畫邊點和標簽
    nx.draw_networkx_nodes(graph,position,node_size=nodesize,alpha=0.4)
    nx.draw_networkx_edges(graph,position,edge_size=edgesize,alpha=0.2)
    nx.draw_networkx_labels(graph,position,font_size=10)
    plt.show()
3、權重的設置與遍歷

設置dges_weights
1、用for遍歷emails.MetadataTo,emails.MetadataFrom,得到temp名字('a','b')
2、將用if not in條件句,將temp導入edges_weights_temp成為defaultdict(<class 'list'>, {('a', 'b'): 1,('c','d'),2……}將所有temp生存隨機字典。
3、將字典的items,dict_items([(('Jake Sullivan', 'Hillary Clinton'), 815),……,利用[key[0],key[1],val轉化為[('Jake Sullivan', 'Hillary Clinton', 815)……list模式

dges_weights_temp=defaultdict()
for row in zip(emails.MetadataTo,emails.MetadataFrom,emails.RawText):
    temp=(row[0],row[1])
    if temp not in edges_weights_temp:
        edges_weights_temp[temp]=1
    else:
        edges_weights_temp[temp]=edges_weights_temp[temp]+1

edges_weight=[(key[0],key[1],val) for key,val in edges_weights_temp.items()]

4、最后利用pagerank和graph設置

先用graph=nx.DiGraph(),再用add代入權重
接著代入nx.pagerank
用nx.set設置參數
運行畫圖函數

graph=nx.DiGraph()
graph.add_weighted_edges_from(edges_weight)
pagerank=nx.pagerank(graph)
nx.set_node_attributes(graph,name='pagerank',values=pagerank)
show_graph(graph)

數據資料:鏈接:https://pan.baidu.com/s/14edTkp4TH9EtkWt5KYhrrg
提取碼:ui5t
復制這段內容后打開百度網盤手機App,操作更方便哦!

最后是完整代碼:

import pandas as pd
import numpy as np
import networkx as nx
import matplotlib.pyplot as plt
from collections import defaultdict

emails = pd.read_csv("E:/數據學習網站/PageRank-master/input/Emails.csv")
file = pd.read_csv("E:/數據學習網站/PageRank-master/input/Aliases.csv")
aliases={}
for index,row in file.iterrows():
    aliases[row['Alias']]=row['PersonId']
# print(aliases)

file = pd.read_csv("E:/數據學習網站/PageRank-master/input/Persons.csv")
persons={}
for index,row in file.iterrows():
    persons[row['Id']]=row['Name']

# 單詞間的轉換
def unify_name(name):
    name =str(name).lower()
    name=name.replace(",","").split("@")[0]
    if name in aliases.keys():
        return persons[aliases[name]]
    return name

def show_graph(graph,layout='spring_layout'):
    if layout=="circular_layout":
        position=nx.circular_layout(graph)
    else:
        position=nx.spring_layout(graph)
    nodesize=[x['pagerank']*20000 for v,x in graph.node(data=True)]
    edgesize=[np.sqrt(e[2]['weight']) for e in graph.edges(data=True)]
    nx.draw_networkx_nodes(graph,position,node_size=nodesize,alpha=0.4)
    nx.draw_networkx_edges(graph,position,edge_size=edgesize,alpha=0.2)
    nx.draw_networkx_labels(graph,position,font_size=10)
    plt.show()

emails.MetadataTo=emails.MetadataTo.apply(unify_name)
emails.MetadataFrom=emails.MetadataFrom.apply(unify_name)

edges_weights_temp=defaultdict()
for row in zip(emails.MetadataTo,emails.MetadataFrom,emails.RawText):
    temp=(row[0],row[1])
    if temp not in edges_weights_temp:
        edges_weights_temp[temp]=1
    else:
        edges_weights_temp[temp]=edges_weights_temp[temp]+1

edges_weight=[(key[0],key[1],val) for key,val in edges_weights_temp.items()]

graph=nx.DiGraph()
graph.add_weighted_edges_from(edges_weight)
pagerank=nx.pagerank(graph)
nx.set_node_attributes(graph,name='pagerank',values=pagerank)
show_graph(graph)

pagerank_threshold=0.005
small_graph=graph.copy()
for n,p_rank in graph.node(data=True):
    if p_rank["pagerank"]<pagerank_threshold:
        small_graph.remove_node(n)
show_graph(small_graph,"circular_layout")

結果:
看不懂?

美事我們設置閥值,節點PR值,將小于0.005的值全部刪除,這樣我們就得到了一張更加清晰的圖。


郵件門關系圖

這位叫nan的伙計,你是干什么的?為什么天天給美國大美女希拉里天天發郵件?

總結:
1、這里運用到了,nexworkx包,G=nx.DiGraph()來畫圖,用add_weighted_edges_from()添加信息為[('A', 'B', 1),('c','d',2)]格式的數據。
然后,nex.pagerank()添加權重信息,最后set進G里。

2、數據清洗,占據了整個代碼中的一半時間。我們要知道如何利用字典的方式把名字統一在一起,x1=x2—{'a':x1}—{'b':x2},我們要做到輸入A找到B。
最后,從字典中取出數據轉化成[('A', 'B', 1),('c','d',2)]格式。

最后編輯于
?著作權歸作者所有,轉載或內容合作請聯系作者
平臺聲明:文章內容(如有圖片或視頻亦包括在內)由作者上傳并發布,文章內容僅代表作者本人觀點,簡書系信息發布平臺,僅提供信息存儲服務。

推薦閱讀更多精彩內容