1 Dataframe Ames

#install.packages("modeldata")
 
dim(ames)
## [1] 2930   74

El conjunto de datos Ames proporcionado por el paquete modeldata contiene información sobre propiedades residenciales vendidas en Ames, Iowa, entre 2006 y 2010. Este conjunto de datos se utiliza comúnmente en la práctica de la ciencia de datos y el análisis predictivo para comprender las características que influyen en el precio de las viviendas y para construir modelos de predicción de precios de viviendas.

El conjunto de datos contiene información sobre 2930 propiedades en Ames, Iowa, incluidas columnas relacionadas con:

  • Características de la casa (dormitorios, garaje, chimenea, piscina, porche, etc.)
  • ubicación (barrio) y latitud-longitud.
  • información del lote (zonificación, forma, tamaño, etc.)
  • calificaciones de condición y calidad
  • precio de venta

Aquí hay una explicación de los atributos más importantes del conjunto de datos “Ames”:

  • Order: Un identificador único para cada fila del conjunto de datos.

  • PID: Un identificador único para cada propiedad.

  • MS_SubClass: La clase de construcción de la propiedad.

  • MS_Zoning: La clasificación de zonificación de la propiedad.

  • Lot_Frontage: La longitud de la línea de frente de la propiedad (en pies).

  • Lot_Area: El tamaño del lote de la propiedad (en pies cuadrados).

  • Street: El tipo de acceso a la propiedad (pavimentado o no pavimentado).

  • Alley: Tipo de acceso a través del callejón (pavimentado, no pavimentado o sin callejón).

  • Lot_Shape: Forma general del lote de la propiedad.

  • Land_Contour: Contorno plano de la propiedad.

  • Utilities: Disponibilidad de servicios públicos (agua, electricidad, gas, etc.).

  • Lot_Configuration: Configuración del lote en relación con la vecindad.

  • Land_Slope: Pendiente de la propiedad.

  • Neighborhood: Vecindario en el que se encuentra la propiedad.

  • Condition_1 y Condition_2: Condiciones de proximidad a diversas características (carreteras, ferrocarriles, etc.).

  • Bldg_Type: Tipo de vivienda (unifamiliar, dúplex, etc.).

  • House_Style: Estilo de la casa (un piso, dos pisos, etc.).

  • Overall_Qual: Calidad general de la casa.

  • Overall_Cond: Condición general de la casa.

  • Year_Built: Año de construcción de la casa.

  • Year_Remod_Add: Año de remodelación de la casa.

  • Roof_Style y Roof_Matl: Estilo y material del techo.

  • Exterior_1st y Exterior_2nd: Revestimiento exterior de la casa.

  • Mas_Vnr_Type: Tipo de revestimiento de mampostería.

  • Mas_Vnr_Area: Área de revestimiento de mampostería en pies cuadrados.

  • Exter_Qual y Exter_Cond: Calidad y condición del material exterior.

  • Foundation: Tipo de cimentación de la casa.

  • Bsmt_Qual y Bsmt_Cond: Calidad y condición del sótano.

  • Bsmt_Exposure: Exposición del sótano.

  • BsmtFin_Type_1 y BsmtFin_Type_2: Calidad del acabado del sótano.

  • BsmtFin_SF_1 y BsmtFin_SF_2: Área acabada del sótano (en pies cuadrados).

  • Bsmt_Unf_SF: Área sin terminar del sótano (en pies cuadrados).

  • Total_Bsmt_SF: Área total del sótano (en pies cuadrados).

  • Heating y Heating_QC: Tipo y calidad del sistema de calefacción.

  • Central_Air: Si la casa tiene aire acondicionado central.

  • Electrical: Tipo de sistema eléctrico.

  • 1st_Flr_SF y 2nd_Flr_SF: Área del primer y segundo piso (en pies cuadrados).

  • Low_Qual_Fin_SF: Área de acabado de baja calidad (en pies cuadrados).

  • Gr_Liv_Area: Área habitable sobre el nivel del suelo (en pies cuadrados).

  • Bsmt_Full_Bath y Bsmt_Half_Bath: Cantidad de baños completos y medios en el sótano.

  • Full_Bath y Half_Bath: Cantidad de baños completos y medios sobre el nivel del suelo.

  • Bedroom_AbvGr: Cantidad de dormitorios sobre el nivel del suelo.

  • Kitchen_AbvGr: Cantidad de cocinas sobre el nivel del suelo.

  • Kitchen_Qual: Calidad de la cocina.

  • TotRms_AbvGrd: Cantidad total de habitaciones sobre el nivel del suelo (excluyendo baños).

  • Functional: Funcionalidad de la casa.

  • Fireplaces: Cantidad de chimeneas.

  • Fireplace_Qu: Calidad de la chimenea.

  • Garage_Type: Ubicación del garaje.

  • Garage_Yr_Blt: Año de construcción del garaje.

  • Garage_Finish: Acabado interior del garaje.

  • Garage_Cars: Capacidad del garaje en términos de número de autos.

  • Garage_Area: Área del garaje (en pies cuadrados).

  • Garage_Qual y Garage_Cond: Calidad y condición del garaje.

  • Paved_Drive: Tipo de entrada pavimentada.

  • Wood_Deck_SF, Open_Porch_SF, Enclosed_Porch, 3Ssn_Porch, y Screen_Porch: Área de porches y terrazas (en pies cuadrados).

  • Pool_Area: Área de la piscina (en pies cuadrados).

  • Pool_QC: Calidad de la piscina.

  • Fence: Calidad de la valla.

  • Misc_Feature: Otras características no cubiertas por los campos anteriores.

  • Misc_Val: Valor de las características adicionales.

  • Mo_Sold y Yr_Sold: Mes y año de venta de la casa.

  • Sale_Type: Tipo de venta.

  • Sale_Condition: Condición de venta.

  • Sale_Price: Precio de venta de la propiedad (variable objetivo).

Este conjunto de datos contiene una amplia gama de características que describen diferentes aspectos de las propiedades residenciales en Ames. Es comúnmente utilizado para la construcción y evaluación de modelos de regresión para predecir el precio de venta de las viviendas.

2 Grafico

2.1 Ploteo básico

plot(ames$Longitude, ames$Latitude, xlab = "Longitud", 
     ylab = "Latitud", 
     main = "Ubicación de propiedades en Ames", 
     col = as.integer(ames$Neighborhood))

2.2 leaflet de un barrio

levels(ames$Neighborhood)
##  [1] "North_Ames"                             
##  [2] "College_Creek"                          
##  [3] "Old_Town"                               
##  [4] "Edwards"                                
##  [5] "Somerset"                               
##  [6] "Northridge_Heights"                     
##  [7] "Gilbert"                                
##  [8] "Sawyer"                                 
##  [9] "Northwest_Ames"                         
## [10] "Sawyer_West"                            
## [11] "Mitchell"                               
## [12] "Brookside"                              
## [13] "Crawford"                               
## [14] "Iowa_DOT_and_Rail_Road"                 
## [15] "Timberland"                             
## [16] "Northridge"                             
## [17] "Stone_Brook"                            
## [18] "South_and_West_of_Iowa_State_University"
## [19] "Clear_Creek"                            
## [20] "Meadow_Village"                         
## [21] "Briardale"                              
## [22] "Bloomington_Heights"                    
## [23] "Veenker"                                
## [24] "Northpark_Villa"                        
## [25] "Blueste"                                
## [26] "Greens"                                 
## [27] "Green_Hills"                            
## [28] "Landmark"                               
## [29] "Hayden_Lake"
boxplot(ames$Lot_Frontage)

summary(ames$Lot_Frontage)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    0.00   43.00   63.00   57.65   78.00  313.00
ames7 <- ames %>%
  filter(Lot_Frontage >= 100)

viviendasPrecioAlto <- ames %>%
  filter(Sale_Price >= 213500)

viviendasPrecioAlto <- viviendasPrecioAlto %>%
  select(Neighborhood, Sale_Price)

tablaPorBarrio <- viviendasPrecioAlto%>%
  group_by(Neighborhood)%>%
  summarise( cantidad = n(),
             Minimo_Valor = min(Sale_Price,na.rm = TRUE),
             Maximo_Valor = max(Sale_Price,na.rm = TRUE),
             Valor_Promedio = mean(Sale_Price,na.rm = TRUE))

