data-raw/create_p97_capturas_historico.R

library(devtools)
original_wd <- getwd()
setwd("data-raw")

catches <- "IEOUPMUEDESTOTSIRENO_2014_2018.TXT"

load("../data/relacion_variables.rda")
load("../data/formato_variables.rda")

# I don't know why but I've to load the function catchesPercentileUP here
# in order to work the scritp... ¿?¿?¿?
catchesPercentileUP <- function(dfs, path = getwd(), per){

  library(dplyr)

  data <- sapmuebase::importRIMCatches(dfs, path)

  percentile <- data %>%
    select(COD_ID, COD_ESP_MUE, ESTRATO_RIM, P_DESEM) %>%
    group_by(COD_ID, COD_ESP_MUE, ESTRATO_RIM) %>%
    summarise('p_desem_total' = sum(P_DESEM)) %>%
    ungroup()%>%
    select(COD_ESP_MUE, ESTRATO_RIM, p_desem_total) %>%
    group_by(COD_ESP_MUE, ESTRATO_RIM) %>%
    summarise('97' = quantile(p_desem_total, probs=per, na.rm = TRUE))

  # convert tibble to dataframe
  percentile <- as.data.frame(percentile)

  return(percentile)

}

p97_capturas_historico <- catchesPercentileUP(catches, path = paste0(getwd(),"/data_source_private"), per = 0.97)

colnames(p97_capturas_historico) <- c("COD_ESP", "ESTRATO_RIM", "P97")

# use_data() create the file in /data. The file created has extension .rda,
# instead of .Rdata.
usethis::use_data(p97_capturas_historico, overwrite = TRUE)

devtools::document()

setwd(original_wd)
rm(original_wd, p97_capturas_historico)
Eucrow/sapmuebase documentation built on April 17, 2024, 10:29 a.m.