Computa o IQA/WQI combinando subindices (Qi) por media geometrica ponderada, conforme a metodologia oficial CETESB e o NSF WQI original (Brown et al., 1970): \(IQA = \prod_{i} Qi_i^{Wi}\).
iqa(
df,
pesos = c(od = 0.17, coliformes = 0.15, dbo = 0.1, nt_total = 0.1, p_total = 0.1,
turbidez = 0.08, tds = 0.08, pH = 0.12, temperatura = 0.1),
method = c("CETESB", "CETESB_equations", "NSF_approx"),
altitude_m = 0,
na_rm = FALSE,
add_status = TRUE,
locale = c("pt", "en"),
...
)O df de entrada com a coluna numerica IQA (0-100) e,
quando add_status = TRUE, a coluna fator IQA_status.
O atributo "iqa_method" e definido no objeto retornado.
Data frame (ou tibble) com as colunas requeridas.
Nomes esperados: od, coliformes, dbo,
nt_total, p_total, turbidez, tds,
ph (ou pH), temperatura (ou temp).
Pesos nomeados para cada parametro. Padroes seguem pratica
CETESB/NSF: od=.17, coliformes=.15, dbo=.10,
nt_total=.10, p_total=.10, turbidez=.08,
tds=.08, pH=.12, temperatura=.10.
Metodo de calculo:
"CETESB" (padrao) — subindices por curvas de interpolacao
+ media geometrica ponderada.
"CETESB_equations" — equacoes polinomiais CETESB com
saturacao de OD dependente de temperatura e altitude +
media geometrica ponderada.
"NSF_approx" — subindices por curvas + media aritmetica
ponderada (metodo legado, mantido para compatibilidade).
Altitude em metros acima do nivel do mar (default 0).
Usado apenas em method = "CETESB_equations" para correcao da
saturacao de oxigenio dissolvido.
Logico; se FALSE (padrao), linhas com Qi ausentes
geram erro. Se TRUE, o IQA e calculado com os parametros
disponiveis e os pesos sao renormalizados por linha.
Logico; se TRUE (padrao), adiciona a coluna
IQA_status com a classificacao qualitativa (0-100).
Idioma de IQA_status: "pt" (padrao) ou
"en".
Reservado para uso futuro.
Metodo de agregacao (correcao em v0.9.0):
O IQA CETESB e o NSF WQI original utilizam media geometrica ponderada.
O metodo "NSF_approx" (media aritmetica) e mantido apenas para
compatibilidade retroativa.
Compatibilidade de nomes de coluna:
ph (minusculo) e aceito como alias de pH.
temp e aceito como alias de temperatura.
CETESB (2021). Qualidade das Aguas Superficiais no Estado de Sao Paulo. CETESB, Sao Paulo.
Brown, R.M. et al. (1970). A Water Quality Index — Do We Dare? Water and Sewage Works, 117, 339-343.
d <- iqa(wq_demo, na_rm = TRUE)
table(d$IQA_status, useNA = "ifany")
# Usando equacoes CETESB com correcao de altitude
# \donttest{
d2 <- iqa(wq_demo, method = "CETESB_equations", altitude_m = 800, na_rm = TRUE)
summary(d2$IQA)
# }
Run the code above in your browser using DataLab