library(readxl)
insurance <- read_excel("C:/Users/gisel/Downloads/insurance.xlsx")
View(insurance)
# install.packages("readxl")
# install.packages("tseries")
library(readxl)
library(tseries)
## Warning: package 'tseries' was built under R version 4.5.3
## Registered S3 method overwritten by 'quantmod':
## method from
## as.zoo.data.frame zoo
modelo_seguro <- lm(charges ~ age + bmi + children + smoker + region, data = insurance)
residuos <- residuals(modelo_seguro)
# A) Prueba de Shapiro-Wilk
shapiro_test <- shapiro.test(residuos)
print(shapiro_test)
##
## Shapiro-Wilk normality test
##
## data: residuos
## W = 0.89909, p-value < 2.2e-16
# B) Prueba de Jarque-Bera
jb_test <- jarque.bera.test(residuos)
print(jb_test)
##
## Jarque Bera Test
##
## data: residuos
## X-squared = 720.52, df = 2, p-value < 2.2e-16
# C) Prueba de Kolmogorov-Smirnov con residuos estandarizados
residuos_estandarizados <- as.vector(scale(residuos))
ks_test <- ks.test(residuos_estandarizados, "pnorm")
## Warning in ks.test.default(residuos_estandarizados, "pnorm"): ties should not
## be present for the one-sample Kolmogorov-Smirnov test
print(ks_test)
##
## Asymptotic one-sample Kolmogorov-Smirnov test
##
## data: residuos_estandarizados
## D = 0.16124, p-value < 2.2e-16
## alternative hypothesis: two-sided
# REPRESENTACIÓN GRÁFICA (Gráfico Q-Q)
qqnorm(residuos,
main = "Gráfico Q-Q de Residuos (Evidencia de No Normalidad)",
xlab = "Cuantiles Teóricos de una Normal",
ylab = "Cuantiles Empíricos de los Residuos",
col = rgb(0.2, 0.6, 0.8, 0.6),
pch = 16)
Conclusion: El comportamiento de los residuos es idéntico: todas las pruebas estadísticas arrojarán un p-valor menor a 0.05, confirmando el rechazo de la hipótesis nula de normalidad. Además, el gráfico Q-Q mostrará una desviación significativa de los puntos respecto a la línea diagonal, evidenciando la falta de normalidad en los residuos del modelo.