Introducción

Como tercera ciudad del país, y una oferta urbana de más del 70%, el análisis habitacional de Cali, capital del Valle del Cauca, brinda una visual de alta relevancia para constructures, comerciales, y agentes inmobiliarios que buscan un entendimiento del mercado de vivienda y sus comportamientos.

Objetivos

Con base en los datos de ofertas de vivienda descargadas del portal Fincaraiz para apartamento de estrato 4 con área construida menor a 200 m2, la inmobiliaria A&C requiere el apoyo de un científico de datos en la construcción de un modelo que lo oriente sobre los precios de inmuebles

Punto 1. Análisis exploratorio de datos

Base de Datos

El insumo obtenido como fuente de información para el presente análisis cuenta con una base de datos de 1.706 registros, representados en 5 distintas columnas (variables), 3 de las cuales son de tipo cuantitativo, mientras las restantes tienen aspecto de calidad. Los datos que soportan los cálculos se pueden encontrar en la siguiente ubicación

library(paqueteMETODOS)
data(vivienda4)

De las variables que componen la base de datos, estas son sus principales características:

viviendas = vivienda4
print(viviendas)
## # A tibble: 1,706 × 5
##    zona       estrato preciom areaconst tipo       
##    <fct>      <fct>     <dbl>     <dbl> <fct>      
##  1 Zona Norte 4          232.        52 Apartamento
##  2 Zona Norte 4          272.       160 Casa       
##  3 Zona Norte 4          255.       108 Apartamento
##  4 Zona Sur   4          258.        96 Apartamento
##  5 Zona Norte 4          250.        82 Apartamento
##  6 Zona Norte 4          261.       117 Casa       
##  7 Zona Norte 4          247.        75 Apartamento
##  8 Zona Norte 4          222.        60 Apartamento
##  9 Zona Norte 4          227.        84 Apartamento
## 10 Zona Norte 4          255.       117 Apartamento
## # ℹ 1,696 more rows
summary(viviendas)
##            zona      estrato     preciom        areaconst     
##  Zona Centro :   8   3:   0   Min.   :207.4   Min.   : 40.00  
##  Zona Norte  : 288   4:1706   1st Qu.:230.7   1st Qu.: 60.00  
##  Zona Oeste  :  60   5:   0   Median :238.8   Median : 75.00  
##  Zona Oriente:   6   6:   0   Mean   :243.7   Mean   : 87.63  
##  Zona Sur    :1344            3rd Qu.:251.5   3rd Qu.: 98.00  
##                               Max.   :309.7   Max.   :200.00  
##           tipo     
##  Apartamento:1363  
##  Casa       : 343  
##                    
##                    
##                    
## 
  • zona - Clasificación en 5 zonas de la ciudad (Centro/Sur/Oeste/Este/Norte)

  • estrato - Estrato socioeconómico del sector de la vivienda, en un rango de 3 a 6

  • preciom - Precio de la viviendia, expresado en millones de pesos, con un rango entre los $207.4 millones y los $309.7 millones

  • areaconst - Área construida en metros cuadrados, con valores que van desde 40 a 200 metros cuadrados

  • tipo - Clasificación del tipo de vivienda, seleccionando entre casa y/o apartamento

table(is.na(viviendas))
## 
## FALSE 
##  8530

En análisis de datos faltantes nos muestra que no existe data que deba ser reemplazada o eliminada. Por tal motivo, no se requiere una depuración de la data en la base de trabajo. Información detallada de clasificación por variables se encuentra en pestaña de análisis

Datos para Análisis

print(viviendas)
## # A tibble: 1,706 × 5
##    zona       estrato preciom areaconst tipo       
##    <fct>      <fct>     <dbl>     <dbl> <fct>      
##  1 Zona Norte 4          232.        52 Apartamento
##  2 Zona Norte 4          272.       160 Casa       
##  3 Zona Norte 4          255.       108 Apartamento
##  4 Zona Sur   4          258.        96 Apartamento
##  5 Zona Norte 4          250.        82 Apartamento
##  6 Zona Norte 4          261.       117 Casa       
##  7 Zona Norte 4          247.        75 Apartamento
##  8 Zona Norte 4          222.        60 Apartamento
##  9 Zona Norte 4          227.        84 Apartamento
## 10 Zona Norte 4          255.       117 Apartamento
## # ℹ 1,696 more rows
summary(viviendas)
##            zona      estrato     preciom        areaconst     
##  Zona Centro :   8   3:   0   Min.   :207.4   Min.   : 40.00  
##  Zona Norte  : 288   4:1706   1st Qu.:230.7   1st Qu.: 60.00  
##  Zona Oeste  :  60   5:   0   Median :238.8   Median : 75.00  
##  Zona Oriente:   6   6:   0   Mean   :243.7   Mean   : 87.63  
##  Zona Sur    :1344            3rd Qu.:251.5   3rd Qu.: 98.00  
##                               Max.   :309.7   Max.   :200.00  
##           tipo     
##  Apartamento:1363  
##  Casa       : 343  
##                    
##                    
##                    
## 

Resultados

Partiendo de la base procesada en la fase previa del presente documento, se dará paso a la construcción de algunos cálculos estadísticos adicionales, que permitan inferir algunos resultados y conclusiones de los datos obtenidos en el análisis

Distribución de Frecuencia

table(is.na(viviendas))
## 
## FALSE 
##  8530
table(viviendas$tipo)
## 
## Apartamento        Casa 
##        1363         343
table(viviendas$estrato)
## 
##    3    4    5    6 
##    0 1706    0    0
table(viviendas$zona)
## 
##  Zona Centro   Zona Norte   Zona Oeste Zona Oriente     Zona Sur 
##            8          288           60            6         1344

Las tres tablas previas nos permiten ir definiendo algunos parámetros de análisis relevantes para la definición de nicho de mercado:

  1. Del total de 8530 registros, encontramos que ninguno tiene algún dato faltante, lo cual lo habíamos concretado en el procesamiento de datos

  2. La mayor cantidad inmuebles ofertados en la ciudad de Cali, 1363 son apartamentos y los 343 restantes son casas, lo cual significa que el 79.89% de la oferta de vivienda está en torres y/o edificios

  3. Para el alcance de la presente base de datos, el 100% de los inmuebles analizados se ubican en estrato 4

  4. El 78.78% de los inmuebles analizados están dentro la zona sur de la ciudad, seguido de la zona norte con el 16.88% y la oeste con el 3.52%


pr.prom_zona = viviendas %>% 
  group_by(zona) %>% 
  summarise(precio_prom = mean(preciom)) #Precio promedio por zona
print(pr.prom_zona)
## # A tibble: 5 × 2
##   zona         precio_prom
##   <fct>              <dbl>
## 1 Zona Centro         257.
## 2 Zona Norte          244.
## 3 Zona Oeste          241.
## 4 Zona Oriente        265.
## 5 Zona Sur            244.
area.prom_zona = viviendas %>% 
  group_by(zona) %>% 
  summarise(area_prom = mean(areaconst)) #Area promedio por zona
print(area.prom_zona)
## # A tibble: 5 × 2
##   zona         area_prom
##   <fct>            <dbl>
## 1 Zona Centro      112. 
## 2 Zona Norte        87.8
## 3 Zona Oeste        81.1
## 4 Zona Oriente     125. 
## 5 Zona Sur          87.6
library(dplyr)
zona_prom = cbind(pr.prom_zona, area.prom_zona[,2])