ames$Latitude
##    [1] 42.05403 42.05301 42.05266 42.05125 42.06090 42.06078 42.06298 42.06073
##    [9] 42.06112 42.05919 42.05848 42.05742 42.05815 42.05727 42.05917 42.06124
##   [17] 42.06035 42.05885 42.05667 42.05445 42.05035 42.05027 42.04930 42.05515
##   [25] 42.05459 42.05523 42.05340 42.05638 42.05532 42.05183 42.05168 42.05168
##   [33] 42.04967 42.05068 42.05068 42.04981 42.06235 42.06211 42.06113 42.06129
##   [41] 42.06087 42.06081 42.06088 42.06030 42.05962 42.05839 42.05860 42.05710
##   [49] 42.05716 42.05825 42.05718 42.06193 42.06257 42.06230 42.06144 42.06098
##   [57] 42.06119 42.06094 42.05710 42.05439 42.05324 42.05172 42.05167 42.05246
##   [65] 42.05244 42.05091 42.05306 42.05183 42.05165 42.05008 42.05235 42.05143
##   [73] 42.03625 42.03611 42.03603 42.03548 42.03666 42.03459 42.03470 42.03663
##   [81] 42.03455 42.03468 42.03608 42.03654 42.03479 42.03569 42.03486 42.03457
##   [89] 42.03533 42.03460 42.04825 42.04682 42.04577 42.04745 42.04809 42.04857
##   [97] 42.04790 42.04788 42.04778 42.04717 42.04615 42.04615 42.04642 42.04611
##  [105] 42.04703 42.04638 42.04257 42.04370 42.03456 42.04861 42.04632 42.04304
##  [113] 42.04331 42.04331 42.04333 42.04382 42.04382 42.04657 42.04221 42.04336
##  [121] 42.03909 42.03897 42.04033 42.04115 42.04039 42.04184 42.03921 42.03921
##  [129] 42.03854 42.03526 42.03548 42.03594 42.03610 42.03742 42.04849 42.04594
##  [137] 42.04951 42.04745 42.04489 42.04398 42.04291 42.04308 42.04485 42.04397
##  [145] 42.04405 42.04418 42.04488 42.04359 42.04231 42.04040 42.04137 42.03851
##  [153] 42.03984 42.03844 42.03852 42.03825 42.03813 42.03587 42.03721 42.04008
##  [161] 42.04112 42.03587 42.03594 42.03476 42.03485 42.03475 42.03485 42.03687
##  [169] 42.03709 42.03331 42.03213 42.03136 42.03138 42.03038 42.03328 42.03364
##  [177] 42.03239 42.03128 42.03031 42.03038 42.03047 42.03039 42.03033 42.02810
##  [185] 42.02915 42.02815 42.02808 42.03005 42.02816 42.02916 42.02817 42.03364
##  [193] 42.03019 42.03362 42.03356 42.03360 42.03366 42.03363 42.03156 42.03134
##  [201] 42.03230 42.02799 42.02889 42.02887 42.02691 42.02069 42.02414 42.02674
##  [209] 42.03345 42.03243 42.02826 42.02675 42.02616 42.02526 42.02282 42.03445
##  [217] 42.03339 42.03339 42.03301 42.03237 42.02987 42.03069 42.02992 42.02973
##  [225] 42.03052 42.03119 42.03136 42.03128 42.03008 42.03009 42.03387 42.03216
##  [233] 42.02655 42.02467 42.02308 42.02823 42.02653 42.02421 42.02517 42.02526
##  [241] 42.02475 42.02413 42.02502 42.02416 42.03374 42.03208 42.03378 42.03068
##  [249] 42.03112 42.02575 42.02461 42.02460 42.02460 42.03060 42.02985 42.02432
##  [257] 42.02132 42.01841 42.01867 42.02219 42.02204 42.02139 42.02113 42.01866
##  [265] 42.01907 42.01612 42.01872 42.01775 42.01628 42.01625 42.01404 42.01409
##  [273] 42.01402 42.02036 42.02125 42.02104 42.02063 42.02116 42.02031 42.01985
##  [281] 42.01963 42.02020 42.02179 42.02158 42.01773 42.02267 42.02270 42.02209
##  [289] 42.01774 42.01780 42.01639 42.01909 42.01909 42.01623 42.01494 42.01468
##  [297] 42.01100 42.01061 42.00935 42.00934 42.01080 42.02277 42.02141 42.02147
##  [305] 42.02146 42.02265 42.02133 42.02152 42.00869 41.99706 41.99705 41.99607
##  [313] 41.99478 41.99401 41.99774 41.99535 41.99528 41.99428 41.99353 41.99549
##  [321] 41.99435 41.99299 41.99321 41.99093 41.99325 41.99100 41.99186 41.99176
##  [329] 41.99188 41.99221 41.99212 41.99170 41.99182 41.99156 41.98873 41.98894
##  [337] 41.98766 41.98665 41.99130 41.98985 41.98918 42.05304 42.05139 42.04932
##  [345] 42.05996 42.06329 42.06090 42.06296 42.06060 42.06290 42.06198 42.06163
##  [353] 42.06118 42.05928 42.05913 42.05916 42.05827 42.05806 42.05965 42.05819
##  [361] 42.05793 42.05817 42.05785 42.05820 42.05875 42.05924 42.06027 42.05982
##  [369] 42.05844 42.05680 42.05353 42.05541 42.05508 42.05557 42.05184 42.05251
##  [377] 42.05210 42.05272 42.05288 42.05291 42.05035 42.04965 42.04939 42.04911
##  [385] 42.05162 42.05038 42.04941 42.05317 42.05235 42.05245 42.05267 42.04891
##  [393] 42.05640 42.05463 42.05523 42.05518 42.05515 42.05337 42.05331 42.05331
##  [401] 42.05644 42.05646 42.05276 42.05234 42.05168 42.05166 42.05182 42.05182
##  [409] 42.05062 42.05071 42.05071 42.05068 42.05053 42.05021 42.05009 42.05013
##  [417] 42.05022 42.05040 42.04889 42.04849 42.06231 42.06330 42.06303 42.06190
##  [425] 42.06128 42.06299 42.06299 42.06299 42.06294 42.06226 42.06226 42.06278
##  [433] 42.06279 42.06240 42.06211 42.06093 42.06035 42.06138 42.06144 42.06058
##  [441] 42.06058 42.06168 42.06043 42.06018 42.05991 42.05990 42.06140 42.05924
##  [449] 42.05878 42.05889 42.05920 42.05918 42.05918 42.05918 42.05834 42.05834
##  [457] 42.05718 42.05851 42.05716 42.05709 42.05702 42.05698 42.05728 42.05747
##  [465] 42.06232 42.06318 42.06206 42.06240 42.06256 42.06256 42.06230 42.06227
##  [473] 42.06208 42.06209 42.06156 42.05938 42.06174 42.06115 42.06111 42.06166
##  [481] 42.06137 42.06328 42.05837 42.05859 42.05948 42.05907 42.05960 42.05784
##  [489] 42.05783 42.05783 42.05710 42.05709 42.05696 42.05867 42.05375 42.05391
##  [497] 42.05563 42.05344 42.05304 42.05267 42.05041 42.04970 42.05161 42.05022
##  [505] 42.04992 42.05097 42.05428 42.05504 42.05430 42.05415 42.05415 42.05321
##  [513] 42.05326 42.05311 42.05180 42.05180 42.05179 42.05165 42.05165 42.05165
##  [521] 42.05041 42.05172 42.05179 42.05234 42.05235 42.05213 42.05142 42.05138
##  [529] 42.05212 42.05141 42.05125 42.05125 42.05089 42.03756 42.03761 42.03766
##  [537] 42.03600 42.03601 42.03528 42.03792 42.03813 42.03776 42.03743 42.03685
##  [545] 42.03607 42.03690 42.03658 42.03453 42.03552 42.03761 42.03597 42.03576
##  [553] 42.03654 42.03529 42.03581 42.03475 42.03462 42.03529 42.03466 42.03533
##  [561] 42.03460 42.03460 42.03458 42.03724 42.04767 42.04763 42.04842 42.04849
##  [569] 42.04734 42.04725 42.04710 42.04614 42.04642 42.04717 42.04378 42.04265
##  [577] 42.04355 42.04302 42.03451 42.04771 42.04801 42.04896 42.04876 42.04876
##  [585] 42.04765 42.04648 42.04648 42.04581 42.04580 42.04324 42.04384 42.04367
##  [593] 42.04469 42.04880 42.04671 42.04664 42.04482 42.04393 42.04315 42.04313
##  [601] 42.04224 42.04484 42.04477 42.04493 42.04488 42.04064 42.03849 42.04073
##  [609] 42.04001 42.03990 42.04089 42.04110 42.03946 42.03860 42.03549 42.03544
##  [617] 42.03464 42.03618 42.03474 42.03745 42.03595 42.03618 42.04931 42.04829
##  [625] 42.04683 42.04861 42.04872 42.04806 42.04579 42.04489 42.04397 42.04413
##  [633] 42.04399 42.04408 42.04547 42.04491 42.04297 42.04223 42.04324 42.04093
##  [641] 42.04086 42.04046 42.03967 42.03853 42.03836 42.03840 42.03712 42.03722
##  [649] 42.03492 42.03469 42.03475 42.03622 42.03583 42.03475 42.03807 42.03768
##  [657] 42.03714 42.03733 42.03601 42.03600 42.03467 42.03457 42.03476 42.03474
##  [665] 42.03526 42.03528 42.04179 42.04015 42.04030 42.04124 42.03851 42.03835
##  [673] 42.03905 42.04077 42.03881 42.03713 42.03571 42.03579 42.03584 42.03723
##  [681] 42.03602 42.03597 42.03572 42.03571 42.03721 42.03475 42.03486 42.03705
##  [689] 42.03895 42.03826 42.03529 42.03422 42.03327 42.03207 42.03233 42.03240
##  [697] 42.03220 42.03211 42.03191 42.03131 42.03047 42.03033 42.03035 42.03028
##  [705] 42.03106 42.03338 42.03358 42.03205 42.03230 42.03223 42.03052 42.03020
##  [713] 42.02967 42.02901 42.02904 42.02896 42.02931 42.02900 42.02811 42.02810
##  [721] 42.02770 42.02856 42.02820 42.02706 42.02720 42.02723 42.02322 42.02322
##  [729] 42.03400 42.03356 42.03338 42.03156 42.03154 42.03341 42.03362 42.03329
##  [737] 42.03335 42.03350 42.03347 42.03140 42.03252 42.03246 42.03236 42.03129
##  [745] 42.03238 42.03235 42.03109 42.03006 42.02676 42.02895 42.02785 42.02800
##  [753] 42.02704 42.02696 42.02531 42.02309 42.02580 42.02603 42.02425 42.02414
##  [761] 42.02283 42.03268 42.02810 42.02810 42.02805 42.02696 42.02282 42.02298
##  [769] 42.03338 42.03363 42.03237 42.03239 42.02986 42.03134 42.03063 42.03114
##  [777] 42.03343 42.03328 42.03230 42.03265 42.03256 42.03250 42.02655 42.02323
##  [785] 42.02403 42.02847 42.02766 42.02751 42.02448 42.02398 42.02392 42.02613
##  [793] 42.02600 42.02631 42.02499 42.02469 42.02472 42.02355 42.02275 42.03418
##  [801] 42.03371 42.03184 42.03254 42.03063 42.03062 42.03068 42.03252 42.03162
##  [809] 42.03183 42.03162 42.03155 42.03147 42.03146 42.03128 42.03057 42.03053
##  [817] 42.03109 42.03093 42.02663 42.02648 42.02641 42.02637 42.02571 42.02529
##  [825] 42.02558 42.02463 42.02460 42.02461 42.02477 42.02433 42.02440 42.02322
##  [833] 42.02977 42.02792 42.02555 42.02117 42.02118 42.02121 42.02082 42.01978
##  [841] 42.01900 42.01972 42.01860 42.01757 42.01690 42.01867 42.01686 42.01640
##  [849] 42.01571 42.01655 42.01613 42.01571 42.02138 42.02009 42.01918 42.01970
##  [857] 42.01955 42.01969 42.02067 42.02099 42.02179 42.02084 42.01711 42.01872
##  [865] 42.01885 42.01877 42.01870 42.01697 42.01698 42.01815 42.01888 42.01887
##  [873] 42.01801 42.01786 42.01786 42.01682 42.01680 42.01627 42.01627 42.01647
##  [881] 42.01396 42.01405 42.01351 42.01339 42.02012 42.02125 42.02122 42.01899
##  [889] 42.01865 42.01967 42.01893 42.01703 42.01896 42.02248 42.01942 42.02032
##  [897] 42.02180 42.02150 42.01900 42.01822 42.01938 42.01859 42.02209 42.02209
##  [905] 42.02209 42.02012 42.01971 42.01756 42.01692 42.01618 42.01909 42.01893
##  [913] 42.01789 42.01783 42.01893 42.01694 42.01684 42.01696 42.01607 42.01913
##  [921] 42.01918 42.01871 42.01758 42.01757 42.01668 42.01585 42.01615 42.01469
##  [929] 42.01317 42.01494 42.01375 42.01250 42.00949 42.00949 42.00949 42.00949
##  [937] 42.01064 42.01141 42.00913 42.02272 42.02143 42.02160 42.02018 42.02000
##  [945] 42.00868 42.01830 41.99689 41.99667 41.99712 41.99539 41.99707 41.99707
##  [953] 41.99692 41.99591 41.99607 42.00125 42.00138 41.99503 41.99506 41.99376
##  [961] 41.99473 41.99470 41.99768 41.99895 41.99817 41.99633 41.99609 41.99536
##  [969] 41.99296 41.99354 41.99219 41.99313 41.99212 41.99179 41.99292 41.99218
##  [977] 41.99282 41.99217 41.99254 41.99179 41.99137 41.98754 41.98878 41.98822
##  [985] 41.98766 41.98650 41.98651 41.99149 41.99095 42.05255 42.05145 42.04959
##  [993] 42.05996 42.06123 42.06077 42.06082 42.06037 42.06064 42.06230 42.06298
## [1001] 42.06137 42.06086 42.06076 42.05941 42.05936 42.05941 42.05844 42.05701
## [1009] 42.05952 42.05964 42.06139 42.06056 42.06024 42.05979 42.05890 42.05867
## [1017] 42.05602 42.05451 42.05456 42.05444 42.05536 42.05450 42.05445 42.05421
## [1025] 42.05543 42.05254 42.05033 42.05147 42.05068 42.04945 42.04986 42.04961
## [1033] 42.05107 42.05631 42.05539 42.05530 42.05433 42.05643 42.05469 42.05279
## [1041] 42.05275 42.05148 42.05233 42.05168 42.05169 42.05166 42.05000 42.05068
## [1049] 42.05026 42.04845 42.06205 42.06327 42.06313 42.06292 42.06226 42.06241
## [1057] 42.06243 42.06113 42.06027 42.06043 42.06077 42.06043 42.05990 42.05848
## [1065] 42.05868 42.05857 42.05865 42.05703 42.05869 42.05716 42.05796 42.05820
## [1073] 42.05677 42.05702 42.05745 42.05718 42.05712 42.05747 42.05709 42.05662
## [1081] 42.06334 42.06230 42.06229 42.06200 42.06330 42.06132 42.05964 42.06154
## [1089] 42.06098 42.06132 42.06131 42.05838 42.05799 42.05859 42.06001 42.05951
## [1097] 42.05695 42.05858 42.05388 42.05367 42.05575 42.05578 42.05223 42.05227
## [1105] 42.05038 42.05162 42.04950 42.04970 42.05069 42.04948 42.05424 42.05516
## [1113] 42.05510 42.05441 42.05437 42.05431 42.05432 42.05552 42.05423 42.05318
## [1121] 42.05165 42.05165 42.05031 42.05228 42.05246 42.05228 42.05205 42.05138
## [1129] 42.05217 42.05122 42.05090 42.03784 42.03527 42.03513 42.03792 42.03781
## [1137] 42.03805 42.03652 42.03726 42.03587 42.03588 42.03588 42.03499 42.03609
## [1145] 42.03453 42.03567 42.03684 42.03650 42.03647 42.03604 42.03518 42.03461
## [1153] 42.03569 42.03570 42.03459 42.03750 42.03453 42.04767 42.04861 42.04747
## [1161] 42.04809 42.04819 42.04781 42.04757 42.04755 42.04710 42.04710 42.04638
## [1169] 42.04638 42.04637 42.04609 42.04608 42.04696 42.04696 42.04665 42.04717
## [1177] 42.04717 42.04717 42.04501 42.04442 42.04464 42.04469 42.03673 42.03758
## [1185] 42.04873 42.04793 42.04925 42.04914 42.04876 42.04640 42.04294 42.04378
## [1193] 42.04585 42.04710 42.04658 42.04754 42.04818 42.04536 42.04426 42.04406
## [1201] 42.04485 42.04314 42.04384 42.04410 42.04240 42.04477 42.04477 42.04123
## [1209] 42.04066 42.03923 42.04011 42.04178 42.04178 42.04181 42.03832 42.04102
## [1217] 42.04187 42.04117 42.04030 42.03558 42.03557 42.03456 42.03574 42.03806
## [1225] 42.03618 42.03610 42.03713 42.03468 42.04895 42.04778 42.04709 42.04708
## [1233] 42.04693 42.04657 42.04237 42.04230 42.04474 42.04399 42.04479 42.04485
## [1241] 42.04430 42.04219 42.04382 42.04202 42.04230 42.04048 42.04030 42.04039
## [1249] 42.04044 42.03953 42.03852 42.03945 42.03720 42.03471 42.03730 42.03746
## [1257] 42.03589 42.03572 42.03594 42.03522 42.04002 42.04014 42.04117 42.03869
## [1265] 42.03999 42.03957 42.03807 42.03713 42.03710 42.03581 42.03716 42.03701
## [1273] 42.03729 42.03481 42.03456 42.03479 42.03469 42.03473 42.03464 42.03480
## [1281] 42.03697 42.03702 42.03572 42.03422 42.03413 42.03338 42.03261 42.03306
## [1289] 42.03239 42.03242 42.03136 42.03043 42.03041 42.03042 42.03048 42.03216
## [1297] 42.03227 42.03041 42.03036 42.03343 42.03357 42.03325 42.03201 42.03209
## [1305] 42.03238 42.03110 42.03048 42.02903 42.02906 42.02928 42.02925 42.02928
## [1313] 42.02829 42.02793 42.02707 42.02703 42.02908 42.02813 42.02883 42.02878
## [1321] 42.02840 42.02921 42.02740 42.02725 42.02774 42.02669 42.02727 42.03361
## [1329] 42.03431 42.03338 42.03047 42.03346 42.03357 42.03370 42.03337 42.03351
## [1337] 42.03332 42.03347 42.03337 42.03349 42.03135 42.03129 42.03149 42.03248
## [1345] 42.03146 42.03150 42.03147 42.03143 42.03105 42.03039 42.03106 42.02940
## [1353] 42.02843 42.02827 42.02707 42.02956 42.02903 42.02904 42.02863 42.02552
## [1361] 42.02411 42.02426 42.02499 42.02426 42.02257 42.02295 42.02648 42.02626
## [1369] 42.02488 42.02423 42.02372 42.03439 42.03441 42.03342 42.03128 42.02393
## [1377] 42.03360 42.03323 42.03236 42.03157 42.03125 42.03112 42.03135 42.03156
## [1385] 42.03055 42.03443 42.03443 42.03345 42.03383 42.03441 42.03249 42.03257
## [1393] 42.03216 42.03392 42.03369 42.03265 42.03356 42.02945 42.02466 42.02430
## [1401] 42.02521 42.02857 42.02860 42.02733 42.02751 42.02755 42.02838 42.02823
## [1409] 42.02819 42.02803 42.02403 42.02572 42.02504 42.02367 42.02446 42.02290
## [1417] 42.02277 42.02283 42.03177 42.03391 42.03085 42.03108 42.03062 42.02541
## [1425] 42.02528 42.02737 42.02719 42.02944 42.03028 42.02492 42.02578 42.02117
## [1433] 42.01915 42.01915 42.01850 42.01879 42.01742 42.01752 42.01791 42.01777
## [1441] 42.01628 42.01597 42.01626 42.01595 42.02214 42.02123 42.01884 42.01977
## [1449] 42.02069 42.02107 42.02128 42.02080 42.02099 42.01873 42.01841 42.01835
## [1457] 42.01871 42.01760 42.01801 42.01727 42.01748 42.01819 42.01698 42.01621
## [1465] 42.01634 42.01443 42.01671 42.01787 42.01628 42.01625 42.01785 42.01626
## [1473] 42.01455 42.01395 42.01391 42.01364 42.01362 42.01360 42.01597 42.01537
## [1481] 42.01383 42.01593 42.01593 42.01612 42.01612 42.01975 42.01994 42.02090
## [1489] 42.02124 42.02000 42.02000 42.02099 42.01986 42.01899 42.01885 42.01891
## [1497] 42.01883 42.01739 42.01680 42.01890 42.02048 42.02029 42.02039 42.01954
## [1505] 42.02011 42.02158 42.01797 42.01796 42.01980 42.01996 42.01854 42.01787
## [1513] 42.02224 42.02224 42.02209 42.02145 42.01999 42.01872 42.01741 42.01708
## [1521] 42.01927 42.01909 42.01900 42.01791 42.01783 42.01769 42.01789 42.01689
## [1529] 42.01687 42.01607 42.01983 42.01962 42.01844 42.01889 42.01877 42.01753
## [1537] 42.01677 42.01601 42.01567 42.01457 42.01264 42.01249 42.00949 42.00949
## [1545] 42.02261 42.02181 42.02280 42.02252 42.02138 42.02128 42.02140 42.02038
## [1553] 42.02139 42.01881 42.02116 42.02151 42.02159 42.02246 42.02263 42.00933
## [1561] 42.00880 42.00881 42.00919 41.99507 41.99712 41.99706 41.99707 41.99692
## [1569] 41.99690 41.99599 41.99854 42.00514 41.99900 41.99514 41.99484 41.99373
## [1577] 41.99514 41.99375 41.99518 41.99981 41.99958 41.99761 41.99855 41.99551
## [1585] 41.99639 41.99423 41.99312 41.99361 41.99338 41.99334 41.99005 41.99306
## [1593] 41.99318 41.99261 41.99137 41.99090 41.99188 41.99178 41.99325 41.99213
## [1601] 41.99193 41.98868 41.98711 41.98893 41.98811 41.98808 41.98729 41.98658
## [1609] 41.98650 41.99189 42.02275 42.05130 42.05136 42.04967 42.05044 42.06127
## [1617] 42.06047 42.06166 42.06133 42.06124 42.06307 42.05959 42.05895 42.05956
## [1625] 42.05952 42.05924 42.05929 42.05701 42.05817 42.05753 42.05702 42.05938
## [1633] 42.05924 42.05931 42.05930 42.06043 42.06014 42.06120 42.06010 42.05974
## [1641] 42.05889 42.05900 42.05900 42.05437 42.05392 42.05555 42.05450 42.05439
## [1649] 42.05494 42.05522 42.05443 42.05622 42.05621 42.05235 42.05327 42.05026
## [1657] 42.04983 42.05007 42.05037 42.04996 42.05104 42.04877 42.05539 42.05526
## [1665] 42.05460 42.05341 42.05461 42.05355 42.05346 42.05350 42.05630 42.05639
## [1673] 42.05467 42.05479 42.05532 42.05267 42.05270 42.05270 42.05229 42.05183
## [1681] 42.05040 42.05057 42.05000 42.04882 42.06257 42.06330 42.06330 42.06327
## [1689] 42.06327 42.06314 42.06195 42.06198 42.06197 42.06128 42.06226 42.06282
## [1697] 42.06241 42.06243 42.06150 42.06113 42.06027 42.06027 42.06113 42.06058
## [1705] 42.06043 42.05956 42.06113 42.05976 42.05968 42.05863 42.05862 42.05921
## [1713] 42.05920 42.05920 42.05919 42.05918 42.05918 42.05918 42.05918 42.05826
## [1721] 42.05828 42.05828 42.05683 42.05712 42.05738 42.05685 42.05721 42.06339
## [1729] 42.06338 42.06248 42.06228 42.06267 42.06208 42.06207 42.06253 42.06289
## [1737] 42.06293 42.06297 42.06302 42.06311 42.06334 42.06258 42.06258 42.06295
## [1745] 42.06330 42.06171 42.06143 42.06148 42.06154 42.06153 42.06132 42.06139
## [1753] 42.06002 42.05886 42.05918 42.05798 42.05840 42.05953 42.05784 42.05890
## [1761] 42.05332 42.05336 42.05575 42.05563 42.05359 42.05344 42.05347 42.05198
## [1769] 42.05236 42.05225 42.05181 42.05048 42.04952 42.04951 42.05242 42.05441
## [1777] 42.05576 42.05432 42.05418 42.05442 42.05435 42.05431 42.05451 42.05415
## [1785] 42.05415 42.05415 42.05416 42.05318 42.05319 42.05323 42.05311 42.05303
## [1793] 42.05307 42.05180 42.05179 42.05165 42.05166 42.05021 42.05030 42.05168
## [1801] 42.05227 42.05212 42.05138 42.05227 42.05215 42.05105 42.03759 42.03764
## [1809] 42.03764 42.03722 42.03603 42.03526 42.03604 42.03507 42.03646 42.03457
## [1817] 42.03568 42.03621 42.03468 42.03595 42.03587 42.03586 42.03453 42.03472
## [1825] 42.03567 42.03465 42.03482 42.03765 42.03750 42.03466 42.03464 42.03467
## [1833] 42.04738 42.04578 42.04762 42.04742 42.04741 42.04771 42.04771 42.04770
## [1841] 42.04769 42.04767 42.04710 42.04710 42.04614 42.04615 42.04615 42.04641
## [1849] 42.04717 42.04632 42.04614 42.04613 42.04641 42.04438 42.04476 42.04499
## [1857] 42.04299 42.04375 42.04305 42.04296 42.03710 42.03497 42.04849 42.04655
## [1865] 42.04776 42.04807 42.04849 42.04762 42.04336 42.04301 42.04548 42.04259
## [1873] 42.04389 42.04217 42.04781 42.04829 42.04650 42.04476 42.04398 42.04393
## [1881] 42.04312 42.04212 42.04227 42.04404 42.04258 42.04248 42.04248 42.04492
## [1889] 42.04361 42.04363 42.03914 42.04014 42.04174 42.03907 42.04139 42.04176
## [1897] 42.04074 42.04125 42.04133 42.03921 42.03920 42.03657 42.03568 42.03602
## [1905] 42.03622 42.03462 42.03620 42.03465 42.03491 42.03475 42.04922 42.04893
## [1913] 42.04671 42.04581 42.04593 42.04688 42.04405 42.04314 42.04474 42.04393
## [1921] 42.04395 42.04231 42.04422 42.04567 42.04395 42.04323 42.04270 42.04273
## [1929] 42.04302 42.04187 42.04124 42.04077 42.04045 42.03981 42.03937 42.03948
## [1937] 42.03846 42.03837 42.03848 42.03854 42.03479 42.03472 42.03481 42.03593
## [1945] 42.03608 42.03603 42.03725 42.03810 42.03570 42.03578 42.03585 42.03488
## [1953] 42.03934 42.03998 42.04015 42.03889 42.03573 42.03592 42.03733 42.03720
## [1961] 42.03595 42.03588 42.03711 42.03586 42.03600 42.03708 42.03469 42.03461
## [1969] 42.03456 42.03695 42.03542 42.03460 42.03503 42.03440 42.03225 42.03147
## [1977] 42.03131 42.03137 42.03111 42.03341 42.03234 42.03339 42.03349 42.03331
## [1985] 42.03228 42.03236 42.03058 42.03180 42.03031 42.03013 42.02998 42.02922
## [1993] 42.02920 42.02922 42.02804 42.02808 42.02885 42.02819 42.02715 42.02820
## [2001] 42.02776 42.02609 42.02681 42.02726 42.02727 42.03331 42.03434 42.03185
## [2009] 42.03082 42.03348 42.03361 42.03334 42.03240 42.03143 42.03144 42.03142
## [2017] 42.03251 42.03239 42.03132 42.03148 42.03149 42.03045 42.02945 42.02691
## [2025] 42.02700 42.02709 42.02997 42.02911 42.02752 42.02758 42.02523 42.02543
## [2033] 42.02158 42.02500 42.02423 42.02270 42.02389 42.02296 42.02384 42.02352
## [2041] 42.02672 42.02589 42.02616 42.02295 42.03352 42.03328 42.02805 42.02800
## [2049] 42.03327 42.03243 42.03445 42.03443 42.03319 42.03238 42.03210 42.03153
## [2057] 42.03120 42.02956 42.03122 42.02968 42.03114 42.02979 42.03059 42.03442
## [2065] 42.03253 42.03442 42.03356 42.03264 42.03373 42.03402 42.03270 42.02868
## [2073] 42.02384 42.02409 42.02618 42.02323 42.02302 42.02842 42.02824 42.02495
## [2081] 42.02582 42.02500 42.02052 42.02019 42.02469 42.02387 42.02394 42.02345
## [2089] 42.02282 42.03445 42.03366 42.03358 42.03155 42.03247 42.03157 42.03080
## [2097] 42.02667 42.02664 42.02567 42.02571 42.02559 42.02526 42.02528 42.02556
## [2105] 42.02460 42.02459 42.02463 42.02461 42.02573 42.02516 42.02464 42.02448
## [2113] 42.02323 42.02969 42.02980 42.02796 42.02797 42.02473 42.02430 42.02223
## [2121] 42.02132 42.02121 42.02097 42.02091 42.01903 42.01908 42.01903 42.01838
## [2129] 42.01791 42.01783 42.01690 42.01683 42.01630 42.01621 42.02199 42.02139
## [2137] 42.02186 42.02224 42.01972 42.02116 42.02106 42.02103 42.02086 42.02174
## [2145] 42.01852 42.01842 42.01885 42.01697 42.01696 42.01621 42.01720 42.01856
## [2153] 42.01886 42.01853 42.01802 42.01786 42.01628 42.01775 42.01597 42.01404
## [2161] 42.01404 42.01482 42.01358 42.01354 42.01343 42.01347 42.01469 42.02103
## [2169] 42.02070 42.02122 42.02119 42.02103 42.02089 42.01993 42.01970 42.01990
## [2177] 42.01917 42.01828 42.01887 42.01892 42.01645 42.01664 42.01897 42.01897
## [2185] 42.02006 42.02078 42.02050 42.02266 42.02144 42.01796 42.01965 42.01894
## [2193] 42.01873 42.01994 42.01902 42.01876 42.01804 42.01804 42.01789 42.01789
## [2201] 42.01892 42.01710 42.01711 42.01628 42.01622 42.01622 42.01622 42.01685
## [2209] 42.01607 42.01995 42.01995 42.01880 42.01757 42.02004 42.01850 42.01894
## [2217] 42.01689 42.01757 42.01696 42.01629 42.01300 42.01455 42.01265 42.01249
## [2225] 42.01151 42.00949 42.01184 42.00932 42.01134 42.01033 42.01008 42.02284
## [2233] 42.02247 42.02266 42.02266 42.02257 42.02269 42.02280 42.02155 42.02128
## [2241] 42.02150 42.02016 42.02008 42.02142 42.02155 42.00940 42.00940 41.99436
## [2249] 41.99421 41.99421 41.99522 41.99597 41.99593 41.99551 41.99419 42.00569
## [2257] 42.00169 41.99854 42.00041 41.99424 41.99414 41.99410 41.99482 41.99509
## [2265] 41.99802 41.99799 41.99736 41.99713 41.99609 41.99551 41.99730 41.99534
## [2273] 41.99399 41.99735 41.99329 41.99298 41.99326 41.99005 41.99232 41.99336
## [2281] 41.99311 41.99279 41.99121 41.99292 41.99284 41.99215 41.99252 41.99248
## [2289] 41.99171 41.99164 41.99187 41.99161 41.99110 41.98958 41.99026 41.98949
## [2297] 41.99114 41.99137 41.98991 41.98761 41.98977 41.98822 41.98817 41.98651
## [2305] 41.98659 42.05341 42.05333 42.05311 42.05252 42.05051 42.05991 42.06141
## [2313] 42.06043 42.06048 42.06074 42.06188 42.06283 42.06281 42.06280 42.06280
## [2321] 42.06289 42.06296 42.06063 42.05910 42.05842 42.05820 42.05960 42.05933
## [2329] 42.05814 42.06046 42.06053 42.06052 42.06148 42.06101 42.06011 42.05996
## [2337] 42.06110 42.05939 42.05670 42.05856 42.05995 42.05882 42.05561 42.05437
## [2345] 42.05434 42.05496 42.05626 42.05327 42.05240 42.05342 42.05051 42.05018
## [2353] 42.04977 42.05108 42.05236 42.05252 42.05236 42.04977 42.04961 42.05513
## [2361] 42.05346 42.05446 42.05439 42.05532 42.05266 42.05247 42.05270 42.05180
## [2369] 42.05184 42.05169 42.05167 42.05167 42.05166 42.05060 42.04895 42.04895
## [2377] 42.05071 42.05044 42.06147 42.06331 42.06190 42.06128 42.06281 42.06211
## [2385] 42.06113 42.06034 42.06027 42.06027 42.06043 42.06083 42.06058 42.06087
## [2393] 42.06042 42.06036 42.06038 42.06154 42.06138 42.05964 42.05982 42.05854
## [2401] 42.05869 42.05922 42.05902 42.05918 42.05716 42.05823 42.05824 42.05716
## [2409] 42.05716 42.05758 42.06330 42.06244 42.06236 42.06191 42.06193 42.06303
## [2417] 42.06258 42.06198 42.06199 42.06257 42.06330 42.06333 42.06318 42.06237
## [2425] 42.06230 42.06330 42.06330 42.06330 42.06283 42.06287 42.06024 42.06029
## [2433] 42.05980 42.06173 42.06084 42.06007 42.05865 42.05853 42.05790 42.05784
## [2441] 42.05696 42.05846 42.05443 42.05384 42.05380 42.05331 42.05560 42.05315
## [2449] 42.05347 42.05178 42.04942 42.05107 42.04969 42.04953 42.05506 42.05432
## [2457] 42.05418 42.05418 42.05318 42.05303 42.05180 42.05226 42.05227 42.05235
## [2465] 42.05136 42.05143 42.05138 42.05125 42.05123 42.05122 42.03784 42.03615
## [2473] 42.03791 42.03791 42.03783 42.03661 42.03718 42.03685 42.03460 42.03674
## [2481] 42.03556 42.03563 42.03748 42.03580 42.03524 42.03605 42.03609 42.03621
## [2489] 42.03530 42.03509 42.03469 42.03570 42.03584 42.03533 42.03554 42.03457
## [2497] 42.03458 42.03466 42.03706 42.04576 42.04861 42.04837 42.04800 42.04806
## [2505] 42.04759 42.04768 42.04767 42.04764 42.04762 42.04760 42.04615 42.04643
## [2513] 42.04613 42.04611 42.04643 42.04642 42.04461 42.04378 42.04237 42.04248
## [2521] 42.04283 42.04335 42.03741 42.04563 42.04820 42.04900 42.04860 42.04647
## [2529] 42.04806 42.04355 42.04467 42.04650 42.04571 42.04573 42.04573 42.04588
## [2537] 42.04588 42.04795 42.04657 42.04568 42.04431 42.04321 42.04498 42.04211
## [2545] 42.04007 42.03843 42.04102 42.04141 42.04131 42.04027 42.04014 42.03919
## [2553] 42.03835 42.03920 42.03584 42.03537 42.03538 42.03567 42.03806 42.03461
## [2561] 42.03602 42.03746 42.03483 42.03599 42.03725 42.04930 42.04903 42.04665
## [2569] 42.04853 42.04230 42.04216 42.04313 42.04474 42.04390 42.04498 42.04498
## [2577] 42.04425 42.04210 42.04363 42.04270 42.04290 42.04182 42.04075 42.04025
## [2585] 42.03945 42.03836 42.03939 42.03860 42.03835 42.03884 42.04028 42.03935
## [2593] 42.03493 42.03462 42.03473 42.03786 42.03718 42.03700 42.03710 42.03704
## [2601] 42.03580 42.03678 42.03586 42.03457 42.03528 42.03476 42.04094 42.04081
## [2609] 42.04010 42.04008 42.04030 42.04015 42.03850 42.03850 42.03830 42.04187
## [2617] 42.04165 42.04077 42.03902 42.03610 42.03578 42.03600 42.03610 42.03605
## [2625] 42.03480 42.03471 42.03466 42.03483 42.03476 42.03471 42.03679 42.03548
## [2633] 42.03462 42.03475 42.03456 42.03327 42.03350 42.03339 42.03214 42.03235
## [2641] 42.03139 42.03151 42.03135 42.03121 42.03040 42.03046 42.03120 42.03362
## [2649] 42.03305 42.03231 42.03027 42.03206 42.03227 42.03196 42.03037 42.03146
## [2657] 42.03114 42.02914 42.02799 42.02914 42.02895 42.02928 42.02795 42.02811
## [2665] 42.02795 42.02896 42.02902 42.02799 42.02727 42.02297 42.03405 42.03179
## [2673] 42.03079 42.03331 42.03337 42.03230 42.03150 42.03240 42.03248 42.03142
## [2681] 42.03144 42.03242 42.03242 42.03108 42.03039 42.02953 42.02686 42.02854
## [2689] 42.02689 42.02689 42.02526 42.02546 42.02495 42.02424 42.02426 42.02312
## [2697] 42.02282 42.02381 42.03441 42.03458 42.03330 42.03315 42.02780 42.02800
## [2705] 42.02800 42.03352 42.03339 42.03330 42.03224 42.03230 42.03238 42.03136
## [2713] 42.03061 42.03222 42.03070 42.03136 42.03135 42.03023 42.03023 42.03328
## [2721] 42.03247 42.03229 42.03231 42.03369 42.03248 42.03164 42.03216 42.03258
## [2729] 42.02945 42.02884 42.02365 42.02470 42.02653 42.02766 42.02667 42.02661
## [2737] 42.02821 42.02819 42.02646 42.02585 42.02552 42.02510 42.02520 42.02631
## [2745] 42.02475 42.02282 42.03208 42.03249 42.03093 42.03152 42.02460 42.02464
## [2753] 42.02455 42.02461 42.02569 42.02522 42.02519 42.02437 42.02448 42.02331
## [2761] 42.02438 42.02331 42.02314 42.02988 42.02817 42.02557 42.02438 42.02459
## [2769] 42.02223 42.02208 42.02121 42.02090 42.01905 42.01804 42.01812 42.01763
## [2777] 42.01687 42.01793 42.01686 42.01622 42.01691 42.02139 42.02014 42.02007
## [2785] 42.01970 42.02181 42.02189 42.02100 42.02179 42.01834 42.01866 42.01785
## [2793] 42.01705 42.01819 42.01449 42.01698 42.01690 42.01869 42.01786 42.01787
## [2801] 42.01786 42.01719 42.01628 42.01625 42.01662 42.01591 42.01609 42.01464
## [2809] 42.01411 42.01579 42.01600 42.01397 42.01399 42.01576 42.01352 42.01459
## [2817] 42.01959 42.02099 42.02034 42.02059 42.02055 42.02111 42.01899 42.02117
## [2825] 42.01865 42.01896 42.01896 42.01896 42.01881 42.01881 42.01881 42.01881
## [2833] 42.01881 42.02011 42.02268 42.02266 42.02264 42.02143 42.01885 42.01773
## [2841] 42.02209 42.02209 42.02197 42.02156 42.02087 42.02029 42.01884 42.01828
## [2849] 42.01924 42.01909 42.01895 42.01881 42.01886 42.01694 42.01658 42.01609
## [2857] 42.01607 42.01761 42.01874 42.01906 42.01976 42.01713 42.01683 42.01638
## [2865] 42.01759 42.01579 42.01612 42.01330 42.01349 42.01494 42.01359 42.01179
## [2873] 42.02279 42.02271 42.02044 42.02085 42.02021 42.02153 42.01902 42.01935
## [2881] 42.01910 42.01910 42.02260 42.00940 42.00838 41.99685 41.99436 41.99599
## [2889] 41.99691 41.99570 41.99606 42.00451 42.00097 41.99984 41.99791 41.99955
## [2897] 41.99802 41.99934 41.99633 41.99435 41.99423 41.99326 41.99305 41.99222
## [2905] 41.99329 41.99317 41.99315 41.99005 41.99099 41.99052 41.99173 41.99306
## [2913] 41.99157 41.99230 41.99216 41.99171 41.99253 41.99252 41.99250 41.99171
## [2921] 41.99153 41.99004 41.99013 41.98874 41.98769 41.98896 41.98831 41.98651
## [2929] 41.99092 41.98927
ames$Longitude
##    [1] -93.61975 -93.61976 -93.61939 -93.61732 -93.63893 -93.63893 -93.63379
##    [8] -93.63383 -93.63285 -93.63907 -93.63695 -93.63595 -93.63865 -93.63463
##   [15] -93.63291 -93.62655 -93.62360 -93.62880 -93.62297 -93.63666 -93.63396
##   [22] -93.63637 -93.63937 -93.62623 -93.62654 -93.62881 -93.62711 -93.62277
##   [29] -93.62437 -93.62727 -93.62754 -93.62728 -93.62723 -93.62592 -93.62597
##   [36] -93.62585 -93.65320 -93.65465 -93.65373 -93.65283 -93.65271 -93.65332
##   [43] -93.65234 -93.65231 -93.65505 -93.65044 -93.65607 -93.65485 -93.65412
##   [50] -93.64945 -93.65101 -93.64436 -93.64163 -93.64225 -93.64408 -93.64149
##   [57] -93.64245 -93.64132 -93.64134 -93.65138 -93.65233 -93.65421 -93.65362
##   [64] -93.65036 -93.65058 -93.65297 -93.64237 -93.64998 -93.65032 -93.64710
##   [71] -93.64277 -93.64222 -93.69231 -93.69241 -93.68981 -93.68642 -93.68769
##   [78] -93.68626 -93.68626 -93.68628 -93.68503 -93.68503 -93.67907 -93.67721
##   [85] -93.67367 -93.67096 -93.67210 -93.67207 -93.66969 -93.66067 -93.64675
##   [92] -93.65177 -93.65277 -93.64799 -93.64614 -93.64560 -93.64489 -93.64489
##   [99] -93.64489 -93.64689 -93.64561 -93.64558 -93.64548 -93.64548 -93.64125
##  [106] -93.64356 -93.64936 -93.64814 -93.65921 -93.63952 -93.63335 -93.63688
##  [113] -93.63746 -93.63753 -93.63779 -93.63143 -93.63255 -93.62614 -93.62593
##  [120] -93.62338 -93.62694 -93.62754 -93.62678 -93.62781 -93.62573 -93.62290
##  [127] -93.62286 -93.62065 -93.62050 -93.62929 -93.62931 -93.62475 -93.62460
##  [134] -93.62143 -93.61822 -93.61842 -93.61431 -93.61424 -93.61769 -93.61755
##  [141] -93.61712 -93.61770 -93.61863 -93.61658 -93.61488 -93.61330 -93.61513
##  [148] -93.61424 -93.61362 -93.61562 -93.61709 -93.61722 -93.61843 -93.61548
##  [155] -93.61164 -93.61879 -93.61920 -93.62034 -93.61168 -93.60728 -93.60690
##  [162] -93.61061 -93.60597 -93.60922 -93.60780 -93.60594 -93.60674 -93.60483
##  [169] -93.60371 -93.61701 -93.61714 -93.62024 -93.62024 -93.61791 -93.61400
##  [176] -93.61543 -93.61545 -93.60895 -93.60908 -93.61047 -93.61047 -93.60679
##  [183] -93.60679 -93.61689 -93.61059 -93.61221 -93.60888 -93.60825 -93.60724
##  [190] -93.60875 -93.60709 -93.62574 -93.62583 -93.62559 -93.62457 -93.62260
##  [197] -93.62352 -93.62245 -93.62348 -93.62148 -93.62041 -93.62338 -93.62167
##  [204] -93.62184 -93.62163 -93.62950 -93.62688 -93.62039 -93.65623 -93.65554
##  [211] -93.65827 -93.65692 -93.65808 -93.65757 -93.65585 -93.67204 -93.67102
##  [218] -93.67108 -93.67340 -93.67238 -93.67814 -93.67666 -93.67799 -93.67665
##  [225] -93.67667 -93.67151 -93.67356 -93.67252 -93.67050 -93.67011 -93.66362
##  [232] -93.66838 -93.66790 -93.67504 -93.67118 -93.66276 -93.66357 -93.66771
##  [239] -93.66526 -93.66610 -93.66589 -93.66580 -93.66337 -93.65981 -93.68500
##  [246] -93.68599 -93.68317 -93.68203 -93.68016 -93.69019 -93.69127 -93.69245
##  [253] -93.69252 -93.68329 -93.68018 -93.68214 -93.69188 -93.69054 -93.68826
##  [260] -93.68773 -93.68781 -93.68775 -93.68373 -93.68786 -93.68552 -93.68776
##  [267] -93.68086 -93.68269 -93.68123 -93.68370 -93.68792 -93.68402 -93.68596
##  [274] -93.67833 -93.67748 -93.67755 -93.67754 -93.67637 -93.67527 -93.67517
##  [281] -93.66634 -93.66751 -93.66484 -93.66017 -93.66194 -93.65576 -93.65565
##  [288] -93.65672 -93.65032 -93.65144 -93.65076 -93.64678 -93.64668 -93.64004
##  [295] -93.63978 -93.64279 -93.63968 -93.64501 -93.64570 -93.64573 -93.64168
##  [302] -93.62529 -93.62950 -93.62669 -93.61527 -93.60659 -93.60827 -93.60775
##  [309] -93.61599 -93.60495 -93.60282 -93.60397 -93.60080 -93.60200 -93.65789
##  [316] -93.64920 -93.64920 -93.64910 -93.64947 -93.64764 -93.64681 -93.65250
##  [323] -93.60743 -93.60837 -93.60762 -93.60822 -93.60359 -93.60432 -93.60443
##  [330] -93.60344 -93.60359 -93.60181 -93.60340 -93.60145 -93.60686 -93.60354
##  [337] -93.60303 -93.60058 -93.59979 -93.60001 -93.60015 -93.61890 -93.61956
##  [344] -93.61552 -93.63939 -93.63992 -93.63640 -93.63584 -93.63744 -93.63426
##  [351] -93.63374 -93.63355 -93.63302 -93.63921 -93.63907 -93.63907 -93.63696
##  [358] -93.63705 -93.63767 -93.63582 -93.63575 -93.63682 -93.63605 -93.63932
##  [365] -93.63464 -93.63325 -93.62775 -93.62864 -93.62936 -93.62398 -93.63734
##  [372] -93.63651 -93.63506 -93.63242 -93.63777 -93.63745 -93.63837 -93.63421
##  [379] -93.63704 -93.63380 -93.63729 -93.63469 -93.63592 -93.63952 -93.63951
##  [386] -93.63923 -93.63937 -93.63291 -93.63218 -93.63025 -93.63143 -93.63179
##  [393] -93.62586 -93.62655 -93.62871 -93.62660 -93.62655 -93.62704 -93.62566
##  [400] -93.62554 -93.62350 -93.62233 -93.62740 -93.62812 -93.62757 -93.62946
##  [407] -93.62985 -93.62988 -93.62724 -93.62473 -93.62475 -93.62595 -93.62592
##  [414] -93.62585 -93.62569 -93.62569 -93.62570 -93.62571 -93.62560 -93.62652
##  [421] -93.65888 -93.65627 -93.65272 -93.65783 -93.65748 -93.65424 -93.65414
##  [428] -93.65395 -93.65286 -93.65444 -93.65457 -93.65822 -93.65783 -93.65705
##  [435] -93.65444 -93.65885 -93.65674 -93.65315 -93.65288 -93.65472 -93.65476
##  [442] -93.65208 -93.65465 -93.65272 -93.65413 -93.65478 -93.64377 -93.65604
##  [449] -93.65486 -93.65452 -93.65048 -93.64978 -93.64974 -93.64972 -93.65026
##  [456] -93.65024 -93.65548 -93.65524 -93.65420 -93.65421 -93.65440 -93.65148
##  [463] -93.65101 -93.65024 -93.64575 -93.64466 -93.64446 -93.64563 -93.64156
##  [470] -93.64151 -93.64209 -93.64260 -93.64291 -93.64291 -93.64275 -93.64331
##  [477] -93.64411 -93.64177 -93.64149 -93.64263 -93.64005 -93.64007 -93.64965
##  [484] -93.64281 -93.64293 -93.64192 -93.64425 -93.64138 -93.64131 -93.64130
##  [491] -93.64137 -93.64140 -93.64131 -93.64088 -93.65429 -93.65716 -93.65096
##  [498] -93.65261 -93.65238 -93.65556 -93.65599 -93.65405 -93.65150 -93.65212
##  [505] -93.65053 -93.65082 -93.64412 -93.64247 -93.64255 -93.64394 -93.64258
##  [512] -93.64236 -93.64125 -93.64130 -93.65032 -93.65040 -93.65044 -93.65042
##  [519] -93.65038 -93.65034 -93.64827 -93.64881 -93.64878 -93.64300 -93.64287
##  [526] -93.64384 -93.64247 -93.64375 -93.64085 -93.64175 -93.64313 -93.64150
##  [533] -93.63983 -93.69153 -93.69202 -93.69247 -93.69143 -93.69151 -93.69171
##  [540] -93.69187 -93.69107 -93.69016 -93.69040 -93.68914 -93.69139 -93.68880
##  [547] -93.68784 -93.68626 -93.68627 -93.68780 -93.68504 -93.68396 -93.67908
##  [554] -93.68020 -93.67636 -93.67701 -93.67405 -93.67572 -93.66955 -93.66954
##  [561] -93.67352 -93.66971 -93.66984 -93.66033 -93.65033 -93.64777 -93.64632
##  [568] -93.64634 -93.64704 -93.64692 -93.64420 -93.64566 -93.64548 -93.64180
##  [575] -93.64592 -93.64936 -93.64806 -93.64555 -93.65921 -93.63534 -93.63463
##  [582] -93.63480 -93.63197 -93.63188 -93.63152 -93.63316 -93.63185 -93.63254
##  [589] -93.63282 -93.63796 -93.63791 -93.63265 -93.63190 -93.62890 -93.62705
##  [596] -93.62461 -93.63071 -93.62742 -93.62578 -93.62299 -93.62247 -93.62349
##  [603] -93.62264 -93.62248 -93.62341 -93.62941 -93.62866 -93.62731 -93.62784
##  [610] -93.62798 -93.62698 -93.62458 -93.62269 -93.62170 -93.62950 -93.62919
##  [617] -93.62563 -93.62568 -93.62371 -93.62461 -93.62264 -93.62249 -93.61967
##  [624] -93.61796 -93.61569 -93.61560 -93.61556 -93.61462 -93.61422 -93.61903
##  [631] -93.61878 -93.61864 -93.61488 -93.61643 -93.61475 -93.61406 -93.61477
##  [638] -93.61588 -93.61454 -93.62002 -93.61848 -93.61655 -93.61855 -93.61904
##  [645] -93.61722 -93.61240 -93.61721 -93.61679 -93.61861 -93.61861 -93.62008
##  [652] -93.62034 -93.61719 -93.61718 -93.61462 -93.61442 -93.61239 -93.61224
##  [659] -93.61405 -93.61547 -93.61405 -93.61390 -93.61174 -93.61203 -93.61225
##  [666] -93.61225 -93.61035 -93.60765 -93.61059 -93.61065 -93.60820 -93.60893
##  [673] -93.60923 -93.60657 -93.60552 -93.60906 -93.60898 -93.61061 -93.60890
##  [680] -93.60891 -93.60689 -93.60689 -93.60781 -93.60597 -93.60679 -93.60994
##  [687] -93.61061 -93.60499 -93.60405 -93.60388 -93.60370 -93.61797 -93.61558
##  [694] -93.61901 -93.61950 -93.62026 -93.61560 -93.61560 -93.61954 -93.61849
##  [701] -93.62007 -93.62022 -93.61764 -93.61538 -93.61291 -93.60887 -93.60778
##  [708] -93.60883 -93.60684 -93.60684 -93.61047 -93.60679 -93.61865 -93.61705
##  [715] -93.61787 -93.61389 -93.61225 -93.61059 -93.61308 -93.61058 -93.61486
##  [722] -93.61004 -93.60765 -93.60870 -93.60871 -93.60706 -93.60521 -93.60542
##  [729] -93.62894 -93.62885 -93.62798 -93.62597 -93.62612 -93.62456 -93.62457
##  [736] -93.62251 -93.62351 -93.62168 -93.62244 -93.62469 -93.62364 -93.62364
##  [743] -93.62364 -93.62363 -93.62256 -93.62135 -93.62136 -93.62726 -93.62536
##  [750] -93.62153 -93.62150 -93.62151 -93.62163 -93.62163 -93.62830 -93.62672
##  [757] -93.62186 -93.62171 -93.62504 -93.62515 -93.59958 -93.65921 -93.65821
##  [764] -93.65817 -93.65711 -93.66019 -93.65724 -93.65838 -93.66984 -93.67540
##  [771] -93.67229 -93.67338 -93.67799 -93.67265 -93.67185 -93.67137 -93.66643
##  [778] -93.66636 -93.66681 -93.66213 -93.65928 -93.66088 -93.66798 -93.67118
##  [785] -93.67216 -93.66624 -93.66600 -93.66599 -93.66778 -93.66821 -93.66721
##  [792] -93.66527 -93.66527 -93.66015 -93.66337 -93.66419 -93.66104 -93.66101
##  [799] -93.65985 -93.68471 -93.68573 -93.68599 -93.68341 -93.68106 -93.68202
##  [806] -93.68218 -93.67979 -93.67978 -93.67914 -93.67911 -93.67978 -93.67953
##  [813] -93.67910 -93.67915 -93.67978 -93.67986 -93.67913 -93.67918 -93.68707
##  [820] -93.68886 -93.68899 -93.68909 -93.69027 -93.69097 -93.69036 -93.68907
##  [827] -93.68907 -93.69143 -93.69157 -93.69011 -93.68907 -93.68908 -93.68512
##  [834] -93.67930 -93.68256 -93.69188 -93.69213 -93.68837 -93.69040 -93.69136
##  [841] -93.69159 -93.69287 -93.69055 -93.69032 -93.69167 -93.68838 -93.68943
##  [848] -93.69118 -93.68944 -93.68994 -93.69167 -93.68863 -93.68692 -93.68721
##  [855] -93.68722 -93.68507 -93.68489 -93.68495 -93.68316 -93.68316 -93.68317
##  [862] -93.68351 -93.68782 -93.68622 -93.68269 -93.68262 -93.68298 -93.68674
##  [869] -93.68567 -93.67915 -93.68071 -93.68111 -93.68134 -93.68141 -93.68078
##  [876] -93.68218 -93.68167 -93.68130 -93.68144 -93.68412 -93.68533 -93.68605
##  [883] -93.68454 -93.68430 -93.67831 -93.67741 -93.67637 -93.67223 -93.67552
##  [890] -93.67397 -93.67525 -93.67661 -93.66995 -93.66478 -93.66631 -93.66469
##  [897] -93.66446 -93.65919 -93.66481 -93.66440 -93.66434 -93.66330 -93.65796
##  [904] -93.65666 -93.65567 -93.65183 -93.65203 -93.65179 -93.64857 -93.65167
##  [911] -93.64656 -93.64841 -93.64665 -93.64843 -93.64594 -93.64688 -93.64842
##  [918] -93.64470 -93.64464 -93.64281 -93.63970 -93.64237 -93.64012 -93.64059
##  [925] -93.64429 -93.64400 -93.63943 -93.64294 -93.64289 -93.64076 -93.64136
##  [932] -93.63941 -93.64574 -93.64574 -93.64574 -93.64574 -93.64486 -93.64354
##  [939] -93.64395 -93.62529 -93.62684 -93.62530 -93.61689 -93.61540 -93.61657
##  [946] -93.58823 -93.60694 -93.60725 -93.60706 -93.60372 -93.60472 -93.60449
##  [953] -93.60452 -93.60236 -93.60242 -93.64978 -93.65212 -93.65315 -93.65183
##  [960] -93.65187 -93.65304 -93.65045 -93.64646 -93.64488 -93.64460 -93.64679
##  [967] -93.64664 -93.64682 -93.65170 -93.60426 -93.60836 -93.60820 -93.60821
##  [974] -93.60425 -93.60422 -93.60348 -93.60064 -93.60311 -93.60197 -93.60176
##  [981] -93.60113 -93.60669 -93.60493 -93.60387 -93.60278 -93.60945 -93.60693
##  [988] -93.59987 -93.60019 -93.61883 -93.61633 -93.61742 -93.63906 -93.63780
##  [995] -93.63877 -93.63878 -93.63748 -93.63672 -93.63695 -93.63388 -93.63378
## [1002] -93.63390 -93.63380 -93.63915 -93.63901 -93.63776 -93.63680 -93.63728
## [1009] -93.63201 -93.63158 -93.62883 -93.62890 -93.62791 -93.62230 -93.62981
## [1016] -93.62257 -93.63915 -93.63971 -93.63783 -93.63781 -93.63651 -93.63650
## [1023] -93.63511 -93.63030 -93.63317 -93.63515 -93.63711 -93.63719 -93.63390
## [1030] -93.63522 -93.63418 -93.63113 -93.63023 -93.62654 -93.62788 -93.62651
## [1037] -93.63021 -93.62426 -93.62487 -93.62897 -93.62763 -93.62711 -93.62946
## [1044] -93.62765 -93.62731 -93.62942 -93.62723 -93.62599 -93.62585 -93.62649
## [1051] -93.65887 -93.65731 -93.65294 -93.65274 -93.65464 -93.65666 -93.65783
## [1058] -93.65789 -93.65719 -93.65680 -93.65369 -93.65473 -93.65444 -93.65696
## [1065] -93.65572 -93.65255 -93.65122 -93.65559 -93.65491 -93.65419 -93.65137
## [1072] -93.64974 -93.65387 -93.65277 -93.65184 -93.65116 -93.65142 -93.65000
## [1079] -93.65061 -93.65089 -93.64018 -93.64096 -93.64244 -93.64304 -93.64233
## [1086] -93.64264 -93.64352 -93.64396 -93.64177 -93.64142 -93.64005 -93.64959
## [1093] -93.64133 -93.64240 -93.64427 -93.64363 -93.64134 -93.64096 -93.65459
## [1100] -93.65412 -93.65332 -93.65093 -93.65156 -93.65528 -93.65604 -93.65393
## [1107] -93.65579 -93.65391 -93.65180 -93.65100 -93.64411 -93.64414 -93.64233
## [1114] -93.64211 -93.64209 -93.64238 -93.64398 -93.64403 -93.64402 -93.64379
## [1121] -93.65040 -93.65036 -93.64826 -93.64384 -93.64119 -93.64092 -93.64304
## [1128] -93.64369 -93.64165 -93.64402 -93.63966 -93.69250 -93.69243 -93.69169
## [1135] -93.69184 -93.69076 -93.69036 -93.69125 -93.69124 -93.68977 -93.68902
## [1142] -93.68898 -93.68980 -93.68880 -93.68782 -93.68783 -93.68785 -93.68643
## [1149] -93.68628 -93.67409 -93.67574 -93.67396 -93.67100 -93.67250 -93.67055
## [1156] -93.66203 -93.66052 -93.65348 -93.65323 -93.64749 -93.64606 -93.64611
## [1163] -93.64489 -93.64489 -93.64489 -93.64417 -93.64412 -93.64698 -93.64699
## [1170] -93.64699 -93.64705 -93.64704 -93.64375 -93.64373 -93.64338 -93.64178
## [1177] -93.64177 -93.64174 -93.64766 -93.64871 -93.64895 -93.64571 -93.65558
## [1184] -93.65832 -93.63923 -93.63510 -93.63433 -93.63498 -93.63153 -93.63178
## [1191] -93.63245 -93.63125 -93.62644 -93.62614 -93.62476 -93.62462 -93.62463
## [1198] -93.63042 -93.62908 -93.63015 -93.63073 -93.62592 -93.62450 -93.62580
## [1205] -93.62134 -93.62267 -93.62266 -93.62926 -93.62955 -93.62774 -93.62988
## [1212] -93.62771 -93.62680 -93.62617 -93.62676 -93.62327 -93.62192 -93.62157
## [1219] -93.62208 -93.62820 -93.62936 -93.62758 -93.62824 -93.62851 -93.62264
## [1226] -93.62264 -93.62356 -93.62049 -93.61670 -93.61655 -93.61786 -93.61699
## [1233] -93.61568 -93.61561 -93.61718 -93.61933 -93.61902 -93.61863 -93.61754
## [1240] -93.61754 -93.61322 -93.61486 -93.61308 -93.61239 -93.61338 -93.61724
## [1247] -93.61708 -93.61709 -93.61940 -93.61562 -93.61707 -93.61226 -93.61652
## [1254] -93.61846 -93.61548 -93.61548 -93.61076 -93.61227 -93.61226 -93.61225
## [1261] -93.60742 -93.61066 -93.60519 -93.60566 -93.60717 -93.60411 -93.60822
## [1268] -93.60796 -93.60891 -93.60688 -93.60597 -93.60597 -93.60679 -93.61060
## [1275] -93.60795 -93.60876 -93.60890 -93.60711 -93.60780 -93.60780 -93.60499
## [1282] -93.60484 -93.60370 -93.61957 -93.61941 -93.61558 -93.62005 -93.61648
## [1289] -93.61560 -93.61699 -93.61711 -93.61894 -93.61958 -93.61985 -93.62022
## [1296] -93.61397 -93.61545 -93.61503 -93.61129 -93.60887 -93.60887 -93.60777
## [1303] -93.60676 -93.60771 -93.60772 -93.60680 -93.60771 -93.61844 -93.62018
## [1310] -93.62018 -93.61830 -93.61830 -93.61739 -93.61688 -93.61700 -93.61700
## [1317] -93.61224 -93.61384 -93.60969 -93.60965 -93.61001 -93.60875 -93.60755
## [1324] -93.60597 -93.60620 -93.60449 -93.60435 -93.62881 -93.62671 -93.62573
## [1331] -93.62593 -93.62471 -93.62558 -93.62559 -93.62367 -93.62457 -93.62351
## [1338] -93.62153 -93.62153 -93.62210 -93.62469 -93.62469 -93.62363 -93.62257
## [1345] -93.62255 -93.62149 -93.62039 -93.62039 -93.62459 -93.62447 -93.62189
## [1352] -93.62700 -93.62556 -93.62554 -93.62543 -93.62562 -93.62168 -93.62107
## [1359] -93.62083 -93.62595 -93.62845 -93.62773 -93.62616 -93.62659 -93.62950
## [1366] -93.62529 -93.62374 -93.62192 -93.62482 -93.62461 -93.62458 -93.65816
## [1373] -93.65698 -93.65796 -93.65703 -93.65839 -93.67555 -93.66981 -93.67224
## [1380] -93.67687 -93.67772 -93.67794 -93.67272 -93.67418 -93.67190 -93.66695
## [1387] -93.66691 -93.66552 -93.66947 -93.66373 -93.66536 -93.66536 -93.66698
## [1394] -93.66213 -93.66067 -93.66081 -93.66052 -93.67812 -93.67622 -93.67399
## [1401] -93.67358 -93.66607 -93.66622 -93.66615 -93.66515 -93.66557 -93.66305
## [1408] -93.66312 -93.66060 -93.66086 -93.66753 -93.66507 -93.66526 -93.66446
## [1415] -93.66188 -93.66292 -93.65853 -93.65995 -93.68322 -93.68330 -93.68203
## [1422] -93.68203 -93.68218 -93.69128 -93.69110 -93.68698 -93.68679 -93.68120
## [1429] -93.68114 -93.68265 -93.67920 -93.69201 -93.69121 -93.69161 -93.69072
## [1436] -93.68977 -93.69028 -93.69011 -93.68900 -93.68902 -93.69171 -93.68910
## [1443] -93.68896 -93.68895 -93.68574 -93.68689 -93.68752 -93.68516 -93.68587
## [1450] -93.68440 -93.68429 -93.68316 -93.68331 -93.68786 -93.68770 -93.68619
## [1457] -93.68423 -93.68440 -93.68272 -93.68557 -93.68557 -93.68522 -93.68664
## [1464] -93.68698 -93.68819 -93.68767 -93.68386 -93.67994 -93.68112 -93.68369
## [1471] -93.67951 -93.68031 -93.68783 -93.68693 -93.68521 -93.68445 -93.68440
## [1478] -93.68434 -93.68343 -93.68343 -93.68384 -93.68313 -93.68299 -93.68144
## [1485] -93.68140 -93.67830 -93.67830 -93.67755 -93.67741 -93.67345 -93.67357
## [1492] -93.67314 -93.67144 -93.67241 -93.67529 -93.67517 -93.67325 -93.67558
## [1499] -93.67490 -93.67063 -93.66504 -93.66499 -93.66538 -93.66556 -93.66190
## [1506] -93.66012 -93.66492 -93.66583 -93.66451 -93.66348 -93.66345 -93.66176
## [1513] -93.65580 -93.65668 -93.65673 -93.65585 -93.64856 -93.65181 -93.65179
## [1520] -93.65163 -93.64702 -93.64675 -93.64841 -93.64644 -93.64644 -93.64471
## [1527] -93.64629 -93.64842 -93.64471 -93.64674 -93.64289 -93.64282 -93.64141
## [1534] -93.63962 -93.63950 -93.63947 -93.64108 -93.64260 -93.64338 -93.64293
## [1541] -93.64031 -93.64029 -93.64574 -93.64574 -93.62841 -93.62746 -93.62529
## [1548] -93.62529 -93.62911 -93.62838 -93.62584 -93.62683 -93.62514 -93.62522
## [1555] -93.61543 -93.61685 -93.61685 -93.60420 -93.60419 -93.61624 -93.61658
## [1562] -93.61598 -93.61592 -93.60734 -93.60712 -93.60489 -93.60455 -93.60455
## [1569] -93.60291 -93.60243 -93.65887 -93.63979 -93.65023 -93.65312 -93.65302
## [1576] -93.65172 -93.64920 -93.64947 -93.65043 -93.64468 -93.64605 -93.64706
## [1583] -93.64427 -93.64673 -93.64665 -93.64783 -93.65322 -93.64721 -93.60827
## [1590] -93.60834 -93.61015 -93.60820 -93.60637 -93.60535 -93.60465 -93.60418
## [1597] -93.60360 -93.60427 -93.60144 -93.60318 -93.60346 -93.61011 -93.60511
## [1604] -93.60476 -93.60392 -93.60392 -93.60389 -93.60829 -93.60372 -93.60007
## [1611] -93.57743 -93.61850 -93.61850 -93.61727 -93.61728 -93.63988 -93.63736
## [1618] -93.63615 -93.63697 -93.63719 -93.63434 -93.63635 -93.64061 -93.63787
## [1625] -93.63788 -93.63907 -93.63785 -93.63730 -93.63792 -93.63583 -93.63333
## [1632] -93.63451 -93.63292 -93.63162 -93.63163 -93.62902 -93.62908 -93.62637
## [1639] -93.62287 -93.62954 -93.63007 -93.62875 -93.62872 -93.63787 -93.63885
## [1646] -93.63788 -93.63782 -93.63716 -93.63466 -93.63397 -93.63153 -93.63267
## [1653] -93.63179 -93.63855 -93.63541 -93.63724 -93.63471 -93.63458 -93.63650
## [1660] -93.63390 -93.63376 -93.63219 -93.62774 -93.62633 -93.62656 -93.62912
## [1667] -93.62672 -93.62711 -93.62805 -93.62805 -93.62475 -93.62165 -93.62480
## [1674] -93.62480 -93.62443 -93.62975 -93.62866 -93.62785 -93.62978 -93.62724
## [1681] -93.62724 -93.62457 -93.62567 -93.62580 -93.65888 -93.65801 -93.65768
## [1688] -93.65698 -93.65664 -93.65424 -93.65704 -93.65666 -93.65627 -93.65797
## [1695] -93.65477 -93.65705 -93.65629 -93.65744 -93.65331 -93.65752 -93.65721
## [1702] -93.65724 -93.65542 -93.65469 -93.65462 -93.65092 -93.65749 -93.65625
## [1709] -93.65654 -93.65699 -93.65452 -93.65060 -93.65050 -93.65045 -93.65041
## [1716] -93.64984 -93.64983 -93.64981 -93.64980 -93.64940 -93.64932 -93.64929
## [1723] -93.65341 -93.65256 -93.65205 -93.65141 -93.65061 -93.64694 -93.64643
## [1730] -93.64588 -93.64573 -93.64479 -93.64441 -93.64430 -93.64573 -93.64438
## [1737] -93.64438 -93.64437 -93.64437 -93.64445 -93.64437 -93.64180 -93.64172
## [1744] -93.64051 -93.64218 -93.64426 -93.64406 -93.64306 -93.64396 -93.64395
## [1751] -93.64178 -93.63998 -93.64021 -93.64229 -93.64334 -93.64128 -93.64284
## [1758] -93.64364 -93.64132 -93.64066 -93.65759 -93.65571 -93.65289 -93.65110
## [1765] -93.65264 -93.65278 -93.64990 -93.65727 -93.65516 -93.65591 -93.65371
## [1772] -93.65371 -93.65334 -93.65295 -93.65386 -93.64412 -93.64377 -93.64207
## [1779] -93.64176 -93.64226 -93.64224 -93.64223 -93.64396 -93.64409 -93.64384
## [1786] -93.64248 -93.64234 -93.64354 -93.64345 -93.64220 -93.64125 -93.64375
## [1793] -93.64227 -93.65034 -93.65042 -93.65030 -93.65028 -93.64758 -93.64793
## [1800] -93.64874 -93.64397 -93.64397 -93.64379 -93.64205 -93.64123 -93.63957
## [1807] -93.69179 -93.69224 -93.68912 -93.69124 -93.68903 -93.68977 -93.68880
## [1814] -93.68945 -93.68643 -93.68641 -93.68519 -93.67911 -93.67701 -93.67448
## [1821] -93.67399 -93.67389 -93.66956 -93.66955 -93.67232 -93.67352 -93.67194
## [1828] -93.66185 -93.66185 -93.66052 -93.66198 -93.66213 -93.65308 -93.65266
## [1835] -93.64797 -93.64733 -93.64717 -93.64665 -93.64664 -93.64661 -93.64658
## [1842] -93.64648 -93.64411 -93.64410 -93.64567 -93.64564 -93.64560 -93.64548
## [1849] -93.64181 -93.64218 -93.63993 -93.63988 -93.64346 -93.64820 -93.64580
## [1856] -93.64559 -93.64577 -93.64817 -93.64555 -93.64554 -93.65823 -93.65750
## [1863] -93.63952 -93.63455 -93.63509 -93.63453 -93.63590 -93.63155 -93.63805
## [1870] -93.63687 -93.63375 -93.63167 -93.63130 -93.62992 -93.62919 -93.62652
## [1877] -93.62590 -93.62476 -93.62459 -93.62454 -93.62475 -93.62364 -93.62461
## [1884] -93.62310 -93.62135 -93.62134 -93.62232 -93.62283 -93.62145 -93.62144
## [1891] -93.62695 -93.63002 -93.62799 -93.62686 -93.62764 -93.62496 -93.62448
## [1898] -93.62326 -93.62137 -93.62280 -93.62150 -93.62966 -93.62954 -93.62865
## [1905] -93.62858 -93.62563 -93.62370 -93.62149 -93.62050 -93.62049 -93.61938
## [1912] -93.61706 -93.61569 -93.61610 -93.61814 -93.61553 -93.61879 -93.61840
## [1919] -93.61935 -93.61752 -93.61528 -93.61687 -93.61305 -93.61170 -93.61235
## [1926] -93.61346 -93.61398 -93.61087 -93.61230 -93.61927 -93.61724 -93.61848
## [1933] -93.61562 -93.61843 -93.61707 -93.61241 -93.61240 -93.61245 -93.61207
## [1940] -93.61225 -93.61861 -93.61975 -93.62034 -93.61719 -93.61844 -93.61562
## [1947] -93.61549 -93.61330 -93.61554 -93.61242 -93.61390 -93.61390 -93.60712
## [1954] -93.60915 -93.61056 -93.60539 -93.61061 -93.61061 -93.61062 -93.60796
## [1961] -93.60796 -93.60796 -93.60693 -93.60597 -93.60674 -93.60678 -93.60594
## [1968] -93.60675 -93.60675 -93.60484 -93.60370 -93.60370 -93.60456 -93.61897
## [1975] -93.61699 -93.61905 -93.62024 -93.61714 -93.61767 -93.61543 -93.61545
## [1982] -93.60861 -93.60887 -93.60686 -93.60684 -93.60772 -93.61047 -93.60826
## [1989] -93.60679 -93.61866 -93.61908 -93.61706 -93.61829 -93.61691 -93.61885
## [1996] -93.61704 -93.61558 -93.61206 -93.61219 -93.60888 -93.60945 -93.60868
## [2003] -93.60476 -93.60418 -93.60512 -93.62819 -93.62627 -93.62761 -93.62749
## [2010] -93.62471 -93.62352 -93.62150 -93.62364 -93.62363 -93.62454 -93.62255
## [2017] -93.62150 -93.62150 -93.62240 -93.62240 -93.62121 -93.62091 -93.62576
## [2024] -93.62543 -93.62529 -93.62528 -93.62416 -93.62033 -93.62215 -93.62200
## [2031] -93.62845 -93.62830 -93.62950 -93.62776 -93.62751 -93.62950 -93.62749
## [2038] -93.62826 -93.62638 -93.62530 -93.62279 -93.62202 -93.62188 -93.62381
## [2045] -93.65570 -93.65612 -93.65709 -93.65592 -93.67555 -93.67556 -93.67221
## [2052] -93.66974 -93.67442 -93.67369 -93.67750 -93.67695 -93.67778 -93.67815
## [2059] -93.67571 -93.67800 -93.67760 -93.67668 -93.67212 -93.66854 -93.66404
## [2066] -93.66224 -93.66067 -93.66069 -93.66198 -93.66199 -93.65936 -93.67627
## [2073] -93.67396 -93.67458 -93.66961 -93.67122 -93.67151 -93.66504 -93.66342
## [2080] -93.66712 -93.66061 -93.66352 -93.66443 -93.66468 -93.66431 -93.65855
## [2087] -93.66057 -93.65966 -93.65854 -93.69315 -93.68410 -93.68410 -93.68347
## [2094] -93.68331 -93.68323 -93.68106 -93.68683 -93.68694 -93.69034 -93.69067
## [2101] -93.69125 -93.69121 -93.69101 -93.69023 -93.68907 -93.68907 -93.69003
## [2108] -93.69135 -93.68894 -93.68892 -93.68988 -93.68892 -93.68893 -93.68514
## [2115] -93.68494 -93.68244 -93.67922 -93.68251 -93.68361 -93.69197 -93.69207
## [2122] -93.68886 -93.69040 -93.69055 -93.69229 -93.69270 -93.69241 -93.69063
## [2129] -93.68895 -93.68892 -93.68934 -93.68952 -93.69147 -93.69205 -93.68496
## [2136] -93.68668 -93.68743 -93.68773 -93.68484 -93.68499 -93.68323 -93.68435
## [2143] -93.68316 -93.68317 -93.68786 -93.68658 -93.68281 -93.68680 -93.68779
## [2150] -93.68675 -93.68381 -93.67958 -93.68120 -93.68008 -93.67981 -93.68127
## [2157] -93.68108 -93.67951 -93.68768 -93.68690 -93.68412 -93.68414 -93.68471
## [2164] -93.68462 -93.68429 -93.68427 -93.68373 -93.67835 -93.67754 -93.67622
## [2171] -93.67145 -93.67145 -93.67145 -93.67320 -93.67470 -93.67525 -93.67433
## [2178] -93.67431 -93.67325 -93.67318 -93.67622 -93.67624 -93.67059 -93.66574
## [2185] -93.66477 -93.66351 -93.66469 -93.66154 -93.66030 -93.66587 -93.66482
## [2192] -93.66346 -93.66345 -93.66333 -93.65016 -93.65181 -93.64662 -93.64669
## [2199] -93.64682 -93.64672 -93.64617 -93.64454 -93.64467 -93.64471 -93.64451
## [2206] -93.64470 -93.64476 -93.64629 -93.64642 -93.64353 -93.64321 -93.64281
## [2213] -93.64195 -93.64443 -93.64132 -93.64236 -93.64285 -93.64050 -93.64116
## [2220] -93.64266 -93.64288 -93.64401 -93.63939 -93.64033 -93.64185 -93.64574
## [2227] -93.64492 -93.64579 -93.64355 -93.64343 -93.64018 -93.62841 -93.62950
## [2234] -93.59657 -93.59700 -93.62687 -93.62826 -93.62531 -93.62950 -93.62824
## [2241] -93.62669 -93.62822 -93.61697 -93.61506 -93.61389 -93.61570 -93.61562
## [2248] -93.60658 -93.60667 -93.60760 -93.60755 -93.60693 -93.60250 -93.60104
## [2255] -93.60178 -93.63966 -93.64740 -93.64976 -93.64945 -93.65059 -93.65057
## [2262] -93.65056 -93.65173 -93.65042 -93.64662 -93.64437 -93.64678 -93.64679
## [2269] -93.64679 -93.64689 -93.64663 -93.64759 -93.64720 -93.66216 -93.65311
## [2276] -93.65243 -93.64757 -93.61010 -93.61012 -93.60819 -93.60551 -93.60728
## [2283] -93.60739 -93.60536 -93.60535 -93.60355 -93.60186 -93.60183 -93.60151
## [2290] -93.60209 -93.60342 -93.60492 -93.60464 -93.60452 -93.60678 -93.60452
## [2297] -93.60256 -93.60152 -93.60711 -93.60486 -93.60512 -93.60415 -93.60308
## [2304] -93.60664 -93.60486 -93.61846 -93.61818 -93.61926 -93.61644 -93.61743
## [2311] -93.63907 -93.63787 -93.63748 -93.63757 -93.63638 -93.63637 -93.63605
## [2318] -93.63597 -93.63594 -93.63584 -93.63580 -93.63593 -93.63160 -93.63572
## [2325] -93.63695 -93.63929 -93.63423 -93.63157 -93.63257 -93.62960 -93.62934
## [2332] -93.62917 -93.62861 -93.62877 -93.62777 -93.62634 -93.62648 -93.62201
## [2339] -93.62517 -93.63007 -93.62406 -93.62827 -93.63803 -93.63787 -93.63266
## [2346] -93.63151 -93.63300 -93.63549 -93.63533 -93.63604 -93.63739 -93.63719
## [2353] -93.63537 -93.63952 -93.63237 -93.63035 -93.63236 -93.63029 -93.63116
## [2360] -93.62626 -93.62537 -93.62798 -93.62785 -93.62441 -93.62991 -93.62923
## [2367] -93.62793 -93.62710 -93.62733 -93.62737 -93.62988 -93.62985 -93.62952
## [2374] -93.62717 -93.62886 -93.62895 -93.62477 -93.62571 -93.65886 -93.65841
## [2381] -93.65823 -93.65799 -93.65665 -93.65477 -93.65798 -93.65879 -93.65712
## [2388] -93.65715 -93.65686 -93.65302 -93.65464 -93.65233 -93.65483 -93.65345
## [2395] -93.65315 -93.64341 -93.64386 -93.65500 -93.65546 -93.65788 -93.65308
## [2402] -93.65063 -93.65341 -93.64971 -93.65414 -93.64955 -93.64952 -93.65412
## [2409] -93.65414 -93.65062 -93.64646 -93.64584 -93.64577 -93.64441 -93.64437
## [2416] -93.64463 -93.64575 -93.64421 -93.64431 -93.64166 -93.64195 -93.64029
## [2423] -93.64023 -93.64132 -93.64190 -93.64226 -93.64224 -93.64212 -93.64152
## [2430] -93.64152 -93.64458 -93.64259 -93.64360 -93.64409 -93.64147 -93.64021
## [2437] -93.64271 -93.64291 -93.64128 -93.64136 -93.64128 -93.64105 -93.65574
## [2444] -93.65452 -93.65444 -93.65785 -93.65320 -93.65227 -93.65006 -93.65391
## [2451] -93.65600 -93.65570 -93.65361 -93.65348 -93.64232 -93.64208 -93.64178
## [2458] -93.64179 -93.64364 -93.64381 -93.65036 -93.64407 -93.64390 -93.64286
## [2465] -93.64231 -93.64237 -93.64357 -93.64317 -93.64386 -93.64398 -93.69258
## [2472] -93.69140 -93.69179 -93.69175 -93.69096 -93.69125 -93.69125 -93.68902
## [2479] -93.68782 -93.68770 -93.68627 -93.68627 -93.68665 -93.68160 -93.68141
## [2486] -93.67417 -93.67596 -93.67712 -93.67572 -93.67701 -93.66955 -93.67104
## [2493] -93.67241 -93.67092 -93.67100 -93.66971 -93.66976 -93.67193 -93.66064
## [2500] -93.65388 -93.65351 -93.65235 -93.65047 -93.65006 -93.64783 -93.64653
## [2507] -93.64489 -93.64489 -93.64489 -93.64489 -93.64562 -93.64548 -93.64548
## [2514] -93.64548 -93.64337 -93.64114 -93.64893 -93.64606 -93.64933 -93.64944
## [2521] -93.64935 -93.64933 -93.65824 -93.63375 -93.63460 -93.63474 -93.63587
## [2528] -93.63347 -93.63026 -93.63247 -93.63310 -93.62793 -93.62598 -93.62610
## [2535] -93.62616 -93.62752 -93.62807 -93.62690 -93.62635 -93.62461 -93.62915
## [2542] -93.62341 -93.62356 -93.62231 -93.63016 -93.62677 -93.62707 -93.62829
## [2549] -93.62162 -93.62136 -93.62202 -93.62357 -93.62265 -93.62140 -93.62857
## [2556] -93.62919 -93.62926 -93.62818 -93.62703 -93.62474 -93.62370 -93.62265
## [2563] -93.62049 -93.62143 -93.62143 -93.61967 -93.61685 -93.61601 -93.61563
## [2570] -93.61812 -93.62036 -93.61839 -93.61933 -93.61750 -93.61862 -93.61289
## [2577] -93.61257 -93.61588 -93.61330 -93.61425 -93.61246 -93.61832 -93.61848
## [2584] -93.61859 -93.61723 -93.61707 -93.61246 -93.61349 -93.61368 -93.61547
## [2591] -93.61547 -93.61226 -93.62035 -93.61718 -93.61687 -93.61435 -93.61337
## [2598] -93.61396 -93.61077 -93.61224 -93.61390 -93.61164 -93.61227 -93.61405
## [2605] -93.61240 -93.61390 -93.60891 -93.60898 -93.60758 -93.60942 -93.61028
## [2612] -93.61050 -93.60726 -93.60842 -93.60907 -93.60651 -93.60538 -93.60674
## [2619] -93.60573 -93.60905 -93.61050 -93.60890 -93.60890 -93.60674 -93.60905
## [2626] -93.60905 -93.60905 -93.60689 -93.60689 -93.60674 -93.60370 -93.60584
## [2633] -93.60469 -93.60457 -93.60507 -93.61854 -93.61955 -93.61716 -93.61852
## [2640] -93.61714 -93.61849 -93.61951 -93.61568 -93.61880 -93.61834 -93.61833
## [2647] -93.61564 -93.61401 -93.61459 -93.61462 -93.61391 -93.60786 -93.60883
## [2654] -93.61035 -93.60894 -93.60774 -93.60680 -93.61829 -93.61689 -93.61389
## [2661] -93.61536 -93.61536 -93.61534 -93.61280 -93.61376 -93.60890 -93.61044
## [2668] -93.60967 -93.60619 -93.60687 -93.62848 -93.62602 -93.62739 -93.62558
## [2675] -93.62043 -93.62364 -93.62454 -93.62455 -93.62455 -93.62348 -93.62148
## [2682] -93.62242 -93.62041 -93.62438 -93.62239 -93.62577 -93.62277 -93.62544
## [2689] -93.62156 -93.62262 -93.62686 -93.62663 -93.62757 -93.62732 -93.62551
## [2696] -93.62842 -93.62950 -93.62659 -93.65702 -93.65565 -93.65815 -93.65783
## [2703] -93.66004 -93.65588 -93.65579 -93.67555 -93.67121 -93.67540 -93.67446
## [2710] -93.67452 -93.67541 -93.67561 -93.67646 -93.67228 -93.67439 -93.67360
## [2717] -93.67265 -93.67020 -93.67051 -93.66557 -93.66550 -93.66563 -93.66708
## [2724] -93.66954 -93.66412 -93.66839 -93.66708 -93.66204 -93.67817 -93.67635
## [2731] -93.67440 -93.67668 -93.66603 -93.66515 -93.66276 -93.66348 -93.66156
## [2738] -93.66066 -93.66023 -93.66626 -93.66719 -93.66609 -93.66610 -93.66019
## [2745] -93.66274 -93.65851 -93.68554 -93.68323 -93.68218 -93.68015 -93.68907
## [2752] -93.68982 -93.69037 -93.69123 -93.68894 -93.68892 -93.68892 -93.68949
## [2759] -93.68907 -93.68908 -93.68879 -93.68893 -93.68893 -93.68216 -93.68135
## [2766] -93.68209 -93.68388 -93.68040 -93.69190 -93.69187 -93.68830 -93.69036
## [2773] -93.69200 -93.69261 -93.69222 -93.69021 -93.69193 -93.68908 -93.68923
## [2780] -93.69124 -93.68858 -93.68663 -93.68682 -93.68716 -93.68498 -93.68347
## [2787] -93.68325 -93.68298 -93.68138 -93.68770 -93.68629 -93.68442 -93.68550
## [2794] -93.68570 -93.68767 -93.68557 -93.68372 -93.67915 -93.68148 -93.67980
## [2801] -93.68269 -93.68189 -93.68119 -93.68373 -93.68378 -93.68768 -93.68769
## [2808] -93.68783 -93.68409 -93.68412 -93.68412 -93.68537 -93.68587 -93.68343
## [2815] -93.68433 -93.68376 -93.67827 -93.67637 -93.67636 -93.67530 -93.67555
## [2822] -93.67145 -93.67238 -93.67130 -93.67161 -93.67002 -93.66998 -93.66993
## [2829] -93.67006 -93.67004 -93.66998 -93.66996 -93.66992 -93.66195 -93.65918
## [2836] -93.65922 -93.65922 -93.66012 -93.66346 -93.66200 -93.65675 -93.65670
## [2843] -93.65564 -93.65904 -93.65563 -93.64856 -93.65181 -93.65129 -93.64694
## [2850] -93.64662 -93.64644 -93.64842 -93.64842 -93.64644 -93.64471 -93.64841
## [2857] -93.64703 -93.64431 -93.64426 -93.64243 -93.64048 -93.64344 -93.64115
## [2864] -93.64232 -93.64100 -93.64167 -93.64061 -93.64288 -93.64420 -93.63962
## [2871] -93.63943 -93.64349 -93.62847 -93.62841 -93.62830 -93.62950 -93.62668
## [2878] -93.62515 -93.61531 -93.61385 -93.61501 -93.61502 -93.60434 -93.61603
## [2885] -93.61624 -93.59997 -93.60674 -93.60367 -93.60317 -93.60119 -93.60385
## [2892] -93.64453 -93.64559 -93.65108 -93.64977 -93.64610 -93.64657 -93.64617
## [2899] -93.64664 -93.64676 -93.64789 -93.64898 -93.64747 -93.61027 -93.60834
## [2906] -93.60819 -93.60819 -93.60714 -93.60439 -93.60392 -93.60352 -93.60428
## [2913] -93.60109 -93.60318 -93.60352 -93.60044 -93.60204 -93.60190 -93.60184
## [2920] -93.60161 -93.60234 -93.60447 -93.60353 -93.60869 -93.60684 -93.60478
## [2927] -93.60268 -93.60685 -93.60019 -93.60000
ames1 <- ames %>%
  filter(Neighborhood == "North_Ames")


