Co to jest test Bartletta jednorodności wariancji?
Test Bartletta (ang. Bartlett’s test of homogeneity of variances) sprawdza, czy wariancje zmiennej w dwóch lub więcej niezależnych grupach są równe. Hipoteza zerowa brzmi σ₁² = σ₂² = … = σₖ², czyli zakłada homogeniczność wariancji. Źródłem testu jest praca Bartletta (1937) o testach dla małych prób. Nie myl go z testem sferyczności Bartletta, który przed analizą czynnikową sprawdza, czy macierz korelacji jest jednostkowa.
Jak działa test?
Statystyka K² porównuje logarytm wariancji połączonej z ważoną sumą logarytmów wariancji w grupach, z poprawką na liczebności. Gdy wariancje są równe, a dane w grupach mają rozkład normalny, K² ma w przybliżeniu rozkład chi-kwadrat z k − 1 stopniami swobody. Właśnie to założenie normalności jest słabym punktem: test myli grube ogony rozkładu z różnicą wariancji. Gastwirth i in. (2009) przypominają, że wiele testów równości wariancji zakłada normalność i nie jest odpornych na jej brak, a w porównaniu Conovera i in. (1981) odporne i zarazem mocne okazało się tylko kilka testów.
Przykład w R
Dla trzech grup po 30 osób o odchyleniach standardowych 8, 10 i 14 test Bartletta dał K²(2) = 17.56, p < .001, a test Levene’a z medianą F(2, 87) = 9.45, p < .001. W symulacji z równymi wariancjami oba testy przy danych normalnych odrzucały hipotezę zerową w około 5% prób (5.3% i 4.9%). Przy rozkładzie t z 5 stopniami swobody test Bartletta dawał fałszywy alarm w 24.8% prób, a test Levene’a w 5.2%.
Jak opisać w pracy
Na co uważać
- Nienormalne dane. Przy skośnych lub grubo-ogonowych rozkładach wybierz test Levene’a z medianą (Browna–Forsythe’a).
- Test wstępny. Nie wybieraj na jego podstawie między ANOVA klasyczną a ANOVA Welcha; Welcha możesz stosować od razu. Oba warianty ze średnich i SD policzysz w kalkulatorze testu t i ANOVA.
- Duże próby. Istotny wynik nie mówi, jak duża jest różnica; podaj odchylenia standardowe w grupach.
Test Bartletta w R
Kod porównuje test Bartletta i test Levene’a z medianą na danych z różnymi wariancjami, a potem w 2000 symulacjach z równymi wariancjami sprawdza, jak często każdy test daje fałszywy alarm przy danych normalnych i grubo-ogonowych.
Test Bartletta a test Levene’a przy danych normalnych i nienormalnych
R
library(car)
set.seed(4361)
g <- factor(rep(c("A", "B", "C"), each = 30))
y <- rnorm(90, mean = 50, sd = rep(c(8, 10, 14), each = 30))
bartlett.test(y ~ g)
leveneTest(y ~ g, center = median)
# Odsetek fałszywych alarmów przy RÓWNYCH wariancjach (α = .05)
alarmy <- function(los) rowMeans(replicate(2000, {
y0 <- los(90)
c(Bartlett = bartlett.test(y0 ~ g)$p.value < .05,
Levene_mediana = leveneTest(y0 ~ g, center = median)[1, 3] < .05)
}))
round(rbind(normalny = alarmy(rnorm),
t_5_df = alarmy(function(n) rt(n, df = 5))), 3)Wynik w R
Bartlett test of homogeneity of variances
data: y by g
Bartlett's K-squared = 17.564, df = 2, p-value = 0.0001535
Levene's Test for Homogeneity of Variance (center = median)
Df F value Pr(>F)
group 2 9.4516 0.000193 ***
87
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Bartlett Levene_mediana
normalny 0.053 0.049
t_5_df 0.248 0.052Bibliografia oprogramowania
Wyniki pochodzą z uruchomienia pokazanego kodu w R 4.5.3 (z pakietem car). W pracy podaj wersje faktycznie użytych pakietów (sessionInfo()).
- R: R Core Team. (2026). R: A language and environment for statistical computing (Version 4.5.3) [Computer software]. R Foundation for Statistical Computing. https://www.R-project.org/
- car (wersja 3.1-5): Fox, J., & Weisberg, S. (2019). An R companion to applied regression (3rd ed.). Sage. https://www.john-fox.ca/Companion/
Najczęstsze pytania
Czym test Bartletta różni się od testu Levene’a?
Test Bartletta porównuje logarytmy wariancji i zakłada normalność, więc grube ogony dają fałszywe alarmy. Test Levene’a analizuje odchylenia od średniej lub mediany; wersja z medianą jest znacznie odporniejsza na nienormalność.
Czy test Bartletta to ten sam test co test sferyczności Bartletta?
Nie. Test jednorodności wariancji porównuje grupy, a test sferyczności sprawdza, czy zmienne w ogóle ze sobą korelują przed analizą czynnikową. Łączy je tylko autor.
Jak wykonać test Bartletta w R?
Funkcją bartlett.test(y ~ grupa) z pakietu stats. Wynik to Bartlett’s K-squared, df = k − 1 i p; w pracy zapisz go jako χ²(df) = …, p = ….
Powiązane hasła
Poradniki, w których użyjesz tego pojęcia
Potrzebujesz wsparcia w analizie danych?
Wybierz ofertę dopasowaną do pracy dyplomowej albo profesjonalnego projektu naukowego.
Wyceń analizęOferta dla studentówBibliografia i dalsza literatura naukowa (APA 7)
- Bartlett, M. S. (1937). Properties of sufficiency and statistical tests. Proceedings of the Royal Society of London. Series A – Mathematical and Physical Sciences, 160(901), 268–282. https://doi.org/10.1098/rspa.1937.0109
- Conover, W. J., Johnson, M. E., & Johnson, M. M. (1981). A comparative study of tests for homogeneity of variances, with applications to the outer continental shelf bidding data. Technometrics, 23(4), 351–361. https://doi.org/10.1080/00401706.1981.10487680
- Gastwirth, J. L., Gel, Y. R., & Miao, W. (2009). The impact of Levene’s test of equality of variances on statistical theory and practice. Statistical Science, 24(3), 343–360. https://doi.org/10.1214/09-sts301 ● Otwarty dostęp
Cytowania BibTeX
@article{bartlett1937properties, author={Bartlett, Maurice Stevenson}, title={{Properties of sufficiency and statistical tests}}, journal={Proceedings of the Royal Society of London. Series A – Mathematical and Physical Sciences}, year={1937}, volume={160}, number={901}, pages={268--282}, doi={10.1098/rspa.1937.0109}}@article{conover1981comparative, author={Conover, W. J. and Johnson, Mark E. and Johnson, Myrle M.}, title={{A comparative study of tests for homogeneity of variances, with applications to the outer continental shelf bidding data}}, journal={Technometrics}, year={1981}, volume={23}, number={4}, pages={351--361}, doi={10.1080/00401706.1981.10487680}}@article{gastwirth2009impact, author={Gastwirth, Joseph L. and Gel, Yulia R. and Miao, Weiwen}, title={{The impact of Levene’s test of equality of variances on statistical theory and practice}}, journal={Statistical Science}, year={2009}, volume={24}, number={3}, pages={343--360}, doi={10.1214/09-sts301}}