ind.prom_zona = viviendas %>% 
  group_by(zona) %>% 
  summarise(ind_prom = mean(preciom)/mean(areaconst)) #Area promedio por zona
print(ind.prom_zona)
## # A tibble: 5 × 2
##   zona         ind_prom
##   <fct>           <dbl>
## 1 Zona Centro      2.29
## 2 Zona Norte       2.78
## 3 Zona Oeste       2.98
## 4 Zona Oriente     2.12
## 5 Zona Sur         2.78
zona_prom = cbind(pr.prom_zona, area.prom_zona[,2], ind.prom_zona[,2])
print(zona_prom)
##           zona precio_prom area_prom ind_prom
## 1  Zona Centro    257.2586 112.12500 2.294391
## 2   Zona Norte    243.6365  87.78819 2.775276
## 3   Zona Oeste    241.4031  81.11667 2.975999
## 4 Zona Oriente    265.2154 125.16667 2.118898
## 5     Zona Sur    243.6433  87.57292 2.782176

La previa construcción de tablas cruzadas entre variables de valor nos permite definir datos de mucha relevancia como: precio promedio por zona, area promedio por zona, y factor de precio por metro cuadrado por zona. De esta manera, podemos concluir que es la zona oeste la que mayor precio por metro cuadrado tiene, con cerca de $2.98 millones, seguida de la zona sur con $2.78 millones y la zona norte con $2.78 millones. Los siguientes dos histogramas apoyan parte de lo concluido en esta sección

hist(viviendas$preciom, main="Histograma de Precios", xlab="Precio (millones)", ylab="Frecuencia",
     axes=TRUE, plot=TRUE, labels=TRUE, col="red", border="yellow", breaks=10, ylim = c(0,700),
     xlim=c(min(viviendas$preciom),max(viviendas$preciom)))

hist(viviendas$areaconst, main="Histograma de Areas", xlab="Area (mts2)", ylab="Frecuencia",
     axes=TRUE, plot=TRUE, labels=TRUE, col="blue", border="green", breaks=10, ylim = c(0,700),
     xlim=c(min(viviendas$areaconst),max(viviendas$areaconst)))

Análisis Descriptivo

summarytools::freq(viviendas)
## Frequencies  
## viviendas$zona  
## Type: Factor  
## 
##                      Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
## ------------------ ------ --------- -------------- --------- --------------
##        Zona Centro      8      0.47           0.47      0.47           0.47
##         Zona Norte    288     16.88          17.35     16.88          17.35
##         Zona Oeste     60      3.52          20.87      3.52          20.87
##       Zona Oriente      6      0.35          21.22      0.35          21.22
##           Zona Sur   1344     78.78         100.00     78.78         100.00
##               <NA>      0                               0.00         100.00
##              Total   1706    100.00         100.00    100.00         100.00
## 
## viviendas$estrato  
## Type: Factor  
## 
##               Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
## ----------- ------ --------- -------------- --------- --------------
##           3      0      0.00           0.00      0.00           0.00
##           4   1706    100.00         100.00    100.00         100.00
##           5      0      0.00         100.00      0.00         100.00
##           6      0      0.00         100.00      0.00         100.00
##        <NA>      0                               0.00         100.00
##       Total   1706    100.00         100.00    100.00         100.00
## 
## viviendas$tipo  
## Type: Factor  
## 
##                     Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
## ----------------- ------ --------- -------------- --------- --------------
##       Apartamento   1363     79.89          79.89     79.89          79.89
##              Casa    343     20.11         100.00     20.11         100.00
##              <NA>      0                               0.00         100.00
##             Total   1706    100.00         100.00    100.00         100.00
summarytools::descr(viviendas)
## Descriptive Statistics  
## viviendas  
## N: 1706  
## 
##                     areaconst   preciom
## ----------------- ----------- ---------
##              Mean       87.63    243.70
##           Std.Dev       36.35     19.56
##               Min       40.00    207.41
##                Q1       60.00    230.73
##            Median       75.00    238.77
##                Q3       98.00    251.51
##               Max      200.00    309.70
##               MAD       22.24     14.19
##               IQR       38.00     20.77
##                CV        0.41      0.08
##          Skewness        1.53      1.26
##       SE.Skewness        0.06      0.06
##          Kurtosis        1.68      1.25
##           N.Valid     1706.00   1706.00
##         Pct.Valid      100.00    100.00

Conclusiones

  1. La mayor concentración de viviendas en Cali están en tipo apartamento, con un 79.89% del total demandado en la ciudad

  2. El mayor precio de metro cuadrado por zona, lo tiene la oeste, con más de $2.98 millones por metro

  3. El estrato 4, tiene el 100% de las viviendas analizadas, lo cual coteja de buena forma con la oferta de la zona oeste, mayor zona de precio por metro cuadrado

  4. El promedio de metro cuadrado construido es de 87.63 metros cuadrados, con un precio promedio de $243.7 millones


Punto 2. Análisis exploratorio bivariado

plot(viviendas$areaconst,viviendas$preciom, xlab = "Area construida", ylab="Precio millones", main="Analisis Exploratorio bivariado Area/Precio" )

El gráfico de dispersión permite ver con claridad la relación positiva de las variables. El cálculo a continuación funcionará como validador de lo visualizado previamente

cor(viviendas$preciom,viviendas$areaconst)
## [1] 0.9309803

Al generar el coeficiente de correlación, y segun la tabla de clasificación de la relación lineal entre las variables, la variable areaconst y preciom tienen una relación positiva fuerte, mayor al 93%. Se puede concluir que si el area construida sube, el precio también lo hará, resultado alineado a la expectativa lógica de encontrar mayores precios ante mayores áreas de construción

Punto 3. Modelo de regresión lineal

Para la construcción del modelo de regresión lineal se utilizará la función “lm”, la cual se utiliza para ajustar modelos lineales, incluidos los multivariados. Se puede utilizar para realizar regresión, análisis de varianza de estrato único y análisis de covarianza

Modelo_lineal = lm(viviendas$preciom ~ viviendas$areaconst)
summary(Modelo_lineal)
## 
## Call:
## lm(formula = viviendas$preciom ~ viviendas$areaconst)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -26.5997  -5.0198  -0.0056   4.6648  24.4010 
## 
## Coefficients:
##                      Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         1.998e+02  4.514e-01   442.7   <2e-16 ***
## viviendas$areaconst 5.009e-01  4.758e-03   105.3   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 7.141 on 1704 degrees of freedom
## Multiple R-squared:  0.8667, Adjusted R-squared:  0.8666 
## F-statistic: 1.108e+04 on 1 and 1704 DF,  p-value: < 2.2e-16
plot(viviendas$areaconst, viviendas$preciom, xlab = "Area construida", ylab="Precio millones", main="Regresión Lineal" )
abline(Modelo_lineal,col="red")

Dados los resultados obtenidos podemos identficar un modelo lineal donde el precio se mueve bajo la siguiente regresión: 199.8 + 0.5*areaconst.

En este caso, tenemos una pendiente positiva en el modelo, es decir, una relación directa en la cual se determina que, en la medida que el area construida va en aumento, igual lo hará el precio del inmueble. Esta relación es de $0.5 millones por cada metro adicional.

Ante la ausencia de incremento en el metraje del inmueble, el precio base del modelo se determinó en $199.8 millones, esto bajo un coeficiente de determinación del 86.6%, y dos variables altamente significativas estadísticamente, con probabilidades que tienden a cero

