getwd()
## [1] "C:/Users/Claro/Desktop"
setwd("C:/Users/Claro/Desktop")
accidentes <- read.csv("Acc17.csv", sep=";", header=TRUE)
accidentes
## GADM Attsv2017
## 1 Quito 8828
## 2 Guayaquil 5048
## 3 Cuenca 1236
## 4 Ambato 923
## 5 Ibarra 781
## 6 Riobamba 621
## 7 Manta 603
## 8 Milagro 569
## 9 Duran 565
## 10 Santo Domingo 548
## 11 Loja 464
## 12 Latacunga 380
## 13 Daule 336
## 14 Mejia 321
## 15 Otavalo 314
## 16 Babahoyo 312
## 17 Naranjal 279
## 18 Machala 270
## 19 Santa Elena 253
## 20 Samborondon 233
## 21 La Libertad 230
## 22 San Jacinto de Yaguachi 229
## 23 Portoviejo 216
## 24 Chone 176
## 25 Quevedo 164
## 26 Salinas 158
## 27 El Empalme 157
## 28 Playas 149
## 29 El Triunfo 140
## 30 Balzar 131
## 31 Antonio Ante 121
## 32 Buena Fe 109
## 33 Quininde 93
## 34 Morona 93
## 35 Naranjito 91
## 36 El Guabo 88
## 37 Ventanas 88
## 38 Jipijapa 88
## 39 Ruminahui 84
## 40 Cotacachi 79
## 41 Pasaje 75
## 42 Guaranda 73
## 43 Esmeraldas 72
## 44 Montufar 71
## 45 Pueblo viejo 70
## 46 Cayambe 70
## 47 San Pedro de Pelileo 68
## 48 Azogues 67
## 49 Santa Rosa 67
## 50 Piñas 64
## 51 Arenillas 62
## 52 Pedro Carbo 60
## 53 Salcedo 59
## 54 Gualaceo 58
## 55 Guano 58
## 56 Santa Isabel 56
## 57 Vinces 54
## 58 Pedro Moncayo 54
## 59 Santa Lucia 50
## 60 Alfredo Baquerizo Moreno (Jujan) 49
## 61 Balao 49
## 62 Cañar 48
## 63 Montecristi 48
## 64 Paute 47
## 65 Zamora 47
## 66 Nobol 45
## 67 Urbina Jado (Salitre) 43
## 68 Lago Agrio 43
## 69 San Lorenzo 42
## 70 General Antonio Elizalde (Bucay) 42
## 71 San Miguel de Bolívar 40
## 72 Tulcan 40
## 73 Giron 38
## 74 Palestina 37
## 75 Tena 37
## 76 Sucua 36
## 77 Colta 35
## 78 Montalvo 35
## 79 Santiago de Pillaro 35
## 80 Pajan 34
## 81 Pastaza 34
## 82 La Concordia 34
## 83 Simón Bolívar 33
## 84 Quero 31
## 85 La Troncal 30
## 86 Camilo Ponce Enríquez 28
## 87 Pujili 28
## 88 Baba 28
## 89 El Carmen 28
## 90 Baños de Agua Santa 28
## 91 Espejo 27
## 92 Valencia 27
## 93 Orellana 26
## 94 Colimes 25
## 95 Quijos 25
## 96 Coronel Marcelino Maridueña 23
## 97 Lomas de Sargentillo 23
## 98 Catamayo 23
## 99 Bolívar_C 20
## 100 Huaquillas 19
## 101 Atacames 18
## 102 Rocafuerte 18
## 103 Archidona 18
## 104 Nabon 17
## 105 Mira 17
## 106 Pallatanga 17
## 107 San Miguel de Urcuqui 17
## 108 Biblian 16
## 109 La Mana 16
## 110 Balsas 16
## 111 Isidro Ayora 16
## 112 Puerto Lopez 16
## 113 Shushufindi 16
## 114 San José de Chimbo 15
## 115 Chambo 15
## 116 Rioverde 15
## 117 Mocha 15
## 118 Yantzaza 14
## 119 Echeandia 13
## 120 Cumanda 13
## 121 Mera 13
## 122 San Pedro de Huaca 12
## 123 Saquisili 12
## 124 Guamote 12
## 125 Saraguro 12
## 126 San Vicente 12
## 127 Caluma 11
## 128 Zaruma 11
## 129 Flavio Alfaro 11
## 130 Santa Ana 11
## 131 Sucre 11
## 132 Huamboya 11
## 133 Tisaleo 11
## 134 Chillanes 10
## 135 Eloy Alfaro 10
## 136 Pimampiro 10
## 137 Calvas 10
## 138 Tosagua 10
## 139 El Tambo 9
## 140 Portovelo 9
## 141 Gualaquiza 9
## 142 Santiago 9
## 143 Penipe 8
## 144 Macara 8
## 145 Paltas 8
## 146 Quinsaloma 8
## 147 El Chaco 8
## 148 Centinela del Condor 8
## 149 Mocache 7
## 150 Junin 7
## 151 Carlos Julio Arosemena Tola 7
## 152 Las Naves 6
## 153 Sigsig 5
## 154 Suscal 5
## 155 Alausi 5
## 156 Chunchi 5
## 157 Bolívar 5
## 158 Jaramijo 5
## 159 Santa Clara 5
## 160 Puerto Quito 5
## 161 Cevallos 5
## 162 Putumayo 5
## 163 Sucumbíos 5
## 164 Chordeleg 4
## 165 Sevilla de Oro 4
## 166 Sigchos 4
## 167 Espindola 4
## 168 Limon Indanza 4
## 169 Palora 4
## 170 Patate 4
## 171 Chinchipe 4
## 172 El Pangui 4
## 173 Cascales 4
## 174 La Joya de los Sachas 4
## 175 Loreto 4
## 176 Las Lajas 3
## 177 24 de Mayo 3
## 178 San Juan Bosco 3
## 179 Tiwintza 3
## 180 Paquisha 3
## 181 Gonzalo Pizarro 3
## 182 Pucara 2
## 183 Celica 2
## 184 Chaguarpamba 2
## 185 Pindal 2
## 186 Pichincha 2
## 187 Logrono 2
## 188 Arajuno 2
## 189 Nangaritza 2
## 190 Santa Cruz 2
## 191 San Fernando 1
## 192 Ona 1
## 193 Deleg 1
## 194 Pangua 1
## 195 Atahualpa 1
## 196 Marcabeli 1
## 197 Puyango 1
## 198 Zapotillo 1
## 199 Palenque 1
## 200 Pedernales 1
## 201 Taisha 1
## 202 Pablo VI 1
## 203 San Miguel de los Bancos 1
## 204 El Pan 0
## 205 Guachapala 0
## 206 Chilla 0
## 207 Muisne 0
## 208 Gonzanama 0
## 209 Sozoranga 0
## 210 Quilanga 0
## 211 Olmedo_L 0
## 212 Urdaneta 0
## 213 Olmedo_M 0
## 214 Jama 0
## 215 Pedro Vicente Maldonado 0
## 216 Yacuambi 0
## 217 Palanda 0
## 218 San Cristóbal 0
## 219 Isabela 0
## 220 Cuyabeno 0
## 221 Aguarico 0
## 222 TOTAL NA
attach(accidentes); Attsv2017
## [1] 8828 5048 1236 923 781 621 603 569 565 548 464 380 336 321
## [15] 314 312 279 270 253 233 230 229 216 176 164 158 157 149
## [29] 140 131 121 109 93 93 91 88 88 88 84 79 75 73
## [43] 72 71 70 70 68 67 67 64 62 60 59 58 58 56
## [57] 54 54 50 49 49 48 48 47 47 45 43 43 42 42
## [71] 40 40 38 37 37 36 35 35 35 34 34 34 33 31
## [85] 30 28 28 28 28 28 27 27 26 25 25 23 23 23
## [99] 20 19 18 18 18 17 17 17 17 16 16 16 16 16
## [113] 16 15 15 15 15 14 13 13 13 12 12 12 12 12
## [127] 11 11 11 11 11 11 11 10 10 10 10 10 9 9
## [141] 9 9 8 8 8 8 8 8 7 7 7 6 5 5
## [155] 5 5 5 5 5 5 5 5 5 4 4 4 4 4
## [169] 4 4 4 4 4 4 4 3 3 3 3 3 3 2
## [183] 2 2 2 2 2 2 2 2 1 1 1 1 1 1
## [197] 1 1 1 1 1 1 1 0 0 0 0 0 0 0
## [211] 0 0 0 0 0 0 0 0 0 0 0 NA
#detach(accidentes)
str(accidentes)
## 'data.frame': 222 obs. of 2 variables:
## $ GADM : Factor w/ 222 levels "24 de Mayo","Aguarico",..: 164 76 46 5 79 165 102 106 51 193 ...
## $ Attsv2017: int 8828 5048 1236 923 781 621 603 569 565 548 ...
View(accidentes)
boxplot(Attsv2017)
y <- boxplot(Attsv2017)
names(y)
## [1] "stats" "n" "conf" "out" "group" "names"
y$out
## [1] 8828 5048 1236 923 781 621 603 569 565 548 464 380 336 321
## [15] 314 312 279 270 253 233 230 229 216 176 164 158 157 149
## [29] 140
y$out %in% Attsv2017 #operacion logical
## [1] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [15] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [29] TRUE
which(Attsv2017 %in% y$out)# devuelve las posiciones de los valores atipicos
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23
## [24] 24 25 26 27 28 29
#OPCION 1
filtro <- Attsv2017 %in% y$out;filtro
## [1] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [12] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [23] TRUE TRUE TRUE TRUE TRUE TRUE TRUE FALSE FALSE FALSE FALSE
## [34] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
## [45] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
## [56] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
## [67] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
## [78] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
## [89] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
## [100] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
## [111] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
## [122] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
## [133] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
## [144] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
## [155] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
## [166] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
## [177] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
## [188] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
## [199] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
## [210] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
## [221] FALSE FALSE
Attsv2017[filtro]
## [1] 8828 5048 1236 923 781 621 603 569 565 548 464 380 336 321
## [15] 314 312 279 270 253 233 230 229 216 176 164 158 157 149
## [29] 140
summary(Attsv2017)
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 0.0 4.0 16.0 131.1 56.0 8828.0 1
#opcion 2
(ind <- which(filtro))
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23
## [24] 24 25 26 27 28 29
Attsv2017[ind]
## [1] 8828 5048 1236 923 781 621 603 569 565 548 464 380 336 321
## [15] 314 312 279 270 253 233 230 229 216 176 164 158 157 149
## [29] 140