| term | estimate | std.error | statistic | p.value |
|---|---|---|---|---|
| (Intercept) | 10816.043255 | 5988.34844 | 1.8061813 | 0.0982963 |
| X2 | -2227.704360 | 920.46572 | -2.4201926 | 0.0339953 |
| X3 | 1251.141202 | 1157.02064 | 1.0813473 | 0.3026793 |
| X4 | 6.282986 | 30.62166 | 0.2051811 | 0.8411781 |
| X5 | -197.399940 | 101.56123 | -1.9436545 | 0.0779553 |
El modelo estimado correspondiente es:
\[\begin{equation} \begin{aligned} \widehat{Y}=10816.0432-2227.7043X_{2t}+1251.1412X_{3t}+6.2829X_{4t}-197.3999X_{5t}. \end{aligned} \tag{1} \end{equation}\]Este modelo presenta una \(R^2\) de \(83.5\%\), lo cual indica que las variables independientes explicarían en buen grado el comportamiento del número de rosas vendidas. Para un nivel de significancia del \(5\%\) la única variable significativa dentro del modelo es el precio de venta al por mayor de las rosas por docena.
| term | estimate | std.error | statistic | p.value |
|---|---|---|---|---|
| (Intercept) | 3.5721557 | 4.6951624 | 0.7608162 | 0.4627645 |
| log(X2) | -1.1707282 | 0.4883244 | -2.3974397 | 0.0353925 |
| log(X3) | 0.7379379 | 0.6528627 | 1.1303110 | 0.2823888 |
| log(X4) | 1.1532130 | 0.9019895 | 1.2785216 | 0.2273823 |
| X5 | -0.0301108 | 0.0164188 | -1.8339251 | 0.0938364 |
La estimación para este modelo es:
\[\begin{equation} \begin{aligned} \widehat{\ln{Y}}=0.6268-1.2735\ln{X_{2t}}+0.9373\ln{X_{3t}}+1.7129\ln{X_{4t}}-0.0301X_{5t}. \end{aligned} \tag{2} \end{equation}\]Esta ecuación permite afirmar que sí el valor al por mayor de la docena de rosas incrementa un \(1\%\) esto propiciaría una reducción del \(1.27\%\) en el número de unidades vendidas. Así mismo, si el valor de los claveles aumenta un \(1\%\) conlleva a que las ventas de rosas aumenten un \(0.94\%\) respecto al ingreso familiar semanal y un \(1\%\) de aumento implica un aumento de unidades de rosas vendidas de \(1.71\%\).
Este modelo tienen una \(R^2 = 77.8\%\), lo cual indica que las variables explican en un grado aceptable la variabilidad del número de rosas vendidas, aunque nuevamente la única variable significativa dentro del modelo corresponde al logaritmo del valor del precio al por mayor de las rosas por docena.
A priori, se espera que si el valor de compra del producto aumenta, la venta del mismo se vea disminuida (oferta y demanda) y por tanto se espera que esta elasticidad sea negativa. Para el precio cruzado, al existir un aumento en el precio de un producto competidor y el precio del producto analizado se mantiene y es inferior a la competencia, se espera que la elasticidad sea positiva. Para el ingreso familiar, un incremento en su valor está directamente relacionado con un aumento de las compras del producto, luego se espera que la elasticidad también sea positiva.
Se encuentra que efectivamente se cumplen las expectativas del comportamiento de las elasticidades, es decir se refleja el efecto de aumento o disminución en términos porcentuales de la variable respuesta frente al comportamiento de cada una de las variables explicativas incluidas en el modelo.
Las elasticidades pueden calcularse mediante la relación:
\[\begin{equation} \begin{aligned} \beta_i\Big(\frac{\overline{X}_i}{\overline{Y}}\Big). \end{aligned} \tag{3} \end{equation}\]Empleando el modelo lineal propuesto inicialmente, se obtiene:
Para deducir la fórmula de la elasticidad en el precio de un modelo lineal se considera la fórmula (\(OpenStax\) , \(s.f.\)):
\[\begin{equation} \begin{aligned} \eta_i=\frac{(\% \Delta X)}{(\% \Delta Y)}=\frac{d X}{d Y}\Bigg(\frac{\overline{X}}{\overline{Y}}\Bigg)=\beta_i\Bigg(\frac{\overline{X}}{\overline{Y}}\Bigg). \end{aligned} \tag{4} \end{equation}\]Para el caso del modelo lineal, el valor de la \(R^2\) indica la proporción de la variación en la variable dependiente explicada por las variables explicativas mientras que en el modelo log log, esta misma medida indica la proporción de la variación en el logaritmo de la variable dependiente explicada por el logaritmo de las variables explicativas y por tanto no son variables que puedan asumirse iguales. Es necesario hacer un cálculo adicional que las haga comparables, para ello se calcula el valor de la \(R^2\) entre \(\ln {Y}\) y \(\widehat{\ln {Y}}\) obteniendo un valor de \(71.2\%\) el cual ya es comparable con la \(R^2\) del modelo log log de \(77.8\%\). Como este valor es inferior, bajo este criterio se escoge el modelo log log para modelar el comportamiento de la ventas de rosas ya que este modelo explica en mayor proporción la variabilidad del precio (teniendo en cuenta que la única variable significativa dentro de éste modelo corresponde a la elasticidad del precio al por mayor de las docenas de rosas).
Los resultados de la prueba MWD se muestran a continuación, donde el modelo log log (M2) muestra la mayor significancia.
| term | Estimate | Std..Error | t.value | Pr…t.. |
|---|---|---|---|---|
| M1 + log(fit(M1))-fit(M2) | -50.7356237 | 6160.1815453 | -0.0082361 | 0.9935906 |
| M2 + fit(M1)-exp(fit(M2)) | 0.0001515 | 0.0001154 | 1.3126663 | 0.2186220 |
## R version 4.2.3 (2023-03-15)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 22.04.2 LTS
##
## Matrix products: default
## BLAS: /usr/lib/x86_64-linux-gnu/blas/libblas.so.3.10.0
## LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.10.0
##
## locale:
## [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
## [3] LC_TIME=es_MX.UTF-8 LC_COLLATE=en_US.UTF-8
## [5] LC_MONETARY=es_MX.UTF-8 LC_MESSAGES=en_US.UTF-8
## [7] LC_PAPER=es_MX.UTF-8 LC_NAME=C
## [9] LC_ADDRESS=C LC_TELEPHONE=C
## [11] LC_MEASUREMENT=es_MX.UTF-8 LC_IDENTIFICATION=C
##
## attached base packages:
## [1] stats graphics grDevices utils datasets methods base
##
## other attached packages:
## [1] lmtest_0.9-40 zoo_1.8-11 knitr_1.42 corrr_0.4.4
## [5] broom_1.0.4 kableExtra_1.3.4 magrittr_2.0.3 lubridate_1.9.2
## [9] forcats_1.0.0 stringr_1.5.0 dplyr_1.1.1 purrr_1.0.1
## [13] readr_2.1.4 tidyr_1.3.0 tibble_3.2.1 ggplot2_3.4.1
## [17] tidyverse_2.0.0 BiocStyle_2.26.0
##
## loaded via a namespace (and not attached):
## [1] tidyselect_1.2.0 xfun_0.38 bslib_0.4.2
## [4] lattice_0.20-45 colorspace_2.1-0 vctrs_0.6.1
## [7] generics_0.1.3 htmltools_0.5.5 viridisLite_0.4.1
## [10] yaml_2.3.7 utf8_1.2.3 rlang_1.1.0
## [13] jquerylib_0.1.4 pillar_1.9.0 glue_1.6.2
## [16] withr_2.5.0 lifecycle_1.0.3 munsell_0.5.0
## [19] gtable_0.3.3 rvest_1.0.3 evaluate_0.20
## [22] tzdb_0.3.0 fastmap_1.1.1 fansi_1.0.4
## [25] highr_0.10 backports_1.4.1 scales_1.2.1
## [28] BiocManager_1.30.20 cachem_1.0.7 webshot_0.5.4
## [31] jsonlite_1.8.4 systemfonts_1.0.4 hms_1.1.3
## [34] digest_0.6.31 stringi_1.7.12 bookdown_0.33
## [37] grid_4.2.3 cli_3.6.1 tools_4.2.3
## [40] sass_0.4.5 pkgconfig_2.0.3 xml2_1.3.3
## [43] timechange_0.2.0 rmarkdown_2.21 svglite_2.1.1
## [46] httr_1.4.5 rstudioapi_0.14 R6_2.5.1
## [49] compiler_4.2.3