par(mfrow=c(2,2))
plot(Modelo_lineal)

Punto 4. Intervalo de confianza al 95%

Para el desarrollo de este punto se utiliza la función confint de R, la cual Calcula intervalos de confianza para uno o más parámetros en un modelo ajustado. Hay un valor predeterminado y un método para los objetos que heredan de la clase “lm”

confint(Modelo_lineal,)
##                           2.5 %      97.5 %
## (Intercept)         198.9248215 200.6954749
## viviendas$areaconst   0.4915592   0.5102243

Del cálculo obtenido se puede determinar un intervalo de confianza para la variable independiente, entre 0.4915 y 0.5102 millones, esto a un 95% de nivel de confianza. El resultado permite confirmar que el coeficiente B1 no es igual a cero. De igual manera, al ser B0 y B1 distintos de cero dentro del modelo, se puede aceptar la prueba de hipótesis.

Punto 5. Indicador de bondad R2

summary(Modelo_lineal)
## 
## Call:
## lm(formula = viviendas$preciom ~ viviendas$areaconst)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -26.5997  -5.0198  -0.0056   4.6648  24.4010 
## 
## Coefficients:
##                      Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         1.998e+02  4.514e-01   442.7   <2e-16 ***
## viviendas$areaconst 5.009e-01  4.758e-03   105.3   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 7.141 on 1704 degrees of freedom
## Multiple R-squared:  0.8667, Adjusted R-squared:  0.8666 
## F-statistic: 1.108e+04 on 1 and 1704 DF,  p-value: < 2.2e-16

Mediante el uso del comando summary de la función lineal, se puede traer el cálculo de R2, encontrando un valor del 86.7%, un nivel de determinación importante y significativo, teniendo en cuenta el nivel de significancia del intercepto y la variable independiente

Punto 6. Precio promedio estimado

Para el cálculo estimado del precio de vivienda, se pueden utilizar diferentes soluciones, que van desde funciones de predicción hasta construcciones de resultados en función de los estimadores ya calculados. Para este ejercicio se realizarán dos cálculos por separado para confrontar resultados

Función de predicción

Se utilizará la función predict() para este primer cálculo, la cual es una función genérica para predicciones a partir de los resultados de varias funciones de ajuste de modelos. La función invoca métodos particulares que dependen de la clase del primer argumento

