def _compute_spect_slope_feat_names(data, **kwargs):
    """Utility function to create feature names compatible with the output of
    :func:`mne_features.univariate.compute_energy_freq_bands`."""
    n_channels = data.shape[0]
    stats = ['intercept', 'slope', 'MSE', 'R2']
    return ['ch%s_%s' % (ch, stat) for ch in range(n_channels)
            for stat in stats]