Mercurial > hg > plosone_underreview

comparison scripts/results.py @ 13:98718fdd8326 branch-tests

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
edits in the core functions
author Maria Panteli <m.x.panteli@gmail.com>
date Tue, 12 Sep 2017 18:03:47 +0100
parents 0f3eba42b425
children
comparison
equal deleted inserted replaced
10:8e897e82af51 13:98718fdd8326
86 # load LDA-transformed frames 86 # load LDA-transformed frames
87 X_list, Y, Yaudio = pickle.load(open('data/lda_data_8.pickle','rb')) 87 X_list, Y, Yaudio = pickle.load(open('data/lda_data_8.pickle','rb'))
88 ddf = load_metadata(Yaudio, metadata_file='data/metadata.csv') 88 ddf = load_metadata(Yaudio, metadata_file='data/metadata.csv')
89 w, data_countries = utils_spatial.get_neighbors_for_countries_in_dataset(Y) 89 w, data_countries = utils_spatial.get_neighbors_for_countries_in_dataset(Y)
90 w_dict = utils_spatial.from_weights_to_dict(w, data_countries) 90 w_dict = utils_spatial.from_weights_to_dict(w, data_countries)
91 Xrhy, Xmel, Xmfc, Xchr = X_list 91 X = np.concatenate(X_list, axis=1)
92 X = np.concatenate((Xrhy, Xmel, Xmfc, Xchr), axis=1)
93 92
94 # global outliers 93 # global outliers
95 df_global, threshold, MD = get_outliers_df(X, Y, chi2thr=0.999) 94 df_global, threshold, MD = get_outliers_df(X, Y, chi2thr=0.999)
96 print_most_least_outliers_topN(df_global, N=10) 95 print_most_least_outliers_topN(df_global, N=10)
97 96