ames.sf <- st_as_sf(ames1, coords = c("Longitude", "Latitude"), crs = 4326) 


library(leaflet)

leaflet() %>%
  addProviderTiles(providers$CartoDB.Positron)%>%
  addMarkers(data = ames.sf)

3 Operaciones

3.1 Cantidad de Viviendas por Barrio

cantidadviviendas <-ames%>%
  group_by(barrio = Neighborhood)%>%
  summarise(cantidad = n() )

ggplot(data=cantidadviviendas)+
  geom_bar(aes(x=barrio,weight=cantidad))

ggplot(data=cantidadviviendas)+
  geom_bar(aes(x=barrio,weight=cantidad))+
  coord_flip()

Ejercicios:

  • ¿Cuantos calles están pavimentadas y cuantas no pavimentadas?

3.2 seleccionamos viviendas de un barrio

barrioGilbert <- ames%>%
  filter(ames$Neighborhood == "Greens")

barrioGilbert <- st_as_sf(barrioGilbert, coords = c("Longitude", "Latitude"), crs = 4326)

library(leaflet)
leaflet() %>%
  addProviderTiles(providers$CartoDB.Positron)%>%
  addMarkers(data = barrioGilbert)

3.3 Transformación de unidad de medidas

Crear dos nuevas columnas a partir de las columnas Lot_Frontage y Lot_Area.