predict(Modelo_lineal, newdata = list(areaconst = 110))
##        1        2        3        4        5        6        7        8 
## 225.8565 279.9528 253.9065 247.8958 240.8833 258.4145 237.3770 229.8637 
##        9       10       11       12       13       14       15       16 
## 241.8851 258.4145 258.9154 229.8637 237.3770 237.8779 227.3592 230.8654 
##       17       18       19       20       21       22       23       24 
## 229.8637 257.4127 229.8637 231.3663 236.8761 228.8619 232.3681 232.3681 
##       25       26       27       28       29       30       31       32 
## 231.3663 229.3628 228.8619 294.9796 298.9867 289.9707 263.4234 268.4323 
##       33       34       35       36       37       38       39       40 
## 239.8815 238.8797 241.3842 224.8547 230.3645 230.3645 227.3592 231.8672 
##       41       42       43       44       45       46       47       48 
## 225.3556 229.8637 234.3717 230.3645 230.3645 230.8654 227.8601 224.8547 
##       49       50       51       52       53       54       55       56 
## 230.3645 226.3574 224.8547 224.8547 227.8601 231.3663 231.3663 226.3574 
##       57       58       59       60       61       62       63       64 
## 229.8637 259.9172 229.8637 230.8654 227.3592 229.8637 238.8797 229.8637 
##       65       66       67       68       69       70       71       72 
## 247.3949 237.3770 241.8851 262.9225 259.9172 257.4127 243.8886 230.8654 
##       73       74       75       76       77       78       79       80 
## 234.3717 229.3628 230.8654 223.3521 229.8637 242.8868 231.8672 229.8637 
##       81       82       83       84       85       86       87       88 
## 241.8851 235.8744 229.8637 234.3717 229.8637 228.3610 232.3681 234.8726 
##       89       90       91       92       93       94       95       96 
## 242.8868 235.8744 232.8690 229.8637 235.8744 229.8637 229.8637 228.8619 
##       97       98       99      100      101      102      103      104 
## 229.8637 231.3663 235.8744 232.3681 244.3895 253.4056 235.8744 239.8815 
##      105      106      107      108      109      110      111      112 
## 230.3645 253.9065 253.9065 253.9065 235.8744 253.9065 244.8904 247.8958 
##      113      114      115      116      117      118      119      120 
## 299.4876 227.3592 229.8637 232.3681 233.8708 234.8726 235.3735 235.3735 
##      121      122      123      124      125      126      127      128 
## 228.8619 231.3663 229.8637 237.3770 244.8904 250.9011 225.3556 233.8708 
##      129      130      131      132      133      134      135      136 
## 252.9047 242.3859 236.3752 242.8868 244.8904 279.9528 299.9885 259.9172 
##      137      138      139      140      141      142      143      144 
## 226.8583 224.8547 244.8904 264.9261 259.9172 236.8761 259.9172 268.4323 
##      145      146      147      148      149      150      151      152 
## 264.9261 244.3895 231.8672 254.9082 237.8779 262.9225 253.9065 261.9207 
##      153      154      155      156      157      158      159      160 
## 299.9885 234.8726 239.8815 224.8547 227.8601 227.3592 237.3770 225.8565 
##      161      162      163      164      165      166      167      168 
## 289.9707 282.9582 233.8708 266.9296 249.8993 230.8654 270.9368 229.8637 
##      169      170      171      172      173      174      175      176 
## 232.8690 284.9617 269.9350 284.9617 229.8637 238.3788 261.9207 248.8975 
##      177      178      179      180      181      182      183      184 
## 266.4287 231.3663 227.8601 242.3859 229.3628 224.3538 226.3574 238.8797 
##      185      186      187      188      189      190      191      192 
## 247.8958 229.8637 229.8637 232.3681 228.8619 225.3556 228.3610 232.3681 
##      193      194      195      196      197      198      199      200 
## 229.8637 229.8637 227.3592 244.3895 226.8583 239.8815 228.8619 242.3859 
##      201      202      203      204      205      206      207      208 
## 242.3859 257.4127 242.3859 290.4716 240.8833 247.3949 225.3556 225.3556 
##      209      210      211      212      213      214      215      216 
## 226.3574 242.8868 294.9796 252.9047 272.4395 241.3842 289.9707 251.9029 
##      217      218      219      220      221      222      223      224 
## 299.9885 258.9154 244.8904 239.8815 247.3949 248.3966 244.8904 247.8958 
##      225      226      227      228      229      230      231      232 
## 237.3770 262.4216 272.4395 237.3770 270.9368 295.9814 287.4662 289.9707 
##      233      234      235      236      237      238      239      240 
## 278.4502 262.9225 299.9885 287.4662 248.8975 246.8940 259.9172 244.8904 
##      241      242      243      244      245      246      247      248 
## 230.8654 279.9528 289.9707 259.9172 244.8904 247.3949 274.9439 294.9796 
##      249      250      251      252      253      254      255      256 
## 294.9796 299.9885 278.9510 281.4555 284.9617 280.9546 289.9707 292.4751 
##      257      258      259      260      261      262      263      264 
## 289.9707 299.9885 289.9707 244.8904 244.8904 232.3681 294.9796 234.8726 
##      265      266      267      268      269      270      271      272 
## 231.8672 234.8726 272.4395 299.9885 289.9707 289.9707 276.9475 234.3717 
##      273      274      275      276      277      278      279      280 
## 232.3681 232.3681 251.4020 244.8904 269.9350 299.9885 289.9707 292.4751 
##      281      282      283      284      285      286      287      288 
## 229.8637 274.9439 244.8904 280.9546 269.9350 264.9261 259.9172 232.3681 
##      289      290      291      292      293      294      295      296 
## 227.8601 283.9600 274.9439 289.9707 295.9814 249.3984 248.8975 235.8744 
##      297      298      299      300      301      302      303      304 
## 226.8583 235.8744 230.3645 229.8637 231.3663 230.8654 231.3663 232.3681 
##      305      306      307      308      309      310      311      312 
## 235.8744 229.8637 231.3663 231.3663 232.8690 229.3628 230.8654 228.8619 
##      313      314      315      316      317      318      319      320 
## 238.3788 239.8815 237.3770 236.3752 238.3788 238.3788 231.3663 229.8637 
##      321      322      323      324      325      326      327      328 
## 233.3699 234.3717 229.8637 242.8868 235.3735 237.3770 256.9118 247.8958 
##      329      330      331      332      333      334      335      336 
## 237.3770 241.8851 248.8975 238.3788 230.8654 237.3770 227.8601 230.3645 
##      337      338      339      340      341      342      343      344 
## 229.8637 274.9439 229.3628 237.3770 238.8797 223.3521 223.3521 299.9885 
##      345      346      347      348      349      350      351      352 
## 289.9707 235.3735 289.9707 229.8637 229.8637 229.8637 229.8637 239.8815 
##      353      354      355      356      357      358      359      360 
## 229.8637 264.9261 299.9885 231.8672 289.9707 260.4180 227.8601 271.4377 
##      361      362      363      364      365      366      367      368 
## 229.8637 227.3592 236.8761 229.8637 234.8726 232.8690 238.8797 236.3752 
##      369      370      371      372      373      374      375      376 
## 299.9885 227.8601 283.4591 238.3788 299.9885 240.8833 268.9332 235.8744 
##      377      378      379      380      381      382      383      384 
## 299.9885 261.9207 232.3681 238.8797 261.4198 262.4216 234.8726 236.3752 
##      385      386      387      388      389      390      391      392 
## 229.8637 259.9172 247.3949 244.8904 260.9189 235.8744 235.8744 239.8815 
##      393      394      395      396      397      398      399      400 
## 264.4252 244.8904 247.8958 282.4573 267.9314 267.9314 280.9546 229.3628 
##      401      402      403      404      405      406      407      408 
## 235.8744 229.8637 232.3681 257.4127 239.8815 233.8708 229.8637 232.3681 
##      409      410      411      412      413      414      415      416 
## 232.3681 226.3574 239.8815 279.9528 281.9564 229.8637 239.8815 228.3610 
##      417      418      419      420      421      422      423      424 
## 227.3592 228.3610 227.3592 228.8619 228.3610 229.8637 246.8940 290.9724 
##      425      426      427      428      429      430      431      432 
## 268.9332 225.8565 229.8637 228.3610 285.9635 284.9617 299.9885 293.9778 
##      433      434      435      436      437      438      439      440 
## 234.8726 241.3842 226.8583 289.9707 244.8904 249.8993 226.8583 226.8583 
##      441      442      443      444      445      446      447      448 
## 279.9528 226.3574 232.3681 284.9617 249.8993 287.4662 274.9439 262.4216 
##      449      450      451      452      453      454      455      456 
## 284.9617 242.8868 299.9885 271.9386 235.8744 229.8637 228.8619 229.8637 
##      457      458      459      460      461      462      463      464 
## 227.3592 227.3592 227.8601 229.8637 259.9172 227.8601 239.8815 227.8601 
##      465      466      467      468      469      470      471      472 
## 227.8601 244.8904 295.9814 244.8904 238.8797 240.8833 274.9439 230.3645 
##      473      474      475      476      477      478      479      480 
## 245.8922 229.8637 226.8583 228.8619 228.8619 298.9867 228.8619 245.8922 
##      481      482      483      484      485      486      487      488 
## 227.3592 226.8583 226.8583 240.8833 240.8833 227.8601 223.8530 224.8547 
##      489      490      491      492      493      494      495      496 
## 226.3574 229.8637 228.8619 239.3806 224.8547 226.8583 247.8958 247.8958 
##      497      498      499      500      501      502      503      504 
## 237.3770 229.8637 245.8922 245.8922 229.8637 222.3503 229.8637 229.8637 
##      505      506      507      508      509      510      511      512 
## 228.3610 229.8637 229.8637 228.8619 229.3628 231.3663 231.3663 228.8619 
##      513      514      515      516      517      518      519      520 
## 242.3859 244.8904 241.8851 274.9439 284.9617 229.8637 232.3681 274.9439 
##      521      522      523      524      525      526      527      528 
## 228.3610 229.8637 231.8672 227.3592 225.8565 229.3628 255.4091 277.9493 
##      529      530      531      532      533      534      535      536 
## 290.4716 229.8637 279.9528 244.8904 229.8637 244.8904 259.9172 244.8904 
##      537      538      539      540      541      542      543      544 
## 237.8779 245.8922 265.9279 284.9617 232.3681 224.3538 230.8654 262.9225 
##      545      546      547      548      549      550      551      552 
## 244.8904 272.4395 298.9867 296.4823 242.3859 298.4858 287.9671 275.4448 
##      553      554      555      556      557      558      559      560 
## 262.4216 260.4180 268.9332 250.4002 297.4840 245.3913 254.9082 279.9528 
##      561      562      563      564      565      566      567      568 
## 234.3717 243.8886 249.8993 254.9082 234.3717 246.3931 235.8744 230.8654 
##      569      570      571      572      573      574      575      576 
## 225.8565 231.3663 232.3681 246.8940 225.8565 231.3663 242.8868 229.8637 
##      577      578      579      580      581      582      583      584 
## 237.3770 244.8904 244.8904 237.3770 231.8672 234.3717 239.3806 227.3592 
##      585      586      587      588      589      590      591      592 
## 244.8904 247.3949 241.3842 252.4038 242.3859 247.3949 245.8922 244.8904 
##      593      594      595      596      597      598      599      600 
## 252.4038 248.3966 237.3770 234.8726 253.9065 237.3770 231.3663 233.8708 
##      601      602      603      604      605      606      607      608 
## 237.3770 242.3859 237.3770 254.9082 271.9386 230.8654 229.8637 242.3859 
##      609      610      611      612      613      614      615      616 
## 242.3859 240.8833 254.9082 249.8993 267.9314 259.9172 236.3752 237.3770 
##      617      618      619      620      621      622      623      624 
## 264.9261 237.3770 245.8922 277.4484 272.4395 238.8797 269.9350 249.8993 
##      625      626      627      628      629      630      631      632 
## 238.8797 299.9885 238.8797 235.8744 239.8815 264.9261 237.3770 237.3770 
##      633      634      635      636      637      638      639      640 
## 233.8708 264.9261 274.9439 256.4109 256.4109 256.4109 264.9261 227.3592 
##      641      642      643      644      645      646      647      648 
## 237.3770 254.9082 244.8904 246.8940 244.8904 234.3717 227.3592 260.9189 
##      649      650      651      652      653      654      655      656 
## 246.8940 246.8940 246.8940 234.8726 251.9029 236.8761 227.8601 237.3770 
##      657      658      659      660      661      662      663      664 
## 236.8761 249.8993 236.8761 242.8868 242.8868 233.8708 238.3788 238.3788 
##      665      666      667      668      669      670      671      672 
## 238.3788 233.8708 254.9082 238.3788 237.8779 239.8815 224.8547 238.8797 
##      673      674      675      676      677      678      679      680 
## 243.3877 275.9457 245.8922 262.4216 259.9172 269.9350 231.3663 238.8797 
##      681      682      683      684      685      686      687      688 
## 248.3966 235.8744 236.3752 247.8958 232.3681 242.3859 229.8637 228.3610 
##      689      690      691      692      693      694      695      696 
## 227.8601 284.9617 234.8726 244.8904 254.9082 232.3681 235.8744 297.9849 
##      697      698      699      700      701      702      703      704 
## 241.8851 242.8868 245.8922 242.3859 238.3788 247.3949 248.8975 242.3859 
##      705      706      707      708      709      710      711      712 
## 253.9065 248.8975 236.3752 262.9225 237.3770 239.8815 237.8779 261.4198 
##      713      714      715      716      717      718      719      720 
## 247.3949 226.8583 229.8637 237.8779 244.8904 249.8993 244.8904 241.3842 
##      721      722      723      724      725      726      727      728 
## 236.3752 232.3681 264.9261 262.4216 294.9796 239.8815 259.9172 229.8637 
##      729      730      731      732      733      734      735      736 
## 229.8637 242.3859 289.9707 229.8637 232.3681 229.8637 242.3859 254.9082 
##      737      738      739      740      741      742      743      744 
## 232.3681 232.3681 232.3681 249.8993 238.8797 249.8993 249.3984 232.3681 
##      745      746      747      748      749      750      751      752 
## 250.4002 227.3592 234.8726 253.4056 229.8637 254.9082 242.3859 235.8744 
##      753      754      755      756      757      758      759      760 
## 298.9867 269.9350 263.9243 249.8993 273.4412 269.9350 299.9885 269.9350 
##      761      762      763      764      765      766      767      768 
## 295.9814 259.9172 277.4484 233.8708 234.8726 249.8993 241.8851 231.8672 
##      769      770      771      772      773      774      775      776 
## 228.8619 289.9707 267.4305 227.8601 254.9082 239.3806 244.8904 254.9082 
##      777      778      779      780      781      782      783      784 
## 242.3859 266.4287 254.9082 299.9885 243.8886 244.8904 239.3806 252.9047 
##      785      786      787      788      789      790      791      792 
## 241.8851 257.4127 299.9885 271.9386 269.9350 236.8761 255.9100 227.3592 
##      793      794      795      796      797      798      799      800 
## 235.8744 232.3681 235.8744 240.8833 259.9172 244.8904 241.3842 259.9172 
##      801      802      803      804      805      806      807      808 
## 226.3574 222.8512 242.3859 237.3770 242.3859 228.8619 228.8619 228.8619 
##      809      810      811      812      813      814      815      816 
## 241.8851 228.8619 231.8672 229.8637 229.8637 229.8637 228.8619 228.8619 
##      817      818      819      820      821      822      823      824 
## 228.8619 260.4180 229.8637 232.3681 236.3752 231.3663 226.3574 246.8940 
##      825      826      827      828      829      830      831      832 
## 227.3592 229.8637 223.8530 284.4609 267.4305 271.4377 241.8851 235.3735 
##      833      834      835      836      837      838      839      840 
## 299.9885 299.9885 246.8940 246.8940 298.9867 289.9707 248.8975 249.8993 
##      841      842      843      844      845      846      847      848 
## 251.4020 232.8690 249.8993 248.8975 252.4038 224.8547 279.9528 272.4395 
##      849      850      851      852      853      854      855      856 
## 259.9172 229.8637 232.8690 238.3788 240.8833 239.8815 264.9261 239.3806 
##      857      858      859      860      861      862      863      864 
## 238.8797 236.8761 238.8797 239.8815 238.3788 237.8779 238.3788 240.8833 
##      865      866      867      868      869      870      871      872 
## 239.8815 239.8815 236.8761 238.8797 259.9172 238.8797 238.3788 239.8815 
##      873      874      875      876      877      878      879      880 
## 249.3984 235.8744 238.3788 252.9047 249.8993 254.9082 249.3984 249.8993 
##      881      882      883      884      885      886      887      888 
## 242.3859 234.8726 235.8744 291.9742 299.9885 245.8922 241.8851 291.4733 
##      889      890      891      892      893      894      895      896 
## 231.8672 279.9528 241.8851 236.3752 235.8744 241.3842 235.8744 237.3770 
##      897      898      899      900      901      902      903      904 
## 249.3984 236.8761 235.8744 235.3735 238.3788 237.3770 244.8904 238.3788 
##      905      906      907      908      909      910      911      912 
## 238.3788 230.3645 236.8761 238.3788 238.3788 229.8637 231.3663 230.3645 
##      913      914      915      916      917      918      919      920 
## 227.3592 229.3628 231.8672 228.3610 230.8654 228.3610 228.3610 229.3628 
##      921      922      923      924      925      926      927      928 
## 233.3699 228.3610 230.8654 233.8708 230.8654 236.3752 234.8726 234.8726 
##      929      930      931      932      933      934      935      936 
## 236.3752 238.3788 233.3699 234.8726 236.8761 230.8654 234.8726 234.8726 
##      937      938      939      940      941      942      943      944 
## 236.3752 230.3645 235.3735 230.8654 232.3681 236.8761 233.3699 236.8761 
##      945      946      947      948      949      950      951      952 
## 233.8708 233.3699 234.8726 234.8726 235.3735 235.3735 235.3735 231.8672 
##      953      954      955      956      957      958      959      960 
## 234.8726 236.8761 229.8637 230.3645 279.9528 295.9814 232.8690 259.9172 
##      961      962      963      964      965      966      967      968 
## 299.9885 239.8815 225.8565 225.8565 256.9118 249.3984 230.8654 296.4823 
##      969      970      971      972      973      974      975      976 
## 244.8904 244.8904 240.8833 240.8833 237.8779 237.3770 228.8619 236.3752 
##      977      978      979      980      981      982      983      984 
## 230.3645 244.8904 232.3681 224.3538 257.9136 229.8637 257.4127 241.8851 
##      985      986      987      988      989      990      991      992 
## 252.4038 229.8637 234.8726 244.3895 254.9082 243.8886 287.4662 246.3931 
##      993      994      995      996      997      998      999     1000 
## 249.8993 256.9118 256.9118 257.9136 271.4377 238.8797 221.8494 294.9796 
##     1001     1002     1003     1004     1005     1006     1007     1008 
## 274.9439 279.9528 279.4519 261.4198 259.9172 271.9386 261.4198 298.9867 
##     1009     1010     1011     1012     1013     1014     1015     1016 
## 247.8958 279.4519 233.3699 291.9742 247.3949 277.4484 299.9885 245.3913 
##     1017     1018     1019     1020     1021     1022     1023     1024 
## 259.9172 259.9172 281.4555 250.9011 251.9029 232.3681 263.9243 239.8815 
##     1025     1026     1027     1028     1029     1030     1031     1032 
## 244.8904 245.8922 259.9172 263.9243 244.8904 230.3645 230.3645 279.9528 
##     1033     1034     1035     1036     1037     1038     1039     1040 
## 233.8708 238.8797 238.8797 243.8886 263.9243 262.9225 242.8868 238.8797 
##     1041     1042     1043     1044     1045     1046     1047     1048 
## 244.3895 259.9172 262.9225 299.9885 256.9118 240.8833 240.8833 240.8833 
##     1049     1050     1051     1052     1053     1054     1055     1056 
## 252.9047 294.9796 243.8886 234.8726 258.9154 249.8993 256.9118 273.4412 
##     1057     1058     1059     1060     1061     1062     1063     1064 
## 240.3824 286.4644 229.8637 229.8637 226.8583 229.3628 223.8530 232.3681 
##     1065     1066     1067     1068     1069     1070     1071     1072 
## 235.8744 254.9082 237.3770 246.8940 274.9439 245.3913 249.8993 268.4323 
##     1073     1074     1075     1076     1077     1078     1079     1080 
## 276.9475 269.9350 230.3645 227.3592 228.8619 235.8744 235.8744 235.8744 
##     1081     1082     1083     1084     1085     1086     1087     1088 
## 240.3824 234.8726 249.8993 260.4180 269.9350 238.3788 238.3788 244.8904 
##     1089     1090     1091     1092     1093     1094     1095     1096 
## 252.4038 246.8940 231.8672 226.3574 229.8637 229.8637 229.8637 236.8761 
##     1097     1098     1099     1100     1101     1102     1103     1104 
## 235.8744 259.9172 248.8975 227.8601 229.8637 227.8601 229.8637 229.8637 
##     1105     1106     1107     1108     1109     1110     1111     1112 
## 239.8815 229.8637 257.4127 229.3628 239.8815 247.8958 239.8815 225.8565 
##     1113     1114     1115     1116     1117     1118     1119     1120 
## 229.8637 228.8619 241.3842 232.3681 224.8547 229.8637 224.8547 224.8547 
##     1121     1122     1123     1124     1125     1126     1127     1128 
## 229.8637 244.8904 229.8637 237.3770 255.9100 229.8637 242.8868 238.8797 
##     1129     1130     1131     1132     1133     1134     1135     1136 
## 244.8904 244.8904 244.8904 227.8601 229.8637 235.3735 236.3752 234.8726 
##     1137     1138     1139     1140     1141     1142     1143     1144 
## 229.8637 247.3949 229.8637 228.3610 221.3485 279.9528 237.3770 247.8958 
##     1145     1146     1147     1148     1149     1150     1151     1152 
## 230.3645 227.3592 228.3610 277.4484 249.8993 237.3770 257.4127 259.9172 
##     1153     1154     1155     1156     1157     1158     1159     1160 
## 230.8654 224.3538 236.8761 250.4002 227.3592 279.9528 224.8547 224.8547 
##     1161     1162     1163     1164     1165     1166     1167     1168 
## 233.8708 227.3592 226.8583 269.9350 235.8744 258.4145 233.8708 230.8654 
##     1169     1170     1171     1172     1173     1174     1175     1176 
## 237.3770 224.8547 229.8637 230.8654 229.3628 228.8619 247.8958 229.3628 
##     1177     1178     1179     1180     1181     1182     1183     1184 
## 229.8637 227.3592 234.8726 247.3949 229.8637 236.3752 226.8583 229.8637 
##     1185     1186     1187     1188     1189     1190     1191     1192 
## 229.8637 235.8744 236.3752 235.8744 237.8779 249.8993 239.3806 234.8726 
##     1193     1194     1195     1196     1197     1198     1199     1200 
## 237.8779 244.8904 254.9082 299.9885 257.4127 251.4020 279.9528 224.8547 
##     1201     1202     1203     1204     1205     1206     1207     1208 
## 231.8672 239.8815 229.8637 239.8815 230.3645 238.8797 229.8637 237.8779 
##     1209     1210     1211     1212     1213     1214     1215     1216 
## 229.8637 230.8654 230.8654 238.8797 238.8797 226.8583 229.8637 242.3859 
##     1217     1218     1219     1220     1221     1222     1223     1224 
## 239.8815 238.8797 229.8637 231.8672 229.8637 235.8744 234.3717 231.8672 
##     1225     1226     1227     1228     1229     1230     1231     1232 
## 228.8619 229.8637 228.8619 229.8637 231.8672 239.8815 261.4198 230.8654 
##     1233     1234     1235     1236     1237     1238     1239     1240 
## 229.8637 231.8672 229.3628 299.9885 228.8619 229.8637 228.8619 229.8637 
##     1241     1242     1243     1244     1245     1246     1247     1248 
## 228.8619 244.8904 236.3752 226.3574 232.3681 230.8654 229.8637 249.8993 
##     1249     1250     1251     1252     1253     1254     1255     1256 
## 234.8726 249.8993 289.9707 261.9207 232.3681 243.3877 248.8975 252.4038 
##     1257     1258     1259     1260     1261     1262     1263     1264 
## 235.8744 229.8637 261.4198 279.9528 229.8637 295.9814 231.3663 233.3699 
##     1265     1266     1267     1268     1269     1270     1271     1272 
## 262.9225 234.8726 230.3645 225.3556 229.8637 227.8601 229.3628 232.3681 
##     1273     1274     1275     1276     1277     1278     1279     1280 
## 239.8815 245.3913 223.8530 229.3628 243.3877 229.8637 239.8815 228.8619 
##     1281     1282     1283     1284     1285     1286     1287     1288 
## 223.3521 229.8637 244.8904 229.8637 224.8547 227.8601 227.3592 224.8547 
##     1289     1290     1291     1292     1293     1294     1295     1296 
## 227.3592 239.8815 236.3752 229.8637 227.3592 229.8637 237.3770 227.3592 
##     1297     1298     1299     1300     1301     1302     1303     1304 
## 230.3645 234.8726 239.8815 234.8726 240.8833 240.3824 240.3824 234.8726 
##     1305     1306     1307     1308     1309     1310     1311     1312 
## 240.3824 240.3824 245.8922 238.8797 229.8637 227.3592 230.8654 224.8547 
##     1313     1314     1315     1316     1317     1318     1319     1320 
## 232.3681 228.3610 229.8637 229.3628 229.8637 232.3681 244.8904 227.8601 
##     1321     1322     1323     1324     1325     1326     1327     1328 
## 244.8904 239.3806 245.3913 225.8565 230.8654 235.8744 254.9082 254.9082 
##     1329     1330     1331     1332     1333     1334     1335     1336 
## 225.3556 228.3610 236.8761 227.8601 229.3628 229.3628 228.3610 233.3699 
##     1337     1338     1339     1340     1341     1342     1343     1344 
## 229.8637 227.8601 229.8637 229.8637 229.8637 225.8565 237.3770 233.8708 
##     1345     1346     1347     1348     1349     1350     1351     1352 
## 227.3592 227.3592 229.8637 229.8637 227.3592 235.8744 227.3592 227.3592 
##     1353     1354     1355     1356     1357     1358     1359     1360 
## 233.3699 228.3610 227.8601 229.3628 234.8726 233.8708 243.3877 235.8744 
##     1361     1362     1363     1364     1365     1366     1367     1368 
## 241.3842 224.8547 241.8851 229.3628 226.8583 230.8654 229.8637 234.3717 
##     1369     1370     1371     1372     1373     1374     1375     1376 
## 229.8637 253.4056 241.8851 244.3895 238.8797 229.8637 222.3503 229.8637 
##     1377     1378     1379     1380     1381     1382     1383     1384 
## 248.8975 244.8904 227.3592 247.8958 230.8654 229.8637 245.3913 239.8815 
##     1385     1386     1387     1388     1389     1390     1391     1392 
## 229.8637 247.3949 239.8815 229.8637 229.3628 238.3788 229.8637 235.8744 
##     1393     1394     1395     1396     1397     1398     1399     1400 
## 237.8779 229.8637 234.8726 230.3645 233.8708 245.8922 227.3592 229.8637 
##     1401     1402     1403     1404     1405     1406     1407     1408 
## 229.8637 233.8708 241.3842 240.3824 227.3592 232.3681 234.8726 229.8637 
##     1409     1410     1411     1412     1413     1414     1415     1416 
## 229.8637 229.8637 229.8637 229.8637 234.8726 235.8744 238.3788 241.8851 
##     1417     1418     1419     1420     1421     1422     1423     1424 
## 229.8637 238.8797 244.8904 238.8797 237.3770 233.8708 239.3806 245.3913 
##     1425     1426     1427     1428     1429     1430     1431     1432 
## 229.8637 237.8779 233.8708 229.8637 234.8726 232.3681 237.3770 231.3663 
##     1433     1434     1435     1436     1437     1438     1439     1440 
## 228.8619 228.8619 238.8797 236.3752 236.8761 234.8726 232.3681 234.8726 
##     1441     1442     1443     1444     1445     1446     1447     1448 
## 229.8637 234.8726 235.8744 239.8815 237.3770 239.8815 240.8833 246.8940 
##     1449     1450     1451     1452     1453     1454     1455     1456 
## 237.3770 234.3717 235.3735 238.8797 244.8904 237.8779 244.3895 244.3895 
##     1457     1458     1459     1460     1461     1462     1463     1464 
## 234.8726 231.3663 234.8726 239.8815 242.8868 240.3824 248.8975 236.3752 
##     1465     1466     1467     1468     1469     1470     1471     1472 
## 242.8868 248.8975 236.3752 238.8797 239.8815 238.8797 236.3752 246.8940 
##     1473     1474     1475     1476     1477     1478     1479     1480 
## 238.8797 248.8975 238.8797 244.3895 238.8797 238.8797 238.8797 243.3877 
##     1481     1482     1483     1484     1485     1486     1487     1488 
## 243.3877 235.8744 235.8744 230.8654 249.8993 238.3788 229.8637 230.3645 
##     1489     1490     1491     1492     1493     1494     1495     1496 
## 229.3628 229.8637 238.8797 229.8637 245.8922 235.8744 230.8654 229.8637 
##     1497     1498     1499     1500     1501     1502     1503     1504 
## 229.8637 229.8637 229.8637 234.8726 234.3717 234.8726 238.3788 238.8797 
##     1505     1506     1507     1508     1509     1510     1511     1512 
## 234.3717 238.8797 244.8904 235.8744 238.8797 229.8637 233.8708 240.8833 
##     1513     1514     1515     1516     1517     1518     1519     1520 
## 240.8833 242.3859 253.4056 238.8797 234.8726 236.8761 239.8815 227.3592 
##     1521     1522     1523     1524     1525     1526     1527     1528 
## 227.3592 227.3592 226.8583 224.8547 227.3592 254.4073 237.3770 234.8726 
##     1529     1530     1531     1532     1533     1534     1535     1536 
## 232.3681 283.9600 249.8993 228.8619 225.8565 229.8637 229.8637 249.8993 
##     1537     1538     1539     1540     1541     1542     1543     1544 
## 233.8708 234.3717 236.8761 227.3592 234.8726 235.8744 226.3574 237.3770 
##     1545     1546     1547     1548     1549     1550     1551     1552 
## 237.3770 229.8637 229.8637 229.8637 231.3663 236.8761 232.3681 236.8761 
##     1553     1554     1555     1556     1557     1558     1559     1560 
## 224.8547 243.3877 234.8726 229.8637 225.8565 229.3628 229.8637 228.8619 
##     1561     1562     1563     1564     1565     1566     1567     1568 
## 229.8637 299.9885 227.8601 223.8530 230.8654 230.3645 260.9189 229.8637 
##     1569     1570     1571     1572     1573     1574     1575     1576 
## 224.8547 234.8726 237.3770 241.3842 235.8744 234.8726 239.8815 230.8654 
##     1577     1578     1579     1580     1581     1582     1583     1584 
## 235.8744 232.3681 234.8726 227.3592 274.9439 250.9011 228.8619 229.8637 
##     1585     1586     1587     1588     1589     1590     1591     1592 
## 250.9011 229.8637 232.3681 229.3628 230.8654 223.8530 246.8940 223.8530 
##     1593     1594     1595     1596     1597     1598     1599     1600 
## 239.8815 244.8904 236.3752 229.8637 229.8637 241.8851 254.9082 245.8922 
##     1601     1602     1603     1604     1605     1606     1607     1608 
## 238.8797 238.3788 240.3824 241.8851 254.9082 246.8940 230.8654 241.3842 
##     1609     1610     1611     1612     1613     1614     1615     1616 
## 229.8637 244.8904 243.3877 240.8833 239.8815 239.8815 259.9172 241.8851 
##     1617     1618     1619     1620     1621     1622     1623     1624 
## 243.3877 240.3824 230.3645 259.9172 235.8744 238.8797 227.8601 244.8904 
##     1625     1626     1627     1628     1629     1630     1631     1632 
## 235.8744 234.8726 240.3824 239.3806 239.8815 240.3824 242.8868 232.8690 
##     1633     1634     1635     1636     1637     1638     1639     1640 
## 237.3770 230.8654 231.3663 249.8993 259.9172 229.8637 289.9707 224.8547 
##     1641     1642     1643     1644     1645     1646     1647     1648 
## 230.8654 231.3663 233.8708 227.3592 219.8458 229.8637 280.9546 226.8583 
##     1649     1650     1651     1652     1653     1654     1655     1656 
## 241.8851 251.4020 242.3859 258.9154 275.9457 236.3752 252.9047 262.9225 
##     1657     1658     1659     1660     1661     1662     1663     1664 
## 251.4020 255.4091 249.8993 236.3752 238.8797 259.9172 259.9172 259.9172 
##     1665     1666     1667     1668     1669     1670     1671     1672 
## 259.9172 237.3770 236.3752 229.8637 238.3788 231.3663 240.8833 236.3752 
##     1673     1674     1675     1676     1677     1678     1679     1680 
## 257.4127 242.8868 285.4626 295.9814 249.8993 278.9510 280.9546 242.3859 
##     1681     1682     1683     1684     1685     1686     1687     1688 
## 226.8583 259.9172 237.3770 249.8993 235.8744 236.8761 247.3949 259.9172 
##     1689     1690     1691     1692     1693     1694     1695     1696 
## 233.3699 224.8547 280.9546 241.8851 244.8904 259.9172 256.9118 264.9261 
##     1697     1698     1699     1700     1701     1702     1703     1704 
## 299.9885 299.9885 263.4234 289.4698 226.8583 231.8672 233.8708 223.3521 
##     1705     1706 
## 229.8637 245.3913

