Stefan Gelissen在他的博客中提供了关于如何确定对数正态分布的中位剩余寿命的优秀概述。此处
然而,当运行他的代码时,我偶然发现了这行代码
(mdrlhat <- uniroot(mdrl, c(0, 1e+9), t=t, Xmu=Xmu, Xsigma=Xsigma)$root)
其中Stefan使用了一个硬编码区间c(0, 1e+9),这对Meeker和Escobar(1998)的例子数据集来说很好。
然而,我不清楚在另一个使用不同数据的情况下,这个区间应该如何确定。
有谁能帮忙吗?谢谢您的帮助。
区间c(0,1e+9)的下界不能低于 "零",上界应该是一个合理的大数,1e+9就是如此,理论上,你甚至可以选择一个更大的数字。但是,目前还没有确定上界的确切方法。