Skip to content

Conversation

@asoplata
Copy link
Collaborator

@asoplata asoplata commented Nov 7, 2025

I believe this adds the last of the code that is extracted from #1168 that is not the actual Duecker model itself, but instead code changes that are used by the default model. These functions are not currently tested at all, and I have not added tests for them for similar reasons to #1185 (comment)

After this, we should be able to make a PR that adds the actual Duecker model itself. However, there will be additional work for that step: we need to update hnn_io.read_network_configuration and the equivalent ...write... such that they support the Duecker model.

Continuing the pattern, this is based on top of #1187.

"""

return zero_val * np.exp(exp_term * x) + offset
# AES: TODO Should the offset be multiplied by zero_val? This was also changed in
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Yes, I double-checked that this is indeed what I wanted it to do. The distribution of the Ih that I need this function for in the new model is:

gbar = zero_val * (-0.8696 + 2.087 * np.exp(x/323))

Taken from Hay et al., 2011.

Actually, maybe we shouldn't call it zero_val, but gbar_proximal or something like that.

Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This checks out from the picture you sent me over email, and I've removed my comment drawing attention to it in the latest commit. I agree about changing the name zero_val now, since zero_val now multiplies both of the offset and the exponential. What do you think about gbar_original? I fear proximal could make a novice become confused and think this relates to our proximal/distal inputs language.

I believe this adds the last of the code that is extracted from jonescompneurolab#1168
that is *not* the actual Duecker model itself, but instead code changes
that are used by the default model. These functions are not currently
tested at all, and I have not added tests for them for similar reasons
to
jonescompneurolab#1185 (comment)

After this, we should be able to make a PR that adds the actual Duecker
model itself. However, there will be additional work for that step: we
need to update `hnn_io.read_network_configuration` and the equivalent
`...write...` such that they support the Duecker model.

Continuing the pattern, this is based on top of jonescompneurolab#1187.
@asoplata
Copy link
Collaborator Author

Hey Katharina, I've rebased this PR so that it's no longer dependent on #1185 or #1187. I've responded to your comments, and you can still see them on the PR page, but you will not see the specific code-location of your comments with respect to the code, due the rebase.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

Status: Backlog

Development

Successfully merging this pull request may close these issues.

2 participants