tikatuwq: basic workflow

tikatuwq authors

Load demo dataset

utils::data("wq_demo", package = "tikatuwq")
head(wq_demo)
#>   ponto       data   ph   od turbidez  dbo coliformes p_total nt_total
#> 1    P1 2024-06-07 7.08 5.33     63.8 5.38        331   0.080     1.30
#> 2    P1 2024-06-27 6.79 8.29     35.7 1.46         98   0.348     3.69
#> 3    P1 2024-01-14 8.43 4.23     60.6 1.22       3712   0.313     4.25
#> 4    P1 2024-07-13 8.30 6.21     50.7 4.09        200   0.447     2.54
#> 5    P1 2024-06-18 7.88 7.99     57.1 2.39         66   0.340     2.00
#> 6    P2 2024-02-19 8.09 4.61      1.0 4.11       3930   0.371     1.31
#>   temperatura tds
#> 1        23.0 249
#> 2        21.2 329
#> 3        25.6 480
#> 4        20.2 731
#> 5        28.4 201
#> 6        27.0 724

Compute IQA and check CONAMA compliance

library(tikatuwq)

df <- wq_demo |>
  validate_wq() |>
  iqa(na_rm = TRUE)

conf <- conama_check(df, classe = "2")
head(conf)
#>   ponto       data   ph   od turbidez  dbo coliformes p_total nt_total
#> 1    P1 2024-06-07 7.08 5.33     63.8 5.38        331   0.080     1.30
#> 2    P1 2024-06-27 6.79 8.29     35.7 1.46         98   0.348     3.69
#> 3    P1 2024-01-14 8.43 4.23     60.6 1.22       3712   0.313     4.25
#> 4    P1 2024-07-13 8.30 6.21     50.7 4.09        200   0.447     2.54
#> 5    P1 2024-06-18 7.88 7.99     57.1 2.39         66   0.340     2.00
#> 6    P2 2024-02-19 8.09 4.61      1.0 4.11       3930   0.371     1.31
#>   temperatura tds      IQA ph_ok ph_status ph__lim_min ph__lim_max ph__delta
#> 1        23.0 249 58.46746  TRUE        ok           6           9         0
#> 2        21.2 329 57.31295  TRUE        ok           6           9         0
#> 3        25.6 480 41.29367  TRUE        ok           6           9         0
#> 4        20.2 731 47.46528  TRUE        ok           6           9         0
#> 5        28.4 201 59.76570  TRUE        ok           6           9         0
#> 6        27.0 724 40.32756  TRUE        ok           6           9         0
#>   od_ok        od_status od__lim_min od__lim_max od__delta turbidez_ok
#> 1  TRUE               ok           5          NA      0.00        TRUE
#> 2  TRUE               ok           5          NA      0.00        TRUE
#> 3 FALSE abaixo_do_minimo           5          NA     -0.77        TRUE
#> 4  TRUE               ok           5          NA      0.00        TRUE
#> 5  TRUE               ok           5          NA      0.00        TRUE
#> 6 FALSE abaixo_do_minimo           5          NA     -0.39        TRUE
#>   turbidez_status turbidez__lim_min turbidez__lim_max turbidez__delta dbo_ok
#> 1              ok                NA               100               0  FALSE
#> 2              ok                NA               100               0   TRUE
#> 3              ok                NA               100               0   TRUE
#> 4              ok                NA               100               0   TRUE
#> 5              ok                NA               100               0   TRUE
#> 6              ok                NA               100               0   TRUE
#>        dbo_status dbo__lim_min dbo__lim_max dbo__delta coliformes_ok
#> 1 acima_do_maximo           NA            5       0.38          TRUE
#> 2              ok           NA            5       0.00          TRUE
#> 3              ok           NA            5       0.00         FALSE
#> 4              ok           NA            5       0.00          TRUE
#> 5              ok           NA            5       0.00          TRUE
#> 6              ok           NA            5       0.00         FALSE
#>   coliformes_status coliformes__lim_min coliformes__lim_max coliformes__delta
#> 1                ok                  NA                1000                 0
#> 2                ok                  NA                1000                 0
#> 3   acima_do_maximo                  NA                1000              2712
#> 4                ok                  NA                1000                 0
#> 5                ok                  NA                1000                 0
#> 6   acima_do_maximo                  NA                1000              2930
#>   p_total_ok  p_total_status p_total__lim_min p_total__lim_max p_total__delta
#> 1      FALSE acima_do_maximo               NA             0.05          0.030
#> 2      FALSE acima_do_maximo               NA             0.05          0.298
#> 3      FALSE acima_do_maximo               NA             0.05          0.263
#> 4      FALSE acima_do_maximo               NA             0.05          0.397
#> 5      FALSE acima_do_maximo               NA             0.05          0.290
#> 6      FALSE acima_do_maximo               NA             0.05          0.321

Summaries and plots

plot_iqa(df)
#> Warning: Removed 16 rows containing missing values or values outside the scale range
#> (`geom_col()`).

conf_long <- conama_summary(df, classe = "2")
head(conf_long)
#> # A tibble: 6 × 7
#>   parametro valor lim_min lim_max status ok    delta
#>   <chr>     <dbl>   <dbl>   <dbl> <chr>  <lgl> <dbl>
#> 1 ph         7.08       6       9 ok     TRUE      0
#> 2 ph         6.79       6       9 ok     TRUE      0
#> 3 ph         8.43       6       9 ok     TRUE      0
#> 4 ph         8.3        6       9 ok     TRUE      0
#> 5 ph         7.88       6       9 ok     TRUE      0
#> 6 ph         8.09       6       9 ok     TRUE      0

Render a report (writes to a temporary directory)

# Avoid running nested rendering during vignette build on CRAN
out_file <- render_report(df)
out_file