MATLAB Answers

0

How to set up estimation method in sbiofit

Asked by Deepa Maheshvare on 26 Apr 2019
Latest activity Edited by 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

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?
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.

Tags

Products

Communities

More Answers in the  SimBiology Community

1 Answer

Answer by 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

Hi Jeremy,
Thanks a lot for the response. I have a few doubts in the first approach.
When I do this ,"You could create a data set with a single time point at tspan" . When I do this ,
would it be incorrect, in the sense, not all species will take the same amount of time to reach steady state. For instance, I have computed transient times of species of all species in the model. The time scale ranges from 0.01 seconds to 1000seconds. When the time scale of equilibration is different for each species, I am not sure how "a data set with a single time point at tspan" can be justified.
I'd be happy to know your views on this.
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.
Thanks a lot for the advice. I'll try it out

Sign in to comment.