Bonferroniho korekce
Pojďme implementovat vícenásobné hypotézové testy pomocí Bonferroniho korekce, o které jsme mluvili v prezentaci. K tomu využiješ importovanou funkci multipletests().
Pracuj s hladinou významnosti pro jednotlivý test .05 a sleduj, jak Bonferroniho korekce ovlivní náš vzorový seznam p-hodnot, který již byl vytvořen.
Toto cvičení je součástí kurzu
Procvičování statistických pohovorových otázek v Pythonu
Pokyny k cvičení
- Vypočítej seznam upravených p-hodnot pomocí Bonferroniho korekce za použití importované funkce
multipletests(). - Vypiš výsledky vícenásobných hypotézových testů vrácené na indexu 0 proměnné
p_adjusted. - Vypiš samotné p-hodnoty vrácené na indexu 1 proměnné
p_adjusted.
Interaktivní cvičení na vyzkoušení si v praxi
Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.
from statsmodels.sandbox.stats.multicomp import multipletests
pvals = [.01, .05, .10, .50, .99]
# Create a list of the adjusted p-values
p_adjusted = multipletests(____, alpha=____, method='bonferroni')
# Print the resulting conclusions
print(____)
# Print the adjusted p-values themselves
print(____)