def print_missing_links_db(data, year, T, log_file_name):
two_way_args = args_for_definition_D(T)
one_way_args = args_for_definition_D(T, mode='one-way')
f = open(OUT_DIR.DEFINITION_D + log_file_name + ".%d.txt" % T, 'w')
for (A, B) in countries.country_pairs(data.countries()):
if definition_D(data, year, A, B, two_way_args) == NO_LINK:
one_way = definition_D(data, year, A, B, one_way_args)
other_way = definition_D(data, year, B, A, one_way_args)
f.write("Y%d,%s,%s,%s,%s\n" % (year, file_safe(A), file_safe(B), one_way, other_way))
if one_way != other_way:
for Y in range(1963, 2001):
f.write("%d,%s,%s,%.4g,%.4g\n" % (
Y, file_safe(A), file_safe(B),
data.export_data_as_percentile(Y, A, B),
data.export_data_as_percentile(Y, B, A)))
f.close()
def print_graph_densities_for_different_thresholds(data, thresholds, years):
f = open(OUT_DIR.DEFINITION_D + 'combinations.txt', 'w')
for T in thresholds:
for year in years:
positive_edges = 0
negative_edges = 0
for (A, B) in countries.country_pairs(data.countries()):
link_sign = definition_D(data, year, A, B, args_for_definition_D(T, log_file=f))
if link_sign == POSITIVE_LINK: positive_edges += 1
if link_sign == NEGATIVE_LINK: negative_edges += 1
N = 203
density = 2.0 * (positive_edges + negative_edges) / (N * (N - 1)) * 100
print "%d,%d,%f,%d,%d,%d" % (T, year, density, positive_edges, negative_edges, N)
f.close()
