r - 确定对数正态分布的中值剩余寿命

Question

Stefan Gelissen 在他的博客中提供了关于如何确定对数正态分布的中值剩余寿命的出色概述

但是，在运行他的代码时，我偶然发现了这一行

(mdrlhat <- uniroot(mdrl, c(0, 1e+9), t=t, Xmu=Xmu, Xsigma=Xsigma)$root)

其中 Stefan 使用硬编码区间 c(0, 1e+9)，这适用于 Meeker 和 Escobar (1998) 的示例数据集。

但是，我不清楚在使用不同数据的另一种情况下应如何确定此间隔。

任何人都可以帮忙吗？谢谢你。

Answer 1

区间 c(0, 1e+9) 的下限不能低于“零”，上限应该是一个相当大的数字，1e+9 就是这种情况。从理论上讲，您甚至可以选择更大的数字。但是，没有确定上限的确切方法。