How to set up estimation method in sbiofit

5 views (last 30 days)
Deepa Maheshvare
Deepa Maheshvare on 26 Apr 2019
Edited: Deepa Maheshvare on 29 Apr 2019
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 Comments
Deepa Maheshvare
Deepa Maheshvare on 26 Apr 2019
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

Sign in to comment.

Answers (1)

Jeremy Huard
Jeremy Huard on 26 Apr 2019
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 Comments
Deepa Maheshvare
Deepa Maheshvare on 29 Apr 2019
Thanks a lot for the advice. I'll try it out

Sign in to comment.

Communities

More Answers in the  SimBiology Community

Tags

Products

Community Treasure Hunt

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

Start Hunting!

Translated by