R motors 电机加速寿命测试

说明

motors 数据帧有 40 行和 3 列。它说明了 10 个电机在四种温度下的加速寿命测试，并且具有相当离散的时间。

用法

motors

格式

该 DataFrame 包含以下列：

temp

测试温度(摄氏度)。

time

故障或审查的时间为 8064 小时(= 336 天)。

cens

死亡的指示变量。

例子

library(survival)
plot(survfit(Surv(time, cens) ~ factor(temp), motors), conf.int = FALSE)
# fit Weibull model
motor.wei <- survreg(Surv(time, cens) ~ temp, motors)
summary(motor.wei)
# and predict at 130C
unlist(predict(motor.wei, data.frame(temp=130), se.fit = TRUE))

motor.cox <- coxph(Surv(time, cens) ~ temp, motors)
summary(motor.cox)
# predict at temperature 200
plot(survfit(motor.cox, newdata = data.frame(temp=200),
     conf.type = "log-log"))
summary( survfit(motor.cox, newdata = data.frame(temp=130)) )

来源

Kalbfleisch, J. D. 和 Prentice, R. L. (1980) 故障时间数据的统计分析。纽约：威利。

取自

Nelson, W. D. 和 Hahn, G. J. (1972) 根据审查数据进行回归关系的线性回归。第 1 部分 - 简单方法及其应用。技术计量学，14, 247-276。

参考

Venables, W. N. and Ripley, B. D. (2002) Modern Applied Statistics with S. Fourth edition. Springer.