diff --git a/ampullary.pdf b/ampullary.pdf index 929c8f9..78369cc 100644 Binary files a/ampullary.pdf and b/ampullary.pdf differ diff --git a/ampullary.png b/ampullary.png index 1a782c2..7be27e2 100644 Binary files a/ampullary.png and b/ampullary.png differ diff --git a/cells_suscept.pdf b/cells_suscept.pdf index abde240..a373dc8 100644 Binary files a/cells_suscept.pdf and b/cells_suscept.pdf differ diff --git a/cells_suscept.png b/cells_suscept.png index fcba576..cb58dba 100644 Binary files a/cells_suscept.png and b/cells_suscept.png differ diff --git a/cells_suscept_high_CV.pdf b/cells_suscept_high_CV.pdf index fc85517..202a9cc 100644 Binary files a/cells_suscept_high_CV.pdf and b/cells_suscept_high_CV.pdf differ diff --git a/cells_suscept_high_CV.png b/cells_suscept_high_CV.png index 8f74682..693f17f 100644 Binary files a/cells_suscept_high_CV.png and b/cells_suscept_high_CV.png differ diff --git a/data/nonlin_regime_afe_0.08_nfft_400000_len_500_psd.npy b/data/nonlin_regime_afe_0.08_nfft_400000_len_500_psd.npy new file mode 100644 index 0000000..9ef739b Binary files /dev/null and b/data/nonlin_regime_afe_0.08_nfft_400000_len_500_psd.npy differ diff --git a/data/nonlin_regime_afe_0.08_nfft_400000_len_500_psdf.npy b/data/nonlin_regime_afe_0.08_nfft_400000_len_500_psdf.npy new file mode 100644 index 0000000..a4fe885 Binary files /dev/null and b/data/nonlin_regime_afe_0.08_nfft_400000_len_500_psdf.npy differ diff --git a/nonlin_regime.pdf b/nonlin_regime.pdf index df34a8f..5803fb6 100644 Binary files a/nonlin_regime.pdf and b/nonlin_regime.pdf differ diff --git a/nonlin_regime.png b/nonlin_regime.png index f7c476f..a1707e9 100644 Binary files a/nonlin_regime.png and b/nonlin_regime.png differ diff --git a/susceptibility1.tex b/susceptibility1.tex index 6f0cc0a..95ac6ad 100644 --- a/susceptibility1.tex +++ b/susceptibility1.tex @@ -450,7 +450,7 @@ Theoretical work on leaky integrate-and-fire and conductance-based models sugges \begin{figure*}[t] \includegraphics[width=\columnwidth]{nonlin_regime.pdf} - \caption{\label{fig:nonlin_regime} Nonlinear response of a model P-unit in response to increasing beat amplitudes in a three-fish or two-beat setting. The model used has the cell identifier 2013-01-08-aa (table~\ref{modelparams}). \figitem{A--D}. Top -- stimulus consisting of two beats. The beat frequencies are 30\,Hz (\bone{}) and 130\,Hz (\btwo{}). \btwo{} is equal to the baseline firing rate \fbase{}. The beat contrasts are equal in the respective example, and increase from \panel{A} to \panel{D}. Middle -- spike response of the model P-unit to the stimulus above. Bottom -- power spectrum of the firing rate of this P-unit. Nonlinear effects at \bsum{} (orange marker) and at \bdiff{} (red marker) increase for intermediate contrasts (\panel{B}), decrease for stronger contrasts (\panel{C}) and again increase for very strong contrasts (\panel{D}). \figitem{E} Amplitude of the linear (\bone{} and \btwo{}) and nonlinear (\bdiff{} and \bsum{}) responses of the model P-units plotted for increasing beat contrasts (contrasts equal for both beats). } + \caption{\label{fig:nonlin_regime} Nonlinear response of a model P-unit to increasing beat amplitudes in a three-fish or two-beat setting. The model used has the cell identifier 2013-01-08-aa (table~\ref{modelparams}). \figitem{A--D}. Top -- stimulus consisting of two beats, with different contrasts. The beat frequencies are 30\,Hz (\bone{}) and 130\,Hz (\btwo{}). \btwo{} is equal to the baseline firing rate \fbase{}. The contrasts of both beats are equal in a panel, and increase from \panel{A} to \panel{D}. Middle -- spike response of the model P-unit to the stimulus above. Bottom -- power spectrum of the firing rate of this P-unit model. Nonlinear effects at \bsum{} (orange marker) and at \bdiff{} (red marker) in the power spectrum of the P-unit firing rate increase for intermediate contrasts (\panel{B}), decrease for stronger contrasts (\panel{C}) and again increase for very strong contrasts (\panel{D}). \figitem{E} Amplitude of the linear (\bone{}, \btwo{}) and nonlinear (\bdiff{}, \bsum{}) responses of the model P-units plotted for increasing beat contrasts (contrasts increase equally for both beats). } \end{figure*} @@ -463,7 +463,7 @@ When stimulating the fish with both frequencies, additional peaks appear in the The beat stimuli in the example were strong and partially revealed saturating nonlinearities of the P-units. For weakly nonlinear responses we need to use stimuli of much lower amplitudes. -The response of a P-unit to varying beat amplitudes has been estimated by leaky-integrate-and-fire (LIF) models, fitted to the baseline firing properties of electrophysiologically measured P-units. In the chosen P-units model nonlinear peaks (red and orange markers) appear for intermediate beat stimuli (\subfigrefb{fig:motivation}{B}), decreases for stronger stimuli (\subfigrefb{fig:motivation}{C}) and again emerges for very strong stimuli (\subfigrefb{fig:motivation}{D}). Thus two regimes of nonlinearity can be observed for intermediate and strong beat amplitudes (\subfigrefb{fig:motivation}{E}). +The response of a P-unit to varying beat amplitudes can be modeled by a leaky-integrate-and-fire (LIF) model, fitted to the baseline firing properties an electrophysiologically measured P-unit. In the chosen P-unit model nonlinear peaks (orange marker) appear for intermediate beat contrasts (\subfigrefb{fig:motivation}{B}), decrease for stronger contrasts (\subfigrefb{fig:motivation}{C}) and again emerges for very strong beat contrasts (\subfigrefb{fig:motivation}{D}). Thus two regimes of nonlinearity can be observed for intermediate and strong beat amplitudes (\subfigrefb{fig:motivation}{E}). For this example, we have chosen two specific stimulus (beat) frequencies. For a full characterization of the nonlinear responses, we need to measure the response of the P-units to many different combinations of stimulus frequencies. diff --git a/susceptibility1.tex.bak b/susceptibility1.tex.bak index 6f0cc0a..95ac6ad 100644 --- a/susceptibility1.tex.bak +++ b/susceptibility1.tex.bak @@ -450,7 +450,7 @@ Theoretical work on leaky integrate-and-fire and conductance-based models sugges \begin{figure*}[t] \includegraphics[width=\columnwidth]{nonlin_regime.pdf} - \caption{\label{fig:nonlin_regime} Nonlinear response of a model P-unit in response to increasing beat amplitudes in a three-fish or two-beat setting. The model used has the cell identifier 2013-01-08-aa (table~\ref{modelparams}). \figitem{A--D}. Top -- stimulus consisting of two beats. The beat frequencies are 30\,Hz (\bone{}) and 130\,Hz (\btwo{}). \btwo{} is equal to the baseline firing rate \fbase{}. The beat contrasts are equal in the respective example, and increase from \panel{A} to \panel{D}. Middle -- spike response of the model P-unit to the stimulus above. Bottom -- power spectrum of the firing rate of this P-unit. Nonlinear effects at \bsum{} (orange marker) and at \bdiff{} (red marker) increase for intermediate contrasts (\panel{B}), decrease for stronger contrasts (\panel{C}) and again increase for very strong contrasts (\panel{D}). \figitem{E} Amplitude of the linear (\bone{} and \btwo{}) and nonlinear (\bdiff{} and \bsum{}) responses of the model P-units plotted for increasing beat contrasts (contrasts equal for both beats). } + \caption{\label{fig:nonlin_regime} Nonlinear response of a model P-unit to increasing beat amplitudes in a three-fish or two-beat setting. The model used has the cell identifier 2013-01-08-aa (table~\ref{modelparams}). \figitem{A--D}. Top -- stimulus consisting of two beats, with different contrasts. The beat frequencies are 30\,Hz (\bone{}) and 130\,Hz (\btwo{}). \btwo{} is equal to the baseline firing rate \fbase{}. The contrasts of both beats are equal in a panel, and increase from \panel{A} to \panel{D}. Middle -- spike response of the model P-unit to the stimulus above. Bottom -- power spectrum of the firing rate of this P-unit model. Nonlinear effects at \bsum{} (orange marker) and at \bdiff{} (red marker) in the power spectrum of the P-unit firing rate increase for intermediate contrasts (\panel{B}), decrease for stronger contrasts (\panel{C}) and again increase for very strong contrasts (\panel{D}). \figitem{E} Amplitude of the linear (\bone{}, \btwo{}) and nonlinear (\bdiff{}, \bsum{}) responses of the model P-units plotted for increasing beat contrasts (contrasts increase equally for both beats). } \end{figure*} @@ -463,7 +463,7 @@ When stimulating the fish with both frequencies, additional peaks appear in the The beat stimuli in the example were strong and partially revealed saturating nonlinearities of the P-units. For weakly nonlinear responses we need to use stimuli of much lower amplitudes. -The response of a P-unit to varying beat amplitudes has been estimated by leaky-integrate-and-fire (LIF) models, fitted to the baseline firing properties of electrophysiologically measured P-units. In the chosen P-units model nonlinear peaks (red and orange markers) appear for intermediate beat stimuli (\subfigrefb{fig:motivation}{B}), decreases for stronger stimuli (\subfigrefb{fig:motivation}{C}) and again emerges for very strong stimuli (\subfigrefb{fig:motivation}{D}). Thus two regimes of nonlinearity can be observed for intermediate and strong beat amplitudes (\subfigrefb{fig:motivation}{E}). +The response of a P-unit to varying beat amplitudes can be modeled by a leaky-integrate-and-fire (LIF) model, fitted to the baseline firing properties an electrophysiologically measured P-unit. In the chosen P-unit model nonlinear peaks (orange marker) appear for intermediate beat contrasts (\subfigrefb{fig:motivation}{B}), decrease for stronger contrasts (\subfigrefb{fig:motivation}{C}) and again emerges for very strong beat contrasts (\subfigrefb{fig:motivation}{D}). Thus two regimes of nonlinearity can be observed for intermediate and strong beat amplitudes (\subfigrefb{fig:motivation}{E}). For this example, we have chosen two specific stimulus (beat) frequencies. For a full characterization of the nonlinear responses, we need to measure the response of the P-units to many different combinations of stimulus frequencies.