ÎncepețiÎncepe gratuit

Estimarea ratei de creștere

Pentru a calcula rata de creștere, poți efectua o regresie liniară a logaritmului ariei totale bacteriene în raport cu timpul. Calculează rata de creștere și obține un interval de încredere de 95% folosind bootstrap pe perechi. Momentele de timp, exprimate în ore, sunt stocate în array-ul numpy t, iar aria bacteriană, exprimată în micrometri pătrați, este stocată în bac_area.

Acest exercițiu face parte din cursul

Studii de caz în gândirea statistică

Vezi cursul

Instrucțiuni pentru exercițiu

  • Calculează logaritmul ariei bacteriene (bac_area) folosind np.log() și stochează rezultatul în variabila log_bac_area.
  • Calculează panta și interceptul curbei de creștere semilogaritmice folosind np.polyfit(). Stochează panta în variabila growth_rate și interceptul în log_a0.
  • Generează 10.000 de replici bootstrap pereche ale ratei de creștere și ale ariei inițiale logaritmice folosind dcst.draw_bs_pairs_linreg(). Stochează rezultatele în growth_rate_bs_reps și log_a0_bs_reps.
  • Folosește np.percentile() pentru a calcula intervalul de încredere de 95% al ratei de creștere (growth_rate_bs_reps).
  • Afișează rata de creștere și intervalul de încredere. Acest pas a fost deja realizat pentru tine, așa că apasă Trimite răspunsul pentru a vedea rezultatele!

Exercițiu interactiv practic

Încearcă acest exercițiu completând acest cod de exemplu.

# 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))
Editează și rulează codul