R中存活分析中MTTF的置信区间

时间:2014-02-13 13:56:49

标签: r confidence-interval survival-analysis weibull

对于不同的测试条件(电压),我必须在存在辐射的情况下评估电气设备的MTTF 对于每个电压,我只能测试几个单位(通常为5个),失败的单位数量范围从0(低电压)到可能5(高电压)。 因此,这是一个生存分析的案例,可能有非常高的“正确”审查(例如,在照射结束时,只有一个单位可能已经失败,其他4个单位可以正常)。我不是统计学家,所以我想征求意见。

我正在使用R和“survival”包来处理数据。根据我的理解,R计算Kaplan-Meier估计量并尝试将得到的经验分布拟合到Weibull故障模型(我使用的是survreg with dist =“weibull”,这是这种设备的故障模型)

我的问题:

我可以计算一个MTTF(我使用Weibull分布的平均值和Weibull拟合产生的参数),但是我有点迷失了MTTF的95%置信区间。
应用于幸存(...)返回的vcov()函数提供了Weibull参数的方差和协方差,但我认为由于Kaplan-Meier估计(参见:Greenwood公式)未集成,存在额外的不确定性对整件事情。

我不知道如何将所有碎片放在一起; R(在生存包或其他地方)是否有设施来实现这一目标?或者,您能否提供有关如何集成所有错误来源的指南?

0 个答案:

没有答案