In jumanbar/siabox: Datos y Funciones Para Trabajar con el SIA

knitr::opts_chunk$set(echo = TRUE, fig.width = 6, fig.height = 4,
                      fig.align='center')

Resumen

vignette('validacion', package = 'siabox')

Un ejemplo para iniciar

library(tidyverse)
library(siabox)
alt <- siabox::datos_sia %>%
  dplyr::mutate(
    limite_deteccion = as.numeric(limpia_num(limite_deteccion)),
    limite_cuantificacion = as.numeric(limpia_num(limite_cuantificacion))
    ) %>%
  dplyr::left_join(
    dplyr::select(siabox::sitio,
                  id_sitio, id_matriz, id_estacion = id_interno),
    by = c("id_matriz", "id_estacion")
  ) %>%
  dplyr::select(id_sitio, codigo_pto, fecha_muestra, id_matriz, id_parametro,
                param, valor, limite_deteccion, limite_cuantificacion) %>%
  v_valor %>%
  filter(menor_minimo | mayor_maximo | valor_menor_ld | valor_menor_lc | ld_mayor_lc)

alt %>%
  dplyr::filter(menor_minimo) %>%
  dplyr::select(codigo_pto, param, fecha_muestra,
                valor_min, valor_min_sitio, vinf, valor)

alt %>%
  dplyr::filter(menor_minimo, !is.na(valor_min_sitio), valor > 0) %>%
  dplyr::select(codigo_pto, param, fecha_muestra,
                valor_min, valor_min_sitio, vinf, valor)

alt %>%
  dplyr::filter(mayor_maximo) %>%
  dplyr::select(codigo_pto, param, fecha_muestra,
                valor_max, valor_max_sitio, vsup, valor)

alt %>%
  dplyr::filter(mayor_maximo, !is.na(valor_max_sitio)) %>%
  dplyr::select(codigo_pto, param, fecha_muestra,
                valor_max, valor_max_sitio, vsup, valor)

alt %>%
  dplyr::filter(valor_menor_ld) %>%
  dplyr::select(codigo_pto, param, fecha_muestra, limite_deteccion, valor)

alt %>%
  dplyr::filter(valor_menor_lc) %>%
  dplyr::select(codigo_pto, param, fecha_muestra, limite_cuantificacion, valor)

alt %>%
  dplyr::filter(ld_mayor_lc) %>%
  dplyr::select(codigo_pto, param, fecha_muestra,
                limite_deteccion, limite_cuantificacion)

alt %>%
  dplyr::filter(ld_mayor_lc, limite_deteccion > limite_cuantificacion) %>%
  dplyr::select(codigo_pto, param, fecha_muestra,
                limite_deteccion, limite_cuantificacion)

Este ejemplo usa la tabla datos_sia, la cual viene incluida en el paquete. Para saber más sobre los datos incluídos, se recomienda la viñeta correspondiente:

vignette('datos-incluidos', package = 'siabox')

# sitios_ru <- filter(sia_estacion, prog_monitoreo == 13L) %>% left_join(sitio, by = c("id" = "id_interno")) %>% pull(id_sitio)
# rango_param_sitio %>% filter(id_sitio %in% sitios_ru) %>% print.data.frame
# d <- listaPD$datos %>% 
#   left_join(filter(sia_estacion, prog_monitoreo == 13L), 
#             by = c("Estacion" = "codigo_pto"))

listaPD$datos$`NO3 (mg NO3-N/L)`[11] <- "7,83"
listaPD$datos$`NO2 (mg NO2-N/L)`[10] <- "0,71"
listaPD$datos$`pH (sin unid)`[11] <- "0,66"
listaPD$datos$`PT (µg P/L)`
listaPD$datos$`PO4 (µg PO4-P/L)`[2:3] <- c("39", "50")
listaPD$datos$`OD (mg/L)`[15] <- "119"
listaPD$datos$`SatO2 (%)`[15] <- "9,94"

datos_largo <-
  largo(listaPD$datos, listaPD$ppd) %>%
  valores_numericos(metodo = "basico") %>%
  # dplyr::filter(id_tipo_dato == 1L) %>% # Esto no cambia nada...
  dplyr::left_join(sia_parametro[c(1, 4)])
v_rel_param(datos_largo, "nfila")

filter(sia_estacion, prog_monitoreo == 13L & tipo_punto_id == 1L) %>% 
  transmute(Estacion = codigo_pto, id_interno = id) %>% 
  right_join(datos_largo, by = "Estacion") %>% 
  mutate(id_matriz = 6L) %>% 
  left_join(sitio, by = c("id_matriz", "id_interno")) %>% 
  v_rango_()