Coleção modelos
Aditividade de Tukey
Teste de não aditividade de Tukey: bloco e tratamento interagem?
O que o bloco faz
models/tukey_additivity Testa aditividade entre bloco e tratamento. A saída é data/test.
Antes de continuar
Termos que esta página usa, em poucas palavras. As definições completas estão no glossário.
- Tratamento
- A condição aplicada a uma unidade experimental: um nível de um fator ou uma combinação de níveis de vários fatores. As comparações do experimento são entre tratamentos.
- Interação
- Há interação quando o efeito de um fator depende do nível de outro: o adubo aumenta a produção numa cultivar e não em outra. Com interação, olhar o efeito médio de cada fator sozinho pode enganar.
- 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.
- Pressuposto
- Uma condição sob a qual o método dá resultados confiáveis, como resíduos normais ou observações independentes. Quando um pressuposto falha, p-valores e intervalos podem ficar errados, mesmo que a conta rode.
- Normalidade
- Pressuposto de que os erros do modelo seguem a distribuição normal (em forma de sino). Em modelos lineares ele se confere nos resíduos, não nos dados brutos; desvios leves pesam pouco em amostras grandes.
Quando usar
Use em DBC para avaliar se tratamento e bloco atuam aditivamente.
Pressupostos
A não aditividade, se existe, tem a forma multiplicativa que o teste procura (efeito proporcional ao produto dos efeitos de bloco e tratamento); outras formas de interação podem passar.
Os erros do modelo aditivo são normais e homocedásticos: o teste é um F.
Verificar:Normalidade dos resíduosLevene
Configuração
Não há parâmetros adicionais; requer ajuste de blocos casualizados.
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/tukey_additivity", from = "ajuste")No DBC simulado, o teste verifica o termo de não aditividade entre hibrido e bloco; o resultado fornece estatística e p-valor para essa hipótese.
Como interpretar
A hipótese nula é ausência de interação bloco × tratamento.
Implementação
trama.models::tr_models_tukey_additivity()versão 0.6.6Cálculo próprio comstats::lmestats::anova: acrescenta o quadrado do valor ajustado ao modelo aditivo (bloco + tratamento, ou linha + coluna + tratamento no DQL) e testa esse termo com 1 gl.
No código O bloco chama trama.models::tr_models_tukey_additivity(), em trama.models/R/pressupostos.R:266.
A conta é implementada pelo próprio trama, sem um pacote R de referência por trás.
Outras funções de pacotes que o bloco usa (6)
stats::as.formula(modelo$formula)stats::as.formula(paste(.tr_models_bt(modelo$resposta), "~", paste(c(.tr_models_bt(controles), ".trat"), collapse = " + ")))stats::lm(f0, data = d)stats::fitted(m0)stats::update(m0, . ~ . + .q, data = d)stats::anova(m0, m1)
Referências
Teoria
- Tukey, J. W. (1949). One degree of freedom for non-additivity. Biometrics, 5(3), 232-242. doi:10.2307/3001938
Livro-texto
- Montgomery, D. C. (2017). Design and Analysis of Experiments. 9. ed. Hoboken: Wiley.