ÎncepețiÎncepe gratuit

Reprezentarea grafică a curbei de creștere

Ai văzut în exercițiul anterior că intervalul de încredere pentru curba de creștere este foarte strâns. Vei explora acest lucru grafic, reprezentând mai multe drepte bootstrap alături de curba de creștere. Vei folosi funcția plt.semilogy() pentru a crea graficul cu axa y pe scară logaritmică. Aceasta înseamnă că va trebui să transformi curba de regresie liniară teoretică înainte de reprezentare, aplicând exponențierea.

Acest exercițiu face parte din cursul

Studii de caz în gândirea statistică

Vezi cursul

Instrucțiuni pentru exercițiu

  • Reprezintă grafic punctele de date folosind plt.semilogy(). Tablourile numpy t și bac_area sunt disponibile în spațiul tău de lucru.
  • Folosește np.array() pentru a genera valorile de timp necesare reprezentării dreptelor bootstrap. Numește acest tablou t_bs. Timpul trebuie să meargă de la 0 la 14 ore.
  • Scrie o buclă for pentru a reprezenta grafic dreptele de regresie corespunzătoare primelor 100 perechi de replici bootstrap. Tablourile numpy growth_rate_bs_reps și log_a0_bs_reps calculate în exercițiul anterior sunt disponibile în spațiul tău de lucru.
    • Calculează curba de creștere prin exponențierea dreptei de regresie liniară, folosind np.exp().
    • Reprezintă grafic dreapta teoretică folosind plt.semilogy() cu argumentele cuvânt cheie linewidth=0.5, alpha=0.05 și color='red'.
  • Etichetează axele și afișează graficul. Etichetele potrivite pentru axele x și y sunt 'time (hr)' și, respectiv, 'area (sq. µm)'.

Exercițiu interactiv practic

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

# Plot data points in a semilog-y plot with axis labeles
_ = ____(____, ____, marker='.', linestyle='none')

# Generate x-values for the bootstrap lines: t_bs
t_bs = ____([____, ____])

# Plot the first 100 bootstrap lines
for i in range(____):
    y = ____(____[i] * ____ + ____[i])
    _ = ____(____, ____, linewidth=____, alpha=____, color=____)
    
# Label axes and show plot
_ = plt.xlabel('____')
_ = plt.ylabel('____')
____
Editează și rulează codul