14 de abril de 2018
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