Testovací statistiky a velikost efektu

Jak lze zkoumat lineární vztahy pomocí bootstrap resamplingu? Zpět na stezku! Každá túra je zobrazena jako jeden bod a vidíme, že mezi celkovou ujetou vzdáleností a uplynulým časem existuje lineární vztah. Pokud budeme vzdálenost považovat za „efekt" uplynulého času, můžeme prozkoumat hlubší spojitost mezi lineární regresí a statistickou inferencí.

V tomto cvičení rozdělíš data do dvou populací neboli „kategorií": rané časy a pozdní časy. Pak se podíváš na rozdíly v celkové ujeté vzdálenosti v rámci každé populace. Tento rozdíl poslouží jako „testovací statistika" a jeho rozdělení otestuje vliv oddělení vzdáleností podle časů.

Toto cvičení je součástí kurzu

Úvod do lineárního modelování v Pythonu

Zobrazit kurz

Pokyny k cvičení

Pomocí "logického indexování" v numpy, např. sample_distances[sample_times < 5], rozděl ukázkové hodnoty distances do populací raných a pozdních časů.
Pomocí np.random.choice() s replacement=True vytvoř resample pro každý ze dvou časových intervalů.
Vypočítej pole test_statistic jako resample_long - resample_short a pomocí np.mean() a np.std() zjisti a vypiš velikost efektu a nejistotu.
Vykresli rozdělení test_statistic pomocí předdefinované funkce fig = plot_test_statistic().

Interaktivní cvičení na vyzkoušení si v praxi

Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.

# Create two poulations, sample_distances for early and late sample_times.
# Then resample with replacement, taking 500 random draws from each population.
group_duration_short = sample_distances[____ < 5]
group_duration_long = sample_distances[____ > 5]
resample_short = np.random.choice(____, size=500, replace=____)
resample_long = np.random.choice(____, size=500, replace=____)

# Difference the resamples to compute a test statistic distribution, then compute its mean and stdev
test_statistic = resample_long - resample_short
effect_size = np.mean(____)
standard_error = np.std(____)

# Print and plot the results
print('Test Statistic: mean={:0.2f}, stdev={:0.2f}'.format(____, ____))
fig = plot_test_statistic(____)

Upravit a spustit kód