docs: def extract_data(neurons, feature, params=None):
      '''Extracts feature from a list of neurons
         and transforms the fitted distribution in the correct format.
         Returns the optimal distribution and corresponding parameters.
         Normal distribution params (mean, std)
         Exponential distribution params (loc, scale)
         Uniform distribution params (min, range)
      '''
      if params is None:
          params = {}
  
      feature_data = [get(FEATURE_MAP[feature], n, **params) for n in neurons]
      feature_data = list(chain(*feature_data))
      return stats.optimal_distribution(feature_data)