Bajo este modelo de predicción, el precio del inmueble de 110 metros cuadrados, debería ser de $253.91 millones

Cálculo de estimadores

intercept <- Modelo_lineal$coefficients[1]
intercept
## (Intercept) 
##    199.8101
beta1 <- Modelo_lineal$coefficients[2]
beta1
## viviendas$areaconst 
##           0.5008917
area <- 110
area
## [1] 110
precio_estimado = intercept + beta1*area
print(precio_estimado)
## (Intercept) 
##    254.9082

Los cálculos arrojan un estimado de precio de $254.91 millones para un inmueble de 110 metros cuadrados. En todo caso, tanto el estimador de la función de predicción, como el construido en esta segunda metodología, llegan a resultados muy similares, los cuales, comparados al precio de $200 millones que propone el ejercicio, como precio de oferta, pareciera muy bajo versus el estimado del mercado, ya que estaría ofertando casi $55 millones por debajo del precio estimado, es decir poco menos de 22% de dicho valor.

Punto 7. Validación de supuestos del modelo

par(mfrow=c(1,2))
plot(Modelo_lineal)

mean(Modelo_lineal$residuals)
## [1] 4.504406e-17
mean(Modelo_lineal$fitted.values)
## [1] 243.7031
summarytools::descr(Modelo_lineal$residuals)
## Descriptive Statistics  
## Modelo_lineal$value  
## N: 1706  
## 
##                                     value
## ----------------- -----------------------
##              Mean                    0.00
##           Std.Dev                    7.14
##               Min                  -26.60
##                Q1                   -5.02
##            Median                   -0.01
##                Q3                    4.67
##               Max                   24.40
##               MAD                    7.19
##               IQR                    9.68
##                CV   158490841253444832.00
##          Skewness                    0.05
##       SE.Skewness                    0.06
##          Kurtosis                   -0.08
##           N.Valid                 1706.00
##         Pct.Valid                  100.00
summarytools::descr(Modelo_lineal$fitted.values)
## Descriptive Statistics  
## Modelo_lineal$value  
## N: 1706  
## 
##                       value
## ----------------- ---------
##              Mean    243.70
##           Std.Dev     18.21
##               Min    219.85
##                Q1    229.86
##            Median    237.38
##                Q3    248.90
##               Max    299.99
##               MAD     11.14
##               IQR     19.03
##                CV      0.07
##          Skewness      1.53
##       SE.Skewness      0.06
##          Kurtosis      1.68
##           N.Valid   1706.00
##         Pct.Valid    100.00

