---
title: "Día 7 #30díasdegráficos (ridgeline). Tamaño de las unidades agropecuarias con superficie agrícola. Perú 2012 (ESP)"
author: "Rubén F. Bustillo"
output:
flexdashboard::flex_dashboard:
orientation: columns
source_code: embed
vertical_layout: fill
theme: simplex
---
```{r setup, include=FALSE}
# PACKAGES / LIBRARIES:
library(flexdashboard)
library(tidyverse)
library(foreign)
library(readxl)
library(ggridges)
library(extrafont)
loadfonts(device = "win")
options(scipen=999)
AMA_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/337-Modulo229/337-Modulo229/01_IVCENAGRO_REC01.dbf")
ANC_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/338-Modulo229/338-Modulo229/01_IVCENAGRO_REC01.dbf")
APU_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/339-Modulo229/339-Modulo229/01_IVCENAGRO_REC01.dbf")
ARE_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/340-Modulo229/340-Modulo229/01_IVCENAGRO_REC01.dbf")
AYA_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/341-Modulo229/341-Modulo229/01_IVCENAGRO_REC01.dbf")
CAJ_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/342-Modulo229/342-Modulo229/01_IVCENAGRO_REC01.dbf")
CAL_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/343-Modulo229/343-Modulo229/01_IVCENAGRO_REC01.dbf")
CUS_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/344-Modulo229/344-Modulo229/01_IVCENAGRO_REC01.dbf")
HUAV_M229 <-read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/345-Modulo229/345-Modulo229/01_IVCENAGRO_REC01.dbf")
HUAC_M229 <-read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/346-Modulo229/346-Modulo229/01_IVCENAGRO_REC01.dbf")
ICA_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/347-Modulo229/347-Modulo229/01_IVCENAGRO_REC01.dbf")
JUN_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/348-Modulo229/348-Modulo229/01_IVCENAGRO_REC01.dbf")
LLI_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/349-Modulo229/349-Modulo229/01_IVCENAGRO_REC01.dbf")
LAM_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/350-Modulo229/350-Modulo229/01_IVCENAGRO_REC01.dbf")
LIM_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/351-Modulo229/351-Modulo229/01_IVCENAGRO_REC01.dbf")
LOR_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/352-Modulo229/352-Modulo229/01_IVCENAGRO_REC01.dbf")
MDD_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/353-Modulo229/353-Modulo229/01_IVCENAGRO_REC01.dbf")
MOQ_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/354-Modulo229/354-Modulo229/01_IVCENAGRO_REC01.dbf")
PAS_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/355-Modulo229/355-Modulo229/01_IVCENAGRO_REC01.dbf")
PIU_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/356-Modulo229/356-Modulo229/01_IVCENAGRO_REC01.dbf")
PUN_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/357-Modulo229/357-Modulo229/01_IVCENAGRO_REC01.dbf")
SMA_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/358-Modulo229/358-Modulo229/01_IVCENAGRO_REC01.dbf")
TAC_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/359-Modulo229/359-Modulo229/01_IVCENAGRO_REC01.dbf")
TUM_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/360-Modulo229/360-Modulo229/01_IVCENAGRO_REC01.dbf")
UCA_M229 <- read.dbf("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/AGRO_PERU/CENSO_NACIONAL_AGROPECUARIO/Modulo229/361-Modulo229/361-Modulo229/01_IVCENAGRO_REC01.dbf")
all_M229 <- rbind(AMA_M229,
ANC_M229,
APU_M229,
ARE_M229,
AYA_M229,
CAJ_M229,
CAL_M229,
CUS_M229,
HUAV_M229,
HUAC_M229,
ICA_M229,
JUN_M229,
LLI_M229,
LAM_M229,
LIM_M229,
LOR_M229,
MDD_M229,
MOQ_M229,
PAS_M229,
PIU_M229,
PUN_M229,
SMA_M229,
TAC_M229,
TUM_M229,
UCA_M229)
#nos quedamos unas cuantas columnas:
all_M229 <- all_M229 %>%
select("TIPO_REC",
"P001",
"P008",
"P020_01",
"WSUP01",
"WSUP02",
"WSUP02A",
"WSUP03",
"WSUP04",
"WSUP05",
"P019",
"WREGION")
# departamentos por ubigeo:
ubigeo_departamentos <- read_excel("C:/Users/Usuario/Desktop/Para_paper_peru/DATOS_PERU/REGIONAL/POR_REGION/UBIGEO/ubigeo_departamentos.xlsx")
departamentos<-ubigeo_departamentos %>%
select(P001=ubigeo, Departamento, id)
all_M229 <- left_join(all_M229, departamentos)
# nueva columna para costa, sierra y selva:
all_M229$WREGION <- as.character(all_M229$WREGION)
all_M229 <- all_M229 %>%
mutate(Geo = if_else(WREGION == "1", "Costa", WREGION)) %>%
mutate(Geo = if_else(WREGION == "2", "Sierra", Geo)) %>%
mutate(Geo = if_else(WREGION == "3", "Selva", Geo))
# Labels en eje y:
all_M229$Departamento <- as.factor(all_M229$Departamento)
y_labels <- paste(levels(all_M229$Departamento),"\n(N=", table(all_M229$Departamento),")",sep="")
```
```{r, out.width="100%", fig.width=10, fig.height=5}
all_M229 %>%
filter(WSUP03 > 0) %>%
ggplot(aes(x= WSUP03, y = Departamento, fill = factor(WREGION))) +
facet_wrap(~Geo, scale = "free_y") +
geom_density_ridges(alpha=0.6, stat="binline", bins=30)+
scale_x_continuous(trans = 'log10', labels = scales::comma) +
theme_minimal() +
scale_fill_manual(values = c("steelblue1", "orange", "limegreen")) +
geom_vline(xintercept = 5, color = "red", size = 1.2) +
labs(title = "Tamaño de las Unidades Agropecuarias con superficie agrícola \nsegún departamento y región natural (2012)",
subtitle = "El número total de UA con superficie agrícola* es de 2,128,087 sumando un total de 7,125,008 hectáreas. \nSe observa una fuerte atomización de la superficie agrícola en tanto que el 79,9% de las UAs tienen una superficie inferior 5 hectáreas (línea roja). \nSin embargo el porcentaje de estas UA de menor tamaño solo suponen el 25,5% de la superficie agrícola total.",
caption = "Fuente: Censo Nacional Agropecuario (CENAGRO) \n* Excluye superficie no agrícola y otra clase de tierras",
y = "",
x = "Superficie Agrícola (has) (escala logarítmica) ") +
theme(legend.position = "none",
plot.title=element_text(size = 11,
family = "Modern No. 20",
hjust = 0),
plot.subtitle = element_text(size = 7,
family = "Modern No. 20",
hjust = 0),
text = element_text(size = 7),
axis.line.x = element_line(size = 1,
colour = "black"),
axis.line.y = element_line(size = 1,
colour = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
panel.background = element_blank(),
axis.text.x=element_text(colour="black",
size = 7),
axis.text.y=element_text(colour="black",
size = 6)) +
scale_y_discrete(labels = y_labels)
```
