##————————————— ## Econometrics in R - Cvičenie 5 ## Autor: ALINA ILNYTSKA ##—————————————

0. Knižnice

library(zoo) library(tseries) library(lmtest) library(sandwich) library(car)

Vyčistiť workspace

rm(list = ls())

1. Nastavenie pracovného adresára

setwd(“/Cloud/project/tyzdne/tyzden5”) # podľa potreby

2. Načítanie dát

udaje <- read.csv(“udaje/Life_Expectancy_Data.csv”, dec = “.”, sep = “,”, header = TRUE)

3. Výber údajov len za rok 2015 a relevantných premenných

udaje.2015 <- udaje[udaje$Year == 2015, c(“Life.expectancy”, “BMI”, “GDP”, “Schooling”)]

4. Imputácia chýbajúcich hodnôt mediánom

column_medians <- sapply(udaje.2015, median, na.rm = TRUE) udaje_imputed <- udaje.2015 for (col in names(udaje.2015)) { udaje_imputed[[col]][is.na(udaje_imputed[[col]])] <- column_medians[col] } udaje.2015 <- udaje_imputed

5. Kontrola nulových a extrémnych hodnôt

summary(udaje.2015) colSums(udaje.2015 == 0)

6. Boxploty pre vizuálnu kontrolu extrémov

par(mfrow = c(2, 2)) par(mar = c(4, 4, 2, 1)) par(oma = c(0, 0, 3, 0)) for (col in names(udaje.2015)) { boxplot(udaje.2015[[col]], main = col, col = “lightblue”) } mtext(“Boxploty jednotlivých premenných”, outer = TRUE, cex = 1.4, font = 2) par(mfrow = c(1, 1)) # reset layout

7. Lineárna regresia

model <- lm(Life.expectancy ~ BMI + GDP + Schooling, data = udaje.2015)

Základný súhrn modelu

summary(model)

8. Diagnostické grafy

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

9. Test normality rezíduí (Jarque-Bera)

residuals <- residuals(model) jb_test <- jarque.bera.test(residuals) jb_test

10. Test odľahlých hodnôt (outlier test)

outlier_test <- outlierTest(model) outlier_test

11. Matica regresorov a leverage

X <- model.matrix(model) hatvalues(model) # diagonálne prvky projekčnej matice