Al analizar el comportamiento de los residuales, encontramos valores medios que tienden a cero, al igual que asimetrías y curtosis casi que nula, por lo cual podemos aceptar la hipótesis de normalidad en los residuos del modelo, e interpretar que los errores son aleatorios e independientes

Como comprobación final, se aplicarán pruebas de normalidad, homocedasticidad y autocorrelación, con el fin de validar la confiabilidad del modelo

Normalidad <- shapiro.test(Modelo_lineal$residuals)
Normalidad
## 
##  Shapiro-Wilk normality test
## 
## data:  Modelo_lineal$residuals
## W = 0.99911, p-value = 0.5907
Homocedasticidad <- lmtest::bptest(Modelo_lineal)
Homocedasticidad
## 
##  studentized Breusch-Pagan test
## 
## data:  Modelo_lineal
## BP = 0.089882, df = 1, p-value = 0.7643
Autocorrelacion <- lmtest::dwtest(Modelo_lineal)
Autocorrelacion
## 
##  Durbin-Watson test
## 
## data:  Modelo_lineal
## DW = 2.0651, p-value = 0.9092
## alternative hypothesis: true autocorrelation is greater than 0

Estos valores de p=0.5907 en la test Shapiro-Wilk permiten confirmar la normalidad de los residuales del modelo. Por otro lado, la prueba de autocorrelación DW=2.0651 indican un nivel de autocorrelación negativa, y totalmente independiente

