tidy script

Author: brendan-ostlund

Python/02-GroupPSDs.ipynb

@@ -225,7 +225,7 @@ m2_freq, m2_spectra = trim_spectrum(freqs, spectra_subsample, m2_PSD_range)
 
 #### Fit models, with different settings
 ```{python, message = FALSE}
-# Fit a model object with model 1 settings
+# Fit model object with model 1 settings
 fg1 = FOOOFGroup(peak_width_limits = m1_peak_width, max_n_peaks = m1_n_peaks, min_peak_height = m1_peak_height)
 fg1.fit(m1_freq, m1_spectra)
 
@@ -237,7 +237,7 @@ for ind in range(len(fg1)):
 ```
 
 ```{python, message = FALSE}
-# Fit a model object with model 2 settings
+# Fit model object with model 2 settings
 fg2 = FOOOFGroup(peak_width_limits = m2_peak_width, max_n_peaks = m2_n_peaks, min_peak_height = m2_peak_height)
 fg2.fit(m2_freq, m2_spectra)
 
@@ -262,6 +262,12 @@ settings2 = FOOOFSettings(peak_width_limits = m2_peak_width, max_n_peaks = m2_n_
                           aperiodic_mode = 'fixed')
 ```
 
+```{python, message = FALSE}
+# Initialize model object for spectral parameterization, with some settings
+fg1 = FOOOFGroup(*settings1)
+fg2 = FOOOFGroup(*settings2)
+```
+
 #### Fit models with group model object
 Note that when fitting power spectra, you can specify a fit range to fit the model to, so you don't have to explicitly trim the spectra. Here we will refit the example data, specifying the fit range, and then save the group reports. 
 ```{python, message = FALSE}

R/03-R_GroupPSDs.Rmd

@@ -39,9 +39,6 @@ from fooof.plts import plot_spectrum
 from fooof.plts.periodic import plot_peak_fits, plot_peak_params
 from fooof.plts.aperiodic import plot_aperiodic_params, plot_aperiodic_fits
 
-# # Import utility functions that manipulate FOOOF objects
-# from fooof.objs.utils import average_fg
-
 # Import functions to examine frequency-by-frequency error of model fits
 from fooof.analysis.error import compute_pointwise_error_fg
 ```