This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. For more details on using R Markdown see http://rmarkdown.rstudio.com.
When you click the Knit button a document will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:
library(dplyr)
##
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
library(tidyr)
## Warning: package 'tidyr' was built under R version 4.5.2
library(ggplot2)
datasets::attenu
## event mag station dist accel
## 1 1 7.0 117 12.0 0.359
## 2 2 7.4 1083 148.0 0.014
## 3 2 7.4 1095 42.0 0.196
## 4 2 7.4 283 85.0 0.135
## 5 2 7.4 135 107.0 0.062
## 6 2 7.4 475 109.0 0.054
## 7 2 7.4 113 156.0 0.014
## 8 2 7.4 1008 224.0 0.018
## 9 2 7.4 1028 293.0 0.010
## 10 2 7.4 2001 359.0 0.004
## 11 2 7.4 117 370.0 0.004
## 12 3 5.3 1117 8.0 0.127
## 13 4 6.1 1438 16.1 0.411
## 14 4 6.1 1083 63.6 0.018
## 15 4 6.1 1013 6.6 0.509
## 16 4 6.1 1014 9.3 0.467
## 17 4 6.1 1015 13.0 0.279
## 18 4 6.1 1016 17.3 0.072
## 19 4 6.1 1095 105.0 0.012
## 20 4 6.1 1011 112.0 0.006
## 21 4 6.1 1028 123.0 0.003
## 22 5 6.6 270 105.0 0.018
## 23 5 6.6 280 122.0 0.048
## 24 5 6.6 116 141.0 0.011
## 25 5 6.6 266 200.0 0.007
## 26 5 6.6 117 45.0 0.142
## 27 5 6.6 113 130.0 0.031
## 28 5 6.6 112 147.0 0.006
## 29 5 6.6 130 187.0 0.010
## 30 5 6.6 475 197.0 0.010
## 31 5 6.6 269 203.0 0.006
## 32 5 6.6 135 211.0 0.013
## 33 6 5.6 1093 62.0 0.005
## 34 7 5.7 1093 62.0 0.003
## 35 8 5.3 111 19.0 0.086
## 36 8 5.3 116 21.0 0.179
## 37 8 5.3 290 13.0 0.205
## 38 8 5.3 112 22.0 0.073
## 39 8 5.3 113 29.0 0.045
## 40 9 6.6 128 17.0 0.374
## 41 9 6.6 126 19.6 0.200
## 42 9 6.6 127 20.2 0.147
## 43 9 6.6 141 21.1 0.188
## 44 9 6.6 266 21.9 0.204
## 45 9 6.6 110 24.2 0.335
## 46 9 6.6 1027 66.0 0.057
## 47 9 6.6 111 87.0 0.021
## 48 9 6.6 125 23.4 0.152
## 49 9 6.6 135 24.6 0.217
## 50 9 6.6 475 25.7 0.114
## 51 9 6.6 262 28.6 0.150
## 52 9 6.6 269 37.4 0.148
## 53 9 6.6 1052 46.7 0.112
## 54 9 6.6 411 56.9 0.043
## 55 9 6.6 290 60.7 0.057
## 56 9 6.6 130 61.4 0.030
## 57 9 6.6 272 62.0 0.027
## 58 9 6.6 1096 64.0 0.028
## 59 9 6.6 1102 82.0 0.034
## 60 9 6.6 112 88.0 0.030
## 61 9 6.6 113 91.0 0.039
## 62 10 5.3 1028 31.0 0.030
## 63 11 7.7 2714 45.0 0.110
## 64 11 7.7 2708 145.0 0.010
## 65 11 7.7 2715 300.0 0.010
## 66 12 6.2 3501 5.0 0.390
## 67 13 5.6 655 50.0 0.031
## 68 13 5.6 272 16.0 0.130
## 69 14 5.2 1032 17.0 0.011
## 70 14 5.2 1377 8.0 0.120
## 71 14 5.2 1028 10.0 0.170
## 72 14 5.2 1250 10.0 0.140
## 73 15 6.0 1051 8.0 0.110
## 74 15 6.0 1293 32.0 0.040
## 75 15 6.0 1291 30.0 0.070
## 76 15 6.0 1292 31.0 0.080
## 77 16 5.1 283 2.9 0.210
## 78 16 5.1 885 3.2 0.390
## 79 16 5.1 <NA> 7.6 0.280
## 80 17 7.6 2734 25.4 0.160
## 81 17 7.6 <NA> 32.9 0.064
## 82 17 7.6 2728 92.2 0.090
## 83 18 5.8 1413 1.2 0.420
## 84 18 5.8 1445 1.6 0.230
## 85 18 5.8 1408 9.1 0.130
## 86 18 5.8 1411 3.7 0.260
## 87 18 5.8 1410 5.3 0.270
## 88 18 5.8 1409 7.4 0.260
## 89 18 5.8 1377 17.9 0.110
## 90 18 5.8 1492 19.2 0.120
## 91 18 5.8 1251 23.4 0.038
## 92 18 5.8 1422 30.0 0.044
## 93 18 5.8 1376 38.9 0.046
## 94 19 6.5 <NA> 23.5 0.170
## 95 19 6.5 286 26.0 0.210
## 96 19 6.5 <NA> 0.5 0.320
## 97 19 6.5 5028 0.6 0.520
## 98 19 6.5 942 1.3 0.720
## 99 19 6.5 <NA> 1.4 0.320
## 100 19 6.5 5054 2.6 0.810
## 101 19 6.5 958 3.8 0.640
## 102 19 6.5 952 4.0 0.560
## 103 19 6.5 5165 5.1 0.510
## 104 19 6.5 117 6.2 0.400
## 105 19 6.5 955 6.8 0.610
## 106 19 6.5 5055 7.5 0.260
## 107 19 6.5 <NA> 7.6 0.240
## 108 19 6.5 <NA> 8.4 0.460
## 109 19 6.5 5060 8.5 0.220
## 110 19 6.5 412 8.5 0.230
## 111 19 6.5 5053 10.6 0.280
## 112 19 6.5 5058 12.6 0.380
## 113 19 6.5 5057 12.7 0.270
## 114 19 6.5 <NA> 12.9 0.310
## 115 19 6.5 5051 14.0 0.200
## 116 19 6.5 <NA> 15.0 0.110
## 117 19 6.5 5115 16.0 0.430
## 118 19 6.5 <NA> 17.7 0.270
## 119 19 6.5 931 18.0 0.150
## 120 19 6.5 5056 22.0 0.150
## 121 19 6.5 5059 22.0 0.150
## 122 19 6.5 5061 23.0 0.130
## 123 19 6.5 <NA> 23.2 0.190
## 124 19 6.5 5062 29.0 0.130
## 125 19 6.5 5052 32.0 0.066
## 126 19 6.5 <NA> 32.7 0.350
## 127 19 6.5 724 36.0 0.100
## 128 19 6.5 <NA> 43.5 0.160
## 129 19 6.5 5066 49.0 0.140
## 130 19 6.5 5050 60.0 0.049
## 131 19 6.5 2316 64.0 0.034
## 132 20 5.0 5055 7.5 0.264
## 133 20 5.0 942 8.8 0.263
## 134 20 5.0 5028 8.9 0.230
## 135 20 5.0 5165 9.4 0.147
## 136 20 5.0 952 9.7 0.286
## 137 20 5.0 958 9.7 0.157
## 138 20 5.0 955 10.5 0.237
## 139 20 5.0 117 10.5 0.133
## 140 20 5.0 412 12.0 0.055
## 141 20 5.0 5053 12.2 0.097
## 142 20 5.0 5054 12.8 0.129
## 143 20 5.0 5058 14.6 0.192
## 144 20 5.0 5057 14.9 0.147
## 145 20 5.0 5115 17.6 0.154
## 146 20 5.0 5056 23.9 0.060
## 147 20 5.0 5060 25.0 0.057
## 148 21 5.8 1030 10.8 0.120
## 149 21 5.8 1418 15.7 0.154
## 150 21 5.8 1383 16.7 0.052
## 151 21 5.8 1308 20.8 0.045
## 152 21 5.8 1298 28.5 0.086
## 153 21 5.8 1299 33.1 0.056
## 154 21 5.8 1219 40.3 0.065
## 155 22 5.5 <NA> 4.0 0.259
## 156 22 5.5 <NA> 10.1 0.267
## 157 22 5.5 1030 11.1 0.071
## 158 22 5.5 1418 17.7 0.275
## 159 22 5.5 1383 22.5 0.058
## 160 22 5.5 <NA> 26.5 0.026
## 161 22 5.5 1299 29.0 0.039
## 162 22 5.5 1308 30.9 0.112
## 163 22 5.5 1219 37.8 0.065
## 164 22 5.5 1456 48.3 0.026
## 165 23 5.3 5045 5.8 0.123
## 166 23 5.3 5044 12.0 0.133
## 167 23 5.3 5160 12.1 0.073
## 168 23 5.3 5043 20.5 0.097
## 169 23 5.3 5047 20.5 0.096
## 170 23 5.3 c168 25.3 0.230
## 171 23 5.3 5068 35.9 0.082
## 172 23 5.3 c118 36.1 0.110
## 173 23 5.3 5042 36.3 0.110
## 174 23 5.3 5067 38.5 0.094
## 175 23 5.3 5049 41.4 0.040
## 176 23 5.3 c204 43.6 0.050
## 177 23 5.3 5070 44.4 0.022
## 178 23 5.3 c266 46.1 0.070
## 179 23 5.3 c203 47.1 0.080
## 180 23 5.3 5069 47.7 0.033
## 181 23 5.3 5073 49.2 0.017
## 182 23 5.3 5072 53.1 0.022
dataset
attenu
## event mag station dist accel
## 1 1 7.0 117 12.0 0.359
## 2 2 7.4 1083 148.0 0.014
## 3 2 7.4 1095 42.0 0.196
## 4 2 7.4 283 85.0 0.135
## 5 2 7.4 135 107.0 0.062
## 6 2 7.4 475 109.0 0.054
## 7 2 7.4 113 156.0 0.014
## 8 2 7.4 1008 224.0 0.018
## 9 2 7.4 1028 293.0 0.010
## 10 2 7.4 2001 359.0 0.004
## 11 2 7.4 117 370.0 0.004
## 12 3 5.3 1117 8.0 0.127
## 13 4 6.1 1438 16.1 0.411
## 14 4 6.1 1083 63.6 0.018
## 15 4 6.1 1013 6.6 0.509
## 16 4 6.1 1014 9.3 0.467
## 17 4 6.1 1015 13.0 0.279
## 18 4 6.1 1016 17.3 0.072
## 19 4 6.1 1095 105.0 0.012
## 20 4 6.1 1011 112.0 0.006
## 21 4 6.1 1028 123.0 0.003
## 22 5 6.6 270 105.0 0.018
## 23 5 6.6 280 122.0 0.048
## 24 5 6.6 116 141.0 0.011
## 25 5 6.6 266 200.0 0.007
## 26 5 6.6 117 45.0 0.142
## 27 5 6.6 113 130.0 0.031
## 28 5 6.6 112 147.0 0.006
## 29 5 6.6 130 187.0 0.010
## 30 5 6.6 475 197.0 0.010
## 31 5 6.6 269 203.0 0.006
## 32 5 6.6 135 211.0 0.013
## 33 6 5.6 1093 62.0 0.005
## 34 7 5.7 1093 62.0 0.003
## 35 8 5.3 111 19.0 0.086
## 36 8 5.3 116 21.0 0.179
## 37 8 5.3 290 13.0 0.205
## 38 8 5.3 112 22.0 0.073
## 39 8 5.3 113 29.0 0.045
## 40 9 6.6 128 17.0 0.374
## 41 9 6.6 126 19.6 0.200
## 42 9 6.6 127 20.2 0.147
## 43 9 6.6 141 21.1 0.188
## 44 9 6.6 266 21.9 0.204
## 45 9 6.6 110 24.2 0.335
## 46 9 6.6 1027 66.0 0.057
## 47 9 6.6 111 87.0 0.021
## 48 9 6.6 125 23.4 0.152
## 49 9 6.6 135 24.6 0.217
## 50 9 6.6 475 25.7 0.114
## 51 9 6.6 262 28.6 0.150
## 52 9 6.6 269 37.4 0.148
## 53 9 6.6 1052 46.7 0.112
## 54 9 6.6 411 56.9 0.043
## 55 9 6.6 290 60.7 0.057
## 56 9 6.6 130 61.4 0.030
## 57 9 6.6 272 62.0 0.027
## 58 9 6.6 1096 64.0 0.028
## 59 9 6.6 1102 82.0 0.034
## 60 9 6.6 112 88.0 0.030
## 61 9 6.6 113 91.0 0.039
## 62 10 5.3 1028 31.0 0.030
## 63 11 7.7 2714 45.0 0.110
## 64 11 7.7 2708 145.0 0.010
## 65 11 7.7 2715 300.0 0.010
## 66 12 6.2 3501 5.0 0.390
## 67 13 5.6 655 50.0 0.031
## 68 13 5.6 272 16.0 0.130
## 69 14 5.2 1032 17.0 0.011
## 70 14 5.2 1377 8.0 0.120
## 71 14 5.2 1028 10.0 0.170
## 72 14 5.2 1250 10.0 0.140
## 73 15 6.0 1051 8.0 0.110
## 74 15 6.0 1293 32.0 0.040
## 75 15 6.0 1291 30.0 0.070
## 76 15 6.0 1292 31.0 0.080
## 77 16 5.1 283 2.9 0.210
## 78 16 5.1 885 3.2 0.390
## 79 16 5.1 <NA> 7.6 0.280
## 80 17 7.6 2734 25.4 0.160
## 81 17 7.6 <NA> 32.9 0.064
## 82 17 7.6 2728 92.2 0.090
## 83 18 5.8 1413 1.2 0.420
## 84 18 5.8 1445 1.6 0.230
## 85 18 5.8 1408 9.1 0.130
## 86 18 5.8 1411 3.7 0.260
## 87 18 5.8 1410 5.3 0.270
## 88 18 5.8 1409 7.4 0.260
## 89 18 5.8 1377 17.9 0.110
## 90 18 5.8 1492 19.2 0.120
## 91 18 5.8 1251 23.4 0.038
## 92 18 5.8 1422 30.0 0.044
## 93 18 5.8 1376 38.9 0.046
## 94 19 6.5 <NA> 23.5 0.170
## 95 19 6.5 286 26.0 0.210
## 96 19 6.5 <NA> 0.5 0.320
## 97 19 6.5 5028 0.6 0.520
## 98 19 6.5 942 1.3 0.720
## 99 19 6.5 <NA> 1.4 0.320
## 100 19 6.5 5054 2.6 0.810
## 101 19 6.5 958 3.8 0.640
## 102 19 6.5 952 4.0 0.560
## 103 19 6.5 5165 5.1 0.510
## 104 19 6.5 117 6.2 0.400
## 105 19 6.5 955 6.8 0.610
## 106 19 6.5 5055 7.5 0.260
## 107 19 6.5 <NA> 7.6 0.240
## 108 19 6.5 <NA> 8.4 0.460
## 109 19 6.5 5060 8.5 0.220
## 110 19 6.5 412 8.5 0.230
## 111 19 6.5 5053 10.6 0.280
## 112 19 6.5 5058 12.6 0.380
## 113 19 6.5 5057 12.7 0.270
## 114 19 6.5 <NA> 12.9 0.310
## 115 19 6.5 5051 14.0 0.200
## 116 19 6.5 <NA> 15.0 0.110
## 117 19 6.5 5115 16.0 0.430
## 118 19 6.5 <NA> 17.7 0.270
## 119 19 6.5 931 18.0 0.150
## 120 19 6.5 5056 22.0 0.150
## 121 19 6.5 5059 22.0 0.150
## 122 19 6.5 5061 23.0 0.130
## 123 19 6.5 <NA> 23.2 0.190
## 124 19 6.5 5062 29.0 0.130
## 125 19 6.5 5052 32.0 0.066
## 126 19 6.5 <NA> 32.7 0.350
## 127 19 6.5 724 36.0 0.100
## 128 19 6.5 <NA> 43.5 0.160
## 129 19 6.5 5066 49.0 0.140
## 130 19 6.5 5050 60.0 0.049
## 131 19 6.5 2316 64.0 0.034
## 132 20 5.0 5055 7.5 0.264
## 133 20 5.0 942 8.8 0.263
## 134 20 5.0 5028 8.9 0.230
## 135 20 5.0 5165 9.4 0.147
## 136 20 5.0 952 9.7 0.286
## 137 20 5.0 958 9.7 0.157
## 138 20 5.0 955 10.5 0.237
## 139 20 5.0 117 10.5 0.133
## 140 20 5.0 412 12.0 0.055
## 141 20 5.0 5053 12.2 0.097
## 142 20 5.0 5054 12.8 0.129
## 143 20 5.0 5058 14.6 0.192
## 144 20 5.0 5057 14.9 0.147
## 145 20 5.0 5115 17.6 0.154
## 146 20 5.0 5056 23.9 0.060
## 147 20 5.0 5060 25.0 0.057
## 148 21 5.8 1030 10.8 0.120
## 149 21 5.8 1418 15.7 0.154
## 150 21 5.8 1383 16.7 0.052
## 151 21 5.8 1308 20.8 0.045
## 152 21 5.8 1298 28.5 0.086
## 153 21 5.8 1299 33.1 0.056
## 154 21 5.8 1219 40.3 0.065
## 155 22 5.5 <NA> 4.0 0.259
## 156 22 5.5 <NA> 10.1 0.267
## 157 22 5.5 1030 11.1 0.071
## 158 22 5.5 1418 17.7 0.275
## 159 22 5.5 1383 22.5 0.058
## 160 22 5.5 <NA> 26.5 0.026
## 161 22 5.5 1299 29.0 0.039
## 162 22 5.5 1308 30.9 0.112
## 163 22 5.5 1219 37.8 0.065
## 164 22 5.5 1456 48.3 0.026
## 165 23 5.3 5045 5.8 0.123
## 166 23 5.3 5044 12.0 0.133
## 167 23 5.3 5160 12.1 0.073
## 168 23 5.3 5043 20.5 0.097
## 169 23 5.3 5047 20.5 0.096
## 170 23 5.3 c168 25.3 0.230
## 171 23 5.3 5068 35.9 0.082
## 172 23 5.3 c118 36.1 0.110
## 173 23 5.3 5042 36.3 0.110
## 174 23 5.3 5067 38.5 0.094
## 175 23 5.3 5049 41.4 0.040
## 176 23 5.3 c204 43.6 0.050
## 177 23 5.3 5070 44.4 0.022
## 178 23 5.3 c266 46.1 0.070
## 179 23 5.3 c203 47.1 0.080
## 180 23 5.3 5069 47.7 0.033
## 181 23 5.3 5073 49.2 0.017
## 182 23 5.3 5072 53.1 0.022
head(attenu)
## event mag station dist accel
## 1 1 7.0 117 12 0.359
## 2 2 7.4 1083 148 0.014
## 3 2 7.4 1095 42 0.196
## 4 2 7.4 283 85 0.135
## 5 2 7.4 135 107 0.062
## 6 2 7.4 475 109 0.054
summary(attenu)
## event mag station dist
## Min. : 1.00 Min. :5.000 117 : 5 Min. : 0.50
## 1st Qu.: 9.00 1st Qu.:5.300 1028 : 4 1st Qu.: 11.32
## Median :18.00 Median :6.100 113 : 4 Median : 23.40
## Mean :14.74 Mean :6.084 112 : 3 Mean : 45.60
## 3rd Qu.:20.00 3rd Qu.:6.600 135 : 3 3rd Qu.: 47.55
## Max. :23.00 Max. :7.700 (Other):147 Max. :370.00
## NA's : 16
## accel
## Min. :0.00300
## 1st Qu.:0.04425
## Median :0.11300
## Mean :0.15422
## 3rd Qu.:0.21925
## Max. :0.81000
##
mendeteksi missing value
is.na(attenu)
## event mag station dist accel
## [1,] FALSE FALSE FALSE FALSE FALSE
## [2,] FALSE FALSE FALSE FALSE FALSE
## [3,] FALSE FALSE FALSE FALSE FALSE
## [4,] FALSE FALSE FALSE FALSE FALSE
## [5,] FALSE FALSE FALSE FALSE FALSE
## [6,] FALSE FALSE FALSE FALSE FALSE
## [7,] FALSE FALSE FALSE FALSE FALSE
## [8,] FALSE FALSE FALSE FALSE FALSE
## [9,] FALSE FALSE FALSE FALSE FALSE
## [10,] FALSE FALSE FALSE FALSE FALSE
## [11,] FALSE FALSE FALSE FALSE FALSE
## [12,] FALSE FALSE FALSE FALSE FALSE
## [13,] FALSE FALSE FALSE FALSE FALSE
## [14,] FALSE FALSE FALSE FALSE FALSE
## [15,] FALSE FALSE FALSE FALSE FALSE
## [16,] FALSE FALSE FALSE FALSE FALSE
## [17,] FALSE FALSE FALSE FALSE FALSE
## [18,] FALSE FALSE FALSE FALSE FALSE
## [19,] FALSE FALSE FALSE FALSE FALSE
## [20,] FALSE FALSE FALSE FALSE FALSE
## [21,] FALSE FALSE FALSE FALSE FALSE
## [22,] FALSE FALSE FALSE FALSE FALSE
## [23,] FALSE FALSE FALSE FALSE FALSE
## [24,] FALSE FALSE FALSE FALSE FALSE
## [25,] FALSE FALSE FALSE FALSE FALSE
## [26,] FALSE FALSE FALSE FALSE FALSE
## [27,] FALSE FALSE FALSE FALSE FALSE
## [28,] FALSE FALSE FALSE FALSE FALSE
## [29,] FALSE FALSE FALSE FALSE FALSE
## [30,] FALSE FALSE FALSE FALSE FALSE
## [31,] FALSE FALSE FALSE FALSE FALSE
## [32,] FALSE FALSE FALSE FALSE FALSE
## [33,] FALSE FALSE FALSE FALSE FALSE
## [34,] FALSE FALSE FALSE FALSE FALSE
## [35,] FALSE FALSE FALSE FALSE FALSE
## [36,] FALSE FALSE FALSE FALSE FALSE
## [37,] FALSE FALSE FALSE FALSE FALSE
## [38,] FALSE FALSE FALSE FALSE FALSE
## [39,] FALSE FALSE FALSE FALSE FALSE
## [40,] FALSE FALSE FALSE FALSE FALSE
## [41,] FALSE FALSE FALSE FALSE FALSE
## [42,] FALSE FALSE FALSE FALSE FALSE
## [43,] FALSE FALSE FALSE FALSE FALSE
## [44,] FALSE FALSE FALSE FALSE FALSE
## [45,] FALSE FALSE FALSE FALSE FALSE
## [46,] FALSE FALSE FALSE FALSE FALSE
## [47,] FALSE FALSE FALSE FALSE FALSE
## [48,] FALSE FALSE FALSE FALSE FALSE
## [49,] FALSE FALSE FALSE FALSE FALSE
## [50,] FALSE FALSE FALSE FALSE FALSE
## [51,] FALSE FALSE FALSE FALSE FALSE
## [52,] FALSE FALSE FALSE FALSE FALSE
## [53,] FALSE FALSE FALSE FALSE FALSE
## [54,] FALSE FALSE FALSE FALSE FALSE
## [55,] FALSE FALSE FALSE FALSE FALSE
## [56,] FALSE FALSE FALSE FALSE FALSE
## [57,] FALSE FALSE FALSE FALSE FALSE
## [58,] FALSE FALSE FALSE FALSE FALSE
## [59,] FALSE FALSE FALSE FALSE FALSE
## [60,] FALSE FALSE FALSE FALSE FALSE
## [61,] FALSE FALSE FALSE FALSE FALSE
## [62,] FALSE FALSE FALSE FALSE FALSE
## [63,] FALSE FALSE FALSE FALSE FALSE
## [64,] FALSE FALSE FALSE FALSE FALSE
## [65,] FALSE FALSE FALSE FALSE FALSE
## [66,] FALSE FALSE FALSE FALSE FALSE
## [67,] FALSE FALSE FALSE FALSE FALSE
## [68,] FALSE FALSE FALSE FALSE FALSE
## [69,] FALSE FALSE FALSE FALSE FALSE
## [70,] FALSE FALSE FALSE FALSE FALSE
## [71,] FALSE FALSE FALSE FALSE FALSE
## [72,] FALSE FALSE FALSE FALSE FALSE
## [73,] FALSE FALSE FALSE FALSE FALSE
## [74,] FALSE FALSE FALSE FALSE FALSE
## [75,] FALSE FALSE FALSE FALSE FALSE
## [76,] FALSE FALSE FALSE FALSE FALSE
## [77,] FALSE FALSE FALSE FALSE FALSE
## [78,] FALSE FALSE FALSE FALSE FALSE
## [79,] FALSE FALSE TRUE FALSE FALSE
## [80,] FALSE FALSE FALSE FALSE FALSE
## [81,] FALSE FALSE TRUE FALSE FALSE
## [82,] FALSE FALSE FALSE FALSE FALSE
## [83,] FALSE FALSE FALSE FALSE FALSE
## [84,] FALSE FALSE FALSE FALSE FALSE
## [85,] FALSE FALSE FALSE FALSE FALSE
## [86,] FALSE FALSE FALSE FALSE FALSE
## [87,] FALSE FALSE FALSE FALSE FALSE
## [88,] FALSE FALSE FALSE FALSE FALSE
## [89,] FALSE FALSE FALSE FALSE FALSE
## [90,] FALSE FALSE FALSE FALSE FALSE
## [91,] FALSE FALSE FALSE FALSE FALSE
## [92,] FALSE FALSE FALSE FALSE FALSE
## [93,] FALSE FALSE FALSE FALSE FALSE
## [94,] FALSE FALSE TRUE FALSE FALSE
## [95,] FALSE FALSE FALSE FALSE FALSE
## [96,] FALSE FALSE TRUE FALSE FALSE
## [97,] FALSE FALSE FALSE FALSE FALSE
## [98,] FALSE FALSE FALSE FALSE FALSE
## [99,] FALSE FALSE TRUE FALSE FALSE
## [100,] FALSE FALSE FALSE FALSE FALSE
## [101,] FALSE FALSE FALSE FALSE FALSE
## [102,] FALSE FALSE FALSE FALSE FALSE
## [103,] FALSE FALSE FALSE FALSE FALSE
## [104,] FALSE FALSE FALSE FALSE FALSE
## [105,] FALSE FALSE FALSE FALSE FALSE
## [106,] FALSE FALSE FALSE FALSE FALSE
## [107,] FALSE FALSE TRUE FALSE FALSE
## [108,] FALSE FALSE TRUE FALSE FALSE
## [109,] FALSE FALSE FALSE FALSE FALSE
## [110,] FALSE FALSE FALSE FALSE FALSE
## [111,] FALSE FALSE FALSE FALSE FALSE
## [112,] FALSE FALSE FALSE FALSE FALSE
## [113,] FALSE FALSE FALSE FALSE FALSE
## [114,] FALSE FALSE TRUE FALSE FALSE
## [115,] FALSE FALSE FALSE FALSE FALSE
## [116,] FALSE FALSE TRUE FALSE FALSE
## [117,] FALSE FALSE FALSE FALSE FALSE
## [118,] FALSE FALSE TRUE FALSE FALSE
## [119,] FALSE FALSE FALSE FALSE FALSE
## [120,] FALSE FALSE FALSE FALSE FALSE
## [121,] FALSE FALSE FALSE FALSE FALSE
## [122,] FALSE FALSE FALSE FALSE FALSE
## [123,] FALSE FALSE TRUE FALSE FALSE
## [124,] FALSE FALSE FALSE FALSE FALSE
## [125,] FALSE FALSE FALSE FALSE FALSE
## [126,] FALSE FALSE TRUE FALSE FALSE
## [127,] FALSE FALSE FALSE FALSE FALSE
## [128,] FALSE FALSE TRUE FALSE FALSE
## [129,] FALSE FALSE FALSE FALSE FALSE
## [130,] FALSE FALSE FALSE FALSE FALSE
## [131,] FALSE FALSE FALSE FALSE FALSE
## [132,] FALSE FALSE FALSE FALSE FALSE
## [133,] FALSE FALSE FALSE FALSE FALSE
## [134,] FALSE FALSE FALSE FALSE FALSE
## [135,] FALSE FALSE FALSE FALSE FALSE
## [136,] FALSE FALSE FALSE FALSE FALSE
## [137,] FALSE FALSE FALSE FALSE FALSE
## [138,] FALSE FALSE FALSE FALSE FALSE
## [139,] FALSE FALSE FALSE FALSE FALSE
## [140,] FALSE FALSE FALSE FALSE FALSE
## [141,] FALSE FALSE FALSE FALSE FALSE
## [142,] FALSE FALSE FALSE FALSE FALSE
## [143,] FALSE FALSE FALSE FALSE FALSE
## [144,] FALSE FALSE FALSE FALSE FALSE
## [145,] FALSE FALSE FALSE FALSE FALSE
## [146,] FALSE FALSE FALSE FALSE FALSE
## [147,] FALSE FALSE FALSE FALSE FALSE
## [148,] FALSE FALSE FALSE FALSE FALSE
## [149,] FALSE FALSE FALSE FALSE FALSE
## [150,] FALSE FALSE FALSE FALSE FALSE
## [151,] FALSE FALSE FALSE FALSE FALSE
## [152,] FALSE FALSE FALSE FALSE FALSE
## [153,] FALSE FALSE FALSE FALSE FALSE
## [154,] FALSE FALSE FALSE FALSE FALSE
## [155,] FALSE FALSE TRUE FALSE FALSE
## [156,] FALSE FALSE TRUE FALSE FALSE
## [157,] FALSE FALSE FALSE FALSE FALSE
## [158,] FALSE FALSE FALSE FALSE FALSE
## [159,] FALSE FALSE FALSE FALSE FALSE
## [160,] FALSE FALSE TRUE FALSE FALSE
## [161,] FALSE FALSE FALSE FALSE FALSE
## [162,] FALSE FALSE FALSE FALSE FALSE
## [163,] FALSE FALSE FALSE FALSE FALSE
## [164,] FALSE FALSE FALSE FALSE FALSE
## [165,] FALSE FALSE FALSE FALSE FALSE
## [166,] FALSE FALSE FALSE FALSE FALSE
## [167,] FALSE FALSE FALSE FALSE FALSE
## [168,] FALSE FALSE FALSE FALSE FALSE
## [169,] FALSE FALSE FALSE FALSE FALSE
## [170,] FALSE FALSE FALSE FALSE FALSE
## [171,] FALSE FALSE FALSE FALSE FALSE
## [172,] FALSE FALSE FALSE FALSE FALSE
## [173,] FALSE FALSE FALSE FALSE FALSE
## [174,] FALSE FALSE FALSE FALSE FALSE
## [175,] FALSE FALSE FALSE FALSE FALSE
## [176,] FALSE FALSE FALSE FALSE FALSE
## [177,] FALSE FALSE FALSE FALSE FALSE
## [178,] FALSE FALSE FALSE FALSE FALSE
## [179,] FALSE FALSE FALSE FALSE FALSE
## [180,] FALSE FALSE FALSE FALSE FALSE
## [181,] FALSE FALSE FALSE FALSE FALSE
## [182,] FALSE FALSE FALSE FALSE FALSE
colSums(is.na(attenu))
## event mag station dist accel
## 0 0 16 0 0
library(VIM)
## Warning: package 'VIM' was built under R version 4.5.2
## Loading required package: colorspace
## Warning: package 'colorspace' was built under R version 4.5.2
## Loading required package: grid
## VIM is ready to use.
## Suggestions and bug-reports can be submitted at: https://github.com/statistikat/VIM/issues
##
## Attaching package: 'VIM'
## The following object is masked from 'package:datasets':
##
## sleep
aggr(attenu, numbers = TRUE, prop = FALSE )
imputasi missing value
attenu$event[is.na(attenu$event)] <- median(attenu$event, na.rm = TRUE)
attenu$mag[is.na(attenu$mag)] <- median(attenu$mag, na.rm = TRUE)
attenu
## event mag station dist accel
## 1 1 7.0 117 12.0 0.359
## 2 2 7.4 1083 148.0 0.014
## 3 2 7.4 1095 42.0 0.196
## 4 2 7.4 283 85.0 0.135
## 5 2 7.4 135 107.0 0.062
## 6 2 7.4 475 109.0 0.054
## 7 2 7.4 113 156.0 0.014
## 8 2 7.4 1008 224.0 0.018
## 9 2 7.4 1028 293.0 0.010
## 10 2 7.4 2001 359.0 0.004
## 11 2 7.4 117 370.0 0.004
## 12 3 5.3 1117 8.0 0.127
## 13 4 6.1 1438 16.1 0.411
## 14 4 6.1 1083 63.6 0.018
## 15 4 6.1 1013 6.6 0.509
## 16 4 6.1 1014 9.3 0.467
## 17 4 6.1 1015 13.0 0.279
## 18 4 6.1 1016 17.3 0.072
## 19 4 6.1 1095 105.0 0.012
## 20 4 6.1 1011 112.0 0.006
## 21 4 6.1 1028 123.0 0.003
## 22 5 6.6 270 105.0 0.018
## 23 5 6.6 280 122.0 0.048
## 24 5 6.6 116 141.0 0.011
## 25 5 6.6 266 200.0 0.007
## 26 5 6.6 117 45.0 0.142
## 27 5 6.6 113 130.0 0.031
## 28 5 6.6 112 147.0 0.006
## 29 5 6.6 130 187.0 0.010
## 30 5 6.6 475 197.0 0.010
## 31 5 6.6 269 203.0 0.006
## 32 5 6.6 135 211.0 0.013
## 33 6 5.6 1093 62.0 0.005
## 34 7 5.7 1093 62.0 0.003
## 35 8 5.3 111 19.0 0.086
## 36 8 5.3 116 21.0 0.179
## 37 8 5.3 290 13.0 0.205
## 38 8 5.3 112 22.0 0.073
## 39 8 5.3 113 29.0 0.045
## 40 9 6.6 128 17.0 0.374
## 41 9 6.6 126 19.6 0.200
## 42 9 6.6 127 20.2 0.147
## 43 9 6.6 141 21.1 0.188
## 44 9 6.6 266 21.9 0.204
## 45 9 6.6 110 24.2 0.335
## 46 9 6.6 1027 66.0 0.057
## 47 9 6.6 111 87.0 0.021
## 48 9 6.6 125 23.4 0.152
## 49 9 6.6 135 24.6 0.217
## 50 9 6.6 475 25.7 0.114
## 51 9 6.6 262 28.6 0.150
## 52 9 6.6 269 37.4 0.148
## 53 9 6.6 1052 46.7 0.112
## 54 9 6.6 411 56.9 0.043
## 55 9 6.6 290 60.7 0.057
## 56 9 6.6 130 61.4 0.030
## 57 9 6.6 272 62.0 0.027
## 58 9 6.6 1096 64.0 0.028
## 59 9 6.6 1102 82.0 0.034
## 60 9 6.6 112 88.0 0.030
## 61 9 6.6 113 91.0 0.039
## 62 10 5.3 1028 31.0 0.030
## 63 11 7.7 2714 45.0 0.110
## 64 11 7.7 2708 145.0 0.010
## 65 11 7.7 2715 300.0 0.010
## 66 12 6.2 3501 5.0 0.390
## 67 13 5.6 655 50.0 0.031
## 68 13 5.6 272 16.0 0.130
## 69 14 5.2 1032 17.0 0.011
## 70 14 5.2 1377 8.0 0.120
## 71 14 5.2 1028 10.0 0.170
## 72 14 5.2 1250 10.0 0.140
## 73 15 6.0 1051 8.0 0.110
## 74 15 6.0 1293 32.0 0.040
## 75 15 6.0 1291 30.0 0.070
## 76 15 6.0 1292 31.0 0.080
## 77 16 5.1 283 2.9 0.210
## 78 16 5.1 885 3.2 0.390
## 79 16 5.1 <NA> 7.6 0.280
## 80 17 7.6 2734 25.4 0.160
## 81 17 7.6 <NA> 32.9 0.064
## 82 17 7.6 2728 92.2 0.090
## 83 18 5.8 1413 1.2 0.420
## 84 18 5.8 1445 1.6 0.230
## 85 18 5.8 1408 9.1 0.130
## 86 18 5.8 1411 3.7 0.260
## 87 18 5.8 1410 5.3 0.270
## 88 18 5.8 1409 7.4 0.260
## 89 18 5.8 1377 17.9 0.110
## 90 18 5.8 1492 19.2 0.120
## 91 18 5.8 1251 23.4 0.038
## 92 18 5.8 1422 30.0 0.044
## 93 18 5.8 1376 38.9 0.046
## 94 19 6.5 <NA> 23.5 0.170
## 95 19 6.5 286 26.0 0.210
## 96 19 6.5 <NA> 0.5 0.320
## 97 19 6.5 5028 0.6 0.520
## 98 19 6.5 942 1.3 0.720
## 99 19 6.5 <NA> 1.4 0.320
## 100 19 6.5 5054 2.6 0.810
## 101 19 6.5 958 3.8 0.640
## 102 19 6.5 952 4.0 0.560
## 103 19 6.5 5165 5.1 0.510
## 104 19 6.5 117 6.2 0.400
## 105 19 6.5 955 6.8 0.610
## 106 19 6.5 5055 7.5 0.260
## 107 19 6.5 <NA> 7.6 0.240
## 108 19 6.5 <NA> 8.4 0.460
## 109 19 6.5 5060 8.5 0.220
## 110 19 6.5 412 8.5 0.230
## 111 19 6.5 5053 10.6 0.280
## 112 19 6.5 5058 12.6 0.380
## 113 19 6.5 5057 12.7 0.270
## 114 19 6.5 <NA> 12.9 0.310
## 115 19 6.5 5051 14.0 0.200
## 116 19 6.5 <NA> 15.0 0.110
## 117 19 6.5 5115 16.0 0.430
## 118 19 6.5 <NA> 17.7 0.270
## 119 19 6.5 931 18.0 0.150
## 120 19 6.5 5056 22.0 0.150
## 121 19 6.5 5059 22.0 0.150
## 122 19 6.5 5061 23.0 0.130
## 123 19 6.5 <NA> 23.2 0.190
## 124 19 6.5 5062 29.0 0.130
## 125 19 6.5 5052 32.0 0.066
## 126 19 6.5 <NA> 32.7 0.350
## 127 19 6.5 724 36.0 0.100
## 128 19 6.5 <NA> 43.5 0.160
## 129 19 6.5 5066 49.0 0.140
## 130 19 6.5 5050 60.0 0.049
## 131 19 6.5 2316 64.0 0.034
## 132 20 5.0 5055 7.5 0.264
## 133 20 5.0 942 8.8 0.263
## 134 20 5.0 5028 8.9 0.230
## 135 20 5.0 5165 9.4 0.147
## 136 20 5.0 952 9.7 0.286
## 137 20 5.0 958 9.7 0.157
## 138 20 5.0 955 10.5 0.237
## 139 20 5.0 117 10.5 0.133
## 140 20 5.0 412 12.0 0.055
## 141 20 5.0 5053 12.2 0.097
## 142 20 5.0 5054 12.8 0.129
## 143 20 5.0 5058 14.6 0.192
## 144 20 5.0 5057 14.9 0.147
## 145 20 5.0 5115 17.6 0.154
## 146 20 5.0 5056 23.9 0.060
## 147 20 5.0 5060 25.0 0.057
## 148 21 5.8 1030 10.8 0.120
## 149 21 5.8 1418 15.7 0.154
## 150 21 5.8 1383 16.7 0.052
## 151 21 5.8 1308 20.8 0.045
## 152 21 5.8 1298 28.5 0.086
## 153 21 5.8 1299 33.1 0.056
## 154 21 5.8 1219 40.3 0.065
## 155 22 5.5 <NA> 4.0 0.259
## 156 22 5.5 <NA> 10.1 0.267
## 157 22 5.5 1030 11.1 0.071
## 158 22 5.5 1418 17.7 0.275
## 159 22 5.5 1383 22.5 0.058
## 160 22 5.5 <NA> 26.5 0.026
## 161 22 5.5 1299 29.0 0.039
## 162 22 5.5 1308 30.9 0.112
## 163 22 5.5 1219 37.8 0.065
## 164 22 5.5 1456 48.3 0.026
## 165 23 5.3 5045 5.8 0.123
## 166 23 5.3 5044 12.0 0.133
## 167 23 5.3 5160 12.1 0.073
## 168 23 5.3 5043 20.5 0.097
## 169 23 5.3 5047 20.5 0.096
## 170 23 5.3 c168 25.3 0.230
## 171 23 5.3 5068 35.9 0.082
## 172 23 5.3 c118 36.1 0.110
## 173 23 5.3 5042 36.3 0.110
## 174 23 5.3 5067 38.5 0.094
## 175 23 5.3 5049 41.4 0.040
## 176 23 5.3 c204 43.6 0.050
## 177 23 5.3 5070 44.4 0.022
## 178 23 5.3 c266 46.1 0.070
## 179 23 5.3 c203 47.1 0.080
## 180 23 5.3 5069 47.7 0.033
## 181 23 5.3 5073 49.2 0.017
## 182 23 5.3 5072 53.1 0.022
colSums(is.na(attenu))
## event mag station dist accel
## 0 0 16 0 0
summary(attenu)
## event mag station dist
## Min. : 1.00 Min. :5.000 117 : 5 Min. : 0.50
## 1st Qu.: 9.00 1st Qu.:5.300 1028 : 4 1st Qu.: 11.32
## Median :18.00 Median :6.100 113 : 4 Median : 23.40
## Mean :14.74 Mean :6.084 112 : 3 Mean : 45.60
## 3rd Qu.:20.00 3rd Qu.:6.600 135 : 3 3rd Qu.: 47.55
## Max. :23.00 Max. :7.700 (Other):147 Max. :370.00
## NA's : 16
## accel
## Min. :0.00300
## 1st Qu.:0.04425
## Median :0.11300
## Mean :0.15422
## 3rd Qu.:0.21925
## Max. :0.81000
##
mendeteksi dan menangani outlier
# Menggunakan metode IQR untuk mendeteksi outlier pada kolom Ozone
Q1 <- quantile(attenu$event, 0.25)
Q3 <- quantile(attenu$event, 0.75)
IQR <- Q3 - Q1
# Batas bawah dan atas
lower_bound <- Q1 - 1.5 * IQR
upper_bound <- Q3 + 1.5 * IQR
# Menandai outlier dengan kondisi apakah nilainya di luar batas bawah atau atas
outliersa <- attenu$event < lower_bound
outliers <- attenu$event > upper_bound
sum(outliersa)
## [1] 0
sum(outliers)
## [1] 0
# Visualisasi boxplot untuk melihat outlier pada kolom Ozone
boxplot(attenu$event, main = "Boxplot event", col = "lightblue")
# Menangani outlier dengan winsorizing (mengganti nilai ekstrem dengan batas)
attenu$event[outliers] <- ifelse(attenu$event[outliers] < lower_bound, lower_bound, upper_bound)
mendeteksi dan menghapus duplikasi
# Cek jumlah duplikasi dalam dataset
sum(duplicated(attenu)) # Menghitung jumlah baris yang duplikat
## [1] 0
# Hapus duplikasi jika ada
attenu <- attenu[!duplicated(attenu), ] # Menyaring hanya baris unik
cek data
# Cek ulang data setelah preprocessing
summary(attenu) # Menampilkan ringkasan statistik setelah preprocessing
## event mag station dist
## Min. : 1.00 Min. :5.000 117 : 5 Min. : 0.50
## 1st Qu.: 9.00 1st Qu.:5.300 1028 : 4 1st Qu.: 11.32
## Median :18.00 Median :6.100 113 : 4 Median : 23.40
## Mean :14.74 Mean :6.084 112 : 3 Mean : 45.60
## 3rd Qu.:20.00 3rd Qu.:6.600 135 : 3 3rd Qu.: 47.55
## Max. :23.00 Max. :7.700 (Other):147 Max. :370.00
## NA's : 16
## accel
## Min. :0.00300
## 1st Qu.:0.04425
## Median :0.11300
## Mean :0.15422
## 3rd Qu.:0.21925
## Max. :0.81000
##
```