Punto 8. Transformación y ajuste del modelo

Modelo_log = lm(log(viviendas$preciom)~log(viviendas$areaconst))
summary(Modelo_log)
## 
## Call:
## lm(formula = log(viviendas$preciom) ~ log(viviendas$areaconst))
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.10654 -0.02199 -0.00015  0.02159  0.09777 
## 
## Coefficients:
##                          Estimate Std. Error t value Pr(>|t|)    
## (Intercept)              4.623576   0.009591  482.07   <2e-16 ***
## log(viviendas$areaconst) 0.197401   0.002171   90.93   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.03181 on 1704 degrees of freedom
## Multiple R-squared:  0.8291, Adjusted R-squared:  0.829 
## F-statistic:  8269 on 1 and 1704 DF,  p-value: < 2.2e-16

En coeficiente de determinación del modelo logarítmico es del 82.91%, menor al que se obtuvo con el modelo lineal inicial. Sin embargo, el error estandar se ubica en 3.18% bastante por debajo en esta ocasión.

par(mfrow=c(2,2))
plot(Modelo_log)

Normalidad <- shapiro.test(Modelo_log$residuals)
Normalidad
## 
##  Shapiro-Wilk normality test
## 
## data:  Modelo_log$residuals
## W = 0.99859, p-value = 0.1742
Homocedasticidad <- lmtest::bptest(Modelo_log)
Homocedasticidad
## 
##  studentized Breusch-Pagan test
## 
## data:  Modelo_log
## BP = 1.2509, df = 1, p-value = 0.2634
Autocorrelacion <- lmtest::dwtest(Modelo_log)
Autocorrelacion
## 
##  Durbin-Watson test
## 
## data:  Modelo_log
## DW = 1.9271, p-value = 0.06473
## alternative hypothesis: true autocorrelation is greater than 0