metros_cuadrados <- 0.3048

ames$FrenteLote <- ames$Lot_Frontage * metros_cuadrados
ames$AreaLote <- ames$Lot_Area * metros_cuadrados * metros_cuadrados

hist(ames$AreaLote)

boxplot(ames$AreaLote)

4 Arbol de Decisión

Los árboles de decisión son un potente algoritmo de aprendizaje automático que se puede utilizar tanto para tareas de clasificación como de regresión. Son fáciles de entender e interpretar y se pueden utilizar para construir modelos complejos sin necesidad de ingeniería de características.

Una vez que haya entrenado un modelo de árbol de decisión, podrá utilizarlo para hacer predicciones sobre nuevos datos. Sin embargo, también puede resultar útil trazar el árbol de decisiones para comprender mejor cómo funciona e identificar posibles problemas.

Te mostrare cómo trazar árboles de decisión en R usando los paquetes rpart y rpart.plot. Empezaremos con un ejemplo utilizando el conjunto de datos del iris y explicaremos los bloques de código en términos simples de usar.

4.1 con Iris

#Load the libraries
library(rpart)
## Warning: package 'rpart' was built under R version 4.3.3
library(rpart.plot)
## Warning: package 'rpart.plot' was built under R version 4.3.3
#Split the data into training and test sets
set.seed(123)
train_index <- sample(1:nrow(iris), size = 0.7 * nrow(iris))
train <- iris[train_index, ]
test <- iris[-train_index, ]

