Coleção modelos
Levene
Levene (Brown-Forsythe): as variâncias dos resíduos são iguais entre os tratamentos?
O que o bloco faz
models/levene Compara variâncias residuais entre grupos de efeitos fixos. 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.
- 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.
- 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.
- 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 avaliar igualdade de variâncias entre grupos de tratamento a partir dos resíduos do ajuste. 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 grupos (as combinações dos tratamentos) têm pelo menos duas observações cada; o bloco recusa senão.
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.Os resíduos vêm de um modelo bem especificado (termos certos, sem tendência sobrando): o teste olha os resíduos, e um modelo errado muda o que eles medem.
Verificar:Diagnóstico dos resíduos
Com bloco (DBC, fatorial em DBC, DQL), o delineamento é equilibrado (cada tratamento o mesmo número de vezes em cada bloco, linha e coluna): a correção de O'Neill & Mathews é um fator do desenho que supõe o equilíbrio; o bloco recusa senão. Ali o centro é o ajuste de mínimos quadrados (a média), qualquer que seja o param Centro.
Se falhar: Leia o painel escala-locação do
models/plot_diagnostics.Com bloco, o desenho não é pequeno demais: a correção acerta a média do F, não a cauda, e em desenho pequeno o teste é conservador (tamanho a 5%, por simulação: DBC 5 × 6 4,6%, 4 × 3 2,9%; DQL 8 × 8 4,7%, 5 × 5 2,9%, 4 × 4 2,3%).
Se falhar: Não rejeitar num desenho pequeno é pouca evidência de variâncias iguais: leia também o painel escala-locação do
models/plot_diagnostics.O modelo não é parcela subdividida: os resíduos de dois estratos de erro não têm correção publicada (a de O'Neill & Mathews supõe um estrato), e o bloco recusa.
Se falhar: Leia o painel escala-locação do
models/plot_diagnostics; com variâncias diferentes, ajuste o misto nomodels/lmer.
Configuração
Centro mediano é Brown–Forsythe (padrão); centro média é Levene original.
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/levene", centro = "mediana", from = "ajuste")No ajuste DIC de PlantGrowth, Levene compara a dispersão dos resíduos entre ctrl, trt1 e trt2, usando a mediana como centro do teste.
Como interpretar
A hipótese nula é igualdade de variâncias. Sem fatores, use Breusch–Pagan para preditores contínuos.
Implementação
car::leveneTest()versão 3.1.5Sem bloco: nos resíduos do modelo, agrupados pelos tratamentos.center = median(padrão do bloco) é a versão de Brown-Forsythe;center = meané o Levene original. Com bloco (DBC, DQL) a conta é própria (O'Neill & Mathews 2002), validada contra ooneilldbcdo ExpDes.pt.
No código O bloco chama trama.models::tr_models_levene(), em trama.models/R/pressupostos.R:160.
A conta é feita nesta chamada:
car::leveneTest(r$residuo, g$g, center = if (centro == "mediana") stats::median else mean)
Outras funções de pacotes que o bloco usa (19)
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)stats::as.formula(fit$formula)stats::model.matrix(stats::reformulate(.tr_models_bt(ctrl)), dd)stats::reformulate(.tr_models_bt(ctrl))stats::model.matrix(stats::reformulate(c(.tr_models_bt(ctrl), ".g")), dd)stats::reformulate(c(.tr_models_bt(ctrl), ".g"))proj(Xc)stats::pf(f, glt, glr, lower.tail = FALSE)
Referências
Teoria
- O'Neill, M. E.; Mathews, K. L. (2002). Levene tests of homogeneity of variance for general block and treatment designs. Biometrics, 58(1), 216-224. doi:10.1111/j.0006-341x.2002.00216.x
- O'Neill, M. E.; Mathews, K. (2000). A weighted least squares approach to Levene's test of homogeneity of variance. Australian & New Zealand Journal of Statistics, 42(1), 81-100. doi:10.1111/1467-842x.00109
- Levene, H. (1960). Robust tests for equality of variances. In: Olkin, I. et al. (ed.). Contributions to Probability and Statistics: Essays in Honor of Harold Hotelling. Stanford: Stanford University Press, p. 278-292.
- Brown, M. B.; Forsythe, A. B. (1974). Robust tests for the equality of variances. Journal of the American Statistical Association, 69(346), 364-367. doi:10.1080/01621459.1974.10482955
Livro-texto
- Montgomery, D. C. (2017). Design and Analysis of Experiments. 9. ed. Hoboken: Wiley.