Las pruebas de normalidad, homocedasticidad y autocorrelación del modelo logarítmico, terminan siendo muy cercanas a las del modelo lineal.

Punto 9. Comparación de modelos

summary(Modelo_lineal)
## 
## Call:
## lm(formula = viviendas$preciom ~ viviendas$areaconst)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -26.5997  -5.0198  -0.0056   4.6648  24.4010 
## 
## Coefficients:
##                      Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         1.998e+02  4.514e-01   442.7   <2e-16 ***
## viviendas$areaconst 5.009e-01  4.758e-03   105.3   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 7.141 on 1704 degrees of freedom
## Multiple R-squared:  0.8667, Adjusted R-squared:  0.8666 
## F-statistic: 1.108e+04 on 1 and 1704 DF,  p-value: < 2.2e-16
summary(Modelo_log)
## 
## Call:
## lm(formula = log(viviendas$preciom) ~ log(viviendas$areaconst))
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.10654 -0.02199 -0.00015  0.02159  0.09777 
## 
## Coefficients:
##                          Estimate Std. Error t value Pr(>|t|)    
## (Intercept)              4.623576   0.009591  482.07   <2e-16 ***
## log(viviendas$areaconst) 0.197401   0.002171   90.93   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.03181 on 1704 degrees of freedom
## Multiple R-squared:  0.8291, Adjusted R-squared:  0.829 
## F-statistic:  8269 on 1 and 1704 DF,  p-value: < 2.2e-16

Conclusiones

Las pruebas de normalidad realizada para los modelos lineal y logarítmico, así como los valores descriptivos principales de ambos modelos, tanto en correlación, como en el comportamiento individual de los estimadores, permiten concluir que no existe motivo suficiente para querer cambiar el modelo lineal en este análisis, ya que terminan explicando en mejor porcentaje la variable dependiente, manteniendo la significancia estadística en el modelo.