Začněte nyníZačněte zdarma

Vyhodnocení rychlosti růstu

Pro výpočet rychlosti růstu můžeš provést lineární regresi logaritmu celkové bakteriální plochy v závislosti na čase. Vypočítej rychlost růstu a získej 95% interval spolehlivosti pomocí párového bootstrapu. Časové body v hodinách jsou uloženy v poli numpy s názvem t a bakteriální plocha v čtverečních mikrometrech je uložena v proměnné bac_area.

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

Případové studie ve statistickém myšlení

Zobrazit kurz

Pokyny k cvičení

  • Vypočítej logaritmus bakteriální plochy (bac_area) pomocí np.log() a výsledek ulož do proměnné log_bac_area.
  • Vypočítej sklon a průsečík semilogaritmické růstové křivky pomocí np.polyfit(). Sklon ulož do proměnné growth_rate a průsečík do log_a0.
  • Vygeneruj 10 000 párových bootstrap replikátů rychlosti růstu a logaritmu počáteční plochy pomocí dcst.draw_bs_pairs_linreg(). Výsledky ulož do growth_rate_bs_reps a log_a0_bs_reps.
  • Pomocí np.percentile() vypočítej 95% interval spolehlivosti rychlosti růstu (growth_rate_bs_reps).
  • Vypiš rychlost růstu a interval spolehlivosti na obrazovku. Tento krok je již připravený, stačí kliknout na Submit Answer a zobrazit výsledky!

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

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

# Compute logarithm of the bacterial area: log_bac_area
log_bac_area = ____

# Compute the slope and intercept: growth_rate, log_a0
____, ____ = ____

# Draw 10,000 pairs bootstrap replicates: growth_rate_bs_reps, log_a0_bs_reps
____, ____ = ____(
    ____, ____, size=____
)
    
# Compute confidence intervals: growth_rate_conf_int
growth_rate_conf_int = ____

# Print the result to the screen
print("""
Growth rate: {0:.4f} 1/hour
95% conf int: [{1:.4f}, {2:.4f}] 1/hour
""".format(growth_rate, *growth_rate_conf_int))
Upravit a spustit kód