Coleção modelos
Bartlett
Bartlett: as variâncias dos resíduos são iguais entre os tratamentos?
O que o bloco faz
models/bartlett Testa igualdade de variâncias residuais entre grupos. 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.
- Variância
- A média dos quadrados das distâncias de cada valor à média: mede o quanto os dados se espalham. O desvio padrão é a raiz quadrada da variância e fica na mesma unidade dos dados.
- Resíduo
- A diferença entre o valor observado e o valor que o modelo ajustado prevê para aquela observação. É nos resíduos que se conferem os pressupostos do modelo: eles devem parecer ruído, sem padrão.
- 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.
- Homocedasticidade
- Pressuposto de que a variância dos erros é a mesma em todos os tratamentos ou ao longo dos valores ajustados. Quando falha (heterocedasticidade), os erros padrão e os testes deixam de ser confiáveis.
Quando usar
Use para testar igualdade de variâncias entre grupos quando a normalidade residual é plausível. Na parcela subdividida o bloco recusa: os resíduos vêm de dois estratos de erro e não há correção publicada; leia o painel escala-locação do models/plot_diagnostics.
Pressupostos
Os resíduos de cada grupo são normais: o teste é muito sensível à falta de normalidade, e rejeita pela cauda, não pela variância.
Verificar:Normalidade dos resíduosDiagnóstico dos resíduos
Se falhar: Use o
models/levene(centro na mediana), robusto à falta de normalidade.Os grupos têm pelo menos duas observações cada.
O modelo não tem bloco (DIC, fatorial em DIC, modelos de fórmula): nos resíduos correlacionados de um DBC, DQL ou parcela subdividida o Bartlett não tem correção publicada, e o bloco recusa.
Se falhar: Use o
models/levene, que no delineamento com bloco aplica a correção de O'Neill & Mathews (2002).As observações de cada grupo são independentes entre si.
Se falhar: Não há teste no trama para isso: é o delineamento (sorteio, uma medida por unidade) que garante. Medidas repetidas na mesma unidade pedem
models/paired_toumodels/lmer.
Configuração
Não há parâmetros adicionais.
Exemplo
library(trama)
reg <- tr_registry()
tr_use("trama.models", registry = reg)
tr_flow(reg) |>
tr_add("dados", "models/example", dataset = "PlantGrowth") |>
tr_add("ajuste", "models/anova_dic", resposta = "weight", tratamento = "group", from = "dados") |>
tr_add("resultado", "models/bartlett", from = "ajuste")O teste recebe os resíduos do DIC de PlantGrowth e compara suas variâncias nos três grupos; o resultado informa qui-quadrado e p-valor.
Como interpretar
A hipótese nula é variâncias iguais. Bartlett é sensível à não normalidade; confira os resíduos.
Implementação
stats::bartlett.test()versão 4.6.0Nos resíduos do modelo, agrupados pelos tratamentos.
No código O bloco chama trama.models::tr_models_bartlett(), em trama.models/R/pressupostos.R:197.
A conta é feita nesta chamada:
stats::bartlett.test(r$residuo, g$g)
Outras funções de pacotes que o bloco usa (12)
stats::fitted(aj)stats::residuals(aj)stats::sigma(aj)stats::fitted(aj)stats::residuals(aj, type = "normalized")stats::fitted(aj)stats::residuals(aj)stats::rstandard(aj)reformulas::nobars(stats::as.formula(fit$formula))stats::as.formula(fit$formula)reformulas::findbars(stats::as.formula(fit$formula))stats::as.formula(fit$formula)
Referências
Teoria
- Bartlett, M. S. (1937). Properties of sufficiency and statistical tests. Proceedings of the Royal Society of London. Series A, 160(901), 268-282. doi:10.1098/rspa.1937.0109
Livro-texto
- Montgomery, D. C. (2017). Design and Analysis of Experiments. 9. ed. Hoboken: Wiley.