fig, ax = plt.subplots(2, 4, figsize=[40, 20], sharex=True, sharey=True) for i in range(2): j = 0 for PERPLEXITY in PERPLEXITIES: SELECTED_RANGES = ([(float(BOTTOMS[i]), float(TOPS[i]))]) print(j, i, PERPLEXITY, SELECTED_RANGES) tsne_data = np.load( ('/Users/pablonavarrobarrachina/Desktop/Results FFTW/' + 'tSNE_results_range_{}_perplexity_{}_SNRof{}.npy').format( SELECTED_RANGES, PERPLEXITY, SNR)) # ---> Initialize DBSCAN dbscan = dbscan_ratio.dbscan_method(tsne_data[:, 0], tsne_data[:, 1], stellar_parameters['binarity'], stellar_parameters, SELECTED_RANGES, SNR, PERPLEXITY) dbscan.normalize_data_tSNE() dbscan.parameter_space = pd.read_csv( ('/Users/pablonavarrobarrachina/Desktop/' + 'Results FFTW/DBSCAN_parameterspace_range_{}_perplexity_' + '{}_SNRof{}_ratio_{}_iterations_{}.csv').format( SELECTED_RANGES, PERPLEXITY, SNR, dbscan.ratio, dbscan.iterations)) dbscan.get_variables_from_imported_data() mask_singles = stellar_parameters['binarity'] == 0 mask_recovered_binaries = ((dbscan.ratio_labels == 1) & (stellar_parameters['binarity'] == 1)) mask_non_recovered_binaries = (
tsne_data = np.load(dirs.tsne_results + 'tSNE_results_range_{}_perplexity_{}_SNRof{}.npy'.format( SELECTED_RANGES, PERPLEXITY, SNR)) stellar_parameters = pd.read_csv( dirs.data + 'stellar_parameters_duchenekrauspopulation.csv') # DBSCAN parameters ____________________________________________________________ minEpsilon = 0.1 maxEpsilon = 0.75 min_minSamples = 25 max_minSamples = 125 # Run the main DBSCAN dbscan = dbscan_ratio.dbscan_method(tsne_data[:, 0], tsne_data[:, 1], stellar_parameters['binarity'], minEpsilon, maxEpsilon, min_minSamples, max_minSamples, stellar_parameters, SELECTED_RANGES, SNR, PERPLEXITY) FILENAME_DBSCAN_RESULT = ( dirs.dbscan_results + 'DBSCAN_parameterspace_range_{}_perplexity_{}_SNRof{}_ratio_{}_iterations_{}_randomnesscheck_{}.csv' .format(SELECTED_RANGES, PERPLEXITY, SNR, dbscan.ratio, dbscan.iterations, NUMBER_RANDOMNESS)) if os.path.isfile(FILENAME_DBSCAN_RESULT) == True: print('File exists already, skipping this analysis..') exit() else: dbscan.normalize_data_tSNE()