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()
data = None
data = ExportData()
data.load_file('../' + WORLD_TRADE_FLOW_DATA_FILE_ORIGINAL, should_read_world_datapoints=True)
#generate_matlab_histogram_code(data, [99], [2000], ['Georgia', 'USA'])
print_histogram_as_text(data, range(1963, 2001), data.countries())
#print_graph_densities_for_different_thresholds(data, thresholds, a_few_years)
print_missing_links_db(data, 2000, 90, 'def_d_histogram_pair_wise')
print_missing_links_db(data, 2000, 99, 'def_d_histogram_pair_wise')
from itertools import combinations
from project.config import WORLD_TRADE_FLOW_DATA_FILE_ORIGINAL
from project.export_data.exportdata import ExportData
from project.signed_networks.definitions import definition_C3, args_for_definition_C
data = ExportData()
data.load_file('../../' + WORLD_TRADE_FLOW_DATA_FILE_ORIGINAL, should_read_world_datapoints=True)
year = 1975
definition = definition_C3
args = args_for_definition_C(10, 5000)
for (A, B) in combinations(data.countries(), 2):
print "%s,%s,%s" % (A, B, definition(data, year, A, B, args))
