Log likelihood is initially high in FOCE but as the iteration proceeds to greater than 3500 the log likelihood reduces and finally ends with error that maximum iteration has been achieved.

2 次查看（过去 30 天）
Praveen Kumar M
Praveen Kumar M 2019-6-16
评论： Praveen Kumar M 2019-6-18
I am trying to develop a population pharmacokinetic model using Simbiology,
My settings are
Estimation method: Estimating mixed effects model with nlmefit
Error model is proportional
Method to approximate non-linear mixed effect model likelihood is FOCE.
The FOCE model fitting progress shows a loglikelihood of around -200 (that is good) till iteration of 3500 but after that suddenly the loglikelihood goes to around -900.
Then it keeps on proceeding, untill finally a error comes that the maximum iteration has been obtained. Now the final model has loglikelihood of around -900.
But initially till 3500 iteration, it had a loglikelihood of -200. So finally should i take loglikelihood as -200 or -900?

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

采纳的回答

Sietse Braakman
Sietse Braakman 2019-6-16
Dear Praveen,
Generally, parameter estimation is a complicated topic that inherently has a lot of ambiguity. This makes it difficult to answer your question. However, I have posted guidance and steps to help you with parameter estimation in our discussion section:
https://www.mathworks.com/matlabcentral/discussions/simbiology/134-guidance-and-troubleshooting-for-parameter-estimation-in-simbiology
I cannot answer your question, other than that aim an optimization is to minimize the value of the objective function. SimBiology's progress plots reports the log-likelihood (log of your objective function) and aims to maximize the log-likelihood (actually, it minimizes the value of the negative log-likelihood. So you are right in saying that -200 would be the function value you are looking for. What might happen during your optimization is that your algorithm gets stuck in a local minimum of the objective function - NLME algorithms are local optimization algorithms and don't guarantee termination in a global minimum.
What I recommend to solve this problem is (if you haven't already done so):
  • that you follow the steps I describe in the link above to make sure your optimization problem is correctly set up
  • start with a regular non-linear regression estimation (no NLME/random effects) and use a global optimization algorithm to find the global minimum
  • Once you have found that minimum (and the associated parameter estimates), you can use these estimates as inital estimates for your NLME algorithm
Hope that helps,
Sietse
  3 个评论
显示 1更早的评论
Praveen Kumar M
Praveen Kumar M 2019-6-17
编辑：Praveen Kumar M 2019-6-17
Dear Sietse Braakman,
I have tried exactly like how you said and it worked well.
That is first I obtained the inital values using the non mixed effects model with lsqnonlin and exponential error. (I tried all other errors also but this was highest log likelihood)
Then I imputed the values into the non linear mixed effect model with exponential error and FOCE method of estimation of likelihood. I got likelihood of -205.0
Perfectly well. But I have got a warning/ error.
Iteration limit exceeded in Laplacian algorithm. Returing result from final iteration.
Please help me overcome that.
Praveen Kumar M
Praveen Kumar M 2019-6-18
On following the same terms, I have got one more query. When i start fitting the data in the model, using non linear mixed effect model with exponential error and FOCE method of estimation of likelihood, the progress bar which shows loglikelihood in realtime is showing somewhere around 820 and it stands long time in that position.
But finally as the maximum of the iteration is achieved and the result of final iteration is shown, the loglikelihood comes around 200. I know that the final loglikelihood only matters. But just want to make sure and understand, why the loglikelihood in progress bar is not reflecting the same

请先登录，再进行评论。

更多回答（0 个）

社区

更多回答在  SimBiology Community

类别

在 Help Center 和 File Exchange 中查找有关 Import Data 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by