#Train a decision tree model
tree <- rpart(Species ~ ., data = train, method = "class")

#Plot the decision tree
rpart.plot(tree, main = "Decision Tree for the Iris Dataset")

4.1.1 Output

La salida de la función rpart.plot() es un diagrama de árbol que muestra las reglas de decisión del modelo. El nodo raíz del árbol está en la parte superior y los nodos de las hojas están en la parte inferior. Cada nodo está etiquetado con la característica que se utiliza para dividir los datos en ese nodo y el valor de la división. Los nodos hoja están etiquetados con la clase prevista para los datos que llegan a ese nodo.

4.1.2 Interpretando el árbol de decisión

Para interpretar el árbol de decisión, comience en el nodo raíz y siga las ramas hasta un nodo hoja. El nodo hoja al que llega es la clase prevista para los datos con los que comenzó.

Por ejemplo, si tiene una nueva flor de iris con una longitud de sépalo de 5,5 cm y una longitud de pétalo de 2,5 cm, comenzaría en el nodo raíz del árbol de decisión. En el nodo raíz, la característica que se utiliza para dividir los datos es la longitud del pétalo. Dado que la longitud de los pétalos de la nueva flor es superior a 2,45 cm, deberás seguir la rama derecha hasta el siguiente nodo. En el siguiente nodo, la característica que se utiliza para dividir los datos es la longitud del sépalo. Dado que la longitud del sépalo de la nueva flor es superior a 5,0 cm, deberá seguir la rama derecha hasta el nudo de la hoja. El nodo de la hoja al que llega está etiquetado como “versicolor”, por lo que la clase prevista para la nueva flor es versicolor.

