Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

GPe fix and params tuning (WIP) #512

Open
wants to merge 7 commits into
base: master
Choose a base branch
from
Open

GPe fix and params tuning (WIP) #512

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
7 commits
Select commit Hold shift + click to select a range
c821577
- fix GPe bug
gabrevaya Dec 23, 2024
0707fc4
add parameter tuning script
gabrevaya Dec 24, 2024
e14cc4a
- make powerspectrum dispatch for EnsembleSolution type-stable
gabrevaya Dec 27, 2024
d2b4ae3
refactor parameter tuning script to adopt a better design pattern and…
gabrevaya Dec 28, 2024
1c6ef90
merge branch 'master' into GPe_fix_and_params_tuning
gabrevaya Dec 30, 2024
c150d71
- use Mason's unroll macro in order to achieve full type stability
gabrevaya Dec 30, 2024
3cdc6e3
Merge branch 'master' of into GPe_fix_and_params_tuning
gabrevaya Dec 30, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
24 changes: 12 additions & 12 deletions docs/src/tutorials/basal_ganglia.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -55,10 +55,10 @@ rasterplot(msn, sol, threshold = -35, title = "Neuron's Spikes - Mean Firing Rat
# Compute and plot the power spectrum of the GABAa current
fig = Figure(size = (1500, 500))

powerspectrumplot(fig[1,1], msn, sol, state = "G",
powerspectrumplot(fig[1,1], msn, sol, state = "I_syn_msn",
title = "FFT with no window")

powerspectrumplot(fig[1,2], msn, sol, state = "G",
powerspectrumplot(fig[1,2], msn, sol, state = "I_syn_msn",
method = welch_pgram, window = hanning,
title = "Welch's method + Hanning window")
fig
Expand Down Expand Up @@ -105,12 +105,12 @@ fig = Figure(size = (1000, 800))
rasterplot(fig[1,1], msn, ens_sol[1], threshold = -35, title = "MSN - Mean Firing Rate: $(round(fr_msn[1], digits=2)) spikes/s")
rasterplot(fig[1,2], fsi, ens_sol[1], threshold = -35, title = "FSI - Mean Firing Rate: $(round(fr_fsi[1], digits=2)) spikes/s")

powerspectrumplot(fig[2,1], msn, ens_sol, state = "G",
powerspectrumplot(fig[2,1], msn, ens_sol, state = "I_syn_msn",
method = welch_pgram, window = hanning,
ylims= (-35, 15),
xlims= (8, 100))

powerspectrumplot(fig[2,2], fsi, ens_sol, state = "G",
powerspectrumplot(fig[2,2], fsi, ens_sol, state = "I_syn_fsi",
method=welch_pgram, window=hanning,
ylims= (-35, 15),
xlims= (8, 100))
Expand Down Expand Up @@ -153,22 +153,22 @@ ens_sol = solve(ens_prob, RKMil(), dt=dt, saveat = dt, trajectories = 3);
# Compute and plot power spectra for all components
fig = Figure(size = (1600, 450))

powerspectrumplot(fig[1,1], msn, ens_sol, state = "G",
powerspectrumplot(fig[1,1], msn, ens_sol, state = "I_syn_msn",
method = welch_pgram, window = hanning,
ylims=(-40, 25),
title = "MSN (Baseline)")

powerspectrumplot(fig[1,2], fsi, ens_sol, state = "G",
powerspectrumplot(fig[1,2], fsi, ens_sol, state = "I_syn_fsi",
method = welch_pgram, window = hanning,
ylims=(-40, 25),
title = "FSI (Baseline)")

powerspectrumplot(fig[1,3], gpe, ens_sol, state = "G",
powerspectrumplot(fig[1,3], gpe, ens_sol, state = "V",
method = welch_pgram, window = hanning,
ylims=(-40, 25),
title = "GPe (Baseline)")

powerspectrumplot(fig[1,4], stn, ens_sol, state = "G",
powerspectrumplot(fig[1,4], stn, ens_sol, state = "V",
method = welch_pgram, window = hanning,
ylims=(-40, 25),
title = "STN (Baseline)")
Expand Down Expand Up @@ -215,22 +215,22 @@ ens_prob = EnsembleProblem(prob)
ens_sol = solve(ens_prob, RKMil(), dt = dt, saveat = dt, trajectories = 3);

# Compute and compare power spectra for all neural populations in Parkinsonian condition against their counterparts in baseline conditions.
powerspectrumplot(fig[2,1], msn, ens_sol, state = "G",
powerspectrumplot(fig[2,1], msn, ens_sol, state = "I_syn_msn",
method = welch_pgram, window = hanning,
ylims=(-40, 25),
title = "MSN (PD)")

powerspectrumplot(fig[2,2], fsi, ens_sol, state = "G",
powerspectrumplot(fig[2,2], fsi, ens_sol, state = "I_syn_fsi",
method = welch_pgram, window = hanning,
ylims=(-40, 25),
title = "FSI (PD)")

powerspectrumplot(fig[2,3], gpe, ens_sol, state = "G",
powerspectrumplot(fig[2,3], gpe, ens_sol, state = "V",
method = welch_pgram, window = hanning,
ylims=(-40, 25),
title = "GPe (PD)")

powerspectrumplot(fig[2,4], stn, ens_sol, state = "G",
powerspectrumplot(fig[2,4], stn, ens_sol, state = "V",
method = welch_pgram, window = hanning,
ylims=(-40, 25),
title = "STN (PD)")
Expand Down
Loading
Loading