Rape_victims_Denmark.R

library(pacman)
p_load(ggplot2, dkstat, stringr, plyr, kirkegaard)

#install kirkegaard and dkstat from github if needed
#library("devtools")
#install_github("deleetdk/kirkegaard")
#install_github("rOpenGov/dkstat")

# load data ---------------------------------------------------------------
#get meta-data
dst_meta = dst_meta(table = "STRAF5", lang = "en")

#time
dst_meta$values$Tid

#fetch data
d_data = dst_get_data(table = "STRAF5", KOEN = "Women", Tid = dst_meta$values$Tid$text, OVERTRÆD = "Rape, etc.", ALDER = dst_meta$values$ALDER[-c(1, 14), 2], 
                       lang = "en")

# transform data --------------------------------------------------------
#sum across years
d_data2 = ddply(.data = d_data, .variables = .(ALDER), .fun = summarize,
                sum = sum(value))

#remove "years"
d_data2$ALDER = str_replace(d_data2$ALDER, pattern = " years", replacement = "")
d_data2$ALDER[12] = "80 and older"

#levels
d_data2$ALDER = ordered(d_data2$ALDER, levels = c(d_data2$ALDER[c(1, 8, 2:7, 9:12)]))

# plots -------------------------------------------------------------------

ggplot(d_data2, aes(ALDER, sum)) + geom_bar(stat = "identity") + 
  xlab("Age group") + ylab("Number of victims (self-reported)") + 
  ggtitle("Rape victims in Denmark " + min(dst_meta$values$Tid$id) + "-" + max(dst_meta$values$Tid$id) + " (self-reported, female-only)")
ggsave("rape_victims_Denmark.png")