# Print out burst statistic information for key, val in burst_stats.items(): print('{:15} \t: {}'.format(key, val)) ################################################################################################### # Burst detection applied to real recordings # ------------------------------------------ # # Next up, we'll load a sample of real neural data, and try out the burst detection. # ################################################################################################### # Download, if needed, and load example data file sig = load_ndsp_data('sample_data_1.npy', folder='data') # Set sampling rate, and create a times vector for plotting fs = 1000 times = create_times(len(sig)/fs, fs) ################################################################################################### # Set the frequency range to look for bursts f_range = (8, 12) # Detect bursts using the dual threshold algorithm bursting = detect_bursts_dual_threshold(sig, fs, (3, 3), f_range) ###################################################################################################
# Import utilities for loading and plotting data from neurodsp.utils import create_times from neurodsp.utils.download import load_ndsp_data from neurodsp.plts.time_series import plot_time_series, plot_instantaneous_measure ################################################################################################### # Load example neural signal # -------------------------- # # First, we will load an example neural signal to use for our time-frequency measures. # ################################################################################################### # Load a neural signal, as well as a filtered version of the same signal sig = load_ndsp_data('sample_data_1.npy', folder='data') sig_filt_true = load_ndsp_data('sample_data_1_filt.npy', folder='data') # Set sampling rate, and create a times vector for plotting fs = 1000 times = create_times(len(sig) / fs, fs) # Set the frequency range to be used f_range = (13, 30) ################################################################################################### # # Throughout this example, we will use # :func:`~neurodsp.plts.time_series.plot_time_series` to plot time series, and # :func:`~neurodsp.plts.time_series.plot_instantaneous_measure` # to plot instantaneous measures.