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.