Analisis de precios de gasolina
Unidad 5: EstadĆstica inferencial
Manuel Arcadio Valenzuela Mata - 22480100
ā Bibliotecas
Aqui se escribiran las bibliotecas necesarias para poder hacer el analisis
carga de bibliotecas
bibliotecas necesarias tidyverse, car, agricolae
library(tidyverse)## āā Attaching core tidyverse packages āāāāāāāāāāāāāāāāāāāāāāāā tidyverse 2.0.0 āā
## ā dplyr 1.1.4 ā readr 2.1.5
## ā forcats 1.0.0 ā stringr 1.5.1
## ā ggplot2 3.5.2 ā tibble 3.2.1
## ā lubridate 1.9.4 ā tidyr 1.3.1
## ā purrr 1.0.4
## āā Conflicts āāāāāāāāāāāāāāāāāāāāāāāāāāāāāāāāāāāāāāāāāā tidyverse_conflicts() āā
## ā dplyr::filter() masks stats::filter()
## ā dplyr::lag() masks stats::lag()
## ā¹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errorslibrary(car)## Loading required package: carData
##
## Attaching package: 'car'
##
## The following object is masked from 'package:dplyr':
##
## recode
##
## The following object is masked from 'package:purrr':
##
## somelibrary(agricolae)
# --- NĆŗmero de decimales
options(digits = 3)se escribe el archivo del cual se hara el analisis, en este caso āgasolina_mx.csvā
gasolina_mx <- read.csv("gasolina_mx.csv", fileEncoding="UTF-8-BOM")
# Despliega encabezado
head(gasolina_mx)ajuste de datos
se crea la variable year como un factor
print(class(gasolina_mx))## [1] "data.frame"# Se transforma la variable Year de integer a factor
gasolina_mx<- mutate(gasolina_mx, Year = as.factor(Year))
# Verifica el tipo de dato de la variable que usaremos como factor
class(gasolina_mx$Year)## [1] "factor"crear variable
se crean designa que los estados de forntera sean baja california, sonora, chihuahua,coahuila, nuevo leon , tamaulipas y seran llamados fronterizos, los que no seran llamados no fronterizos
frontera = c("Baja California", "Sonora", "Chihuahua", "Coahuila", "Nuevo LeĆ³n", "Tamaulipas")
gasolina_mx <- mutate(gasolina_mx, Region = as.factor(ifelse(State %in% frontera, "Fronteriza", "No fronteriza")))
head(gasolina_mx )analisis exploratorio
se realiza un analisis por aƱo
# HINT: agrupe los datos por Region para obtener el resumen estadĆstico
gasolina_mx %>% group_by(Year) %>%
summarise(
N = n(),
Mean = mean(Price, na.rm = TRUE),
Sd = sd(Price, na.rm = TRUE),
Median = median(Price, na.rm = TRUE),
Min = min(Price, na.rm = TRUE),
Max = max(Price, na.rm = TRUE),
Range = Max - Min,
.groups = 'drop'
)se realiza un analisis por estado
# HINT: agrupe los datos por Region para obtener el resumen estadĆstico
gasolina_mx %>% group_by(State) %>%
summarise(
N = n(),
Mean = mean(Price, na.rm = TRUE),
Sd = sd(Price, na.rm = TRUE),
Median = median(Price, na.rm = TRUE),
Min = min(Price, na.rm = TRUE),
Max = max(Price, na.rm = TRUE),
Range = Max - Min,
.groups = 'drop'
)#prueba t-student
t.test(Price ~ Region, data = gasolina_mx, alternative = "two.sided", var.equal = FALSE)##
## Welch Two Sample t-test
##
## data: Price by Region
## t = -3, df = 62, p-value = 0.01
## alternative hypothesis: true difference in means between group Fronteriza and group No fronteriza is not equal to 0
## 95 percent confidence interval:
## -1.960 -0.276
## sample estimates:
## mean in group Fronteriza mean in group No fronteriza
## 19.1 20.2conclusion
dados los resultados,se puede concluir que existe una diferencia significativa entre los precios de las regiones āFronterizaā y āNo fronterizaā (p value = 0.01, menor a a=0.05) el grupo āFronterizaā tiene un precio promedio significativamente menor que el grupo āNo fronterizaā.
gasolina_mx %>%
ggplot(aes(Year, Price, fill = Region)) +
geom_boxplot(color = "#5C5C5CFF", alpha = 0.6, show.legend = TRUE) +
scale_fill_manual(values = c("lightblue","gray" )) +
# HINT: en la siguiente lĆena coloque un tĆtulo adecuado para la grĆ”fica
labs(title = "comparacion de foroneriza y no fronteriza ",
x = "AƱo",
y = "Precio (MXN)") +
theme_bw() + theme(text = element_text(size=12, face = "bold"))
validacion de supuestos
se realiza una grafica quantile-quantile
qqPlot(gasolina_mx$Price, main = "QQ Plot de Price")
## [1] 60 279shapiro.test(gasolina_mx$Price)##
## Shapiro-Wilk normality test
##
## data: gasolina_mx$Price
## W = 1, p-value = 1e-05#ANOVA se realizara un analisis de anova con bloqueo y sin bloqueo
#CON BLOQUEO
# HINT: en la siguiente lĆnea escriba el factor y el bloque adecuados para el anĆ”lisis de varianza
ANOVA_bloqueo <- aov(Price ~ State , data = gasolina_mx)
summary(ANOVA_bloqueo)## Df Sum Sq Mean Sq F value Pr(>F)
## State 31 129 4.15 0.56 0.97
## Residuals 256 1901 7.43# Prueba de Tukey
TukeyHSD(ANOVA_bloqueo)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Price ~ State, data = gasolina_mx)
##
## $State
## diff lwr upr p adj
## Baja California-Aguascalientes -0.97444 -5.89 3.94 1.000
## Baja California Sur-Aguascalientes 0.46111 -4.45 5.38 1.000
## Campeche-Aguascalientes 0.24222 -4.67 5.16 1.000
## Chiapas-Aguascalientes -0.30444 -5.22 4.61 1.000
## Chihuahua-Aguascalientes -1.85778 -6.77 3.06 1.000
## Ciudad de MĆ©xico-Aguascalientes 0.23889 -4.68 5.15 1.000
## Coahuila de Zaragoza-Aguascalientes -0.50000 -5.41 4.41 1.000
## Colima-Aguascalientes -0.08111 -5.00 4.83 1.000
## Durango-Aguascalientes 0.16778 -4.75 5.08 1.000
## Guanajuato-Aguascalientes 0.05444 -4.86 4.97 1.000
## Guerrero-Aguascalientes 0.55111 -4.36 5.47 1.000
## Hidalgo-Aguascalientes -0.49889 -5.41 4.42 1.000
## Jalisco-Aguascalientes 0.35333 -4.56 5.27 1.000
## MĆ©xico-Aguascalientes -0.18778 -5.10 4.73 1.000
## MichoacƔn de Ocampo-Aguascalientes 0.35889 -4.56 5.27 1.000
## Morelos-Aguascalientes -0.29889 -5.21 4.62 1.000
## Nayarit-Aguascalientes 0.47000 -4.44 5.38 1.000
## Nuevo LeĆ³n-Aguascalientes 0.31000 -4.60 5.22 1.000
## Oaxaca-Aguascalientes 0.36778 -4.55 5.28 1.000
## Puebla-Aguascalientes -0.51111 -5.43 4.40 1.000
## QuerƩtaro-Aguascalientes -0.36889 -5.28 4.55 1.000
## Quintana Roo-Aguascalientes 0.48889 -4.43 5.40 1.000
## San Luis PotosĆ-Aguascalientes -0.06778 -4.98 4.85 1.000
## Sinaloa-Aguascalientes 0.17333 -4.74 5.09 1.000
## Sonora-Aguascalientes -0.58444 -5.50 4.33 1.000
## Tabasco-Aguascalientes -0.77556 -5.69 4.14 1.000
## Tamaulipas-Aguascalientes -2.60889 -7.52 2.31 0.978
## Tlaxcala-Aguascalientes -0.55000 -5.46 4.36 1.000
## Veracruz de Ignacio de la Llave-Aguascalientes -0.56333 -5.48 4.35 1.000
## YucatƔn-Aguascalientes -0.04444 -4.96 4.87 1.000
## Zacatecas-Aguascalientes 0.15000 -4.76 5.06 1.000
## Baja California Sur-Baja California 1.43556 -3.48 6.35 1.000
## Campeche-Baja California 1.21667 -3.70 6.13 1.000
## Chiapas-Baja California 0.67000 -4.24 5.58 1.000
## Chihuahua-Baja California -0.88333 -5.80 4.03 1.000
## Ciudad de MĆ©xico-Baja California 1.21333 -3.70 6.13 1.000
## Coahuila de Zaragoza-Baja California 0.47444 -4.44 5.39 1.000
## Colima-Baja California 0.89333 -4.02 5.81 1.000
## Durango-Baja California 1.14222 -3.77 6.06 1.000
## Guanajuato-Baja California 1.02889 -3.89 5.94 1.000
## Guerrero-Baja California 1.52556 -3.39 6.44 1.000
## Hidalgo-Baja California 0.47556 -4.44 5.39 1.000
## Jalisco-Baja California 1.32778 -3.59 6.24 1.000
## MĆ©xico-Baja California 0.78667 -4.13 5.70 1.000
## MichoacƔn de Ocampo-Baja California 1.33333 -3.58 6.25 1.000
## Morelos-Baja California 0.67556 -4.24 5.59 1.000
## Nayarit-Baja California 1.44444 -3.47 6.36 1.000
## Nuevo LeĆ³n-Baja California 1.28444 -3.63 6.20 1.000
## Oaxaca-Baja California 1.34222 -3.57 6.26 1.000
## Puebla-Baja California 0.46333 -4.45 5.38 1.000
## QuerƩtaro-Baja California 0.60556 -4.31 5.52 1.000
## Quintana Roo-Baja California 1.46333 -3.45 6.38 1.000
## San Luis PotosĆ-Baja California 0.90667 -4.01 5.82 1.000
## Sinaloa-Baja California 1.14778 -3.77 6.06 1.000
## Sonora-Baja California 0.39000 -4.52 5.30 1.000
## Tabasco-Baja California 0.19889 -4.72 5.11 1.000
## Tamaulipas-Baja California -1.63444 -6.55 3.28 1.000
## Tlaxcala-Baja California 0.42444 -4.49 5.34 1.000
## Veracruz de Ignacio de la Llave-Baja California 0.41111 -4.50 5.33 1.000
## YucatƔn-Baja California 0.93000 -3.98 5.84 1.000
## Zacatecas-Baja California 1.12444 -3.79 6.04 1.000
## Campeche-Baja California Sur -0.21889 -5.13 4.70 1.000
## Chiapas-Baja California Sur -0.76556 -5.68 4.15 1.000
## Chihuahua-Baja California Sur -2.31889 -7.23 2.60 0.996
## Ciudad de MĆ©xico-Baja California Sur -0.22222 -5.14 4.69 1.000
## Coahuila de Zaragoza-Baja California Sur -0.96111 -5.88 3.95 1.000
## Colima-Baja California Sur -0.54222 -5.46 4.37 1.000
## Durango-Baja California Sur -0.29333 -5.21 4.62 1.000
## Guanajuato-Baja California Sur -0.40667 -5.32 4.51 1.000
## Guerrero-Baja California Sur 0.09000 -4.82 5.00 1.000
## Hidalgo-Baja California Sur -0.96000 -5.87 3.95 1.000
## Jalisco-Baja California Sur -0.10778 -5.02 4.81 1.000
## MĆ©xico-Baja California Sur -0.64889 -5.56 4.27 1.000
## MichoacƔn de Ocampo-Baja California Sur -0.10222 -5.02 4.81 1.000
## Morelos-Baja California Sur -0.76000 -5.67 4.15 1.000
## Nayarit-Baja California Sur 0.00889 -4.91 4.92 1.000
## Nuevo LeĆ³n-Baja California Sur -0.15111 -5.07 4.76 1.000
## Oaxaca-Baja California Sur -0.09333 -5.01 4.82 1.000
## Puebla-Baja California Sur -0.97222 -5.89 3.94 1.000
## QuerƩtaro-Baja California Sur -0.83000 -5.74 4.08 1.000
## Quintana Roo-Baja California Sur 0.02778 -4.89 4.94 1.000
## San Luis PotosĆ-Baja California Sur -0.52889 -5.44 4.39 1.000
## Sinaloa-Baja California Sur -0.28778 -5.20 4.63 1.000
## Sonora-Baja California Sur -1.04556 -5.96 3.87 1.000
## Tabasco-Baja California Sur -1.23667 -6.15 3.68 1.000
## Tamaulipas-Baja California Sur -3.07000 -7.98 1.84 0.867
## Tlaxcala-Baja California Sur -1.01111 -5.93 3.90 1.000
## Veracruz de Ignacio de la Llave-Baja California Sur -1.02444 -5.94 3.89 1.000
## YucatƔn-Baja California Sur -0.50556 -5.42 4.41 1.000
## Zacatecas-Baja California Sur -0.31111 -5.23 4.60 1.000
## Chiapas-Campeche -0.54667 -5.46 4.37 1.000
## Chihuahua-Campeche -2.10000 -7.01 2.81 0.999
## Ciudad de MĆ©xico-Campeche -0.00333 -4.92 4.91 1.000
## Coahuila de Zaragoza-Campeche -0.74222 -5.66 4.17 1.000
## Colima-Campeche -0.32333 -5.24 4.59 1.000
## Durango-Campeche -0.07444 -4.99 4.84 1.000
## Guanajuato-Campeche -0.18778 -5.10 4.73 1.000
## Guerrero-Campeche 0.30889 -4.61 5.22 1.000
## Hidalgo-Campeche -0.74111 -5.66 4.17 1.000
## Jalisco-Campeche 0.11111 -4.80 5.03 1.000
## MĆ©xico-Campeche -0.43000 -5.34 4.48 1.000
## MichoacƔn de Ocampo-Campeche 0.11667 -4.80 5.03 1.000
## Morelos-Campeche -0.54111 -5.46 4.37 1.000
## Nayarit-Campeche 0.22778 -4.69 5.14 1.000
## Nuevo LeĆ³n-Campeche 0.06778 -4.85 4.98 1.000
## Oaxaca-Campeche 0.12556 -4.79 5.04 1.000
## Puebla-Campeche -0.75333 -5.67 4.16 1.000
## QuerƩtaro-Campeche -0.61111 -5.53 4.30 1.000
## Quintana Roo-Campeche 0.24667 -4.67 5.16 1.000
## San Luis PotosĆ-Campeche -0.31000 -5.22 4.60 1.000
## Sinaloa-Campeche -0.06889 -4.98 4.85 1.000
## Sonora-Campeche -0.82667 -5.74 4.09 1.000
## Tabasco-Campeche -1.01778 -5.93 3.90 1.000
## Tamaulipas-Campeche -2.85111 -7.77 2.06 0.937
## Tlaxcala-Campeche -0.79222 -5.71 4.12 1.000
## Veracruz de Ignacio de la Llave-Campeche -0.80556 -5.72 4.11 1.000
## YucatƔn-Campeche -0.28667 -5.20 4.63 1.000
## Zacatecas-Campeche -0.09222 -5.01 4.82 1.000
## Chihuahua-Chiapas -1.55333 -6.47 3.36 1.000
## Ciudad de MĆ©xico-Chiapas 0.54333 -4.37 5.46 1.000
## Coahuila de Zaragoza-Chiapas -0.19556 -5.11 4.72 1.000
## Colima-Chiapas 0.22333 -4.69 5.14 1.000
## Durango-Chiapas 0.47222 -4.44 5.39 1.000
## Guanajuato-Chiapas 0.35889 -4.56 5.27 1.000
## Guerrero-Chiapas 0.85556 -4.06 5.77 1.000
## Hidalgo-Chiapas -0.19444 -5.11 4.72 1.000
## Jalisco-Chiapas 0.65778 -4.26 5.57 1.000
## MĆ©xico-Chiapas 0.11667 -4.80 5.03 1.000
## MichoacƔn de Ocampo-Chiapas 0.66333 -4.25 5.58 1.000
## Morelos-Chiapas 0.00556 -4.91 4.92 1.000
## Nayarit-Chiapas 0.77444 -4.14 5.69 1.000
## Nuevo LeĆ³n-Chiapas 0.61444 -4.30 5.53 1.000
## Oaxaca-Chiapas 0.67222 -4.24 5.59 1.000
## Puebla-Chiapas -0.20667 -5.12 4.71 1.000
## QuerƩtaro-Chiapas -0.06444 -4.98 4.85 1.000
## Quintana Roo-Chiapas 0.79333 -4.12 5.71 1.000
## San Luis PotosĆ-Chiapas 0.23667 -4.68 5.15 1.000
## Sinaloa-Chiapas 0.47778 -4.44 5.39 1.000
## Sonora-Chiapas -0.28000 -5.19 4.63 1.000
## Tabasco-Chiapas -0.47111 -5.39 4.44 1.000
## Tamaulipas-Chiapas -2.30444 -7.22 2.61 0.996
## Tlaxcala-Chiapas -0.24556 -5.16 4.67 1.000
## Veracruz de Ignacio de la Llave-Chiapas -0.25889 -5.17 4.66 1.000
## YucatƔn-Chiapas 0.26000 -4.65 5.17 1.000
## Zacatecas-Chiapas 0.45444 -4.46 5.37 1.000
## Ciudad de MĆ©xico-Chihuahua 2.09667 -2.82 7.01 0.999
## Coahuila de Zaragoza-Chihuahua 1.35778 -3.56 6.27 1.000
## Colima-Chihuahua 1.77667 -3.14 6.69 1.000
## Durango-Chihuahua 2.02556 -2.89 6.94 1.000
## Guanajuato-Chihuahua 1.91222 -3.00 6.83 1.000
## Guerrero-Chihuahua 2.40889 -2.51 7.32 0.993
## Hidalgo-Chihuahua 1.35889 -3.56 6.27 1.000
## Jalisco-Chihuahua 2.21111 -2.70 7.13 0.998
## MĆ©xico-Chihuahua 1.67000 -3.24 6.58 1.000
## MichoacƔn de Ocampo-Chihuahua 2.21667 -2.70 7.13 0.998
## Morelos-Chihuahua 1.55889 -3.36 6.47 1.000
## Nayarit-Chihuahua 2.32778 -2.59 7.24 0.996
## Nuevo LeĆ³n-Chihuahua 2.16778 -2.75 7.08 0.999
## Oaxaca-Chihuahua 2.22556 -2.69 7.14 0.998
## Puebla-Chihuahua 1.34667 -3.57 6.26 1.000
## QuerƩtaro-Chihuahua 1.48889 -3.43 6.40 1.000
## Quintana Roo-Chihuahua 2.34667 -2.57 7.26 0.995
## San Luis PotosĆ-Chihuahua 1.79000 -3.12 6.70 1.000
## Sinaloa-Chihuahua 2.03111 -2.88 6.95 1.000
## Sonora-Chihuahua 1.27333 -3.64 6.19 1.000
## Tabasco-Chihuahua 1.08222 -3.83 6.00 1.000
## Tamaulipas-Chihuahua -0.75111 -5.67 4.16 1.000
## Tlaxcala-Chihuahua 1.30778 -3.61 6.22 1.000
## Veracruz de Ignacio de la Llave-Chihuahua 1.29444 -3.62 6.21 1.000
## YucatƔn-Chihuahua 1.81333 -3.10 6.73 1.000
## Zacatecas-Chihuahua 2.00778 -2.91 6.92 1.000
## Coahuila de Zaragoza-Ciudad de MĆ©xico -0.73889 -5.65 4.18 1.000
## Colima-Ciudad de MĆ©xico -0.32000 -5.23 4.59 1.000
## Durango-Ciudad de MĆ©xico -0.07111 -4.99 4.84 1.000
## Guanajuato-Ciudad de MĆ©xico -0.18444 -5.10 4.73 1.000
## Guerrero-Ciudad de MĆ©xico 0.31222 -4.60 5.23 1.000
## Hidalgo-Ciudad de MĆ©xico -0.73778 -5.65 4.18 1.000
## Jalisco-Ciudad de MĆ©xico 0.11444 -4.80 5.03 1.000
## MĆ©xico-Ciudad de MĆ©xico -0.42667 -5.34 4.49 1.000
## MichoacƔn de Ocampo-Ciudad de MƩxico 0.12000 -4.79 5.03 1.000
## Morelos-Ciudad de MĆ©xico -0.53778 -5.45 4.38 1.000
## Nayarit-Ciudad de MĆ©xico 0.23111 -4.68 5.15 1.000
## Nuevo LeĆ³n-Ciudad de MĆ©xico 0.07111 -4.84 4.99 1.000
## Oaxaca-Ciudad de MĆ©xico 0.12889 -4.79 5.04 1.000
## Puebla-Ciudad de MĆ©xico -0.75000 -5.66 4.16 1.000
## QuerƩtaro-Ciudad de MƩxico -0.60778 -5.52 4.31 1.000
## Quintana Roo-Ciudad de MĆ©xico 0.25000 -4.66 5.16 1.000
## San Luis PotosĆ-Ciudad de MĆ©xico -0.30667 -5.22 4.61 1.000
## Sinaloa-Ciudad de MĆ©xico -0.06556 -4.98 4.85 1.000
## Sonora-Ciudad de MĆ©xico -0.82333 -5.74 4.09 1.000
## Tabasco-Ciudad de MĆ©xico -1.01444 -5.93 3.90 1.000
## Tamaulipas-Ciudad de MĆ©xico -2.84778 -7.76 2.07 0.938
## Tlaxcala-Ciudad de MĆ©xico -0.78889 -5.70 4.13 1.000
## Veracruz de Ignacio de la Llave-Ciudad de MĆ©xico -0.80222 -5.72 4.11 1.000
## YucatƔn-Ciudad de MƩxico -0.28333 -5.20 4.63 1.000
## Zacatecas-Ciudad de MĆ©xico -0.08889 -5.00 4.83 1.000
## Colima-Coahuila de Zaragoza 0.41889 -4.50 5.33 1.000
## Durango-Coahuila de Zaragoza 0.66778 -4.25 5.58 1.000
## Guanajuato-Coahuila de Zaragoza 0.55444 -4.36 5.47 1.000
## Guerrero-Coahuila de Zaragoza 1.05111 -3.86 5.97 1.000
## Hidalgo-Coahuila de Zaragoza 0.00111 -4.91 4.92 1.000
## Jalisco-Coahuila de Zaragoza 0.85333 -4.06 5.77 1.000
## MĆ©xico-Coahuila de Zaragoza 0.31222 -4.60 5.23 1.000
## MichoacƔn de Ocampo-Coahuila de Zaragoza 0.85889 -4.06 5.77 1.000
## Morelos-Coahuila de Zaragoza 0.20111 -4.71 5.12 1.000
## Nayarit-Coahuila de Zaragoza 0.97000 -3.94 5.88 1.000
## Nuevo LeĆ³n-Coahuila de Zaragoza 0.81000 -4.10 5.72 1.000
## Oaxaca-Coahuila de Zaragoza 0.86778 -4.05 5.78 1.000
## Puebla-Coahuila de Zaragoza -0.01111 -4.93 4.90 1.000
## QuerƩtaro-Coahuila de Zaragoza 0.13111 -4.78 5.05 1.000
## Quintana Roo-Coahuila de Zaragoza 0.98889 -3.93 5.90 1.000
## San Luis PotosĆ-Coahuila de Zaragoza 0.43222 -4.48 5.35 1.000
## Sinaloa-Coahuila de Zaragoza 0.67333 -4.24 5.59 1.000
## Sonora-Coahuila de Zaragoza -0.08444 -5.00 4.83 1.000
## Tabasco-Coahuila de Zaragoza -0.27556 -5.19 4.64 1.000
## Tamaulipas-Coahuila de Zaragoza -2.10889 -7.02 2.81 0.999
## Tlaxcala-Coahuila de Zaragoza -0.05000 -4.96 4.86 1.000
## Veracruz de Ignacio de la Llave-Coahuila de Zaragoza -0.06333 -4.98 4.85 1.000
## YucatƔn-Coahuila de Zaragoza 0.45556 -4.46 5.37 1.000
## Zacatecas-Coahuila de Zaragoza 0.65000 -4.26 5.56 1.000
## Durango-Colima 0.24889 -4.67 5.16 1.000
## Guanajuato-Colima 0.13556 -4.78 5.05 1.000
## Guerrero-Colima 0.63222 -4.28 5.55 1.000
## Hidalgo-Colima -0.41778 -5.33 4.50 1.000
## Jalisco-Colima 0.43444 -4.48 5.35 1.000
## MĆ©xico-Colima -0.10667 -5.02 4.81 1.000
## MichoacƔn de Ocampo-Colima 0.44000 -4.47 5.35 1.000
## Morelos-Colima -0.21778 -5.13 4.70 1.000
## Nayarit-Colima 0.55111 -4.36 5.47 1.000
## Nuevo LeĆ³n-Colima 0.39111 -4.52 5.31 1.000
## Oaxaca-Colima 0.44889 -4.47 5.36 1.000
## Puebla-Colima -0.43000 -5.34 4.48 1.000
## QuerƩtaro-Colima -0.28778 -5.20 4.63 1.000
## Quintana Roo-Colima 0.57000 -4.34 5.48 1.000
## San Luis PotosĆ-Colima 0.01333 -4.90 4.93 1.000
## Sinaloa-Colima 0.25444 -4.66 5.17 1.000
## Sonora-Colima -0.50333 -5.42 4.41 1.000
## Tabasco-Colima -0.69444 -5.61 4.22 1.000
## Tamaulipas-Colima -2.52778 -7.44 2.39 0.986
## Tlaxcala-Colima -0.46889 -5.38 4.45 1.000
## Veracruz de Ignacio de la Llave-Colima -0.48222 -5.40 4.43 1.000
## YucatƔn-Colima 0.03667 -4.88 4.95 1.000
## Zacatecas-Colima 0.23111 -4.68 5.15 1.000
## Guanajuato-Durango -0.11333 -5.03 4.80 1.000
## Guerrero-Durango 0.38333 -4.53 5.30 1.000
## Hidalgo-Durango -0.66667 -5.58 4.25 1.000
## Jalisco-Durango 0.18556 -4.73 5.10 1.000
## MĆ©xico-Durango -0.35556 -5.27 4.56 1.000
## MichoacƔn de Ocampo-Durango 0.19111 -4.72 5.11 1.000
## Morelos-Durango -0.46667 -5.38 4.45 1.000
## Nayarit-Durango 0.30222 -4.61 5.22 1.000
## Nuevo LeĆ³n-Durango 0.14222 -4.77 5.06 1.000
## Oaxaca-Durango 0.20000 -4.71 5.11 1.000
## Puebla-Durango -0.67889 -5.59 4.24 1.000
## QuerƩtaro-Durango -0.53667 -5.45 4.38 1.000
## Quintana Roo-Durango 0.32111 -4.59 5.24 1.000
## San Luis PotosĆ-Durango -0.23556 -5.15 4.68 1.000
## Sinaloa-Durango 0.00556 -4.91 4.92 1.000
## Sonora-Durango -0.75222 -5.67 4.16 1.000
## Tabasco-Durango -0.94333 -5.86 3.97 1.000
## Tamaulipas-Durango -2.77667 -7.69 2.14 0.953
## Tlaxcala-Durango -0.71778 -5.63 4.20 1.000
## Veracruz de Ignacio de la Llave-Durango -0.73111 -5.65 4.18 1.000
## YucatƔn-Durango -0.21222 -5.13 4.70 1.000
## Zacatecas-Durango -0.01778 -4.93 4.90 1.000
## Guerrero-Guanajuato 0.49667 -4.42 5.41 1.000
## Hidalgo-Guanajuato -0.55333 -5.47 4.36 1.000
## Jalisco-Guanajuato 0.29889 -4.62 5.21 1.000
## MĆ©xico-Guanajuato -0.24222 -5.16 4.67 1.000
## MichoacƔn de Ocampo-Guanajuato 0.30444 -4.61 5.22 1.000
## Morelos-Guanajuato -0.35333 -5.27 4.56 1.000
## Nayarit-Guanajuato 0.41556 -4.50 5.33 1.000
## Nuevo LeĆ³n-Guanajuato 0.25556 -4.66 5.17 1.000
## Oaxaca-Guanajuato 0.31333 -4.60 5.23 1.000
## Puebla-Guanajuato -0.56556 -5.48 4.35 1.000
## QuerƩtaro-Guanajuato -0.42333 -5.34 4.49 1.000
## Quintana Roo-Guanajuato 0.43444 -4.48 5.35 1.000
## San Luis PotosĆ-Guanajuato -0.12222 -5.04 4.79 1.000
## Sinaloa-Guanajuato 0.11889 -4.80 5.03 1.000
## Sonora-Guanajuato -0.63889 -5.55 4.28 1.000
## Tabasco-Guanajuato -0.83000 -5.74 4.08 1.000
## Tamaulipas-Guanajuato -2.66333 -7.58 2.25 0.972
## Tlaxcala-Guanajuato -0.60444 -5.52 4.31 1.000
## Veracruz de Ignacio de la Llave-Guanajuato -0.61778 -5.53 4.30 1.000
## YucatƔn-Guanajuato -0.09889 -5.01 4.82 1.000
## Zacatecas-Guanajuato 0.09556 -4.82 5.01 1.000
## Hidalgo-Guerrero -1.05000 -5.96 3.86 1.000
## Jalisco-Guerrero -0.19778 -5.11 4.72 1.000
## MĆ©xico-Guerrero -0.73889 -5.65 4.18 1.000
## MichoacƔn de Ocampo-Guerrero -0.19222 -5.11 4.72 1.000
## Morelos-Guerrero -0.85000 -5.76 4.06 1.000
## Nayarit-Guerrero -0.08111 -5.00 4.83 1.000
## Nuevo LeĆ³n-Guerrero -0.24111 -5.16 4.67 1.000
## Oaxaca-Guerrero -0.18333 -5.10 4.73 1.000
## Puebla-Guerrero -1.06222 -5.98 3.85 1.000
## QuerƩtaro-Guerrero -0.92000 -5.83 3.99 1.000
## Quintana Roo-Guerrero -0.06222 -4.98 4.85 1.000
## San Luis PotosĆ-Guerrero -0.61889 -5.53 4.30 1.000
## Sinaloa-Guerrero -0.37778 -5.29 4.54 1.000
## Sonora-Guerrero -1.13556 -6.05 3.78 1.000
## Tabasco-Guerrero -1.32667 -6.24 3.59 1.000
## Tamaulipas-Guerrero -3.16000 -8.07 1.75 0.828
## Tlaxcala-Guerrero -1.10111 -6.02 3.81 1.000
## Veracruz de Ignacio de la Llave-Guerrero -1.11444 -6.03 3.80 1.000
## YucatƔn-Guerrero -0.59556 -5.51 4.32 1.000
## Zacatecas-Guerrero -0.40111 -5.32 4.51 1.000
## Jalisco-Hidalgo 0.85222 -4.06 5.77 1.000
## MĆ©xico-Hidalgo 0.31111 -4.60 5.23 1.000
## MichoacƔn de Ocampo-Hidalgo 0.85778 -4.06 5.77 1.000
## Morelos-Hidalgo 0.20000 -4.71 5.11 1.000
## Nayarit-Hidalgo 0.96889 -3.95 5.88 1.000
## Nuevo LeĆ³n-Hidalgo 0.80889 -4.11 5.72 1.000
## Oaxaca-Hidalgo 0.86667 -4.05 5.78 1.000
## Puebla-Hidalgo -0.01222 -4.93 4.90 1.000
## QuerƩtaro-Hidalgo 0.13000 -4.78 5.04 1.000
## Quintana Roo-Hidalgo 0.98778 -3.93 5.90 1.000
## San Luis PotosĆ-Hidalgo 0.43111 -4.48 5.35 1.000
## Sinaloa-Hidalgo 0.67222 -4.24 5.59 1.000
## Sonora-Hidalgo -0.08556 -5.00 4.83 1.000
## Tabasco-Hidalgo -0.27667 -5.19 4.64 1.000
## Tamaulipas-Hidalgo -2.11000 -7.02 2.80 0.999
## Tlaxcala-Hidalgo -0.05111 -4.97 4.86 1.000
## Veracruz de Ignacio de la Llave-Hidalgo -0.06444 -4.98 4.85 1.000
## YucatƔn-Hidalgo 0.45444 -4.46 5.37 1.000
## Zacatecas-Hidalgo 0.64889 -4.27 5.56 1.000
## MĆ©xico-Jalisco -0.54111 -5.46 4.37 1.000
## MichoacƔn de Ocampo-Jalisco 0.00556 -4.91 4.92 1.000
## Morelos-Jalisco -0.65222 -5.57 4.26 1.000
## Nayarit-Jalisco 0.11667 -4.80 5.03 1.000
## Nuevo LeĆ³n-Jalisco -0.04333 -4.96 4.87 1.000
## Oaxaca-Jalisco 0.01444 -4.90 4.93 1.000
## Puebla-Jalisco -0.86444 -5.78 4.05 1.000
## QuerƩtaro-Jalisco -0.72222 -5.64 4.19 1.000
## Quintana Roo-Jalisco 0.13556 -4.78 5.05 1.000
## San Luis PotosĆ-Jalisco -0.42111 -5.34 4.49 1.000
## Sinaloa-Jalisco -0.18000 -5.09 4.73 1.000
## Sonora-Jalisco -0.93778 -5.85 3.98 1.000
## Tabasco-Jalisco -1.12889 -6.04 3.79 1.000
## Tamaulipas-Jalisco -2.96222 -7.88 1.95 0.906
## Tlaxcala-Jalisco -0.90333 -5.82 4.01 1.000
## Veracruz de Ignacio de la Llave-Jalisco -0.91667 -5.83 4.00 1.000
## YucatƔn-Jalisco -0.39778 -5.31 4.52 1.000
## Zacatecas-Jalisco -0.20333 -5.12 4.71 1.000
## MichoacƔn de Ocampo-MƩxico 0.54667 -4.37 5.46 1.000
## Morelos-MĆ©xico -0.11111 -5.03 4.80 1.000
## Nayarit-MĆ©xico 0.65778 -4.26 5.57 1.000
## Nuevo LeĆ³n-MĆ©xico 0.49778 -4.42 5.41 1.000
## Oaxaca-MĆ©xico 0.55556 -4.36 5.47 1.000
## Puebla-MĆ©xico -0.32333 -5.24 4.59 1.000
## QuerƩtaro-MƩxico -0.18111 -5.10 4.73 1.000
## Quintana Roo-MĆ©xico 0.67667 -4.24 5.59 1.000
## San Luis PotosĆ-MĆ©xico 0.12000 -4.79 5.03 1.000
## Sinaloa-MĆ©xico 0.36111 -4.55 5.28 1.000
## Sonora-MĆ©xico -0.39667 -5.31 4.52 1.000
## Tabasco-MĆ©xico -0.58778 -5.50 4.33 1.000
## Tamaulipas-MĆ©xico -2.42111 -7.34 2.49 0.992
## Tlaxcala-MĆ©xico -0.36222 -5.28 4.55 1.000
## Veracruz de Ignacio de la Llave-MĆ©xico -0.37556 -5.29 4.54 1.000
## YucatƔn-MƩxico 0.14333 -4.77 5.06 1.000
## Zacatecas-MĆ©xico 0.33778 -4.58 5.25 1.000
## Morelos-MichoacƔn de Ocampo -0.65778 -5.57 4.26 1.000
## Nayarit-MichoacƔn de Ocampo 0.11111 -4.80 5.03 1.000
## Nuevo LeĆ³n-MichoacĆ”n de Ocampo -0.04889 -4.96 4.87 1.000
## Oaxaca-MichoacƔn de Ocampo 0.00889 -4.91 4.92 1.000
## Puebla-MichoacƔn de Ocampo -0.87000 -5.78 4.04 1.000
## QuerƩtaro-MichoacƔn de Ocampo -0.72778 -5.64 4.19 1.000
## Quintana Roo-MichoacƔn de Ocampo 0.13000 -4.78 5.04 1.000
## San Luis PotosĆ-MichoacĆ”n de Ocampo -0.42667 -5.34 4.49 1.000
## Sinaloa-MichoacƔn de Ocampo -0.18556 -5.10 4.73 1.000
## Sonora-MichoacƔn de Ocampo -0.94333 -5.86 3.97 1.000
## Tabasco-MichoacƔn de Ocampo -1.13444 -6.05 3.78 1.000
## Tamaulipas-MichoacƔn de Ocampo -2.96778 -7.88 1.95 0.904
## Tlaxcala-MichoacƔn de Ocampo -0.90889 -5.82 4.01 1.000
## Veracruz de Ignacio de la Llave-MichoacƔn de Ocampo -0.92222 -5.84 3.99 1.000
## YucatƔn-MichoacƔn de Ocampo -0.40333 -5.32 4.51 1.000
## Zacatecas-MichoacƔn de Ocampo -0.20889 -5.12 4.71 1.000
## Nayarit-Morelos 0.76889 -4.15 5.68 1.000
## Nuevo LeĆ³n-Morelos 0.60889 -4.31 5.52 1.000
## Oaxaca-Morelos 0.66667 -4.25 5.58 1.000
## Puebla-Morelos -0.21222 -5.13 4.70 1.000
## QuerƩtaro-Morelos -0.07000 -4.98 4.84 1.000
## Quintana Roo-Morelos 0.78778 -4.13 5.70 1.000
## San Luis PotosĆ-Morelos 0.23111 -4.68 5.15 1.000
## Sinaloa-Morelos 0.47222 -4.44 5.39 1.000
## Sonora-Morelos -0.28556 -5.20 4.63 1.000
## Tabasco-Morelos -0.47667 -5.39 4.44 1.000
## Tamaulipas-Morelos -2.31000 -7.22 2.60 0.996
## Tlaxcala-Morelos -0.25111 -5.17 4.66 1.000
## Veracruz de Ignacio de la Llave-Morelos -0.26444 -5.18 4.65 1.000
## YucatƔn-Morelos 0.25444 -4.66 5.17 1.000
## Zacatecas-Morelos 0.44889 -4.47 5.36 1.000
## Nuevo LeĆ³n-Nayarit -0.16000 -5.07 4.75 1.000
## Oaxaca-Nayarit -0.10222 -5.02 4.81 1.000
## Puebla-Nayarit -0.98111 -5.90 3.93 1.000
## QuerƩtaro-Nayarit -0.83889 -5.75 4.08 1.000
## Quintana Roo-Nayarit 0.01889 -4.90 4.93 1.000
## San Luis PotosĆ-Nayarit -0.53778 -5.45 4.38 1.000
## Sinaloa-Nayarit -0.29667 -5.21 4.62 1.000
## Sonora-Nayarit -1.05444 -5.97 3.86 1.000
## Tabasco-Nayarit -1.24556 -6.16 3.67 1.000
## Tamaulipas-Nayarit -3.07889 -7.99 1.84 0.863
## Tlaxcala-Nayarit -1.02000 -5.93 3.89 1.000
## Veracruz de Ignacio de la Llave-Nayarit -1.03333 -5.95 3.88 1.000
## YucatƔn-Nayarit -0.51444 -5.43 4.40 1.000
## Zacatecas-Nayarit -0.32000 -5.23 4.59 1.000
## Oaxaca-Nuevo LeĆ³n 0.05778 -4.86 4.97 1.000
## Puebla-Nuevo LeĆ³n -0.82111 -5.74 4.09 1.000
## QuerĆ©taro-Nuevo LeĆ³n -0.67889 -5.59 4.24 1.000
## Quintana Roo-Nuevo LeĆ³n 0.17889 -4.74 5.09 1.000
## San Luis PotosĆ-Nuevo LeĆ³n -0.37778 -5.29 4.54 1.000
## Sinaloa-Nuevo LeĆ³n -0.13667 -5.05 4.78 1.000
## Sonora-Nuevo LeĆ³n -0.89444 -5.81 4.02 1.000
## Tabasco-Nuevo LeĆ³n -1.08556 -6.00 3.83 1.000
## Tamaulipas-Nuevo LeĆ³n -2.91889 -7.83 2.00 0.919
## Tlaxcala-Nuevo LeĆ³n -0.86000 -5.77 4.05 1.000
## Veracruz de Ignacio de la Llave-Nuevo LeĆ³n -0.87333 -5.79 4.04 1.000
## YucatĆ”n-Nuevo LeĆ³n -0.35444 -5.27 4.56 1.000
## Zacatecas-Nuevo LeĆ³n -0.16000 -5.07 4.75 1.000
## Puebla-Oaxaca -0.87889 -5.79 4.04 1.000
## QuerƩtaro-Oaxaca -0.73667 -5.65 4.18 1.000
## Quintana Roo-Oaxaca 0.12111 -4.79 5.04 1.000
## San Luis PotosĆ-Oaxaca -0.43556 -5.35 4.48 1.000
## Sinaloa-Oaxaca -0.19444 -5.11 4.72 1.000
## Sonora-Oaxaca -0.95222 -5.87 3.96 1.000
## Tabasco-Oaxaca -1.14333 -6.06 3.77 1.000
## Tamaulipas-Oaxaca -2.97667 -7.89 1.94 0.901
## Tlaxcala-Oaxaca -0.91778 -5.83 4.00 1.000
## Veracruz de Ignacio de la Llave-Oaxaca -0.93111 -5.85 3.98 1.000
## YucatƔn-Oaxaca -0.41222 -5.33 4.50 1.000
## Zacatecas-Oaxaca -0.21778 -5.13 4.70 1.000
## QuerƩtaro-Puebla 0.14222 -4.77 5.06 1.000
## Quintana Roo-Puebla 1.00000 -3.91 5.91 1.000
## San Luis PotosĆ-Puebla 0.44333 -4.47 5.36 1.000
## Sinaloa-Puebla 0.68444 -4.23 5.60 1.000
## Sonora-Puebla -0.07333 -4.99 4.84 1.000
## Tabasco-Puebla -0.26444 -5.18 4.65 1.000
## Tamaulipas-Puebla -2.09778 -7.01 2.82 0.999
## Tlaxcala-Puebla -0.03889 -4.95 4.88 1.000
## Veracruz de Ignacio de la Llave-Puebla -0.05222 -4.97 4.86 1.000
## YucatƔn-Puebla 0.46667 -4.45 5.38 1.000
## Zacatecas-Puebla 0.66111 -4.25 5.58 1.000
## Quintana Roo-QuerƩtaro 0.85778 -4.06 5.77 1.000
## San Luis PotosĆ-QuerĆ©taro 0.30111 -4.61 5.22 1.000
## Sinaloa-QuerƩtaro 0.54222 -4.37 5.46 1.000
## Sonora-QuerƩtaro -0.21556 -5.13 4.70 1.000
## Tabasco-QuerƩtaro -0.40667 -5.32 4.51 1.000
## Tamaulipas-QuerƩtaro -2.24000 -7.15 2.67 0.998
## Tlaxcala-QuerƩtaro -0.18111 -5.10 4.73 1.000
## Veracruz de Ignacio de la Llave-QuerƩtaro -0.19444 -5.11 4.72 1.000
## YucatƔn-QuerƩtaro 0.32444 -4.59 5.24 1.000
## Zacatecas-QuerƩtaro 0.51889 -4.40 5.43 1.000
## San Luis PotosĆ-Quintana Roo -0.55667 -5.47 4.36 1.000
## Sinaloa-Quintana Roo -0.31556 -5.23 4.60 1.000
## Sonora-Quintana Roo -1.07333 -5.99 3.84 1.000
## Tabasco-Quintana Roo -1.26444 -6.18 3.65 1.000
## Tamaulipas-Quintana Roo -3.09778 -8.01 1.82 0.856
## Tlaxcala-Quintana Roo -1.03889 -5.95 3.88 1.000
## Veracruz de Ignacio de la Llave-Quintana Roo -1.05222 -5.97 3.86 1.000
## YucatƔn-Quintana Roo -0.53333 -5.45 4.38 1.000
## Zacatecas-Quintana Roo -0.33889 -5.25 4.58 1.000
## Sinaloa-San Luis PotosĆ 0.24111 -4.67 5.16 1.000
## Sonora-San Luis PotosĆ -0.51667 -5.43 4.40 1.000
## Tabasco-San Luis PotosĆ -0.70778 -5.62 4.21 1.000
## Tamaulipas-San Luis PotosĆ -2.54111 -7.46 2.37 0.985
## Tlaxcala-San Luis PotosĆ -0.48222 -5.40 4.43 1.000
## Veracruz de Ignacio de la Llave-San Luis PotosĆ -0.49556 -5.41 4.42 1.000
## YucatĆ”n-San Luis PotosĆ 0.02333 -4.89 4.94 1.000
## Zacatecas-San Luis PotosĆ 0.21778 -4.70 5.13 1.000
## Sonora-Sinaloa -0.75778 -5.67 4.16 1.000
## Tabasco-Sinaloa -0.94889 -5.86 3.97 1.000
## Tamaulipas-Sinaloa -2.78222 -7.70 2.13 0.952
## Tlaxcala-Sinaloa -0.72333 -5.64 4.19 1.000
## Veracruz de Ignacio de la Llave-Sinaloa -0.73667 -5.65 4.18 1.000
## YucatƔn-Sinaloa -0.21778 -5.13 4.70 1.000
## Zacatecas-Sinaloa -0.02333 -4.94 4.89 1.000
## Tabasco-Sonora -0.19111 -5.11 4.72 1.000
## Tamaulipas-Sonora -2.02444 -6.94 2.89 1.000
## Tlaxcala-Sonora 0.03444 -4.88 4.95 1.000
## Veracruz de Ignacio de la Llave-Sonora 0.02111 -4.89 4.94 1.000
## YucatƔn-Sonora 0.54000 -4.37 5.45 1.000
## Zacatecas-Sonora 0.73444 -4.18 5.65 1.000
## Tamaulipas-Tabasco -1.83333 -6.75 3.08 1.000
## Tlaxcala-Tabasco 0.22556 -4.69 5.14 1.000
## Veracruz de Ignacio de la Llave-Tabasco 0.21222 -4.70 5.13 1.000
## YucatƔn-Tabasco 0.73111 -4.18 5.65 1.000
## Zacatecas-Tabasco 0.92556 -3.99 5.84 1.000
## Tlaxcala-Tamaulipas 2.05889 -2.86 6.97 0.999
## Veracruz de Ignacio de la Llave-Tamaulipas 2.04556 -2.87 6.96 1.000
## YucatƔn-Tamaulipas 2.56444 -2.35 7.48 0.983
## Zacatecas-Tamaulipas 2.75889 -2.16 7.67 0.957
## Veracruz de Ignacio de la Llave-Tlaxcala -0.01333 -4.93 4.90 1.000
## YucatƔn-Tlaxcala 0.50556 -4.41 5.42 1.000
## Zacatecas-Tlaxcala 0.70000 -4.21 5.61 1.000
## YucatƔn-Veracruz de Ignacio de la Llave 0.51889 -4.40 5.43 1.000
## Zacatecas-Veracruz de Ignacio de la Llave 0.71333 -4.20 5.63 1.000
## Zacatecas-YucatƔn 0.19444 -4.72 5.11 1.000residuos <- residuals(ANOVA_bloqueo)
plot(residuos)
# GrƔfica quantile - quantile
qqPlot(residuos, main = "Residuos con bloqueo", xlab = "quantile teĆ³rico", ylab = "quantile real")
## [1] 266 274library(car)
# AsegĆŗrate de que Year sea un factor
gasolina_mx$Year <- as.factor(gasolina_mx$Year)
# ANOVA con bloqueo
ANOVA_bloqueo <- aov(Price ~ Year, data = gasolina_mx)
summary(ANOVA_bloqueo)## Df Sum Sq Mean Sq F value Pr(>F)
## Year 8 1869 233.6 405 <2e-16 ***
## Residuals 279 161 0.6
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1# Prueba de Tukey
if (nlevels(gasolina_mx$Year) > 2) {
print(TukeyHSD(ANOVA_bloqueo))
} else {
message("TukeyHSD no se puede aplicar con menos de 3 niveles en Year.")
}## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Price ~ Year, data = gasolina_mx)
##
## $Year
## diff lwr upr p adj
## 2018-2017 0.777 0.18281 1.3703 0.002
## 2019-2017 3.346 2.75249 3.9400 0.000
## 2020-2017 3.617 3.02343 4.2109 0.000
## 2021-2017 3.090 2.49624 3.6838 0.000
## 2022-2017 4.779 4.18499 5.3725 0.000
## 2023-2017 5.919 5.32531 6.5128 0.000
## 2024-2017 6.512 5.91781 7.1053 0.000
## 2025-2017 8.477 7.88281 9.0703 0.000
## 2019-2018 2.570 1.97593 3.1634 0.000
## 2020-2018 2.841 2.24687 3.4344 0.000
## 2021-2018 2.313 1.71968 2.9072 0.000
## 2022-2018 4.002 3.40843 4.5959 0.000
## 2023-2018 5.142 4.54874 5.7363 0.000
## 2024-2018 5.735 5.14124 6.3288 0.000
## 2025-2018 7.700 7.10624 8.2938 0.000
## 2020-2019 0.271 -0.32282 0.8647 0.887
## 2021-2019 -0.256 -0.85001 0.3375 0.915
## 2022-2019 1.433 0.83874 2.0263 0.000
## 2023-2019 2.573 1.97906 3.1666 0.000
## 2024-2019 3.165 2.57156 3.7591 0.000
## 2025-2019 5.130 4.53656 5.7241 0.000
## 2021-2020 -0.527 -1.12094 0.0666 0.127
## 2022-2020 1.162 0.56781 1.7553 0.000
## 2023-2020 2.302 1.70812 2.8956 0.000
## 2024-2020 2.894 2.30062 3.4881 0.000
## 2025-2020 4.859 4.26562 5.4531 0.000
## 2022-2021 1.689 1.09499 2.2825 0.000
## 2023-2021 2.829 2.23531 3.4228 0.000
## 2024-2021 3.422 2.82781 4.0153 0.000
## 2025-2021 5.387 4.79281 5.9803 0.000
## 2023-2022 1.140 0.54656 1.7341 0.000
## 2024-2022 1.733 1.13906 2.3266 0.000
## 2025-2022 3.698 3.10406 4.2916 0.000
## 2024-2023 0.592 -0.00126 1.1863 0.051
## 2025-2023 2.558 1.96374 3.1513 0.000
## 2025-2024 1.965 1.37124 2.5588 0.000# Residuos
residuos <- residuals(ANOVA_bloqueo)
plot(residuos, main = "Residuos del modelo")
# QQ Plot
qqPlot(residuos, main = "QQ Plot de los residuos", xlab = "Quantil teĆ³rico", ylab = "Quantil real")
## [1] 220 124#ANOVA sin bloqueo
ANOVA_sin_bloqueo <- aov(Price ~ State , data = gasolina_mx)
summary(ANOVA_sin_bloqueo)## Df Sum Sq Mean Sq F value Pr(>F)
## State 31 129 4.15 0.56 0.97
## Residuals 256 1901 7.43# Prueba de Tukey
TukeyHSD(ANOVA_sin_bloqueo)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Price ~ State, data = gasolina_mx)
##
## $State
## diff lwr upr p adj
## Baja California-Aguascalientes -0.97444 -5.89 3.94 1.000
## Baja California Sur-Aguascalientes 0.46111 -4.45 5.38 1.000
## Campeche-Aguascalientes 0.24222 -4.67 5.16 1.000
## Chiapas-Aguascalientes -0.30444 -5.22 4.61 1.000
## Chihuahua-Aguascalientes -1.85778 -6.77 3.06 1.000
## Ciudad de MĆ©xico-Aguascalientes 0.23889 -4.68 5.15 1.000
## Coahuila de Zaragoza-Aguascalientes -0.50000 -5.41 4.41 1.000
## Colima-Aguascalientes -0.08111 -5.00 4.83 1.000
## Durango-Aguascalientes 0.16778 -4.75 5.08 1.000
## Guanajuato-Aguascalientes 0.05444 -4.86 4.97 1.000
## Guerrero-Aguascalientes 0.55111 -4.36 5.47 1.000
## Hidalgo-Aguascalientes -0.49889 -5.41 4.42 1.000
## Jalisco-Aguascalientes 0.35333 -4.56 5.27 1.000
## MĆ©xico-Aguascalientes -0.18778 -5.10 4.73 1.000
## MichoacƔn de Ocampo-Aguascalientes 0.35889 -4.56 5.27 1.000
## Morelos-Aguascalientes -0.29889 -5.21 4.62 1.000
## Nayarit-Aguascalientes 0.47000 -4.44 5.38 1.000
## Nuevo LeĆ³n-Aguascalientes 0.31000 -4.60 5.22 1.000
## Oaxaca-Aguascalientes 0.36778 -4.55 5.28 1.000
## Puebla-Aguascalientes -0.51111 -5.43 4.40 1.000
## QuerƩtaro-Aguascalientes -0.36889 -5.28 4.55 1.000
## Quintana Roo-Aguascalientes 0.48889 -4.43 5.40 1.000
## San Luis PotosĆ-Aguascalientes -0.06778 -4.98 4.85 1.000
## Sinaloa-Aguascalientes 0.17333 -4.74 5.09 1.000
## Sonora-Aguascalientes -0.58444 -5.50 4.33 1.000
## Tabasco-Aguascalientes -0.77556 -5.69 4.14 1.000
## Tamaulipas-Aguascalientes -2.60889 -7.52 2.31 0.978
## Tlaxcala-Aguascalientes -0.55000 -5.46 4.36 1.000
## Veracruz de Ignacio de la Llave-Aguascalientes -0.56333 -5.48 4.35 1.000
## YucatƔn-Aguascalientes -0.04444 -4.96 4.87 1.000
## Zacatecas-Aguascalientes 0.15000 -4.76 5.06 1.000
## Baja California Sur-Baja California 1.43556 -3.48 6.35 1.000
## Campeche-Baja California 1.21667 -3.70 6.13 1.000
## Chiapas-Baja California 0.67000 -4.24 5.58 1.000
## Chihuahua-Baja California -0.88333 -5.80 4.03 1.000
## Ciudad de MĆ©xico-Baja California 1.21333 -3.70 6.13 1.000
## Coahuila de Zaragoza-Baja California 0.47444 -4.44 5.39 1.000
## Colima-Baja California 0.89333 -4.02 5.81 1.000
## Durango-Baja California 1.14222 -3.77 6.06 1.000
## Guanajuato-Baja California 1.02889 -3.89 5.94 1.000
## Guerrero-Baja California 1.52556 -3.39 6.44 1.000
## Hidalgo-Baja California 0.47556 -4.44 5.39 1.000
## Jalisco-Baja California 1.32778 -3.59 6.24 1.000
## MĆ©xico-Baja California 0.78667 -4.13 5.70 1.000
## MichoacƔn de Ocampo-Baja California 1.33333 -3.58 6.25 1.000
## Morelos-Baja California 0.67556 -4.24 5.59 1.000
## Nayarit-Baja California 1.44444 -3.47 6.36 1.000
## Nuevo LeĆ³n-Baja California 1.28444 -3.63 6.20 1.000
## Oaxaca-Baja California 1.34222 -3.57 6.26 1.000
## Puebla-Baja California 0.46333 -4.45 5.38 1.000
## QuerƩtaro-Baja California 0.60556 -4.31 5.52 1.000
## Quintana Roo-Baja California 1.46333 -3.45 6.38 1.000
## San Luis PotosĆ-Baja California 0.90667 -4.01 5.82 1.000
## Sinaloa-Baja California 1.14778 -3.77 6.06 1.000
## Sonora-Baja California 0.39000 -4.52 5.30 1.000
## Tabasco-Baja California 0.19889 -4.72 5.11 1.000
## Tamaulipas-Baja California -1.63444 -6.55 3.28 1.000
## Tlaxcala-Baja California 0.42444 -4.49 5.34 1.000
## Veracruz de Ignacio de la Llave-Baja California 0.41111 -4.50 5.33 1.000
## YucatƔn-Baja California 0.93000 -3.98 5.84 1.000
## Zacatecas-Baja California 1.12444 -3.79 6.04 1.000
## Campeche-Baja California Sur -0.21889 -5.13 4.70 1.000
## Chiapas-Baja California Sur -0.76556 -5.68 4.15 1.000
## Chihuahua-Baja California Sur -2.31889 -7.23 2.60 0.996
## Ciudad de MĆ©xico-Baja California Sur -0.22222 -5.14 4.69 1.000
## Coahuila de Zaragoza-Baja California Sur -0.96111 -5.88 3.95 1.000
## Colima-Baja California Sur -0.54222 -5.46 4.37 1.000
## Durango-Baja California Sur -0.29333 -5.21 4.62 1.000
## Guanajuato-Baja California Sur -0.40667 -5.32 4.51 1.000
## Guerrero-Baja California Sur 0.09000 -4.82 5.00 1.000
## Hidalgo-Baja California Sur -0.96000 -5.87 3.95 1.000
## Jalisco-Baja California Sur -0.10778 -5.02 4.81 1.000
## MĆ©xico-Baja California Sur -0.64889 -5.56 4.27 1.000
## MichoacƔn de Ocampo-Baja California Sur -0.10222 -5.02 4.81 1.000
## Morelos-Baja California Sur -0.76000 -5.67 4.15 1.000
## Nayarit-Baja California Sur 0.00889 -4.91 4.92 1.000
## Nuevo LeĆ³n-Baja California Sur -0.15111 -5.07 4.76 1.000
## Oaxaca-Baja California Sur -0.09333 -5.01 4.82 1.000
## Puebla-Baja California Sur -0.97222 -5.89 3.94 1.000
## QuerƩtaro-Baja California Sur -0.83000 -5.74 4.08 1.000
## Quintana Roo-Baja California Sur 0.02778 -4.89 4.94 1.000
## San Luis PotosĆ-Baja California Sur -0.52889 -5.44 4.39 1.000
## Sinaloa-Baja California Sur -0.28778 -5.20 4.63 1.000
## Sonora-Baja California Sur -1.04556 -5.96 3.87 1.000
## Tabasco-Baja California Sur -1.23667 -6.15 3.68 1.000
## Tamaulipas-Baja California Sur -3.07000 -7.98 1.84 0.867
## Tlaxcala-Baja California Sur -1.01111 -5.93 3.90 1.000
## Veracruz de Ignacio de la Llave-Baja California Sur -1.02444 -5.94 3.89 1.000
## YucatƔn-Baja California Sur -0.50556 -5.42 4.41 1.000
## Zacatecas-Baja California Sur -0.31111 -5.23 4.60 1.000
## Chiapas-Campeche -0.54667 -5.46 4.37 1.000
## Chihuahua-Campeche -2.10000 -7.01 2.81 0.999
## Ciudad de MĆ©xico-Campeche -0.00333 -4.92 4.91 1.000
## Coahuila de Zaragoza-Campeche -0.74222 -5.66 4.17 1.000
## Colima-Campeche -0.32333 -5.24 4.59 1.000
## Durango-Campeche -0.07444 -4.99 4.84 1.000
## Guanajuato-Campeche -0.18778 -5.10 4.73 1.000
## Guerrero-Campeche 0.30889 -4.61 5.22 1.000
## Hidalgo-Campeche -0.74111 -5.66 4.17 1.000
## Jalisco-Campeche 0.11111 -4.80 5.03 1.000
## MĆ©xico-Campeche -0.43000 -5.34 4.48 1.000
## MichoacƔn de Ocampo-Campeche 0.11667 -4.80 5.03 1.000
## Morelos-Campeche -0.54111 -5.46 4.37 1.000
## Nayarit-Campeche 0.22778 -4.69 5.14 1.000
## Nuevo LeĆ³n-Campeche 0.06778 -4.85 4.98 1.000
## Oaxaca-Campeche 0.12556 -4.79 5.04 1.000
## Puebla-Campeche -0.75333 -5.67 4.16 1.000
## QuerƩtaro-Campeche -0.61111 -5.53 4.30 1.000
## Quintana Roo-Campeche 0.24667 -4.67 5.16 1.000
## San Luis PotosĆ-Campeche -0.31000 -5.22 4.60 1.000
## Sinaloa-Campeche -0.06889 -4.98 4.85 1.000
## Sonora-Campeche -0.82667 -5.74 4.09 1.000
## Tabasco-Campeche -1.01778 -5.93 3.90 1.000
## Tamaulipas-Campeche -2.85111 -7.77 2.06 0.937
## Tlaxcala-Campeche -0.79222 -5.71 4.12 1.000
## Veracruz de Ignacio de la Llave-Campeche -0.80556 -5.72 4.11 1.000
## YucatƔn-Campeche -0.28667 -5.20 4.63 1.000
## Zacatecas-Campeche -0.09222 -5.01 4.82 1.000
## Chihuahua-Chiapas -1.55333 -6.47 3.36 1.000
## Ciudad de MĆ©xico-Chiapas 0.54333 -4.37 5.46 1.000
## Coahuila de Zaragoza-Chiapas -0.19556 -5.11 4.72 1.000
## Colima-Chiapas 0.22333 -4.69 5.14 1.000
## Durango-Chiapas 0.47222 -4.44 5.39 1.000
## Guanajuato-Chiapas 0.35889 -4.56 5.27 1.000
## Guerrero-Chiapas 0.85556 -4.06 5.77 1.000
## Hidalgo-Chiapas -0.19444 -5.11 4.72 1.000
## Jalisco-Chiapas 0.65778 -4.26 5.57 1.000
## MĆ©xico-Chiapas 0.11667 -4.80 5.03 1.000
## MichoacƔn de Ocampo-Chiapas 0.66333 -4.25 5.58 1.000
## Morelos-Chiapas 0.00556 -4.91 4.92 1.000
## Nayarit-Chiapas 0.77444 -4.14 5.69 1.000
## Nuevo LeĆ³n-Chiapas 0.61444 -4.30 5.53 1.000
## Oaxaca-Chiapas 0.67222 -4.24 5.59 1.000
## Puebla-Chiapas -0.20667 -5.12 4.71 1.000
## QuerƩtaro-Chiapas -0.06444 -4.98 4.85 1.000
## Quintana Roo-Chiapas 0.79333 -4.12 5.71 1.000
## San Luis PotosĆ-Chiapas 0.23667 -4.68 5.15 1.000
## Sinaloa-Chiapas 0.47778 -4.44 5.39 1.000
## Sonora-Chiapas -0.28000 -5.19 4.63 1.000
## Tabasco-Chiapas -0.47111 -5.39 4.44 1.000
## Tamaulipas-Chiapas -2.30444 -7.22 2.61 0.996
## Tlaxcala-Chiapas -0.24556 -5.16 4.67 1.000
## Veracruz de Ignacio de la Llave-Chiapas -0.25889 -5.17 4.66 1.000
## YucatƔn-Chiapas 0.26000 -4.65 5.17 1.000
## Zacatecas-Chiapas 0.45444 -4.46 5.37 1.000
## Ciudad de MĆ©xico-Chihuahua 2.09667 -2.82 7.01 0.999
## Coahuila de Zaragoza-Chihuahua 1.35778 -3.56 6.27 1.000
## Colima-Chihuahua 1.77667 -3.14 6.69 1.000
## Durango-Chihuahua 2.02556 -2.89 6.94 1.000
## Guanajuato-Chihuahua 1.91222 -3.00 6.83 1.000
## Guerrero-Chihuahua 2.40889 -2.51 7.32 0.993
## Hidalgo-Chihuahua 1.35889 -3.56 6.27 1.000
## Jalisco-Chihuahua 2.21111 -2.70 7.13 0.998
## MĆ©xico-Chihuahua 1.67000 -3.24 6.58 1.000
## MichoacƔn de Ocampo-Chihuahua 2.21667 -2.70 7.13 0.998
## Morelos-Chihuahua 1.55889 -3.36 6.47 1.000
## Nayarit-Chihuahua 2.32778 -2.59 7.24 0.996
## Nuevo LeĆ³n-Chihuahua 2.16778 -2.75 7.08 0.999
## Oaxaca-Chihuahua 2.22556 -2.69 7.14 0.998
## Puebla-Chihuahua 1.34667 -3.57 6.26 1.000
## QuerƩtaro-Chihuahua 1.48889 -3.43 6.40 1.000
## Quintana Roo-Chihuahua 2.34667 -2.57 7.26 0.995
## San Luis PotosĆ-Chihuahua 1.79000 -3.12 6.70 1.000
## Sinaloa-Chihuahua 2.03111 -2.88 6.95 1.000
## Sonora-Chihuahua 1.27333 -3.64 6.19 1.000
## Tabasco-Chihuahua 1.08222 -3.83 6.00 1.000
## Tamaulipas-Chihuahua -0.75111 -5.67 4.16 1.000
## Tlaxcala-Chihuahua 1.30778 -3.61 6.22 1.000
## Veracruz de Ignacio de la Llave-Chihuahua 1.29444 -3.62 6.21 1.000
## YucatƔn-Chihuahua 1.81333 -3.10 6.73 1.000
## Zacatecas-Chihuahua 2.00778 -2.91 6.92 1.000
## Coahuila de Zaragoza-Ciudad de MĆ©xico -0.73889 -5.65 4.18 1.000
## Colima-Ciudad de MĆ©xico -0.32000 -5.23 4.59 1.000
## Durango-Ciudad de MĆ©xico -0.07111 -4.99 4.84 1.000
## Guanajuato-Ciudad de MĆ©xico -0.18444 -5.10 4.73 1.000
## Guerrero-Ciudad de MĆ©xico 0.31222 -4.60 5.23 1.000
## Hidalgo-Ciudad de MĆ©xico -0.73778 -5.65 4.18 1.000
## Jalisco-Ciudad de MĆ©xico 0.11444 -4.80 5.03 1.000
## MĆ©xico-Ciudad de MĆ©xico -0.42667 -5.34 4.49 1.000
## MichoacƔn de Ocampo-Ciudad de MƩxico 0.12000 -4.79 5.03 1.000
## Morelos-Ciudad de MĆ©xico -0.53778 -5.45 4.38 1.000
## Nayarit-Ciudad de MĆ©xico 0.23111 -4.68 5.15 1.000
## Nuevo LeĆ³n-Ciudad de MĆ©xico 0.07111 -4.84 4.99 1.000
## Oaxaca-Ciudad de MĆ©xico 0.12889 -4.79 5.04 1.000
## Puebla-Ciudad de MĆ©xico -0.75000 -5.66 4.16 1.000
## QuerƩtaro-Ciudad de MƩxico -0.60778 -5.52 4.31 1.000
## Quintana Roo-Ciudad de MĆ©xico 0.25000 -4.66 5.16 1.000
## San Luis PotosĆ-Ciudad de MĆ©xico -0.30667 -5.22 4.61 1.000
## Sinaloa-Ciudad de MĆ©xico -0.06556 -4.98 4.85 1.000
## Sonora-Ciudad de MĆ©xico -0.82333 -5.74 4.09 1.000
## Tabasco-Ciudad de MĆ©xico -1.01444 -5.93 3.90 1.000
## Tamaulipas-Ciudad de MĆ©xico -2.84778 -7.76 2.07 0.938
## Tlaxcala-Ciudad de MĆ©xico -0.78889 -5.70 4.13 1.000
## Veracruz de Ignacio de la Llave-Ciudad de MĆ©xico -0.80222 -5.72 4.11 1.000
## YucatƔn-Ciudad de MƩxico -0.28333 -5.20 4.63 1.000
## Zacatecas-Ciudad de MĆ©xico -0.08889 -5.00 4.83 1.000
## Colima-Coahuila de Zaragoza 0.41889 -4.50 5.33 1.000
## Durango-Coahuila de Zaragoza 0.66778 -4.25 5.58 1.000
## Guanajuato-Coahuila de Zaragoza 0.55444 -4.36 5.47 1.000
## Guerrero-Coahuila de Zaragoza 1.05111 -3.86 5.97 1.000
## Hidalgo-Coahuila de Zaragoza 0.00111 -4.91 4.92 1.000
## Jalisco-Coahuila de Zaragoza 0.85333 -4.06 5.77 1.000
## MĆ©xico-Coahuila de Zaragoza 0.31222 -4.60 5.23 1.000
## MichoacƔn de Ocampo-Coahuila de Zaragoza 0.85889 -4.06 5.77 1.000
## Morelos-Coahuila de Zaragoza 0.20111 -4.71 5.12 1.000
## Nayarit-Coahuila de Zaragoza 0.97000 -3.94 5.88 1.000
## Nuevo LeĆ³n-Coahuila de Zaragoza 0.81000 -4.10 5.72 1.000
## Oaxaca-Coahuila de Zaragoza 0.86778 -4.05 5.78 1.000
## Puebla-Coahuila de Zaragoza -0.01111 -4.93 4.90 1.000
## QuerƩtaro-Coahuila de Zaragoza 0.13111 -4.78 5.05 1.000
## Quintana Roo-Coahuila de Zaragoza 0.98889 -3.93 5.90 1.000
## San Luis PotosĆ-Coahuila de Zaragoza 0.43222 -4.48 5.35 1.000
## Sinaloa-Coahuila de Zaragoza 0.67333 -4.24 5.59 1.000
## Sonora-Coahuila de Zaragoza -0.08444 -5.00 4.83 1.000
## Tabasco-Coahuila de Zaragoza -0.27556 -5.19 4.64 1.000
## Tamaulipas-Coahuila de Zaragoza -2.10889 -7.02 2.81 0.999
## Tlaxcala-Coahuila de Zaragoza -0.05000 -4.96 4.86 1.000
## Veracruz de Ignacio de la Llave-Coahuila de Zaragoza -0.06333 -4.98 4.85 1.000
## YucatƔn-Coahuila de Zaragoza 0.45556 -4.46 5.37 1.000
## Zacatecas-Coahuila de Zaragoza 0.65000 -4.26 5.56 1.000
## Durango-Colima 0.24889 -4.67 5.16 1.000
## Guanajuato-Colima 0.13556 -4.78 5.05 1.000
## Guerrero-Colima 0.63222 -4.28 5.55 1.000
## Hidalgo-Colima -0.41778 -5.33 4.50 1.000
## Jalisco-Colima 0.43444 -4.48 5.35 1.000
## MĆ©xico-Colima -0.10667 -5.02 4.81 1.000
## MichoacƔn de Ocampo-Colima 0.44000 -4.47 5.35 1.000
## Morelos-Colima -0.21778 -5.13 4.70 1.000
## Nayarit-Colima 0.55111 -4.36 5.47 1.000
## Nuevo LeĆ³n-Colima 0.39111 -4.52 5.31 1.000
## Oaxaca-Colima 0.44889 -4.47 5.36 1.000
## Puebla-Colima -0.43000 -5.34 4.48 1.000
## QuerƩtaro-Colima -0.28778 -5.20 4.63 1.000
## Quintana Roo-Colima 0.57000 -4.34 5.48 1.000
## San Luis PotosĆ-Colima 0.01333 -4.90 4.93 1.000
## Sinaloa-Colima 0.25444 -4.66 5.17 1.000
## Sonora-Colima -0.50333 -5.42 4.41 1.000
## Tabasco-Colima -0.69444 -5.61 4.22 1.000
## Tamaulipas-Colima -2.52778 -7.44 2.39 0.986
## Tlaxcala-Colima -0.46889 -5.38 4.45 1.000
## Veracruz de Ignacio de la Llave-Colima -0.48222 -5.40 4.43 1.000
## YucatƔn-Colima 0.03667 -4.88 4.95 1.000
## Zacatecas-Colima 0.23111 -4.68 5.15 1.000
## Guanajuato-Durango -0.11333 -5.03 4.80 1.000
## Guerrero-Durango 0.38333 -4.53 5.30 1.000
## Hidalgo-Durango -0.66667 -5.58 4.25 1.000
## Jalisco-Durango 0.18556 -4.73 5.10 1.000
## MĆ©xico-Durango -0.35556 -5.27 4.56 1.000
## MichoacƔn de Ocampo-Durango 0.19111 -4.72 5.11 1.000
## Morelos-Durango -0.46667 -5.38 4.45 1.000
## Nayarit-Durango 0.30222 -4.61 5.22 1.000
## Nuevo LeĆ³n-Durango 0.14222 -4.77 5.06 1.000
## Oaxaca-Durango 0.20000 -4.71 5.11 1.000
## Puebla-Durango -0.67889 -5.59 4.24 1.000
## QuerƩtaro-Durango -0.53667 -5.45 4.38 1.000
## Quintana Roo-Durango 0.32111 -4.59 5.24 1.000
## San Luis PotosĆ-Durango -0.23556 -5.15 4.68 1.000
## Sinaloa-Durango 0.00556 -4.91 4.92 1.000
## Sonora-Durango -0.75222 -5.67 4.16 1.000
## Tabasco-Durango -0.94333 -5.86 3.97 1.000
## Tamaulipas-Durango -2.77667 -7.69 2.14 0.953
## Tlaxcala-Durango -0.71778 -5.63 4.20 1.000
## Veracruz de Ignacio de la Llave-Durango -0.73111 -5.65 4.18 1.000
## YucatƔn-Durango -0.21222 -5.13 4.70 1.000
## Zacatecas-Durango -0.01778 -4.93 4.90 1.000
## Guerrero-Guanajuato 0.49667 -4.42 5.41 1.000
## Hidalgo-Guanajuato -0.55333 -5.47 4.36 1.000
## Jalisco-Guanajuato 0.29889 -4.62 5.21 1.000
## MĆ©xico-Guanajuato -0.24222 -5.16 4.67 1.000
## MichoacƔn de Ocampo-Guanajuato 0.30444 -4.61 5.22 1.000
## Morelos-Guanajuato -0.35333 -5.27 4.56 1.000
## Nayarit-Guanajuato 0.41556 -4.50 5.33 1.000
## Nuevo LeĆ³n-Guanajuato 0.25556 -4.66 5.17 1.000
## Oaxaca-Guanajuato 0.31333 -4.60 5.23 1.000
## Puebla-Guanajuato -0.56556 -5.48 4.35 1.000
## QuerƩtaro-Guanajuato -0.42333 -5.34 4.49 1.000
## Quintana Roo-Guanajuato 0.43444 -4.48 5.35 1.000
## San Luis PotosĆ-Guanajuato -0.12222 -5.04 4.79 1.000
## Sinaloa-Guanajuato 0.11889 -4.80 5.03 1.000
## Sonora-Guanajuato -0.63889 -5.55 4.28 1.000
## Tabasco-Guanajuato -0.83000 -5.74 4.08 1.000
## Tamaulipas-Guanajuato -2.66333 -7.58 2.25 0.972
## Tlaxcala-Guanajuato -0.60444 -5.52 4.31 1.000
## Veracruz de Ignacio de la Llave-Guanajuato -0.61778 -5.53 4.30 1.000
## YucatƔn-Guanajuato -0.09889 -5.01 4.82 1.000
## Zacatecas-Guanajuato 0.09556 -4.82 5.01 1.000
## Hidalgo-Guerrero -1.05000 -5.96 3.86 1.000
## Jalisco-Guerrero -0.19778 -5.11 4.72 1.000
## MĆ©xico-Guerrero -0.73889 -5.65 4.18 1.000
## MichoacƔn de Ocampo-Guerrero -0.19222 -5.11 4.72 1.000
## Morelos-Guerrero -0.85000 -5.76 4.06 1.000
## Nayarit-Guerrero -0.08111 -5.00 4.83 1.000
## Nuevo LeĆ³n-Guerrero -0.24111 -5.16 4.67 1.000
## Oaxaca-Guerrero -0.18333 -5.10 4.73 1.000
## Puebla-Guerrero -1.06222 -5.98 3.85 1.000
## QuerƩtaro-Guerrero -0.92000 -5.83 3.99 1.000
## Quintana Roo-Guerrero -0.06222 -4.98 4.85 1.000
## San Luis PotosĆ-Guerrero -0.61889 -5.53 4.30 1.000
## Sinaloa-Guerrero -0.37778 -5.29 4.54 1.000
## Sonora-Guerrero -1.13556 -6.05 3.78 1.000
## Tabasco-Guerrero -1.32667 -6.24 3.59 1.000
## Tamaulipas-Guerrero -3.16000 -8.07 1.75 0.828
## Tlaxcala-Guerrero -1.10111 -6.02 3.81 1.000
## Veracruz de Ignacio de la Llave-Guerrero -1.11444 -6.03 3.80 1.000
## YucatƔn-Guerrero -0.59556 -5.51 4.32 1.000
## Zacatecas-Guerrero -0.40111 -5.32 4.51 1.000
## Jalisco-Hidalgo 0.85222 -4.06 5.77 1.000
## MĆ©xico-Hidalgo 0.31111 -4.60 5.23 1.000
## MichoacƔn de Ocampo-Hidalgo 0.85778 -4.06 5.77 1.000
## Morelos-Hidalgo 0.20000 -4.71 5.11 1.000
## Nayarit-Hidalgo 0.96889 -3.95 5.88 1.000
## Nuevo LeĆ³n-Hidalgo 0.80889 -4.11 5.72 1.000
## Oaxaca-Hidalgo 0.86667 -4.05 5.78 1.000
## Puebla-Hidalgo -0.01222 -4.93 4.90 1.000
## QuerƩtaro-Hidalgo 0.13000 -4.78 5.04 1.000
## Quintana Roo-Hidalgo 0.98778 -3.93 5.90 1.000
## San Luis PotosĆ-Hidalgo 0.43111 -4.48 5.35 1.000
## Sinaloa-Hidalgo 0.67222 -4.24 5.59 1.000
## Sonora-Hidalgo -0.08556 -5.00 4.83 1.000
## Tabasco-Hidalgo -0.27667 -5.19 4.64 1.000
## Tamaulipas-Hidalgo -2.11000 -7.02 2.80 0.999
## Tlaxcala-Hidalgo -0.05111 -4.97 4.86 1.000
## Veracruz de Ignacio de la Llave-Hidalgo -0.06444 -4.98 4.85 1.000
## YucatƔn-Hidalgo 0.45444 -4.46 5.37 1.000
## Zacatecas-Hidalgo 0.64889 -4.27 5.56 1.000
## MĆ©xico-Jalisco -0.54111 -5.46 4.37 1.000
## MichoacƔn de Ocampo-Jalisco 0.00556 -4.91 4.92 1.000
## Morelos-Jalisco -0.65222 -5.57 4.26 1.000
## Nayarit-Jalisco 0.11667 -4.80 5.03 1.000
## Nuevo LeĆ³n-Jalisco -0.04333 -4.96 4.87 1.000
## Oaxaca-Jalisco 0.01444 -4.90 4.93 1.000
## Puebla-Jalisco -0.86444 -5.78 4.05 1.000
## QuerƩtaro-Jalisco -0.72222 -5.64 4.19 1.000
## Quintana Roo-Jalisco 0.13556 -4.78 5.05 1.000
## San Luis PotosĆ-Jalisco -0.42111 -5.34 4.49 1.000
## Sinaloa-Jalisco -0.18000 -5.09 4.73 1.000
## Sonora-Jalisco -0.93778 -5.85 3.98 1.000
## Tabasco-Jalisco -1.12889 -6.04 3.79 1.000
## Tamaulipas-Jalisco -2.96222 -7.88 1.95 0.906
## Tlaxcala-Jalisco -0.90333 -5.82 4.01 1.000
## Veracruz de Ignacio de la Llave-Jalisco -0.91667 -5.83 4.00 1.000
## YucatƔn-Jalisco -0.39778 -5.31 4.52 1.000
## Zacatecas-Jalisco -0.20333 -5.12 4.71 1.000
## MichoacƔn de Ocampo-MƩxico 0.54667 -4.37 5.46 1.000
## Morelos-MĆ©xico -0.11111 -5.03 4.80 1.000
## Nayarit-MĆ©xico 0.65778 -4.26 5.57 1.000
## Nuevo LeĆ³n-MĆ©xico 0.49778 -4.42 5.41 1.000
## Oaxaca-MĆ©xico 0.55556 -4.36 5.47 1.000
## Puebla-MĆ©xico -0.32333 -5.24 4.59 1.000
## QuerƩtaro-MƩxico -0.18111 -5.10 4.73 1.000
## Quintana Roo-MĆ©xico 0.67667 -4.24 5.59 1.000
## San Luis PotosĆ-MĆ©xico 0.12000 -4.79 5.03 1.000
## Sinaloa-MĆ©xico 0.36111 -4.55 5.28 1.000
## Sonora-MĆ©xico -0.39667 -5.31 4.52 1.000
## Tabasco-MĆ©xico -0.58778 -5.50 4.33 1.000
## Tamaulipas-MĆ©xico -2.42111 -7.34 2.49 0.992
## Tlaxcala-MĆ©xico -0.36222 -5.28 4.55 1.000
## Veracruz de Ignacio de la Llave-MĆ©xico -0.37556 -5.29 4.54 1.000
## YucatƔn-MƩxico 0.14333 -4.77 5.06 1.000
## Zacatecas-MĆ©xico 0.33778 -4.58 5.25 1.000
## Morelos-MichoacƔn de Ocampo -0.65778 -5.57 4.26 1.000
## Nayarit-MichoacƔn de Ocampo 0.11111 -4.80 5.03 1.000
## Nuevo LeĆ³n-MichoacĆ”n de Ocampo -0.04889 -4.96 4.87 1.000
## Oaxaca-MichoacƔn de Ocampo 0.00889 -4.91 4.92 1.000
## Puebla-MichoacƔn de Ocampo -0.87000 -5.78 4.04 1.000
## QuerƩtaro-MichoacƔn de Ocampo -0.72778 -5.64 4.19 1.000
## Quintana Roo-MichoacƔn de Ocampo 0.13000 -4.78 5.04 1.000
## San Luis PotosĆ-MichoacĆ”n de Ocampo -0.42667 -5.34 4.49 1.000
## Sinaloa-MichoacƔn de Ocampo -0.18556 -5.10 4.73 1.000
## Sonora-MichoacƔn de Ocampo -0.94333 -5.86 3.97 1.000
## Tabasco-MichoacƔn de Ocampo -1.13444 -6.05 3.78 1.000
## Tamaulipas-MichoacƔn de Ocampo -2.96778 -7.88 1.95 0.904
## Tlaxcala-MichoacƔn de Ocampo -0.90889 -5.82 4.01 1.000
## Veracruz de Ignacio de la Llave-MichoacƔn de Ocampo -0.92222 -5.84 3.99 1.000
## YucatƔn-MichoacƔn de Ocampo -0.40333 -5.32 4.51 1.000
## Zacatecas-MichoacƔn de Ocampo -0.20889 -5.12 4.71 1.000
## Nayarit-Morelos 0.76889 -4.15 5.68 1.000
## Nuevo LeĆ³n-Morelos 0.60889 -4.31 5.52 1.000
## Oaxaca-Morelos 0.66667 -4.25 5.58 1.000
## Puebla-Morelos -0.21222 -5.13 4.70 1.000
## QuerƩtaro-Morelos -0.07000 -4.98 4.84 1.000
## Quintana Roo-Morelos 0.78778 -4.13 5.70 1.000
## San Luis PotosĆ-Morelos 0.23111 -4.68 5.15 1.000
## Sinaloa-Morelos 0.47222 -4.44 5.39 1.000
## Sonora-Morelos -0.28556 -5.20 4.63 1.000
## Tabasco-Morelos -0.47667 -5.39 4.44 1.000
## Tamaulipas-Morelos -2.31000 -7.22 2.60 0.996
## Tlaxcala-Morelos -0.25111 -5.17 4.66 1.000
## Veracruz de Ignacio de la Llave-Morelos -0.26444 -5.18 4.65 1.000
## YucatƔn-Morelos 0.25444 -4.66 5.17 1.000
## Zacatecas-Morelos 0.44889 -4.47 5.36 1.000
## Nuevo LeĆ³n-Nayarit -0.16000 -5.07 4.75 1.000
## Oaxaca-Nayarit -0.10222 -5.02 4.81 1.000
## Puebla-Nayarit -0.98111 -5.90 3.93 1.000
## QuerƩtaro-Nayarit -0.83889 -5.75 4.08 1.000
## Quintana Roo-Nayarit 0.01889 -4.90 4.93 1.000
## San Luis PotosĆ-Nayarit -0.53778 -5.45 4.38 1.000
## Sinaloa-Nayarit -0.29667 -5.21 4.62 1.000
## Sonora-Nayarit -1.05444 -5.97 3.86 1.000
## Tabasco-Nayarit -1.24556 -6.16 3.67 1.000
## Tamaulipas-Nayarit -3.07889 -7.99 1.84 0.863
## Tlaxcala-Nayarit -1.02000 -5.93 3.89 1.000
## Veracruz de Ignacio de la Llave-Nayarit -1.03333 -5.95 3.88 1.000
## YucatƔn-Nayarit -0.51444 -5.43 4.40 1.000
## Zacatecas-Nayarit -0.32000 -5.23 4.59 1.000
## Oaxaca-Nuevo LeĆ³n 0.05778 -4.86 4.97 1.000
## Puebla-Nuevo LeĆ³n -0.82111 -5.74 4.09 1.000
## QuerĆ©taro-Nuevo LeĆ³n -0.67889 -5.59 4.24 1.000
## Quintana Roo-Nuevo LeĆ³n 0.17889 -4.74 5.09 1.000
## San Luis PotosĆ-Nuevo LeĆ³n -0.37778 -5.29 4.54 1.000
## Sinaloa-Nuevo LeĆ³n -0.13667 -5.05 4.78 1.000
## Sonora-Nuevo LeĆ³n -0.89444 -5.81 4.02 1.000
## Tabasco-Nuevo LeĆ³n -1.08556 -6.00 3.83 1.000
## Tamaulipas-Nuevo LeĆ³n -2.91889 -7.83 2.00 0.919
## Tlaxcala-Nuevo LeĆ³n -0.86000 -5.77 4.05 1.000
## Veracruz de Ignacio de la Llave-Nuevo LeĆ³n -0.87333 -5.79 4.04 1.000
## YucatĆ”n-Nuevo LeĆ³n -0.35444 -5.27 4.56 1.000
## Zacatecas-Nuevo LeĆ³n -0.16000 -5.07 4.75 1.000
## Puebla-Oaxaca -0.87889 -5.79 4.04 1.000
## QuerƩtaro-Oaxaca -0.73667 -5.65 4.18 1.000
## Quintana Roo-Oaxaca 0.12111 -4.79 5.04 1.000
## San Luis PotosĆ-Oaxaca -0.43556 -5.35 4.48 1.000
## Sinaloa-Oaxaca -0.19444 -5.11 4.72 1.000
## Sonora-Oaxaca -0.95222 -5.87 3.96 1.000
## Tabasco-Oaxaca -1.14333 -6.06 3.77 1.000
## Tamaulipas-Oaxaca -2.97667 -7.89 1.94 0.901
## Tlaxcala-Oaxaca -0.91778 -5.83 4.00 1.000
## Veracruz de Ignacio de la Llave-Oaxaca -0.93111 -5.85 3.98 1.000
## YucatƔn-Oaxaca -0.41222 -5.33 4.50 1.000
## Zacatecas-Oaxaca -0.21778 -5.13 4.70 1.000
## QuerƩtaro-Puebla 0.14222 -4.77 5.06 1.000
## Quintana Roo-Puebla 1.00000 -3.91 5.91 1.000
## San Luis PotosĆ-Puebla 0.44333 -4.47 5.36 1.000
## Sinaloa-Puebla 0.68444 -4.23 5.60 1.000
## Sonora-Puebla -0.07333 -4.99 4.84 1.000
## Tabasco-Puebla -0.26444 -5.18 4.65 1.000
## Tamaulipas-Puebla -2.09778 -7.01 2.82 0.999
## Tlaxcala-Puebla -0.03889 -4.95 4.88 1.000
## Veracruz de Ignacio de la Llave-Puebla -0.05222 -4.97 4.86 1.000
## YucatƔn-Puebla 0.46667 -4.45 5.38 1.000
## Zacatecas-Puebla 0.66111 -4.25 5.58 1.000
## Quintana Roo-QuerƩtaro 0.85778 -4.06 5.77 1.000
## San Luis PotosĆ-QuerĆ©taro 0.30111 -4.61 5.22 1.000
## Sinaloa-QuerƩtaro 0.54222 -4.37 5.46 1.000
## Sonora-QuerƩtaro -0.21556 -5.13 4.70 1.000
## Tabasco-QuerƩtaro -0.40667 -5.32 4.51 1.000
## Tamaulipas-QuerƩtaro -2.24000 -7.15 2.67 0.998
## Tlaxcala-QuerƩtaro -0.18111 -5.10 4.73 1.000
## Veracruz de Ignacio de la Llave-QuerƩtaro -0.19444 -5.11 4.72 1.000
## YucatƔn-QuerƩtaro 0.32444 -4.59 5.24 1.000
## Zacatecas-QuerƩtaro 0.51889 -4.40 5.43 1.000
## San Luis PotosĆ-Quintana Roo -0.55667 -5.47 4.36 1.000
## Sinaloa-Quintana Roo -0.31556 -5.23 4.60 1.000
## Sonora-Quintana Roo -1.07333 -5.99 3.84 1.000
## Tabasco-Quintana Roo -1.26444 -6.18 3.65 1.000
## Tamaulipas-Quintana Roo -3.09778 -8.01 1.82 0.856
## Tlaxcala-Quintana Roo -1.03889 -5.95 3.88 1.000
## Veracruz de Ignacio de la Llave-Quintana Roo -1.05222 -5.97 3.86 1.000
## YucatƔn-Quintana Roo -0.53333 -5.45 4.38 1.000
## Zacatecas-Quintana Roo -0.33889 -5.25 4.58 1.000
## Sinaloa-San Luis PotosĆ 0.24111 -4.67 5.16 1.000
## Sonora-San Luis PotosĆ -0.51667 -5.43 4.40 1.000
## Tabasco-San Luis PotosĆ -0.70778 -5.62 4.21 1.000
## Tamaulipas-San Luis PotosĆ -2.54111 -7.46 2.37 0.985
## Tlaxcala-San Luis PotosĆ -0.48222 -5.40 4.43 1.000
## Veracruz de Ignacio de la Llave-San Luis PotosĆ -0.49556 -5.41 4.42 1.000
## YucatĆ”n-San Luis PotosĆ 0.02333 -4.89 4.94 1.000
## Zacatecas-San Luis PotosĆ 0.21778 -4.70 5.13 1.000
## Sonora-Sinaloa -0.75778 -5.67 4.16 1.000
## Tabasco-Sinaloa -0.94889 -5.86 3.97 1.000
## Tamaulipas-Sinaloa -2.78222 -7.70 2.13 0.952
## Tlaxcala-Sinaloa -0.72333 -5.64 4.19 1.000
## Veracruz de Ignacio de la Llave-Sinaloa -0.73667 -5.65 4.18 1.000
## YucatƔn-Sinaloa -0.21778 -5.13 4.70 1.000
## Zacatecas-Sinaloa -0.02333 -4.94 4.89 1.000
## Tabasco-Sonora -0.19111 -5.11 4.72 1.000
## Tamaulipas-Sonora -2.02444 -6.94 2.89 1.000
## Tlaxcala-Sonora 0.03444 -4.88 4.95 1.000
## Veracruz de Ignacio de la Llave-Sonora 0.02111 -4.89 4.94 1.000
## YucatƔn-Sonora 0.54000 -4.37 5.45 1.000
## Zacatecas-Sonora 0.73444 -4.18 5.65 1.000
## Tamaulipas-Tabasco -1.83333 -6.75 3.08 1.000
## Tlaxcala-Tabasco 0.22556 -4.69 5.14 1.000
## Veracruz de Ignacio de la Llave-Tabasco 0.21222 -4.70 5.13 1.000
## YucatƔn-Tabasco 0.73111 -4.18 5.65 1.000
## Zacatecas-Tabasco 0.92556 -3.99 5.84 1.000
## Tlaxcala-Tamaulipas 2.05889 -2.86 6.97 0.999
## Veracruz de Ignacio de la Llave-Tamaulipas 2.04556 -2.87 6.96 1.000
## YucatƔn-Tamaulipas 2.56444 -2.35 7.48 0.983
## Zacatecas-Tamaulipas 2.75889 -2.16 7.67 0.957
## Veracruz de Ignacio de la Llave-Tlaxcala -0.01333 -4.93 4.90 1.000
## YucatƔn-Tlaxcala 0.50556 -4.41 5.42 1.000
## Zacatecas-Tlaxcala 0.70000 -4.21 5.61 1.000
## YucatƔn-Veracruz de Ignacio de la Llave 0.51889 -4.40 5.43 1.000
## Zacatecas-Veracruz de Ignacio de la Llave 0.71333 -4.20 5.63 1.000
## Zacatecas-YucatƔn 0.19444 -4.72 5.11 1.000residuos <- residuals(ANOVA_sin_bloqueo)
plot(residuos)
# GrƔfica quantile - quantile
qqPlot(residuos, main = "Residuos sin bloqueo", xlab = "quantile teĆ³rico", ylab = "quantile real")
## [1] 266 274library(car)
# AsegĆŗrate de que Year sea un factor
gasolina_mx$Year <- as.factor(gasolina_mx$Year)
# ANOVA sin bloqueo
ANOVA_sin_bloqueo <- aov(Price ~ Year, data = gasolina_mx)
summary(ANOVA_sin_bloqueo)## Df Sum Sq Mean Sq F value Pr(>F)
## Year 8 1869 233.6 405 <2e-16 ***
## Residuals 279 161 0.6
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1# Prueba de Tukey
if (nlevels(gasolina_mx$Year) > 2) {
print(TukeyHSD(ANOVA_sin_bloqueo))
} else {
message("TukeyHSD no se puede aplicar con menos de 3 niveles en Year.")
}## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = Price ~ Year, data = gasolina_mx)
##
## $Year
## diff lwr upr p adj
## 2018-2017 0.777 0.18281 1.3703 0.002
## 2019-2017 3.346 2.75249 3.9400 0.000
## 2020-2017 3.617 3.02343 4.2109 0.000
## 2021-2017 3.090 2.49624 3.6838 0.000
## 2022-2017 4.779 4.18499 5.3725 0.000
## 2023-2017 5.919 5.32531 6.5128 0.000
## 2024-2017 6.512 5.91781 7.1053 0.000
## 2025-2017 8.477 7.88281 9.0703 0.000
## 2019-2018 2.570 1.97593 3.1634 0.000
## 2020-2018 2.841 2.24687 3.4344 0.000
## 2021-2018 2.313 1.71968 2.9072 0.000
## 2022-2018 4.002 3.40843 4.5959 0.000
## 2023-2018 5.142 4.54874 5.7363 0.000
## 2024-2018 5.735 5.14124 6.3288 0.000
## 2025-2018 7.700 7.10624 8.2938 0.000
## 2020-2019 0.271 -0.32282 0.8647 0.887
## 2021-2019 -0.256 -0.85001 0.3375 0.915
## 2022-2019 1.433 0.83874 2.0263 0.000
## 2023-2019 2.573 1.97906 3.1666 0.000
## 2024-2019 3.165 2.57156 3.7591 0.000
## 2025-2019 5.130 4.53656 5.7241 0.000
## 2021-2020 -0.527 -1.12094 0.0666 0.127
## 2022-2020 1.162 0.56781 1.7553 0.000
## 2023-2020 2.302 1.70812 2.8956 0.000
## 2024-2020 2.894 2.30062 3.4881 0.000
## 2025-2020 4.859 4.26562 5.4531 0.000
## 2022-2021 1.689 1.09499 2.2825 0.000
## 2023-2021 2.829 2.23531 3.4228 0.000
## 2024-2021 3.422 2.82781 4.0153 0.000
## 2025-2021 5.387 4.79281 5.9803 0.000
## 2023-2022 1.140 0.54656 1.7341 0.000
## 2024-2022 1.733 1.13906 2.3266 0.000
## 2025-2022 3.698 3.10406 4.2916 0.000
## 2024-2023 0.592 -0.00126 1.1863 0.051
## 2025-2023 2.558 1.96374 3.1513 0.000
## 2025-2024 1.965 1.37124 2.5588 0.000# Residuos
residuos <- residuals(ANOVA_sin_bloqueo)
plot(residuos, main = "Residuos del modelo")
# QQ Plot
qqPlot(residuos, main = "QQ Plot de los residuos", xlab = "Quantil teĆ³rico", ylab = "Quantil real")
## [1] 220 124conclusion
dados los resultados del anova se puede concluir que hay evidencia estadistica que demuestre que existe cambio entre los precios de la gasolina en los aƱos por region, ya que el p- value= 0-97, mayor a a=0.05