Przedział ufności
![Wykres przedstawiający sto przedziałów ufności; na osi X oszacowanie. Zdecydowana większość przedziałów pokrywa prawdziwej średniej, ale różnią się pomiędzy sobą szerokością i położeniem.](http://upload.wikimedia.org/wikipedia/commons/thumb/3/31/Confidence_intervall_normal_dist.svg/220px-Confidence_intervall_normal_dist.svg.png)
Przedział ufności – podstawowe narzędzie estymacji przedziałowej. Pojęcie to zostało wprowadzone do statystyki przez polskiego matematyka Jerzego Spławę-Neymana. Występuje w wielu wariantach, w klasycznym wąskim rozumieniu opiera się na błędzie standardowym. Szczególny przypadek przedziału ufności w badaniach ankietowych jest zwyczajowo określany marginesem błędu[1].
Definicja
Niech cecha ma rozkład w populacji z nieznanym parametrem Z populacji wybieramy próbę losową Przedziałem ufności o współczynniku ufności nazywamy taki przedział który spełnia warunek:
gdzie i są funkcjami wyznaczonymi na podstawie próby losowej.
Podobnie jak w przypadku estymatorów definicja pozwala na dowolność wyboru funkcji z próby, jednak tutaj kryterium wyboru najlepszych funkcji narzuca się automatycznie – zazwyczaj będziemy poszukiwać przedziałów najkrótszych.
Współczynnik ufności jest wielkością, którą można interpretować w następujący sposób: jest to prawdopodobieństwo wyznaczenia takiego przedziału, że rzeczywista wartość parametru w populacji znajdzie się w tym przedziale. Im większa wartość tego współczynnika, tym szerszy przedział ufności, a więc mniejsza dokładność estymacji parametru. Im mniejsza wartość tym większa dokładność estymacji, ale jednocześnie tym większe prawdopodobieństwo popełnienia błędu. Wybór odpowiedniego współczynnika jest więc kompromisem pomiędzy dokładnością estymacji a ryzykiem błędu. W praktyce przyjmuje się zazwyczaj wartości: 0,99; 0,95 lub 0,90, zależnie od parametru.
Interpretacja
![Wykres punktowy przedstawiający oszacowania prędkości światła z przedziałami ufności (na osi Y) w kolejnych latach między 1874 a 1979 (na osi X). Duża liczba przedziałów nie pokrywa przyjętej po 1983, prawdziwej z definicji wartości.](http://upload.wikimedia.org/wikipedia/commons/thumb/2/21/Historical_estimates_of_the_speed_of_light_in_vacuum_with_standard_error_bars.svg/220px-Historical_estimates_of_the_speed_of_light_in_vacuum_with_standard_error_bars.svg.png)
Przedziały ufności konstruowane w ramach paradygmatu częstościowego weryfikacji hipotez statystycznych niosą ze sobą trudności interpretacyjne typowe dla tego podejścia. Metodologia ta nie była tworzona jako narzędzie służące do jednorazowego szacowania parametrów, lecz długoterminowej kontroli błędów, co jest z reguły kontrintuicyjne dla użytkowników. Zdaniem niektórych statystyków preferujących paradygmat bayesowski, jak Hoekstra i in., lub Gelman, prawidłowa interpretacja przedziałów ufności nie uprawnia w ścisłym sensie do przyjęcia, że z daną „ufnością” prawdziwa wartość parametru znajduje się w zaobserwowanym przedziale. Przy wielokrotnym powtarzaniu eksperymentu, 95% przedziałów ufności na poziomie 95% powinno w idealnych warunkach zawierać rzeczywistą wartość parametru, ale nie przekłada się to na równoważną pewność co do żadnego konkretnego przedziału, ani tym bardziej – średniej na której jest on skonstruowany. Do wniosków takich w pełni uprawniają na przykład metody bayesowskie[3][4][5][6].
Dla ilustracji, osoba, która kierowałaby się przedziałami ufności w grze hazardowej, w przypadku wielu rodzajów danych przegrywałaby statystycznie częściej, niż sugerują nominalne właściwości tych przedziałów[7]. Według badań z 2004 r. większość naukowców jest przekonana, że średnie arytmetyczne uzyskane dla kolejnych prób losowych będą się zawierały w skonstruowanym 95% przedziale ufności z częstotliwością 95%. W rzeczywistości kolejne średnie z próbek zawierają się w skonstruowanym przedziale z prawdopodobieństwem ok. 83,4% dla próbek z rozkładu normalnego[8]. W wielu przypadkach przedziały skonstruowane z użyciem statystyk odpornościowych są bardziej niezawodne, ale i w tym przypadku wykryto wyjątki[9].
Z przyjęciem tej nieintuicyjnej i zniuansowanej interpretacji ma jednak problem większość osób, nawet badaczy-praktyków[3][6]. Dodatkowo, w praktyce dobór prób i pomiary nie są nigdy realizowane w idealnej zgodności z założeniami modeli statystycznych. Ogólnie rzecz biorąc, zaleca się traktowanie nominalnego poziomu ufności jako przeszacowującego epistemologiczną pewność wyników[10], oraz używanie statystyk bayesowskich do zastosowań ściśle epistemologicznych[3].
Przykłady przedziałów ufności
Ponieważ szukamy jak najkrótszych przedziałów ufności, dlatego przy wyznaczaniu przedziału staramy się wykorzystać jak najwięcej dostępnych informacji o rozkładzie cechy w populacji. Jeśli np. cecha ma rozkład normalny z odchyleniem standardowym to zastosowanie wzoru na przedział ufności dla nieznanego również da poprawny wynik, jednak przedział otrzymany tą metodą będzie szerszy, czyli mniej dokładny. Z kolei wzory ogólniejsze, np. dla nieznanego rozkładu, często korzystają z rozkładów granicznych estymatorów i dlatego wymagają dużej liczebności próby.
Przedział ufności dla średniej
Znane odchylenie standardowe
Cecha ma w populacji rozkład normalny przy czym odchylenie standardowe jest znane. Przedział ufności dla parametru tego rozkładu ma postać:
lub równoważnie
gdzie:
![]() | Zobacz w Wikiźródłach tablicę rozkładu normalnego |
- – liczebność próby losowej,
- – średnia z próby losowej,
- – odchylenie standardowe populacji,
- to kwantyl rzędu standaryzowanego rozkładu normalnego
Nieznane odchylenie standardowe
Cecha ma w populacji rozkład normalny przy czym odchylenie standardowe jest nieznane. Przedział ufności dla parametru tego rozkładu ma postać:
gdzie:
![]() | Zobacz w Wikiźródłach tablicę rozkładu Studenta |
- – liczebność próby losowej,
- – średnia z próby losowej,
- – odchylenie standardowe z próby,
- to kwantyl rzędu rozkładu Studenta z stopniami swobody.
Zwykle stosuje się ten wzór dla małej próby Tak naprawdę działa on dla każdej wielkości próby, jednak dla dużych prób można przybliżyć rozkład t Studenta rozkładem normalnym, co jest łatwiejsze do wyliczenia a dające niemal takie same wartości (patrz niżej).
Nieznane odchylenie standardowe – duża próba (n > 30)
Cecha ma w populacji rozkład normalny przy czym odchylenie standardowe jest nieznane, a próba jest duża (n>30). Granica 30 jest czysto umowna, im jest większe, tym wzór dokładniejszy. Przedział ufności dla parametru m tego rozkładu ma postać:
gdzie:
- – liczebność próby losowej,
- – średnia z próby losowej,
- – odchylenie standardowe z próby,
- – kwantyl rzędu standaryzowanego rozkładu normalnego
Przedział ufności dla wariancji
Poniższy wzór pozwala wyznaczyć przedział ufności dla wariancji w populacji o rozkładzie normalnym
gdzie:
- – liczebność próby losowej,
- – odchylenie standardowe z próby,
- i – liczby będące kwantylami rozkładu chi-kwadrat spełniające odpowiednio nierówności:
![]() | Zobacz w Wikiźródłach tablicę rozkładu chi-kwadrat |
gdzie ma rozkład chi-kwadrat z stopniami swobody.
Podobnie jak poprzednio, zwykle stosuje się ten wzór dla małej próby (n < 30), choć również działa on dla każdej wielkości próby.
Duża próba (n > 30)
Poniższy wzór pozwala wyznaczyć przedział ufności dla wariancji w populacji o rozkładzie normalnym dla dużej próby, czyli umownie dla
gdzie:
- – liczebność próby losowej,
- – odchylenie standardowe z próby,
- – statystyka spełniająca warunek:
gdzie jest zmienną losową o rozkładzie normalnym
Przedział ufności dla odsetka (proporcji, wskaźnika struktury)
Poniższy wzór pozwala wyznaczyć przedział ufności dla wskaźnika struktury w populacji o rozkładzie normalnym
gdzie:
- – liczebność próby losowej,
- – liczebność wybranej grupy z próby,
- – statystyka spełniającą warunek:
- gdzie jest zmienną losową o rozkładzie normalnym
Przedział ufności dla współczynnika korelacji
Poniższy wzór pozwala wyznaczyć przedział ufności dla współczynnika korelacji w populacji o rozkładzie normalnym Tak jak poprzednio, działa on dla dowolnej próby, choć jest zwykle stosowany tylko dla prób małych
gdzie:
- – liczebność próby losowej,
- – statystyka spełniająca warunek:
- gdzie jest zmienną losową o rozkładzie normalnym
- – współczynnik korelacji.
Duża próba (n > 30)
Poniższy wzór pozwala wyznaczyć przedział ufności dla współczynnika korelacji w populacji o rozkładzie normalnym
gdzie:
- – liczebność próby losowej,
- jest statystyką spełniającą warunek:
- gdzie jest zmienną losową o rozkładzie normalnym
- – współczynnik korelacji.
Przedział ufności dla współczynnika
Poniższy wzór pozwala wyznaczyć przedział ufności dla współczynnika w populacji o rozkładzie normalnym
gdzie:
- – wartość z próby losowej,
- – średnia z próby losowej,
- ma rozkład Studenta z stopniami swobody.
Minimalna liczebność próby
Jeśli chcemy oszacować parametr z określoną dokładnością możemy – po odpowiednich przekształceniach wzorów na przedziały ufności – wyznaczyć liczebność próby losowej potrzebną do osiągnięcia zakładanej dokładności.
Przykład: Wiemy, że wzrost wikipedystów ma rozkład normalny z odchyleniem standardowym 25,28 cm (dane chyba nieprawdziwe). Obliczmy, ilu wikipedystów wystarczy zmierzyć, aby z prawdopodobieństwem 95% wyznaczyć średni wzrost wikipedysty z dokładnością do 5 cm.
Jeśli chcemy uzyskać dokładność 5 cm, należy zadbać o to, aby połowa długości przedziału ufności była mniejsza lub równa niż 5 cm. Ze wzoru na przedział ufności dla rozkładu normalnego o znanym odchyleniu standardowym wynika, że dokładność estymacji powinna spełniać zależność:
Przekształcamy podaną nierówność, uzyskując pożądany wzór na liczebność próby:
Podstawiając do wzoru wartości = 25,28; d = 5 cm; = 1,96 (wartość obliczona na podstawie tablic rozkładu normalnego), uzyskujemy minimalną wielkość próby na poziomie 99 wikipedystów.
Zobacz też
Przypisy
- ↑ DennisD. Gilliland DennisD., VinceV. Melfi VinceV., A Note on Confidence Interval Estimation and Margin of Error, „Journal of Statistics Education”, 18 (1), 2010, DOI: 10.1080/10691898.2010.11889474, ISSN 1069-1898 [dostęp 2019-03-31] (ang.).
- ↑ MaxM. Henrion MaxM., BaruchB. Fischhoff BaruchB., Assessing uncertainty in physical constants, „American Journal of Physics”, 54 (9), 1986, s. 791–798, DOI: 10.1119/1.14447, ISSN 0002-9505 [dostęp 2019-03-30] (ang.).
- ↑ a b c RinkR. Hoekstra RinkR. i inni, Robust misinterpretation of confidence intervals, „Psychonomic Bulletin & Review”, 5, 2014, s. 1157–1164, DOI: 10.3758/s13423-013-0572-3, ISSN 1531-5320, PMID: 24420726 [dostęp 2017-01-06] .
- ↑ Problematic interpretations of confidence intervals - Statistical Modeling, Causal Inference, and Social Science, „Statistical Modeling, Causal Inference, and Social Science”, 15 marca 2014 [dostęp 2017-01-06] (ang.).
- ↑ Abraham Lincoln and confidence intervals - Statistical Modeling, Causal Inference, and Social Science, „Statistical Modeling, Causal Inference, and Social Science”, 23 listopada 2016 [dostęp 2017-01-06] (ang.).
- ↑ a b SanderS. Greenland SanderS. i inni, Statistical tests, P values, confidence intervals, and power: a guide to misinterpretations, „European Journal of Epidemiology”, 31 (4), 2016, s. 337–350, DOI: 10.1007/s10654-016-0149-3, ISSN 0393-2990, PMID: 27209009, PMCID: PMC4877414 [dostęp 2017-02-09] (ang.).
- ↑ Ulrich K.U.K. Müller Ulrich K.U.K., AndriyA. Norets AndriyA., Credibility of Confidence Sets in Nonstandard Econometric Problems, „Econometrica”, 6, 2016, s. 2183–2213, DOI: 10.3982/ECTA14023, ISSN 1468-0262 [dostęp 2017-01-06] (ang.).
- ↑ GeoffG. Cumming GeoffG., JenniferJ. Williams JenniferJ., FionaF. Fidler FionaF., Replication and Researchers’ Understanding of Confidence Intervals and Standard Error Bars, „Understanding Statistics”, 3 (4), 2004, s. 299–311, DOI: 10.1207/s15328031us0304_5, ISSN 1534-844X [dostęp 2017-04-01] .
- ↑ GeoffG. Cumming GeoffG., RobertR. Maillardet RobertR., Confidence intervals and replication: Where will the next mean fall?, „Psychological Methods”, 11 (3), s. 217–227, DOI: 10.1037/1082-989x.11.3.217 [dostęp 2017-02-27] .
- ↑ MaxM. Henrion MaxM., BaruchB. Fischhoff BaruchB., Assessing uncertainty in physical constants, „American Journal of Physics”, 9, 1986, s. 791–798, DOI: 10.1119/1.14447, ISSN 0002-9505 [dostęp 2017-01-06] .
- Britannica: science/confidence-interval
- NE.se: konfidensintervall
- SNL: konfidensintervall
- Catalana: 0226011
- DSDE: konfidensinterval