日志最初可能是高FOCE但随着迭代的进行,大于3500的日志可能性降低,最后结尾错误,取得了最大迭代。

1视图(30天)
我想开发一个人口使用Simbiology药代动力学模型,
我的设置是
评估方法:估计与nlmefit混合效应模型
误差模型是成比例的
方法来近似非线性混合效应模型是FOCE可能性。
FOCE模型拟合的进步显示了loglikelihood约-200(这很好)直到3500年但之后突然loglikelihood迭代到-900左右。
然后继续进行,直到最后一个错误的最大迭代获得。现在的最终模式loglikelihood -900左右。
但一开始直到3500年迭代,loglikelihood -200。所以最后我应该以loglikelihood为-200年或-900年吗?

接受的答案

Sietse Braakman
Sietse Braakman 2019年6月16日
亲爱的Praveen,
一般来说,参数估计本身是一个复杂的话题,有很多歧义。这个很难回答你的问题。然而,我发布了指导和步骤来帮助你与参数估计在我们的讨论:
我不能回答你的问题,除此之外,一个优化目标是最小化目标函数的值。SimBiology的 情节进展 报告对数似(你的目标函数的日志),旨在最大化对数似(实际上,它最大限度地减少的价值负对数似。所以你就在说,-200年将是你正在寻找的函数值。会发生什么在你优化你的算法陷入局部最小值的目标函数——NLME算法局部优化算法和不保证终止在全球最低。
我建议来解决这个问题是什么(如果您还没有这样做):
  • 你跟我描述的步骤在上面的链接,以确保你的优化问题是正确设置
  • 从常规非线性回归估计(没有NLME /随机效应)和使用一个全局优化算法寻找全球最低
  • 一旦你发现最低(和相关的参数估计),您可以使用这些估计作为NLME初估计算法
希望有所帮助,
Sietse
3评论
Praveen Kumar米
Praveen Kumar米 2019年6月18日
在相同的条款后,我有一个查询。当我开始拟合模型中的数据,使用 非线性混合效应模型与指数的错误 FOCE的似然估计方法, 在实时显示进度条显示loglikelihood 820左右 长时间站在那个位置。
但最终的最大迭代实现和最后一次迭代的结果显示, loglikelihood是200左右。 我只知道最后loglikelihood至关重要。但只是想确定和理解, 为什么loglikelihood在进度条反映一样吗

登录置评。

更多的答案(0)

社区

更多的答案SimBiology社区

类别

找到更多的在非线性Mixed-Effects建模帮助中心文件交换

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!