forked from vincentarelbundock/Rdatasets
-
Notifications
You must be signed in to change notification settings - Fork 0
/
motors.html
101 lines (77 loc) · 2.72 KB
/
motors.html
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
<!DOCTYPE html><html><head><title>R: Accelerated Life Testing of Motorettes</title>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
<meta name="viewport" content="width=device-width, initial-scale=1.0, user-scalable=yes" />
<link rel="stylesheet" href="https://cdn.jsdelivr.net/npm/[email protected]/dist/katex.min.css">
<script type="text/javascript">
const macros = { "\\R": "\\textsf{R}", "\\code": "\\texttt"};
function processMathHTML() {
var l = document.getElementsByClassName('reqn');
for (let e of l) { katex.render(e.textContent, e, { throwOnError: false, macros }); }
return;
}</script>
<script defer src="https://cdn.jsdelivr.net/npm/[email protected]/dist/katex.min.js"
onload="processMathHTML();"></script>
<link rel="stylesheet" type="text/css" href="R.css" />
</head><body><div class="container">
<table style="width: 100%;"><tr><td>motors</td><td style="text-align: right;">R Documentation</td></tr></table>
<h2>
Accelerated Life Testing of Motorettes
</h2>
<h3>Description</h3>
<p>The <code>motors</code> data frame has 40 rows and 3 columns. It describes an
accelerated life test at each of four temperatures of 10 motorettes,
and has rather discrete times.
</p>
<h3>Usage</h3>
<pre><code class='language-R'>motors
</code></pre>
<h3>Format</h3>
<p>This data frame contains the following columns:
</p>
<dl>
<dt><code>temp</code></dt><dd>
<p>the temperature (degrees C) of the test.
</p>
</dd>
<dt><code>time</code></dt><dd>
<p>the time in hours to failure or censoring at 8064 hours (= 336 days).
</p>
</dd>
<dt><code>cens</code></dt><dd>
<p>an indicator variable for death.
</p>
</dd>
</dl>
<h3>Source</h3>
<p>Kalbfleisch, J. D. and Prentice, R. L. (1980)
<em>The Statistical Analysis of Failure Time Data.</em>
New York: Wiley.
</p>
<p>taken from
</p>
<p>Nelson, W. D. and Hahn, G. J. (1972)
Linear regression of a regression relationship from censored data.
Part 1 – simple methods and their application.
<em>Technometrics</em>, <b>14</b>, 247–276.
</p>
<h3>References</h3>
<p>Venables, W. N. and Ripley, B. D. (2002)
<em>Modern Applied Statistics with S.</em> Fourth edition. Springer.
</p>
<h3>Examples</h3>
<pre><code class='language-R'>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)) )
</code></pre>
</div>
</body></html>