Jeg udvider min seneste analyse af udviklingen i tysk transportvolumen ved at sammenligne forskellige indlandsgodskategorier i et ggplot-diagram. Dataene forespurgte jeg igen ved hjælp af OECD-pakken i R.
library(OECD)
Fra det forrige indlæg kender vi allerede ID-nøglen til det transportrelaterede datasæt af interesse. Ved hjælp af get_dataset (R-funktion fra OECD-databasen) trækker jeg dataene via OECD-grænsefladen i R .
data_df <- as.data.frame(get_dataset(dataset = "ITF_GOODS_TRANSPORT"))
Ved hjælp af dplyr bortfiltrerer jeg dataindtastningerne af interesse:
library(dplyr)
colnames(data_df) <- c("country","variable","timeformat","unit","powercode","obsTime","obsValue","obsStatus")
data_df <- dplyr::filter(data_df,country=="DEU")
data_df <- dplyr::filter(data_df,timeformat=="P1Y")
data_df <- dplyr::filter(data_df,unit=="TONNEKM")
data_df <- data_df[is.na(data_df$obsStatus),]
Efter filtrering er poster stadig kendetegnet ved variable indikatorer. For at jeg skal kunne fortolke dem trækker jeg datastrukturen ved at bruge get_data_structure-funktionen fra OECD-pakken i R :
data_struct <- get_data_structure("ITF_GOODS_TRANSPORT")
data_struct$VARIABLE
## id label
## 1 T-GOODS-TOT-INLD Total inland freight transport
## 2 T-GOODS-RL-TOT Rail freight transport
## 3 T-GOODS-RD-TOT Road freight transport
## 4 T-GOODS-RD-REW Road freight transport for hire and reward
## 5 T-GOODS-RD-OWN Road freight transport on own account
## 6 T-GOODS-IW-TOT Inland waterways freight transport
## 7 T-GOODS-PP-TOT Pipelines transport
## 8 T-SEA-CAB Coastal shipping (national transport)
## 9 T-SEA Maritime transport
## 10 T-CONT-RL-TEU Rail containers transport (TEU)
## 11 T-CONT Containers transport
## 12 T-CONT-RL-TON Rail containers transport (weight)
## 13 T-CONT-SEA-TEU Maritime containers transport (TEU)
## 14 T-CONT-SEA-TON Maritime containers transport (weight)
Jeg kan nu oprette et ggplot-stiplot, der sammenligner følgende kategorier af interesse for mig: – Indlandsvejsgods – Indenrigsgodstransport – Indlandsrørledningstransport – Transport ad indre vandveje
library(ggplot2)
ggplot(data_df[data_df$variable == c("T-GOODS-IW-TOT",
"T-GOODS-RD-TOT",
"T-GOODS-RL-TOT"),]) +
geom_path(mapping = aes(x=as.numeric(obsTime),y=obsValue, color=variable)) +
ggtitle("German inland freight development by considered category") +
xlab("year") +
ylab("in millions of TONNEKM")
Industriingeniør som gerne beskæftiger sig med optimering, simulation og matematisk modellering i R, SQL, VBA og Python
Leave a Reply