Pomocou simulácií riešte nasledujúce úlohy:
Predpokladajme, že úspešnosť študenta na písomke má rovnomerné rozdelenie na intervale \((0,1)\) (po vynásobení 100 dostaneme úspešnosť v percentách). Výsledky jednotlivých študentov sú nezávislé.
Písomku písali aj Dlhý, Široký a Bystrozraký. Dlhý sa dozvedel neoficiálnu informáciu, že jeho úspešnosť bola viac ako \(0.6\), presný výsledok, ani výsledky svojich kamarátov ešte nevie. Aká je pravdepodobnosť, že má viac bodov ako oni?
Na modelovanie hodnôt medzi 0 a 1 sa používa často beta rozdelenie (rovnomerné rozdelenie je jeho špeciálnym prípadom). Jeho hustota je nenulová na intervale \((0,1)\) a môže mať rôzny tvar v závislosti od parametrov \(a\), \(b\). Jeho hustota je úmerná \(x^{a-1} (1-x)^{b-1}\).
V R-ku:
shape1 (zodpovedá parametru \(a\) z definície) a shape2
(zodpovedá parametru \(b\))dbeta, distribučná funkcia pbeta, náhodné
čísla rbetaUkážky :
curve(dbeta(x, shape1 = 50, shape2 = 20), from = 0, to = 1, col = "red", ylab = "hustota", n = 1000)
curve(dbeta(x, shape1 = 5, shape2 = 2), from = 0, to = 1, col = "green", add = TRUE)
curve(dbeta(x, shape1 = 2, shape2 = 5), from = 0, to = 1, col = "blue", add = TRUE)
curve(dbeta(x, shape1 = 2, shape2 = 2), from = 0, to = 1, col = "black", add = TRUE)
Predpokladajme teraz, že úspešnosť študenta na písomke má beta rozdelenie. Výsledky jednotlivých študentov sú nezávislé.
Zadanie 1.
Písomku písali aj Dlhý, Široký a Bystrozraký. Dlhý sa znovu dozvedel neoficiálnu informáciu, že jeho úspešnosť bola viac ako \(0.6\), presný výsledok, ani výsledky svojich kamarátov ešte nevie. Aká je pravdepodobnosť, že má viac bodov ako oni, ak parametre beta rozdelenia úspešnosti sú
Na základe grafov hustôt vysvetlite rozdiely.
Zadanie 2.
Aká je pravdepodobnosť, že študent s parametrami \(a=3\), \(b=6\) napíše písomku lepšie ako študent s parametrami \(a=8\), \(b=4\)? Nakreslite hustoty a skomentujte výsledok.
Zadanie 3.
Uvažujme postupnosť nezávislých pokusov, v ktorých má pravdepodobnosť úspechu beta rozdelenie so zvolenými parametrami. Stredná hodnota beta rozdelenia je \(a/(a+b)\).
Zvoľte si parametre a vygenerujte počty úspechov v takýchto pokusoch.
Graficky porovnajte relatívne početnosti jednotlivých výsledkov s pravdepodobnosťami binomického rozdelenia, ktoré uvažuje rovnaký počet pokusov a pravdepodobnosť úspechu rovnú strednej hodnote beta rozdelenia.
Ktoré rozdelenie má väčšiu disperziu a aké je intuitívne vysvetlenie? Porovnajte simulačne disperziu počtu úspechov