Este cuaderno ilustra algunas estadisticas agricolas del año 2021 en los municipios del departamento Casanare usando como fuente datos la evaluación agrícola municipal elaborada por la Unidad de Planeación Rural Agropecuaria (UPRA).
El archivo de UPRA se puede descargar de este enlace.
# setup
# install.packages("leaflet")
# install.packages('tidyverse')
# install.packages('readxl')
# cargar librerias
library(leaflet)
library(tidyverse)
## ── Attaching packages ─────────────────────────────────────── tidyverse 1.3.2 ──
## ✔ ggplot2 3.3.6 ✔ purrr 0.3.4
## ✔ tibble 3.1.8 ✔ dplyr 1.0.10
## ✔ tidyr 1.2.1 ✔ stringr 1.4.1
## ✔ readr 2.1.2 ✔ forcats 0.5.2
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag() masks stats::lag()
library(readxl)
#
#(archivos = list.files(pattern='xls'))
list.files(path='./inputs', pattern='xls')
## [1] "BaseEVA_Agrícola2021.xlsx" "EVA2021.xlsx"
# como saber cuales hojas hay en el archivo
(hojas = readxl::excel_sheets("./inputs/EVA2021.xlsx"))
## [1] "AgricolaDesagregados"
# ahora si lectura
eva2021 = readxl::read_excel("./inputs/EVA2021.xlsx", sheet = "AgricolaDesagregados")
# despliegue de los datos almacenados en la variable EVA2021
eva2021
El Departamento está ubicado en la parte oriental del pais, al noroccidente de la Orinoquia, abarcando una amplia faja de la cordillera Oriental de los Andes.
El departamento del Casanare esta dividido en 19 municipios, 11 corregimientos, 106 inspecciones de policia, asi como, numerosos caserios y sitios poblados.
-Superficie 44.640 km2 -Poblacion 325.389 Habitantes -Densidad 7,29 Habitantes por km2 -Capital Yopal: 88.124 Habitantes
casanare2021 = dplyr::filter(eva2021, DEPTO == "Casanare")
casanare2021
casanare2021.tmp <- casanare2021 %>% select ('COD_MUN':'CICLO')
casanare2021.tmp
casanare2021 %>% mutate (AREA_SEMB = as.numeric(AREA_SEMB),
AREA_COS = as.numeric(AREA_COS),
PROD = as.numeric(PROD),
REND = as.numeric(REND)) -> casanare2021
casanare2021
Encontramos los cultivos mas representativos en el departamento del Casanare.
casanare2021 %>%
filter(PROD > 0) %>%
group_by(CULTIVO) %>%
summarize(total_produccion = sum(PROD)) %>%
arrange(desc(total_produccion))
Ahora tenemos los municipios con mayor produccion para cada cultivo:
casanare2021 %>%
group_by(CULTIVO, MUN) %>%
summarize(max_prod = max(PROD, na.rm = TRUE)) %>%
slice(which.max(max_prod)) %>%
arrange(desc(max_prod))
## `summarise()` has grouped output by 'CULTIVO'. You can override using the
## `.groups` argument.
casanare2021 %>%
group_by(GRUPO,MUN) %>%
summarize(max_prod = max(PROD, na.rm = TRUE)) %>%
slice(which.max(max_prod)) %>%
arrange(desc(max_prod))
## `summarise()` has grouped output by 'GRUPO'. You can override using the
## `.groups` argument.
casanare2021 %>%
group_by(COD_MUN, MUN, GRUPO) %>%
filter(GRUPO=='Frutales') %>%
summarize(max_prod = max(PROD, na.rm = TRUE)) %>%
arrange(desc(max_prod)) -> FRUTALES2021
## `summarise()` has grouped output by 'COD_MUN', 'MUN'. You can override using
## the `.groups` argument.
FRUTALES2021
casanare2021 %>%
group_by(COD_MUN, MUN, GRUPO) %>%
filter(GRUPO=='Oleaginosas') %>%
summarize(max_prod = max(PROD, na.rm = TRUE)) %>%
arrange(desc(max_prod)) -> OLEAGINOSAS2021
## `summarise()` has grouped output by 'COD_MUN', 'MUN'. You can override using
## the `.groups` argument.
OLEAGINOSAS2021
# write_csv (OLEAGINOSAS2021, "./CASANARE_oleag_2021.csv")
Unidad de Planificacion Rural Agropecuaria (UPRA) este enlace
Gobernacion de Casanare este enlace
Lizarazo, I. Reading and processing municipal agricultural statistics for 2021. Available at este enlace