Przedział ufności – Wikipedia, wolna encyklopedia

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.
Przedziały ufności (95%) oparte na stu symulacjach prób z tego samego rozkładu normalnego. Szerokość i położenie każdego przedziału to parametry oparte na próbie, w związku z czym cechują się one zmiennością. Można jednak oczekiwać, że w idealnej sytuacji na sto przypadków tylko około pięć nie będzie zawierać średniej oryginalnego rozkładu.

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ę o błąd standardowy. Szczególny przypadek przedziału ufności w badaniach ankietowych jest zwyczajowo określany marginesem błędu[1].

Definicja[edytuj | edytuj kod]

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[edytuj | edytuj kod]

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.
55 wybranych historycznych oszacowań prędkości światła w próżni wraz z przybliżonymi przedziałami ufności (86%; 1.48 × szacowany w badaniu błąd; na podstawie pracy Henriona i Fischoffa[2]). Gdyby eksperymenty nie były obarczone błędami systematycznymi, można byłoby oczekiwać aby tylko około 8 przedziałów nie zawierało poprawnej wartości. W tym przypadku aż 24 nie są trafne (zaznaczone na żółto).

Przedziały ufności konstruowane w oparciu o paradygmat częstościowy 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[edytuj | edytuj kod]

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[edytuj | edytuj kod]

Znane odchylenie standardowe[edytuj | edytuj kod]

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:

  • – liczebność próby losowej,
  • – średnia z próby losowej,
  • – odchylenie standardowe populacji,
  • to kwantyl rzędów rozkładu

Nieznane odchylenie standardowe[edytuj | edytuj kod]

Cecha ma w populacji rozkład normalny przy czym odchylenie standardowe jest nieznane. Przedział ufności dla parametru tego rozkładu ma postać:

gdzie:

  • – liczebność próby losowej,
  • – średnia z próby losowej,
  • – odchylenie standardowe z próby,
  • to kwantyl rzędu rozkład 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)[edytuj | edytuj kod]

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,
  • – statystyka ze zmienną losową o rozkładzie normalnym

Przedział ufności dla wariancji[edytuj | edytuj kod]

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 – statystyki spełniające odpowiednio nierówności:


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)[edytuj | edytuj kod]

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 (wskaźnik struktury)[edytuj | edytuj kod]

Poniższy wzór pozwala wyznaczyć przedział ufności dla odsetka 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[edytuj | edytuj kod]

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)[edytuj | edytuj kod]

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 [edytuj | edytuj kod]

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[edytuj | edytuj kod]

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.

Przypisy[edytuj | edytuj kod]

  1. Dennis Gilliland, Vince Melfi, A Note on Confidence Interval Estimation and Margin of Error, „Journal of Statistics Education”, 18 (1), 2010, DOI10.1080/10691898.2010.11889474, ISSN 1069-1898 [dostęp 2019-03-31] (ang.).
  2. Max Henrion, Baruch Fischhoff, Assessing uncertainty in physical constants, „American Journal of Physics”, 54 (9), 1986, s. 791–798, DOI10.1119/1.14447, ISSN 0002-9505 [dostęp 2019-03-30] (ang.).
  3. a b c Rink Hoekstra i inni, Robust misinterpretation of confidence intervals, „Psychonomic Bulletin & Review”, 5, 2014, s. 1157–1164, DOI10.3758/s13423-013-0572-3, ISSN 1531-5320, PMID24420726 [dostęp 2017-01-06].
  4. 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.).
  5. 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.).
  6. a b Sander Greenland i inni, Statistical tests, P values, confidence intervals, and power: a guide to misinterpretations, „European Journal of Epidemiology”, 31 (4), 2016, s. 337–350, DOI10.1007/s10654-016-0149-3, ISSN 0393-2990, PMID27209009, PMCIDPMC4877414 [dostęp 2017-02-09] (ang.).
  7. Ulrich K. Müller, Andriy Norets, Credibility of Confidence Sets in Nonstandard Econometric Problems, „Econometrica”, 6, 2016, s. 2183–2213, DOI10.3982/ECTA14023, ISSN 1468-0262 [dostęp 2017-01-06] (ang.).
  8. Geoff Cumming, Jennifer Williams, Fiona Fidler, Replication and Researchers’ Understanding of Confidence Intervals and Standard Error Bars, „Understanding Statistics”, 3 (4), 2004, s. 299–311, DOI10.1207/s15328031us0304_5, ISSN 1534-844X [dostęp 2017-04-01].
  9. Geoff Cumming, Robert Maillardet, Confidence intervals and replication: Where will the next mean fall?, „Psychological Methods”, 11 (3), s. 217–227, DOI10.1037/1082-989x.11.3.217 [dostęp 2017-02-27].
  10. Max Henrion, Baruch Fischhoff, Assessing uncertainty in physical constants, „American Journal of Physics”, 9, 1986, s. 791–798, DOI10.1119/1.14447, ISSN 0002-9505 [dostęp 2017-01-06].