4.1.3 Algoritmo de alto nivel

  • Para comenzar, cargue las bibliotecas rpart y rpart.plot y cargue su conjunto de datos.

  • Luego, divida los datos en conjuntos de entrenamiento y prueba. + Entrene un modelo de árbol de decisión utilizando la función rpart().

  • Finalmente, trace el árbol de decisión usando la función rpart.plot().

  • Una vez que haya trazado el árbol de decisión, tómese un tiempo para interpretarlo.

  • Intente comprender cómo el modelo hace predicciones e identifique cualquier problema potencial.

  • También puede intentar mejorar el modelo agregando o eliminando funciones o cambiando los hiperparámetros de la función rpart().

4.2 Arbol de decision para Ames

El precio de la vivienda puede depender de Bldg_Type y de Gr_Liv_Area

4.2.1 Gr_Liv_Area

El campo Gr_Liv_Area en el conjunto de datos de AmesHousing del paquete modeldata representa el área habitable sobre el nivel del suelo (en pies cuadrados) de cada propiedad. Es decir, incluye el área de todas las habitaciones habitables (como salas de estar, dormitorios, comedores, etc.) que están sobre el nivel del suelo, excluyendo sótanos y áreas no habitables como garajes o porches.

Este atributo es importante en el análisis de precios de viviendas, ya que el área habitable es uno de los factores clave que influyen en el precio de una propiedad. Las propiedades con áreas habitables más grandes tienden a tener precios más altos en comparación con propiedades con áreas habitables más pequeñas, siempre y cuando otros factores se mantengan constantes.

