【发布时间】:2017-12-19 19:10:58
【问题描述】:
我正在使用以下代码分析加权化学反应网络:
import networkx as nx
import pandas as pd
edge_data = pd.read_table('VULCAN 800.dat', sep=',')
edge_list=[]
edge_data = edge_data.fillna(value=0.1)
col1=edge_data['Column_1']
col2=edge_data['Column_2']
col3=edge_data['Column_3']
col4=edge_data['Column_4']
col5=edge_data['Column_5']
col6=edge_data['Column_6']
for i in range(1,560):
edge_list.append((col1.iloc[i],col2.iloc[i],col5.iloc[i]))
if pd.isnull(col3.iloc[i]) != True:
edge_list.append((col1.iloc[i],col3.iloc[i],col5.iloc[i]))
edge_list.append((col2.iloc[i],col3.iloc[i],col5.iloc[i]))
if pd.isnull(col4.iloc[i]) != True:
edge_list.append((col1.iloc[i],col4.iloc[i],col5.iloc[i]))
edge_list.append((col2.iloc[i],col4.iloc[i],col5.iloc[i]))
if pd.isnull(col3.iloc[i]) != True and pd.isnull(col4.iloc[i]) != True:
edge_list.append((col3.iloc[i],col4.iloc[i],col5.iloc[i]))
G=nx.Graph()
G.add_weighted_edges_from(edge_list)
G.remove_node(0.1)
nx.write_graphml(G, '/home/tessa/Git/Network_biosignatures/VULCAN 800k.graphml')
我想将结果图导出为 graphml 文件,以便在 Cytoscape 中对其进行可视化,但由于某种原因,它会抛出错误消息:
Traceback (most recent call last):
File "topo_measure_pn_hot jupiter_weighted.py", line 196, in <module>
nx.write_graphml(G, '/home/tessa/Git/Network_biosignatures/VULCAN 800k.graphml')
File "<decorator-gen-202>", line 2, in write_graphml
File "/usr/lib/python2.7/dist-packages/networkx/utils/decorators.py", line 220, in _open_file
result = func(*new_args, **kwargs)
File "/usr/lib/python2.7/dist-packages/networkx/readwrite/graphml.py", line 82, in write_graphml
writer.add_graph_element(G)
File "/usr/lib/python2.7/dist-packages/networkx/readwrite/graphml.py", line 351, in add_graph_element
self.add_edges(G,graph_element)
File "/usr/lib/python2.7/dist-packages/networkx/readwrite/graphml.py", line 325, in add_edges
self.add_attributes("edge", edge_element, data, default)
File "/usr/lib/python2.7/dist-packages/networkx/readwrite/graphml.py", line 300, in add_attributes
scope=scope, default=default_value)
File "/usr/lib/python2.7/dist-packages/networkx/readwrite/graphml.py", line 288, in add_data
'%s as data values.'%element_type)
networkx.exception.NetworkXError: GraphML writer does not support <type 'numpy.float64'> as data values.
显然,graphml 编写器不能很好地处理除了简单的图形属性之外的任何事情,我猜是它的权重引发了错误。但是,我不确定 在哪里将权重值分配为 numpy.float64,或者如何将其从该分配中剥离。
正在读入的数据格式如下:
1,H,H2O,OH,H2,1.75116588E-16,,,,,,,,,,,,,,,,,,,,,,,,,,,
2,H,H2O,OH,H2,4.00975292E-13,,,,,,,,,,,,,,,,,,,,,,,,,,,
3,O,H2,OH,H,9.25180896E-14,,,,,,,,,,,,,,,,,,,,,,,,,,,
4,O,H2,OH,H,1.04176774E-13,,,,,,,,,,,,,,,,,,,,,,,,,,,
5,O,H2O,OH,OH,1.04560994E-15,,,,,,,,,,,,,,,,,,,,,,,,,,,
6,O,H2O,OH,OH,2.6959031E-12,,,,,,,,,,,,,,,,,,,,,,,,,,,
感谢任何和所有建议。谢谢!
【问题讨论】:
-
大概你的重量是浮动的?如果没有访问示例数据集,就无法确定,因此很难调试。
-
哦,对了!这是一个示例:1,H,H2O,OH,H2,1.75116588E-16,,,,,,,,,,,,,,,,,,,,,,,,,,,, 1,H,H2O ,OH,H2,4.00975292E-13,,,,,,,,,,,,,,,,,,,,,,,,,,,, 3,O,H2,OH,H,9.25180896E-14 ,,,,,,,,,,,,,,,,,,,,,,,,,,, 3,O,H2,OH,H,1.04176774E-13,,,,,,,,,, ,,,,,,,,,,,,,,,,,, 5,O,H2O,OH,OH,1.04560994E-15,,,,,,,,,,,,,,,,,,, ,,,,,,,,, 5,O,H2O,OH,OH,2.6959031E-12,,,,,,,,,,,,,,,,,,,,,,,,,,,,
-
if 语句的缩进级别是否符合预期?
标签: python numpy networkx graphml