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.

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

Can the output be a ratio between state variables like: y(t) = x(t)/(S(t)+E(t)+x(t))? #392

Closed
qchenmath opened this issue Mar 8, 2025 · 4 comments
Labels
question Further information is requested

Comments

@qchenmath
Copy link

For all the examples, the output(s) are always simple function of the state variables. Will it work if the output is as given in the title? Thanks.

@qchenmath qchenmath added the question Further information is requested label Mar 8, 2025
@sumiya11
Copy link
Collaborator

sumiya11 commented Mar 8, 2025

Hi @qchenmath . Yes, ratios in the output should work.

For example,

using StructuralIdentifiability
ode = @ODEmodel(
           x'(t) = a*x,
           y(t) = x^2 // (x + a + 1)
       )
assess_identifiability(ode)

Note the double slash. (edit: perhaps / works as well)

@pogudingleb
Copy link
Collaborator

Thanks, @sumiya11 !
Indeed, single slash for division is fine.

@qchenmath
Copy link
Author

Thanks. Guess the model is somewhat complicated. It didn't finish after 5 days, and I had to kill it.

@pogudingleb
Copy link
Collaborator

@qchenmath
Yes, this may happen, especially with fractions. Feel free to send me the model to look at, I might be able to propose something to speed up the calculations.
Note that you can also always get at least local identifiability info with assess_local_identifiability function.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
question Further information is requested
Projects
None yet
Development

No branches or pull requests

3 participants