İlkÖdev

Dünya Bankasından Veri Seti İndirme

library(WDI)
library(ggplot2)

data <- WDI(country = "TR", 
            indicator = c("NY.GDP.PCAP.CD", "SP.DYN.LE00.IN"))

head(data)

##   country iso2c iso3c year NY.GDP.PCAP.CD SP.DYN.LE00.IN
## 1 Turkiye    TR   TUR 1960       275.0417         50.740
## 2 Turkiye    TR   TUR 1961       282.7425         51.550
## 3 Turkiye    TR   TUR 1962       307.3063         52.382
## 4 Turkiye    TR   TUR 1963       347.1771         53.173
## 5 Turkiye    TR   TUR 1964       365.1339         53.714
## 6 Turkiye    TR   TUR 1965       381.4133         54.258

Kişi başına düşen GSYİH’nin yıllara göre grafiği:

gdp_plot <- ggplot(data, aes(x = year, y = NY.GDP.PCAP.CD)) +
  geom_line(color = "blue", size = 1) +
  geom_point(color = "darkblue") +
  labs(title = "Kişi Başına Düşen GSYİH (Türkiye)",
       x = "Yıl", 
       y = "Kişi Başına GSYİH (USD)") +
  theme_minimal()

## Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
## ℹ Please use `linewidth` instead.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.

gdp_plot

Yaşam beklentisinin yıllara göre grafiği:

life_exp_plot <- ggplot(data, aes(x = year, y = SP.DYN.LE00.IN)) +
  geom_line(color = "green", size = 1) +
  geom_point(color = "darkgreen") +
  labs(title = "Yaşam Beklentisi (Türkiye)",
       x = "Yıl", 
       y = "Yaşam Beklentisi (Yıl)") +
  theme_minimal()

life_exp_plot

## Warning: Removed 1 row containing missing values or values outside the scale range
## (`geom_line()`).

## Warning: Removed 1 row containing missing values or values outside the scale range
## (`geom_point()`).

## Regresyon Modelini Oluşturma:

regres_model <- lm(SP.DYN.LE00.IN ~ NY.GDP.PCAP.CD, data = data)

Regresyon modelini Özetleme:

summary(regres_model)

## 
## Call:
## lm(formula = SP.DYN.LE00.IN ~ NY.GDP.PCAP.CD, data = data)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -9.1525 -3.0832 -0.2042  3.8385  7.4614 
## 
## Coefficients:
##                 Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    59.417686   0.766005   77.57   <2e-16 ***
## NY.GDP.PCAP.CD  0.001726   0.000131   13.17   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 4.192 on 61 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.7399, Adjusted R-squared:  0.7356 
## F-statistic: 173.5 on 1 and 61 DF,  p-value: < 2.2e-16

GSYİH’nin varyansı

var_gdp <- var(data$NY.GDP.PCAP.CD, na.rm = TRUE)

Yaşam Beklentisi’nin varyansı

var_life_exp <- var(data$SP.DYN.LE00.IN, na.rm = TRUE)

var_gdp

## [1] 17436234

var_life_exp

## [1] 66.45478

GSYİH ile Yaşam Beklentisi’nin kovaryansı

cov_gdp_life_exp <- cov(data$NY.GDP.PCAP.CD, data$SP.DYN.LE00.IN, use = "complete.obs")

cov_gdp_life_exp

## [1] 28484.47

Regresyon Yorumu

GSYİH ve Yaşam Beklentisi Arasındaki İlişki: Eğer eğim pozitifse (örneğin 0.0005), GSYİH arttıkça yaşam beklentisinin de arttığı anlamına gelir. Bu, zengin ülkelerde yaşam beklentisinin daha uzun olduğunu gösteren bir bulgu olabilir. R-Kare Değeri: Modelin, yaşam beklentisini ne kadar iyi açıkladığını gösterir. Yüksek bir R² değeri, bağımsız değişkenin bağımlı değişken üzerinde güçlü bir etkisi olduğunu gösterir. P-Değeri: Eğer p-değeri 0.05’ten küçükse, bu ilişkinin istatistiksel olarak anlamlı olduğu ve modelin geçerli olduğu anlamına gelir.