From Łukasz Graczykowski
(Difference between revisions)
|
|
| (12 intermediate revisions not shown) |
| Line 1: |
Line 1: |
| | + | <!-- |
| | == Zadanie == | | == Zadanie == |
| | '''Weryfikacja hipotez statystycznych''' (5 pkt.) | | '''Weryfikacja hipotez statystycznych''' (5 pkt.) |
| | | | |
| - | * Przeprowadzono eksperyment naświetlania wodorowej komory pęcherzykowej wiązką fotonów w celu badania oddziaływań fotonów z protonami. Fotony powodują powstawanie par elektron-pozyton, które mogą być wykorzystane do monitorowania wiązki fotonów. Częstość występowania zdjęć z 0,1,2,... parami elektron-pozyton powinna podlegać rozkładowi Poissona. Należy wczytać dane z [http://hirg.if.pw.edu.pl/~mszymanski/kadd/dane.txt pliku] (w pierwszej kolumnie znajduje się liczba par elektronowych na zdjęciu <code>k</code>, a w drugiej liczba zdjęć zawierających <code>k</code> par elektronowych). Widzimy, że rozkład te przypomina rozkład Poissona - próbujemy (1 pkt.) | + | * Przeprowadzono eksperyment naświetlania wodorowej komory pęcherzykowej wiązką fotonów w celu badania oddziaływań fotonów z protonami. Fotony powodują powstawanie par elektron-pozyton, które mogą być wykorzystane do monitorowania wiązki fotonów. Częstość występowania zdjęć z 0,1,2,... parami elektron-pozyton powinna podlegać rozkładowi Poissona. Należy wczytać dane z [http://www.if.pw.edu.pl/~lgraczyk/KADD2016/lab10/dane10.txt pliku] (w pierwszej kolumnie znajduje się liczba par elektronowych na zdjęciu <code>k</code>, a w drugiej liczba zdjęć zawierających <code>k</code> par elektronowych). Widzimy, że rozkład ten przypomina rozkład Poissona - próbujemy zatem obliczyć estymator największej wiarygodności dla parametry rozkładu Poissona (patrz [http://www.if.pw.edu.pl/~lgraczyk/KADD2016/Wyklad10-2016.pdf Wykład 10] slajd 13) (1 pkt.) |
| | | | |
| - | * Narysować na jednym wykresie punkty pomiarowe i dopasowanie. (1 pkt.) | + | * Narysować na jednym wykresie punkty pomiarowe i dopasowanie (metodą estymatora największej wiarygodności i funkcją Fit z ROOT'a użytą z parametrami "LR" - dopasowanie metodą największej wiarygodności). Funkcja TF1 do rysowania (i dopasowania ROOT'em) to TMath::PoissonI (1 pkt.) |
| | | | |
| | * Sprawdzić jakość dopasowania za pomocą testu χ2. W tym celu należy zaimplementować funkcję obliczającą statystykę testową χ2 zgodnie z wzorem [[File:wzor.png]] | | * Sprawdzić jakość dopasowania za pomocą testu χ2. W tym celu należy zaimplementować funkcję obliczającą statystykę testową χ2 zgodnie z wzorem [[File:wzor.png]] |
| Line 28: |
Line 29: |
| | == Wynik == | | == Wynik == |
| | | | |
| - | [[File:lab10.png]] | + | [[File:Lab10_2.png]] |
| | | | |
| | Output: | | Output: |
| Line 38: |
Line 39: |
| | 2 p1 2.33737e+00 8.17264e-02 1.40382e-04 -2.26405e-03 | | 2 p1 2.33737e+00 8.17264e-02 1.40382e-04 -2.26405e-03 |
| | ERR DEF= 0.5 | | ERR DEF= 0.5 |
| - | | + | |
| - | | + | |
| | Lambda najwiekszej wiarygodnosci: 2.33239 | | Lambda najwiekszej wiarygodnosci: 2.33239 |
| | + | Lambda (ROOT Fit): 2.33737 |
| | chi2 (wartosc statystyki testowej T): 10.5336 | | chi2 (wartosc statystyki testowej T): 10.5336 |
| | chi2/NDF: 1.7556 | | chi2/NDF: 1.7556 |
| Line 47: |
Line 48: |
| | Poziom istotnosci alpha: 0.01 | | Poziom istotnosci alpha: 0.01 |
| | Wynik testu: nie ma podstaw do odrzucenia hipotezy | | Wynik testu: nie ma podstaw do odrzucenia hipotezy |
| | + | --> |
Latest revision as of 15:20, 24 February 2017
