How to set up estimation method in sbiofit

2 次查看（过去 30 天）
Deepa Maheshvare
Deepa Maheshvare 2019-4-26
编辑： Deepa Maheshvare 2019-4-29
Hi All,
I'm trying to estimate parametrs in my model using parameter estimation task in SimBiology toolbox.In the documentation available here it is mentioned that
"The default estimation method that sbiofit uses will change depending on which toolboxes are available."
From the list of methods available here, I want to select fmincon .
Could someone explain how to set up fmincon in the following syntax?
sbiofit(model,gData,responseMap,estimatedParam,'fmincon')
I want to define an objective function and pass the following arguments to fmincon
fmincon(@objective,p0,A,b,Aeq,beq,lb,ub,nlcon)
I have a confusion here,
p0 is the initial value of parameters, A,b,Aeq,beq are linear equality and inequality constraints, which is (A=[ ],b=[ ],Aeq=[ ],beq=[ ]) in my case.lb and ub are the bounds of paramaters and ncoln is the nonlinear constraints that are the differential equations taht will be generated from my model present in xml file.
I've already specified p0,lb,ub in estimatedParam array with properties Name, InitialValue and Bounds
estimatedParams = estimatedInfo(paramsToEstimate,'InitialValue',ParameterInitialValues,'Bounds',ParameterBounds)
It's not to me how to define
fmincon(@objective,p0,A,b,Aeq,beq,lb,ub,nlcon)
since p0,lb,ub is already defined in estimatedParam array. However, I'd like to know
how to pass the other arguments objective ,A,b,Aeq,beq, nlcon in fmincon.
nlcon function is defined like this
Any suggestions will be highly appreciated.
  2 个评论
Jeremy Huard
Jeremy Huard 2019-4-26
The estimatedInfo object will let you define the parameters initial values and their lower and upper bounds.
But can you explain what you would like to achieve with nlcon?
Deepa Maheshvare
Deepa Maheshvare 2019-4-26
编辑：Deepa Maheshvare 2019-4-26
Hi, Thanks you very much for the response. I want to pass ode's as constraints to the minimization problem as explained here. I'd like to implement it like this. The lb and ub in the inputs of the nonlincon(ncoln) function defined here are the bounds of the variables(Please note: lb and ub in nlcon are not the bounds of parameters that are already specified in estimatedParam) in my model. I'm trying to restrict the values of the solution of the variables at all times to be in a certain range[lb ub] . The range[lb ub] is the only experimental input that I have. [lb ub] is the range of the values of transient variables obtained from different experimental observations . i.e if there are 20 variables in my model, there are 20 differential equations generated form the sbml/xml file and the [lb ub] is available for all 20 variables.These bounds are steady state measurements from experiments.
In addition to this, I have the initial values of the parameters to be estimated and the [lb ub] of the parameters.
I don't have time varying experimental input to assign to gData,responseMap in sbiofit.
I'm not sure how to implemet this in SimBiology Toolbox. If there is a way to generate the ode's from the modelObject generated after reading the sml file using sbmlimport , I would be able to implement the above easily.
Any suggestions will be highly appreciated

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Jeremy Huard
Jeremy Huard 2019-4-26
Hi Deepa,
let me summarize my understanding of your problem and please correct me if I'm wrong:
  • you only have steady state data for several experimental conditions
  • there is no other data to fit the model to
  • you also have parameter bounds
Is this correct?
If yes, I suggest to create a dataset with these steady state data as explained here. If this is experimental data, you will know when this data was measured so you can build a table with time points and measurements.
Then you can use sbiofit to fit your model to those data points. Again, the App will be of great help.
  7 个评论
显示 5更早的评论
Jeremy Huard
Jeremy Huard 2019-4-29
This is what you are currently doing by calling ode15s for [0,tspan] and taking the last data points.
And this is fine in my opinion. All you know from your data is that it was taken at steady state but you don't know exactly when. So all you need to do is take tspan large for everything in your system to be at steady state.
The idea with the data set with a single time point is to reproduce what you actually have, which is a single time point.
Deepa Maheshvare
Deepa Maheshvare 2019-4-29
编辑：Deepa Maheshvare 2019-4-29
Thanks a lot for the advice. I'll try it out

请先登录，再进行评论。

社区

更多回答在  SimBiology Community

类别

在 Help Center 和 File Exchange 中查找有关 Perform Sensitivity Analysis 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by