How can I calculate the Hessian matrix for fit results with sbiofit (SimBiology)?
5 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
Hi all,
With sbiofit (SimBiology) the fit results give me standard errors obtained from the Jacobian calculation. I would like to compute the standard errors with the Hessiana matrix, but I don't know how to get it analytically from my object function and the fit parameters.
Thank you in advance,
Valentina
0 comentarios
Respuestas (1)
Arthur Goldsipe
el 11 de En. de 2023
I can't think of an easy way to calculate the Hessian matrix. The objective function is not analytical (since it involves approximate numerical solutions of differential equations). It's possible to estimate the Hessian numerically for situations like this, but I worry that it would be difficult to get an accurate estimate. (We already have challenges calculating the Jacobian accurately.)
Have you looked at the different options available for calculating confidence intervals using sbioparameterci or sbiopredictionci? If none of those are satisfactory, can you tell me more about what you're looking for? Can you point me to any references that describe the calculations you're trying to do? I'd like record this in our enhancment request database, and the more details the better.
-Arthur
5 comentarios
Arthur Goldsipe
el 13 de En. de 2023
I'm sorry, but I still don't fully understand what you're trying to do. I think we're at the point where I need to see sample code before I can recommend anything more specific. If you're interested in that, please feel free to share code here or to reach out to me via my profile if you want to share code privately.
I'm familiar with Numerical Recipes in C, but I don't recall anything in it about Bayesian approaches. What I was hoping to get was some sort of reference that describes whatever algorithm or approach you're trying to implement, in sufficient detail that I could think about how to implement it in SimBiology.
Comunidades de usuarios
Más respuestas en SimBiology Community
Ver también
Categorías
Más información sobre Import Data en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!