I've created a graph in networkx
and get its bfs tree.
G = nx.Graph()
# build a graph
tree = nx.bfs_tree(G, '1')
Now I would like to save the tree in newick format to a file. What's the best way to do this?
Inspired by this answer, I do something like this:
import networkx as nx
import matplotlib.pyplot as plt
def recursive_search(dict, key):
if key in dict:
return dict[key]
for k, v in dict.items():
item = recursive_search(v, key)
if item is not None:
return item
def bfs_edge_lst(graph, n):
return list(nx.bfs_edges(graph, n))
def load_graph(filename):
G = nx.Graph()
# build the graph
return G
def tree_from_edge_lst(elst, n):
tree = {n: {}}
for src, dst in elst:
subt = recursive_search(tree, src)
subt[dst] = {}
return tree
def tree_to_newick(tree):
items = []
for k in tree.keys():
s = ''
if len(tree[k].keys()) > 0:
subt = tree_to_newick(tree[k])
if subt != '':
s += '(' + subt + ')'
s += k
items.append(s)
return ','.join(items)
g = load_graph('dataset.txt')
elst = bfs_edge_lst(g, '1') #'1' being the root node of the graph
tree = tree_from_edge_lst(elst, '1')
newick = tree_to_newick(tree) + ';'