docs: def extract_data(data_path, feature):
      '''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.
      '''
      population = nm.load_neurons(data_path)
  
      feature_data = [nm.get(feature, n) for n in population]
      feature_data = list(chain(*feature_data))
  
      return stats.optimal_distribution(feature_data)