return nodes def maxCC(ccs): max = 0 for cc in ccs: if len(cc) > max: max = len(cc) return max graph = None if opts.network.endswith(".sif"): graph = convert.read_sif(open(opts.network)) else: graph = convert.read_spf(open(opts.network)) heats = parseHeats(opts.heats) if opts.cut_graph and not opts.node_attributes: cut_val = None try: cutoff = float(opts.cut_graph) except: raise Exception( "Error: value supplied to cut_graph must be a positive number") cutG = cutGraph(graph, heats, cutoff) convert.write_sif(cutG, sys.stdout) sys.exit(0)
current = update return current if __name__ == "__main__": parser = argparse.ArgumentParser() parser.add_argument("network", default=None) parser.add_argument("values", default=None) parser.add_argument("--cycles", type=int, default=10) parser.add_argument("--timestep", type=float, default=0.0001) parser.add_argument("--header", action="store_true", default=False) args = parser.parse_args() handle = open(args.network) gr = networkx.Graph(convert.read_sif(handle)) handle.close() val = read_vec_file(args.values) current = {} for n in gr.node: current[n] = val.get(n, None) found = True while found: update = {} found = False for n in current: if current[n] is not None: update[n] = current[n]
nodes.add(target) return nodes def maxCC(ccs): max = 0 for cc in ccs: if len(cc) > max: max = len(cc) return max graph = None if opts.network.endswith(".sif"): graph = convert.read_sif(open(opts.network)) else: graph = convert.read_spf(open(opts.network)) heats = parseHeats(opts.heats) if opts.cut_graph and not opts.node_attributes: cut_val = None try: cutoff = float(opts.cut_graph) except: raise Exception("Error: value supplied to cut_graph must be a positive number") cutG = cutGraph(graph, heats, cutoff) convert.write_sif(cutG, sys.stdout) sys.exit(0)