def write_all_correlation_files(data, definition, def_args):
window_size = 4
for year in data.all_years:
write_correlation_list('adj-y%d.txt' % year, adjacency_matrix(data, definition, def_args, year), 0)
if 1963 + window_size <= year:
write_correlation_list('pn-y%d-to-y%d.txt' % (year - window_size + 1, year),
positives_and_negatives_matrix(data, definition, def_args,
range(year - window_size + 1, year + 1)), 0.5)
def matlab_code_for_rcm_ordered_corr_coef_for_sliding_window_degree_matrix(data, definition, def_args,
normalize_row_or_column):
f = open(OUT_DIR.RCM_MATRIC + 'iraniraqpn.m', 'w')
all_countries = DEFAULT_COUNTRIES_LIST
allowed_countries = iran_iraq_countries
# all_countries = ['USA', 'UK', 'Australia', 'Greece']
# allowed_countries = ['Australia', 'Greece']
window_size = 5
for window_start_year in data.all_years:
sliding_window = range(window_start_year, window_start_year + window_size)
window_end_year = sliding_window[-1:][0]
if window_end_year > 2000: break
pn_matrix = positives_and_negatives_matrix(data, definition, def_args, sliding_window, all_countries,
normalize_row_or_column)
f.write(matrix_py_matlab_with_name('pnmatrix', pn_matrix))
data_matrix = corrcoef(pn_matrix)
for threshold in [0]:
for line in adjacency_rcm_ordered(data_matrix, threshold, [all_countries, allowed_countries],
'%s-%s' % (window_start_year, window_end_year)):
f.write("%s\n" % line)
f.close()
