-
Notifications
You must be signed in to change notification settings - Fork 0
/
Module 9
369 lines (361 loc) · 9.79 KB
/
Module 9
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
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
File chosen for Visualization Assignment:
Package: datasets
Item: faithful
Title Old Faithful Geyser Data
Rows: 272
Cols: 2
has_logical: FALSE
has_binary: FALSE
has_numberic: TRUE
has_character: FALSE
datasets faithful Old Faithful Geyser Data 272 2 FALSE FALSE TRUE FALSE
> install.packages("dataset")
> faithful
eruptions waiting
1 3.600 79
2 1.800 54
3 3.333 74
4 2.283 62
5 4.533 85
6 2.883 55
7 4.700 88
8 3.600 85
9 1.950 51
10 4.350 85
11 1.833 54
12 3.917 84
13 4.200 78
14 1.750 47
15 4.700 83
16 2.167 52
17 1.750 62
18 4.800 84
19 1.600 52
20 4.250 79
21 1.800 51
22 1.750 47
23 3.450 78
24 3.067 69
25 4.533 74
26 3.600 83
27 1.967 55
28 4.083 76
29 3.850 78
30 4.433 79
31 4.300 73
32 4.467 77
33 3.367 66
34 4.033 80
35 3.833 74
36 2.017 52
37 1.867 48
38 4.833 80
39 1.833 59
40 4.783 90
41 4.350 80
42 1.883 58
43 4.567 84
44 1.750 58
45 4.533 73
46 3.317 83
47 3.833 64
48 2.100 53
49 4.633 82
50 2.000 59
51 4.800 75
52 4.716 90
53 1.833 54
54 4.833 80
55 1.733 54
56 4.883 83
57 3.717 71
58 1.667 64
59 4.567 77
60 4.317 81
61 2.233 59
62 4.500 84
63 1.750 48
64 4.800 82
65 1.817 60
66 4.400 92
67 4.167 78
68 4.700 78
69 2.067 65
70 4.700 73
71 4.033 82
72 1.967 56
73 4.500 79
74 4.000 71
75 1.983 62
76 5.067 76
77 2.017 60
78 4.567 78
79 3.883 76
80 3.600 83
81 4.133 75
82 4.333 82
83 4.100 70
84 2.633 65
85 4.067 73
86 4.933 88
87 3.950 76
88 4.517 80
89 2.167 48
90 4.000 86
91 2.200 60
92 4.333 90
93 1.867 50
94 4.817 78
95 1.833 63
96 4.300 72
97 4.667 84
98 3.750 75
99 1.867 51
100 4.900 82
101 2.483 62
102 4.367 88
103 2.100 49
104 4.500 83
105 4.050 81
106 1.867 47
107 4.700 84
108 1.783 52
109 4.850 86
110 3.683 81
111 4.733 75
112 2.300 59
113 4.900 89
114 4.417 79
115 1.700 59
116 4.633 81
117 2.317 50
118 4.600 85
119 1.817 59
120 4.417 87
121 2.617 53
122 4.067 69
123 4.250 77
124 1.967 56
125 4.600 88
126 3.767 81
127 1.917 45
128 4.500 82
129 2.267 55
130 4.650 90
131 1.867 45
132 4.167 83
133 2.800 56
134 4.333 89
135 1.833 46
136 4.383 82
137 1.883 51
138 4.933 86
139 2.033 53
140 3.733 79
141 4.233 81
142 2.233 60
143 4.533 82
144 4.817 77
145 4.333 76
146 1.983 59
147 4.633 80
148 2.017 49
149 5.100 96
150 1.800 53
151 5.033 77
152 4.000 77
153 2.400 65
154 4.600 81
155 3.567 71
156 4.000 70
157 4.500 81
158 4.083 93
159 1.800 53
160 3.967 89
161 2.200 45
162 4.150 86
163 2.000 58
164 3.833 78
165 3.500 66
166 4.583 76
167 2.367 63
168 5.000 88
169 1.933 52
170 4.617 93
171 1.917 49
172 2.083 57
173 4.583 77
174 3.333 68
175 4.167 81
176 4.333 81
177 4.500 73
178 2.417 50
179 4.000 85
180 4.167 74
181 1.883 55
182 4.583 77
183 4.250 83
184 3.767 83
185 2.033 51
186 4.433 78
187 4.083 84
188 1.833 46
189 4.417 83
190 2.183 55
191 4.800 81
192 1.833 57
193 4.800 76
194 4.100 84
195 3.966 77
196 4.233 81
197 3.500 87
198 4.366 77
199 2.250 51
200 4.667 78
201 2.100 60
202 4.350 82
203 4.133 91
204 1.867 53
205 4.600 78
206 1.783 46
207 4.367 77
208 3.850 84
209 1.933 49
210 4.500 83
211 2.383 71
212 4.700 80
213 1.867 49
214 3.833 75
215 3.417 64
216 4.233 76
217 2.400 53
218 4.800 94
219 2.000 55
220 4.150 76
221 1.867 50
222 4.267 82
223 1.750 54
224 4.483 75
225 4.000 78
226 4.117 79
227 4.083 78
228 4.267 78
229 3.917 70
230 4.550 79
231 4.083 70
232 2.417 54
233 4.183 86
234 2.217 50
235 4.450 90
236 1.883 54
237 1.850 54
238 4.283 77
239 3.950 79
240 2.333 64
241 4.150 75
242 2.350 47
243 4.933 86
244 2.900 63
245 4.583 85
246 3.833 82
247 2.083 57
248 4.367 82
249 2.133 67
250 4.350 74
251 2.200 54
252 4.450 83
253 3.567 73
254 4.500 73
255 4.150 88
256 3.817 80
257 3.917 71
258 4.450 83
259 2.000 56
260 4.283 79
261 4.767 78
262 4.533 84
263 1.850 58
264 4.250 83
265 1.983 43
266 2.250 60
267 4.750 75
268 4.117 81
269 2.150 46
270 4.417 90
271 1.817 46
272 4.467 74
Given Example from Dataset Description, which is basic visualtions within R
> require(stats); require(graphics)
> f.tit <- "faithful data: Eruptions of Old Faithful"
> ne60 <- round(e60 <-60 * faithful$eruptions)
> all.equal(e60, ne60)
[1] "Mean relative difference: 9.515332e-05"
> all.equal(e60, ne60) #relative diff. ~1/10000
[1] "Mean relative difference: 9.515332e-05"
> f.tit
[1] "faithful data: Eruptions of Old Faithful"
> ne60
[1] 216 108 200 137 272 173 282 216 117 261 110 235 252 105 282 130
[17] 105 288 96 255 108 105 207 184 272 216 118 245 231 266 258 268
[33] 202 242 230 121 112 290 110 287 261 113 274 105 272 199 230 126
[49] 278 120 288 283 110 290 104 293 223 100 274 259 134 270 105 288
[65] 109 264 250 282 124 282 242 118 270 240 119 304 121 274 233 216
[81] 248 260 246 158 244 296 237 271 130 240 132 260 112 289 110 258
[97] 280 225 112 294 149 262 126 270 243 112 282 107 291 221 284 138
[113] 294 265 102 278 139 276 109 265 157 244 255 118 276 226 115 270
[129] 136 279 112 250 168 260 110 263 113 296 122 224 254 134 272 289
[145] 260 119 278 121 306 108 302 240 144 276 214 240 270 245 108 238
[161] 132 249 120 230 210 275 142 300 116 277 115 125 275 200 250 260
[177] 270 145 240 250 113 275 255 226 122 266 245 110 265 131 288 110
[193] 288 246 238 254 210 262 135 280 126 261 248 112 276 107 262 231
[209] 116 270 143 282 112 230 205 254 144 288 120 249 112 256 105 269
[225] 240 247 245 256 235 273 245 145 251 133 267 113 111 257 237 140
[241] 249 141 296 174 275 230 125 262 128 261 132 267 214 270 249 229
[257] 235 267 120 257 286 272 111 255 119 135 285 247 129 265 109 268
> table(zapsmall(abs(e60 -ne60))) # 0, 0.02, or 004
0 0.02 0.04
106 163 3
> te <- table(ne60)
> te
ne60
96 100 102 104 105 107 108 109 110 111 112 113 115 116 117 118 119
1 1 1 1 6 2 4 3 7 2 8 4 2 2 1 3 3
120 121 122 124 125 126 128 129 130 131 132 133 134 135 136 137 138
4 3 2 1 2 3 1 1 2 1 3 1 2 2 1 1 1
139 140 141 142 143 144 145 149 157 158 168 173 174 184 199 200 202
1 1 1 1 1 2 2 1 1 1 1 1 1 1 1 2 1
205 207 210 214 216 221 223 224 225 226 229 230 231 233 235 237 238
1 1 2 2 4 1 1 1 1 2 1 5 2 1 3 2 2
240 242 243 244 245 246 247 248 249 250 251 252 254 255 256 257 258
6 2 1 2 5 2 2 2 4 4 1 1 3 4 2 2 2
259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275
1 5 4 4 1 1 4 2 3 2 1 8 1 5 1 3 4
276 277 278 279 280 282 283 284 285 286 287 288 289 290 291 293 294
4 1 3 1 2 6 1 1 1 1 1 6 2 2 1 1 2
296 300 302 304 306
3 1 1 1 1
> te[te>=4]
ne60
105 108 110 112 113 120 216 230 240 245 249 250 255 260 261 262 265
6 4 7 8 4 4 4 5 6 5 4 4 4 5 4 4 4
270 272 275 276 282 288
8 5 4 4 6 6
> te[te>=4] # (too) many multiples of 5!
ne60
105 108 110 112 113 120 216 230 240 245 249 250 255 260 261 262 265
6 4 7 8 4 4 4 5 6 5 4 4 4 5 4 4 4
270 272 275 276 282 288
8 5 4 4 6 6
> plot(names(te), te, type = "h", main=f.tit, xlab="Eruption time (sec)")
> plot(faithful[,-3], main=f.tit)
> plot(faithful[, -3], main = f.tit,
+ xlab = "Eruption time (min)",
+ ylab = "Waiting time to next eruption (min)")
> lines(lowess(faithful$eruptions, faithful$waiting, f = 2/3, iter = 3),
+ col = "red")
TRY USING LATTICE
> library(lattice)
> xyplot(eruptions ~ waiting, faithful)
TRYP USING GGPLOT2
> library(ggplot2)
> ggplot(faithful, aes(eruptions, waiting)) + geom_point()