Đánh giá tốc độ tăng trưởng

Để tính tốc độ tăng trưởng, bạn có thể hồi quy tuyến tính của logarit diện tích vi khuẩn tổng theo thời gian. Hãy tính tốc độ tăng trưởng và lấy khoảng tin cậy 95% bằng pairs bootstrap. Các mốc thời gian, tính theo giờ, được lưu trong mảng numpy t và diện tích vi khuẩn, tính theo micromet vuông, được lưu trong bac_area.

Bài tập này là một phần của khóa học

Các nghiên cứu tình huống về tư duy thống kê

Xem khóa học

Hướng dẫn bài tập

Tính logarit của diện tích vi khuẩn (bac_area) bằng np.log() và lưu kết quả vào biến log_bac_area.
Tính hệ số góc và hệ số chặn của đường cong tăng trưởng semilog bằng np.polyfit(). Lưu hệ số góc vào biến growth_rate và hệ số chặn vào log_a0.
Vẽ 10.000 bản sao bootstrap theo cặp cho tốc độ tăng trưởng và log diện tích ban đầu bằng dcst.draw_bs_pairs_linreg(). Lưu kết quả vào growth_rate_bs_reps và log_a0_bs_reps.
Dùng np.percentile() để tính khoảng tin cậy 95% của tốc độ tăng trưởng (growth_rate_bs_reps).
In tốc độ tăng trưởng và khoảng tin cậy ra màn hình. Phần này đã được làm sẵn, nên hãy bấm "Gửi câu trả lời" để xem kết quả!

Bài tập tương tác thực hành trực tiếp

Hãy thử làm bài tập này bằng cách hoàn thành đoạn mã mẫu này.

# 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))

Chỉnh sửa và Chạy Mã