1. Introduction

This is the third notebook

library (tidyverse)
## -- Attaching packages --------------------------------------- tidyverse 1.3.0 --
## v ggplot2 3.3.3     v purrr   0.3.4
## v tibble  3.0.6     v dplyr   1.0.5
## v tidyr   1.1.3     v stringr 1.4.0
## v readr   1.4.0     v forcats 0.5.1
## -- Conflicts ------------------------------------------ tidyverse_conflicts() --
## x dplyr::filter() masks stats::filter()
## x dplyr::lag()    masks stats::lag()
library(dplyr)
library (ggplot2)
eva_santander <- read_csv(file="C:/users/morea/Documents/UNAL/2021 - 1S/GeomaticaBasica/eva_santander_c.csv")
## 
## -- Column specification --------------------------------------------------------
## cols(
##   DPTO = col_character(),
##   COD_MUN = col_double(),
##   MUNICIPIO = col_character(),
##   GRUPO = col_character(),
##   SUBGRUPO = col_character(),
##   CULTIVO = col_character(),
##   DESAGREG = col_character(),
##   YEAR = col_double(),
##   PERIODO = col_character(),
##   HA_SIEMBRA = col_double(),
##   HA_COSECHA = col_double(),
##   TON_PROD = col_double(),
##   RENDIM = col_double(),
##   CICLO = col_character()
## )
head(eva_santander)
class(eva_santander)
## [1] "spec_tbl_df" "tbl_df"      "tbl"         "data.frame"
evita  <- select(eva_santander, MUNICIPIO, CULTIVO, YEAR, TON_PROD, RENDIM)
evita
evita_cacao <- filter(evita, CULTIVO=="CACAO")
evita_cacao
evita_carmen <- filter(evita, MUNICIPIO=="EL CARMEN DE CHUCURI")
evita_carmen
evita_carmen_2013 <- evita %>%
  filter(MUNICIPIO == "EL CARMEN DE CHUCURI") %>%
  filter(YEAR == 2013) %>%
  select(CULTIVO, RENDIM)
evita_carmen_2013
eva_santander %>%
  select(MUNICIPIO, CULTIVO, HA_SIEMBRA, TON_PROD, RENDIM) %>%
  filter(HA_SIEMBRA!=0) %>%
  mutate(RENDIM_SIEMBRA = TON_PROD/HA_SIEMBRA)
evita %>%
  group_by(CULTIVO) %>%
  summarize(mean_rend = mean(RENDIM, na.rm = TRUE))
eva_santander %>%
  group_by(CULTIVO, MUNICIPIO) %>%
  summarize(max_rend = max(RENDIM, na.rm = TRUE)) %>%
  slice(which.max(max_rend))
## Warning in max(RENDIM, na.rm = TRUE): ningun argumento finito para max;
## retornando -Inf

## Warning in max(RENDIM, na.rm = TRUE): ningun argumento finito para max;
## retornando -Inf

## Warning in max(RENDIM, na.rm = TRUE): ningun argumento finito para max;
## retornando -Inf

## Warning in max(RENDIM, na.rm = TRUE): ningun argumento finito para max;
## retornando -Inf

## Warning in max(RENDIM, na.rm = TRUE): ningun argumento finito para max;
## retornando -Inf

## Warning in max(RENDIM, na.rm = TRUE): ningun argumento finito para max;
## retornando -Inf

## Warning in max(RENDIM, na.rm = TRUE): ningun argumento finito para max;
## retornando -Inf

## Warning in max(RENDIM, na.rm = TRUE): ningun argumento finito para max;
## retornando -Inf

## Warning in max(RENDIM, na.rm = TRUE): ningun argumento finito para max;
## retornando -Inf

## Warning in max(RENDIM, na.rm = TRUE): ningun argumento finito para max;
## retornando -Inf

## Warning in max(RENDIM, na.rm = TRUE): ningun argumento finito para max;
## retornando -Inf

## Warning in max(RENDIM, na.rm = TRUE): ningun argumento finito para max;
## retornando -Inf

## Warning in max(RENDIM, na.rm = TRUE): ningun argumento finito para max;
## retornando -Inf

## Warning in max(RENDIM, na.rm = TRUE): ningun argumento finito para max;
## retornando -Inf

## Warning in max(RENDIM, na.rm = TRUE): ningun argumento finito para max;
## retornando -Inf

## Warning in max(RENDIM, na.rm = TRUE): ningun argumento finito para max;
## retornando -Inf

## Warning in max(RENDIM, na.rm = TRUE): ningun argumento finito para max;
## retornando -Inf
## `summarise()` has grouped output by 'CULTIVO'. You can override using the `.groups` argument.
eva_santander %>%
  group_by(GRUPO, MUNICIPIO) %>%
  summarize(max_rend = max(RENDIM, na.rm = TRUE)) %>%
  slice(which.max(max_rend))
## Warning in max(RENDIM, na.rm = TRUE): ningun argumento finito para max;
## retornando -Inf

## Warning in max(RENDIM, na.rm = TRUE): ningun argumento finito para max;
## retornando -Inf

## Warning in max(RENDIM, na.rm = TRUE): ningun argumento finito para max;
## retornando -Inf
## `summarise()` has grouped output by 'GRUPO'. You can override using the `.groups` argument.
eva_santander %>% 
  filter(YEAR==2013) %>% 
  group_by(GRUPO, MUNICIPIO) %>% 
  summarize(max_area_cosecha = max(HA_COSECHA, na.rm = TRUE)) %>%
  slice(which.max(max_area_cosecha)) %>%
  
    arrange(desc(max_area_cosecha)) -> area_cosecha_max
## `summarise()` has grouped output by 'GRUPO'. You can override using the `.groups` argument.
area_cosecha_max
eva_santander %>% 
  group_by(GRUPO, MUNICIPIO, YEAR) %>%
  summarize(max_prod = max(TON_PROD, na.rm = TRUE)) %>%
    slice(which.max(max_prod)) %>%
    arrange(desc(max_prod)) -> ton_prod_max
## `summarise()` has grouped output by 'GRUPO', 'MUNICIPIO'. You can override using the `.groups` argument.
ton_prod_max
eva_santander %>% 
  filter(MUNICIPIO=="EL CARMEN DE CHUCURI" & CULTIVO=="CACAO") %>% 
  group_by(YEAR, CULTIVO) %>%
  select(MUNICIPIO, CULTIVO, TON_PROD, YEAR) ->  carmen_cacao

carmen_cacao
# we use the ggplot 2 library
g <- ggplot(aes(x=YEAR, y=TON_PROD/1000), data = carmen_cacao) + geom_bar(stat='identity') + labs(y='Produccion de cacao [Ton x 1000]')
g + ggtitle("Evolution of Cacao Crop Production in El Carmen de Chucuri from 2007 to 2018") + labs(caption= "Based on EVA data (Minagricultura, 2020)")

# Prepare data
data <- eva_santander %>%
       filter(CULTIVO=="CACAO" & YEAR==2018 & TON_PROD>35.000) %>% 
       select(MUNICIPIO, TON_PROD)
# Basic piechart
ggplot(data, aes(x="", y= TON_PROD/1000, fill= MUNICIPIO)) +
    geom_bar (stat="identity", width=1) +
    coord_polar("y", start=0)

`