Coleção modelos
Duncan
Teste de Duncan (amplitude múltipla): letras de agrupamento das médias.
O que o bloco faz
models/duncan Aplica o teste de Duncan às médias do tratamento e retorna grupos por letras. A saída é models/emm.
Antes de continuar
Termos que esta página usa, em poucas palavras. As definições completas estão no glossário.
- Intervalo de confiança
- Uma faixa de valores calculada a partir da amostra. Com 95% de confiança, o procedimento produz intervalos que contêm o valor verdadeiro em 95% das amostras possíveis; não quer dizer que um intervalo específico tenha 95% de probabilidade de acertar.
- Modelo linear generalizado (GLM)
- Uma extensão do modelo linear para respostas que não são normais, como contagens (Poisson) e proporções ou sim/não (binomial). Uma função de ligação conecta a média da resposta à parte linear do modelo.
- Fator
- Uma variável explicativa categórica cujo efeito se quer estudar, como adubo, cultivar ou dose. Num experimento fatorial, dois ou mais fatores são combinados.
- Nível
- Cada valor que um fator assume no estudo: as doses 0, 50 e 100 kg/ha são três níveis do fator dose.
- Comparações múltiplas
- Comparar muitos pares de tratamentos aumenta a chance de achar alguma diferença só por acaso. Os procedimentos de comparações múltiplas (Tukey, Dunnett, Scott-Knott...) controlam essa chance de formas diferentes, e por isso podem discordar entre si.
Quando usar
Use em experimentos de tratamento quando o plano analítico especifica o procedimento de Duncan para agrupar médias.
Pressupostos
O modelo que entra atende aos seus pressupostos (erros normais, variância constante, independência): as comparações usam o erro padrão dele.
Verificar:Normalidade dos resíduosLeveneDiagnóstico dos resíduos
Se falhar: Acerte o modelo antes: transformação,
models/glm, oumodels/lmerpara erros correlacionados.As médias vêm de uma ANOVA com um erro só por fator (o bloco escolhe o erro (a) ou (b) na parcela subdividida), e os dados são balanceados: as médias usadas são as da tabela.
Se falhar: No desbalanceado, ou em GLM e misto, use o
models/emmeans.Sem interação relevante com outro fator; havendo, a média de um fator somada sobre o outro esconde o que aconteceu.
Verificar:Quadro da ANOVA
Se falhar: Compare dentro de cada nível do outro fator (Por no
models/emmeans).Aceita-se uma taxa de erro por família maior que a nominal: o nível de proteção cai com o número de médias entre as comparadas.
Se falhar: Para controlar o erro da família, use o Tukey do
models/emmeans.
Configuração
Tratamento escolhe o fator; Confiança (padrão 0,95) define o nível de decisão, com alfa = 1 − confiança.
Exemplo
library(trama)
reg <- tr_registry()
tr_use("trama.models", registry = reg)
tr_flow(reg) |>
tr_add("dados", "models/example", dataset = "milho_dbc") |>
tr_add("ajuste", "models/anova_dbc", resposta = "producao", tratamento = "hibrido", bloco = "bloco", from = "dados") |>
tr_add("resultado", "models/duncan", tratamento = "hibrido", from = "ajuste")O DBC de milho fornece médias de produção dos híbridos; as letras de Duncan indicam quais híbridos o procedimento não separa ao nível de confiança usado.
Como interpretar
Letras comuns indicam médias não separadas pelo procedimento. Duncan é menos conservador que Tukey no controle de erro familiar.
Implementação
agricolae::duncan.test()versão 1.3.7Com o QM e os gl do erro do quadro (DFerror,MSerror) egroup = TRUE.
No código O bloco chama trama.models::tr_models_duncan(), em trama.models/R/agricolae.R:92.
A conta é feita nesta chamada:
agricolae::duncan.test(b$y, b$trt, DFerror = b$gl, MSerror = b$qm, alpha = alfa, group = TRUE, console = FALSE)
Outras funções de pacotes que o bloco usa (34)
stats::anova(aj)reformulas::nobars(stats::formula(aj))stats::formula(aj)stats::setNames(rep(list("contr.sum"), length(usados)), usados)MASS::glm.nb(stats::formula(aj), data = fit$dados, contrasts = if (length(ctr)) ctr else NULL)stats::formula(aj)stats::glm(stats::formula(aj), family = stats::family(aj), data = fit$dados, contrasts = if (length(ctr)) ctr else NULL)stats::formula(aj)stats::family(aj)lme4::glmer(stats::formula(aj), data = as.data.frame(fit$dados), family = stats::family(aj), contrasts = if (length(ctr)) ctr else NULL)stats::formula(aj)stats::family(aj)stats::contrasts(d[[v]])stats::contr.sum(nlevels(d[[v]]))stats::lm(stats::formula(aj), data = fit$dados, contrasts = if (length(ctr)) ctr else NULL)car::Anova(m, type = if (tipo == "II") 2 else 3)stats::family(aj)stats::anova(aj, test = if (usa_f) "F" else "Chisq")car::Anova(m, type = if (tipo == "II") 2 else 3, test.statistic = if (usa_f) "F" else "LR")stats::anova(modelo$ajuste, type = tipo)stats::anova(m, type = if (tipo == "I") "sequential" else "marginal")car::Anova(m, type = tipo)stats::anova(fit$ajuste)stats::pf(tab$F[i], tab$gl[i], gl_a, lower.tail = FALSE)car::Anova(modelo$ajuste, type = 2)reformulas::nobars(stats::as.formula(fit$formula))stats::as.formula(fit$formula)stats::terms(fixa)stats::as.formula(paste("~", t))stats::anova(aj)car::Anova(aj, type = 2)stats::residuals(fit$ajuste)stats::df.residual(fit$ajuste)stats::qt(1 - alfa/2, b$gl)
Referências
Teoria
- Duncan, D. B. (1955). Multiple range and multiple F tests. Biometrics, 11(1), 1-42. doi:10.2307/3001478
Livro-texto
- Banzatto, D. A.; Kronka, S. N. (2006). Experimentação agrícola. 4. ed. Jaboticabal: Funep.