def extract_data(data_path, feature):
population = nm.load_neurons(data_path)<EOL>feature_data = [nm.get(feature, n) for n in population]<EOL>feature_data = list(chain(*feature_data))<EOL>return stats.optimal_distribution(feature_data)<EOL>
Loads a list of neurons, extracts feature
       and transforms the fitted distribution in the correct format.
       Returns the optimal distribution, corresponding parameters,
       minimun and maximum values.