Pruebas de Normalidad para los residuos, Guia practica 3
Gabriela Lisbeth Palencia Flores PF13012 Econometria GT 02
26 de abril de 2019
CARGA DE DATOS
## Warning: package 'dplyr' was built under R version 3.5.3
## Warning: package 'readr' was built under R version 3.5.3
library(stargazer)
ejemplo_regresion <- read.csv("F:/ejercicioregresion.csv")
head(ejemplo_regresion, n=6)
## X1 X2 Y
## 1 3.92 7298 0.75
## 2 3.61 6855 0.71
## 3 3.32 6636 0.66
## 4 3.07 6506 0.61
## 5 3.06 6450 0.70
## 6 3.11 6402 0.72
ESTIMANDO EL MODELO:
library(stargazer)
options(scipen = 9999)
modelo_lineal <- lm(formula = Y~X1+X2, data = ejemplo_regresion)
stargazer(modelo_lineal, title = "ejemplo de Regresion Lineal Multiple", type = "text", digits = 8)
##
## ejemplo de Regresion Lineal Multiple
## ===============================================
## Dependent variable:
## ---------------------------
## Y
## -----------------------------------------------
## X1 0.23719750***
## (0.05555937)
##
## X2 -0.00024908***
## (0.00003205)
##
## Constant 1.56449700***
## (0.07939598)
##
## -----------------------------------------------
## Observations 25
## R2 0.86529610
## Adjusted R2 0.85305030
## Residual Std. Error 0.05330222 (df = 22)
## F Statistic 70.66057000*** (df = 2; 22)
## ===============================================
## Note: *p<0.1; **p<0.05; ***p<0.01
AJUSTE DE LOS RESIDUOS DE LA DISTRIBUCION NORMAL
## Warning: package 'fitdistrplus' was built under R version 3.5.3
library(stargazer)
fit_normal<- fitdist(data = modelo_lineal$residuals, distr = "norm")
plot(fit_normal)
## Fitting of the distribution ' norm ' by maximum likelihood
## Parameters :
## estimate Std. Error
## mean 0.000000000000000007770748 0.010000382
## sd 0.050001911895951975384200 0.007058615
## Loglikelihood: 39.41889 AIC: -74.83778 BIC: -72.40002
## Correlation matrix:
## mean sd
## mean 1 0
## sd 0 1
PRUEBA DE NORMALIDAD DE JARQUE-BERA
library(normtest) #carga los comandos para las pruebas de normalidad
jb.norm.test(modelo_lineal$residuals) #ejecuta la prueba de Jarque Bera
##
## Jarque-Bera test for normality
##
## data: modelo_lineal$residuals
## JB = 0.93032, p-value = 0.48
## Al realizar la prueba de JARQUE-BERA dio como resultado el valor p=0.49 siendo este mayor a alfa(0.05), por lo que se concluye que los residuos provienen de una distribucion normal, no se rechaza Ho.
qqnorm(modelo_lineal$residuals)
qqline(modelo_lineal$residuals)
hist(modelo_lineal$residuals, main = "Histograma de residuos", xlab = "Residuos", ylab = "Frecuencia")
PRUEBA DE NORMALIDAD DE KOLMOGOROV-SMIRNOV
library(nortest) #carga los comandos para las pruebas de normalidad
lillie.test(modelo_lineal$residuals) #ejecuta la prueba KS, con la corrección de Lilliefors
##
## Lilliefors (Kolmogorov-Smirnov) normality test
##
## data: modelo_lineal$residuals
## D = 0.082345, p-value = 0.9328
## Al realizar la prueba de KOLMOGOROV-SMIRNOV dio como resultado el valor p=0.9328 siendo este mayor a alfa(0.05), por lo que se concluye que los residuos provienen de una distribucion normal, no se rechaza Ho.
qqnorm(modelo_lineal$residuals)
qqline(modelo_lineal$residuals)
hist(modelo_lineal$residuals, main = "Histograma de Residuos", xlab = "Residuos", ylab = "Frecuencia")
PRUEBA DE NORMALIDAD DE SHAPIRO-WILK
shapiro.test(modelo_lineal$residuals)
##
## Shapiro-Wilk normality test
##
## data: modelo_lineal$residuals
## W = 0.97001, p-value = 0.6453
## Al realizar la prueba de SHAPIRO-WILK dio como resultado el valor p=0.6453 siendo este mayor a alfa(0.05), por lo que se concluye que los residuos provienen de una distribucion normal, no se rechaza Ho.
CARGA DE DATOS
library (dplyr)
library (readr)
library(stargazer)
ejemplo_regresion2<- read.csv("F:/ejemplo_regresion2.csv")
head(ejemplo_regresion2, n=6)
## X1 X2 Y
## 1 50 7.4 320
## 2 53 5.1 450
## 3 60 4.2 370
## 4 63 3.9 470
## 5 69 1.4 420
## 6 82 2.2 500
ESTIMANDO EL MODELO:
library(stargazer)
options(scipen = 9999)
modelo_lineal2<- lm(formula = Y~X1+X2, data = ejemplo_regresion2)
stargazer(modelo_lineal2, title = "ejemplo de Regresion Lineal Multiple", type = "text", digits = 8)
##
## ejemplo de Regresion Lineal Multiple
## ================================================
## Dependent variable:
## ----------------------------
## Y
## ------------------------------------------------
## X1 3.94065500***
## (0.26658770)
##
## X2 14.80694000**
## (5.51843400)
##
## Constant 85.57848000
## (54.61834000)
##
## ------------------------------------------------
## Observations 20
## R2 0.93167600
## Adjusted R2 0.92363790
## Residual Std. Error 89.72397000 (df = 17)
## F Statistic 115.90730000*** (df = 2; 17)
## ================================================
## Note: *p<0.1; **p<0.05; ***p<0.01
AJUSTE DE LOS RESIDUOS DE LA DISTRIBUCION NORMAL
library(fitdistrplus)
library(stargazer)
fit_normal2<- fitdist(data = modelo_lineal2$residuals, distr = "norm")
plot(fit_normal2)
## Fitting of the distribution ' norm ' by maximum likelihood
## Parameters :
## estimate Std. Error
## mean -0.000000000000002488114 18.49709
## sd 82.721409061040219512506 13.07940
## Loglikelihood: -116.6883 AIC: 237.3767 BIC: 239.3681
## Correlation matrix:
## mean sd
## mean 1.000000000000 -0.000001289266
## sd -0.000001289266 1.000000000000
PRUEBA DE NORMALIDAD JARQUE-BERA
library(normtest) #carga los comandos para las pruebas de normalidad
jb.norm.test(modelo_lineal2$residuals) #ejecuta la prueba de Jarque Bera
##
## Jarque-Bera test for normality
##
## data: modelo_lineal2$residuals
## JB = 2.7402, p-value = 0.077
## Al realizar la prueba de JARQUE-BERA dio como resultado el valor p=0.0805 siendo este mayor a alfa(0.05), por lo que se concluye que los residuos provienen de una distribucion normal, no se rechaza Ho.
qqnorm(modelo_lineal2$residuals)
qqline(modelo_lineal2$residuals)
hist(modelo_lineal2$residuals, main = "Histograma de residuos", xlab = "Residuos", ylab = "Frecuencia")
PRUEBA DE NORMALIDAD DE KOLMOGOROV-SMIRNOV
library(nortest) #carga los comandos para las pruebas de normalidad
lillie.test(modelo_lineal2$residuals) #ejecuta la prueba KS, con la corrección de Lilliefors
##
## Lilliefors (Kolmogorov-Smirnov) normality test
##
## data: modelo_lineal2$residuals
## D = 0.12278, p-value = 0.5968
## Al realizar la prueba de KOLMOGOROV-SMIRNOV dio como resultado el valor p=0.5968 siendo este mayor a alfa(0.05), por lo que se concluye que los residuos provienen de una distribucion normal, no se rechaza Ho.
qqnorm(modelo_lineal2$residuals)
qqline(modelo_lineal2$residuals)
hist(modelo_lineal2$residuals, main = "Histograma de Residuos", xlab = "Residuos", ylab = "Frecuencia")
PRUEBA DE NORMALIDAD DE SHAPIRO-WILK
shapiro.test(modelo_lineal2$residuals)
##
## Shapiro-Wilk normality test
##
## data: modelo_lineal2$residuals
## W = 0.94463, p-value = 0.2929
## Al realizar la prueba de SHAPIRO-WILK dio como resultado el valor p=0.2929 siendo este mayor a alfa(0.05), por lo que se concluye que los residuos provienen de una distribucion normal, no se rechaza Ho.