4.2.2 Bldg_Type

Cada valor de Bldg_Type indica un tipo específico de estructura de vivienda. Aquí está el significado general de cada tipo de Bldg_Type:

  • 1Fam: Casa unifamiliar.
  • 2FmCon: Casa con dos unidades residenciales (dúplex).
  • Duplx: Casa con dos unidades residenciales (dúplex).
  • TwnhsE: Casa adosada de estilo europeo (townhouse end unit).
  • TwnhsI: Casa adosada de estilo interior (townhouse inside unit).

4.2.3 Ploteo

# Función para convertir pies cuadrados a metros cuadrados
convertir_pies_cuadrados_a_metros_cuadrados <- function(pies_cuadrados) {
  conversion_factor <- 0.092903
  metros_cuadrados <- pies_cuadrados * conversion_factor
  return(metros_cuadrados)
}

# Crear una nueva columna en el dataframe 'ames' con 'Gr_Liv_Area' en metros cuadrados
ames$Gr_Liv_Area_m2 <- convertir_pies_cuadrados_a_metros_cuadrados(ames$Gr_Liv_Area)


ggplot(data = ames, aes(x = Gr_Liv_Area_m2, y = Sale_Price, color = Bldg_Type)) +
  geom_point() +
  labs(x = "Above Grade Living Area (sqft)", y = "Sale Price", 
       title = "Sale Price vs. Above Grade Living Area by Building Type")

