Ders 5

library(WDI)

## Warning: package 'WDI' was built under R version 4.2.2

Türkiye, Almanya, Fransa ve Japonya’nın bu konu hakkındaki verileri 1990-2020

df = WDI(indicator='SL.UEM.TOTL.NE.ZS', country=c('TR','DE','FR','JP'), start=1990, end=2020)

head(df)

##   country iso2c iso3c year SL.UEM.TOTL.NE.ZS
## 1 Germany    DE   DEU 2020              3.86
## 2 Germany    DE   DEU 2019              3.14
## 3 Germany    DE   DEU 2018              3.38
## 4 Germany    DE   DEU 2017              3.75
## 5 Germany    DE   DEU 2016              4.12
## 6 Germany    DE   DEU 2015              4.62

library(dplyr)

## Warning: package 'dplyr' was built under R version 4.2.2

## 
## Attaching package: 'dplyr'

## The following objects are masked from 'package:stats':
## 
##     filter, lag

## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

df <- df %>%
  rename(ulke = 1,
         ulkekoduiso2 = 2,
         ulkekoduiso3 = 3,
         sene = 4,
         İşsizlik = 5)

library(reshape2)

## Warning: package 'reshape2' was built under R version 4.2.2

data_genis <- dcast(df, sene ~ ulke, value.var="İşsizlik")

tail(data_genis)

##    sene France Germany Japan Turkiye
## 26 2015  10.35    4.62  3.39   10.24
## 27 2016  10.05    4.12  3.13   10.84
## 28 2017   9.41    3.75  2.82   10.82
## 29 2018   9.02    3.38  2.47   10.89
## 30 2019   8.41    3.14  2.35   13.67
## 31 2020   8.01    3.86  2.80   13.11

df2 <- WDI(
  country = "all",
  indicator = "SL.UEM.TOTL.NE.ZS",
  start = 1990,
  end = 2020,
  extra = TRUE,
  cache = NULL,
  latest = NULL,
  language = "tr")

## Warning in WDI(country = "all", indicator = "SL.UEM.TOTL.NE.ZS", start = 1990, :
## This language is only supported partially.

df3<- df2 %>%
  select(iso2c, year, SL.UEM.TOTL.NE.ZS)

data_genis<-  data_genis %>% mutate(ulkekoduiso2=NULL, ulkekoduiso3=NULL)

tsveri <- ts(data_genis, start=1990, frequency=1)

Data_genis veri setimizi, bir zaman serisine çevirelim. R’da bu tür zaman serilerine, ts denilir.

tail(tsveri)

##       sene France Germany Japan Turkiye
## [26,] 2015  10.35    4.62  3.39   10.24
## [27,] 2016  10.05    4.12  3.13   10.84
## [28,] 2017   9.41    3.75  2.82   10.82
## [29,] 2018   9.02    3.38  2.47   10.89
## [30,] 2019   8.41    3.14  2.35   13.67
## [31,] 2020   8.01    3.86  2.80   13.11

frequency, seçeneği aylık verilerde 12, çeyreklik verilerde 4, haftalık verilerde 52 olarak girilir. Türkiye’nin verisini, zaman serisi grafiği olarak çizmek istersek.

library(ggplot2)

## Warning: package 'ggplot2' was built under R version 4.2.2

library(ggfortify)

## Warning: package 'ggfortify' was built under R version 4.2.2

Türkiye, Almanya, Japonya ve Fransa’nın verileri, ikinci, üçüncü, dördüncü ve beşinci sütundadır.

plot(tsveri[,2:5])

plot(tsveri[,2:5], plot.type = "single")

plot(tsveri[,2],
     type = "l",
     col = 2,
     ylim = c(0, 25),
     xlab = "Sene",
     ylab = "İşsizlik")
lines(tsveri[,3],
      type = "l",
      col = 3)
lines(tsveri[,4],
      type = "l",
      col = 4)
lines(tsveri[,5],
      type = "l",
      col = 5)
legend("topright",                           
       c("Türkiye", "Almanya", "Fransa","Japonya"),
       lty = 1,
       col = 2:5)

Veri setinde sene isimli değişkene göre, veri setini uzun formatına geri dönmesi

data_uzun <- melt(data_genis, id.vars = "sene") 

data_uzun veri seti artık variable ve value isimli iki değişkene sahip, variable değişkeni dört ülkeyi sırasıyla gösteriyor, value değişkeni ise bu ülkelerdeki GSYİH değerlerini veriyor.

ggplot(data_uzun,                           
       aes(x = sene,
           y = value,
           col = variable)) +
  geom_line()