Hello,
I have a model built in simbiology that I have confirmed to be stable and running with a set of nominal parameters (10+ parameters). Now I am trying to fit one of these parameters from the model with sbiofit using the instructions here, but am running into some issues. Unfortunately I cannot share the model publicly here at the moment, which I realize makes it a bit difficult to help, but I thought I'd ask anyways to see if I'm doing anything that pops out as immediately incorrect to someone.
The general structure of my code is as follows:
T = table(time,conc); % time and conc are equally sized vectors of time-concentration data points
gData = groupedData(T); responseMap = {'compartment.species= conc'}; paramsToEstimate = {'log(p1)'}; % keep p1 positive
p1_i = 0.003; p1_limits = [0.0001 3]; estimatedParams = estimatedInfo(paramsToEstimate,'InitialValue',[p1_i],'Bounds',[p1_limits]); fitConst = sbiofit(model,gData,responseMap,estimatedParams,dose); s.Labels.XLabel = 'Time (hour)'; s.Labels.YLabel = 'Amount (nanomole)'; plot(fitConst,'AxesStyle',s);
This code runs successfully, but the estimate never deviates from my initial value guess specified for p1 (0.003) and the StandardError that gets spit out is consistently infinity ('Inf').
I realize this may be difficult to debug without my exact model posted here, but I was wondering if anyone has any pointers as to what might be causing the fitter to behave like this? I've searched around and haven't been able to find anyone who has had this problem on these forums or the rest of the web as far as I can tell.
Thanks in advance!
Best Answer