Test de ipoteză bootstrap cu un singur eșantion
A mai fost studiată o broască juvenilă, Broasca C, și vrei să verifici dacă Broasca B și Broasca C au forțe de impact similare. Din păcate, nu ai la dispoziție forțele de impact ale Broaștei C, dar știi că au o medie de 0,55 N. Deoarece nu ai datele originale, nu poți efectua un test de permutare și nu poți evalua ipoteza că forțele celor două broaște provin din aceeași distribuție. Vei testa, prin urmare, o ipoteză mai puțin restrictivă: forța medie de lovire a Broaștei B este egală cu cea a Broaștei C.
Pentru a configura testul de ipoteză bootstrap, vei folosi media ca statistică de test. Reține că scopul tău este să calculezi probabilitatea de a obține o forță de impact medie mai mică sau egală cu cea observată la Broasca B, în condițiile în care ipoteza că media reală a forțelor Broaștei B este egală cu cea a Broaștei C este adevărată. Mai întâi, translezi toate datele Broaștei B astfel încât media să fie 0,55 N. Aceasta presupune adăugarea mediei forței Broaștei C și scăderea mediei forței Broaștei B din fiecare măsurătoare a Broaștei B. Astfel, alte proprietăți ale distribuției Broaștei B, cum ar fi varianța, rămân nemodificate.
Acest exercițiu face parte din cursul
Gândire statistică în Python (Partea 2)
Instrucțiuni pentru exercițiu
- Translează forțele de impact ale Broaștei B astfel încât media lor să fie 0,55 N.
- Folosește funcția
draw_bs_reps()pentru a genera 10.000 de replicate bootstrap ale mediei forțelor translatate. - Calculează valoarea p găsind fracția de replicate bootstrap care sunt mai mici decât media observată a forțelor de impact ale Broaștei B. Reține că variabila de interes aici este
force_b. - Afișează valoarea p.
Exercițiu interactiv practic
Încearcă acest exercițiu completând acest cod de exemplu.
# Make an array of translated impact forces: translated_force_b
translated_force_b = ____
# Take bootstrap replicates of Frog B's translated impact forces: bs_replicates
bs_replicates = draw_bs_reps(____, ____, 10000)
# Compute fraction of replicates that are less than the observed Frog B force: p
p = np.sum(____ <= np.mean(____)) / 10000
# Print the p-value
print('p = ', ____)