Evolución de la edad media de la población por provincias

Estadística, R

Abundo en la entrada de ayer. Lo hago para mostrar

En el gráfico anterior se muestra la evolución de la edad media de la población de las provincias españolas como diferencia con respecto a una evolución media calculada como la regresión lineal de todas las edades medias con respecto al año. Es decir, algo así como evolución relativa.

Se aprecian claramente los rejuvenecimientos relativos de Guadalajara y, en menor medida, Toledo. Especialmente acusados durante este siglo.

Los datos son del INE y descargados en formato PC-Axis. El código,

library(pxR)
library(ggplot2)

pob <- read.px("3199.px")
pob <- as.data.frame(pob)
pob$Sexo <- NULL

pob$Periodo <- as.numeric(as.character(pob$Periodo))

# ggplot(pob, aes(x = Periodo, y = value)) +
#   geom_line() + facet_wrap(~ Provincias)

modelo <- lm(value ~ Periodo, dat = pob)

pob$pred <- predict(modelo, pob)
pob$diff <- pob$value - pob$pred

tmp <- pob[pob$Periodo == max(pob$Periodo),]
tmp$Provincias <- reorder(tmp$Provincias, tmp$diff, max)
pob$Provincias <- factor(pob$Provincias, levels = levels(tmp$Provincias))
#pob$Provincias <- reorder(pob$Provincias, pob$diff, max)

pob$fecha <- as.Date(paste0(pob$Periodo, "-12-31"))

ggplot(pob, aes(x = fecha, y = diff)) +
  geom_hline(aes(yintercept = 0, col = "red", alpha = 0.5)) +
  geom_line() + facet_wrap(~ Provincias, ncol = 4) +
  theme(legend.position="none")