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()