确定对数正态分布的中值剩余寿命
Determine Median Residual Lifetime of Lognormal Distribution
Stefan Gelissen 在他的博客中提供了关于如何确定对数正态分布的中值剩余寿命的出色概述 here
然而,当运行他的代码时,我偶然发现了这一行
(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) 的下界不能小于 "zero" 并且上界应该是一个相当大的数,1e+9 就是这种情况。理论上,您甚至可以选择更大的数字。但是,没有确定上限的确切方法。
Stefan Gelissen 在他的博客中提供了关于如何确定对数正态分布的中值剩余寿命的出色概述 here
然而,当运行他的代码时,我偶然发现了这一行
(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) 的下界不能小于 "zero" 并且上界应该是一个相当大的数,1e+9 就是这种情况。理论上,您甚至可以选择更大的数字。但是,没有确定上限的确切方法。