### con facetas

# Create the plot with facets
ggplot(data = ames, aes(x = Gr_Liv_Area_m2, y = Sale_Price)) +
  geom_point() +
  facet_wrap(~ Bldg_Type, ncol = 2) +
  labs(x = "Above Grade Living Area (sqft)", y = "Sale Price", 
       title = "Sale Price vs. Above Grade Living Area by Building Type")

otro grafico con colores

colors <- c("#FF5733", "#33FF57", "#3357FF", "#FF33ED", "#FF5733")

# Crea el gráfico con facetas y colores personalizados
ggplot(data = ames, aes(x = Gr_Liv_Area_m2, y = Sale_Price)) +
  geom_point(aes(color = Bldg_Type), size = 3) +
  facet_wrap(~ Bldg_Type, ncol = 2) +
  labs(x = "Above Grade Living Area (sqft)", y = "Sale Price", 
       title = "Sale Price vs. Above Grade Living Area by Building Type") +
  scale_color_manual(values = colors) +
  theme(legend.position = "none")  # Opcional: quitar la leyenda

### Dolar format

# Crea el gráfico con facetas y colores personalizados
ggplot(data = ames, aes(x = Gr_Liv_Area, y = Sale_Price)) +
  geom_point(aes(color = Bldg_Type), size = 3) +
  facet_wrap(~ Bldg_Type, ncol = 2) +
  labs(x = "Above Grade Living Area (sqft)", y = "Sale Price (USD)", 
       title = "Sale Price vs. Above Grade Living Area by Building Type") +
  scale_color_manual(values = colors) +
  theme(legend.position = "none") +  # Opcional: quitar la leyenda
  scale_y_continuous(labels = scales::dollar_format())

4.2.4 Arbol

# https://cs.calvin.edu/courses/data/202/21fa/ex/ex07/ex07-modeling-inst.html
set.seed(1234)

# Split our data randomly
library(rsample)
## Warning: package 'rsample' was built under R version 4.3.3
ames_split <- initial_split(ames, prop = 2/3) 
ames_train <- training(ames_split)
ames_test <- testing(ames_split)

#install.packages("rminer")
library(rminer)
## Warning: package 'rminer' was built under R version 4.3.3
# Instala y carga el paquete rpart
#install.packages("rpart")
library(rpart)
library(rpart.plot)
# Ajusta el modelo de árbol de decisión
decision_tree_fit <- rpart(Sale_Price ~ Gr_Liv_Area_m2 + Bldg_Type, data = ames_train, method = "anova", control = rpart.control(maxdepth = 3))


# Grafica el árbol de decisión
rpart.plot(decision_tree_fit, roundint = FALSE)