Munucipality age development 1987-2019

mypath <- "C://Users//np83zg//OneDrive - Aalborg Universitet//Skrivebord//kommunedata//data"
setwd(mypath)
library(data.table)
library(tmap)
library(sf)

## Linking to GEOS 3.8.0, GDAL 3.0.4, PROJ 6.3.1

library(cartogram)

## Warning: package 'cartogram' was built under R version 4.0.3

library(gganimate)

## Warning: package 'gganimate' was built under R version 4.0.3

## Loading required package: ggplot2

library(ggplot2)


mydata <- as.data.table(read.table("bef_tab2.txt"))
setkey(mydata,residence,year)
T <- length(table(mydata$year))
mydata[,dm_age:=m_age-m_age[1],by=.(residence)]
mydata[,index:=-dm_age[T],by=.(residence)]
setkey(mydata,index,year)
mydata[,id:=rep(1:98,each=33)]


theme_set(theme_bw())

p <- ggplot(
  mydata, 
  aes(x =id, y=dm_age, colour = as.factor(residence))
  ) +
  geom_point(show.legend = FALSE, alpha = 0.7) +
  scale_color_viridis_d() +
  scale_size(range = c(2, 12)) +
  labs(x = "Municipality ID", y = "Change in mean age relative to year 1987") +
  transition_time(year) +
  labs(title = "Year: {frame_time}") +
  shadow_mark(alpha = 0.3, size = 0.5)

mypath <- "C://Users//np83zg//OneDrive - Aalborg Universitet//Skrivebord//kommunedata//scripts"
setwd(mypath)
anim_save("p.gif",p)

mypath <- "C://Users//np83zg//OneDrive - Aalborg Universitet//Skrivebord//kommunedata//data"
setwd(mypath)
library(data.table)
library(tmap)
library(sf)
library(cartogram)
library(gganimate)
library(ggplot2)


mydata <- as.data.table(read.table("bef_tab3.txt"))
setkey(mydata,residence,year)
T <- length(table(mydata$year))
cities <- c(101,751,461,851,561)


mydata[,s:=age_16+age_16_60+age_60]
mydata[,s1:=age_16/s]
mydata[,s2:=age_16_60/s]
mydata[,s3:=age_60/s]
mydata[,ds2:=s2-s2[1],by=.(residence)]


theme_set(theme_bw())
p <- ggplot(
  mydata, 
  aes(x =residence, y=ds2,size=s, colour = as.factor(residence))
  ) +
  geom_point(show.legend = FALSE, alpha = 0.9) +
  scale_color_viridis_d() +
  scale_size(range = c(2, 12)) +
  labs(x = "Municipality ID", y = "Change in share in agegroup 16-60 relative 1987") +
  transition_time(year) +
  labs(title = "Year: {frame_time}") +
  shadow_mark(alpha = 0.3, size = 0.5)

mypath <- "C://Users//np83zg//OneDrive - Aalborg Universitet//Skrivebord//kommunedata//scripts"
setwd(mypath)
anim_save("p2.gif",p)

theme_set(theme_bw())
p <- ggplot(
  mydata, 
  aes(x =residence, y=s2,size=s, colour = as.factor(residence))
  ) +
  geom_point(show.legend = FALSE, alpha = 0.9) +
  scale_color_viridis_d() +
  scale_size(range = c(2, 12)) +
  labs(x = "Municipality ID", y = "Share in agegroup 16-60 relative 1987") +
  transition_time(year) +
  labs(title = "Year: {frame_time}") +
  shadow_mark(alpha = 0.3, size = 0.5)

mypath <- "C://Users//np83zg//OneDrive - Aalborg Universitet//Skrivebord//kommunedata//scripts"
setwd(mypath)
anim_save("p3.gif",p)