def extract_data(neurons, feature, params=None):
if params is None:<EOL><INDENT>params = {}<EOL><DEDENT>feature_data = [get(FEATURE_MAP[feature], n, **params) for n in neurons]<EOL>feature_data = list(chain(*feature_data))<EOL>return stats.optimal_distribution(feature_data)<EOL>
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)