Bắt đầu ngayBắt đầu miễn phí

Mô hình hóa tăng trưởng đậu nành với GAM

Trong bài tập này, bạn sẽ mô hình hóa trọng lượng lá trung bình trên một cây đậu nành theo thời gian (sau khi gieo trồng). Bạn sẽ thấy cây đậu nành không lớn lên với tốc độ ổn định mà có một “đợt bứt phá” rồi dần chậm lại. Vì vậy, mô hình tuyến tính không mô tả tốt trọng lượng lá.

Hãy nhớ rằng bạn có thể chỉ định biến muốn mô hình hóa phi tuyến trong công thức bằng hàm s() (docs):

y ~ s(x)

Cũng lưu ý rằng gam() (docs) từ gói mgcv có giao diện gọi hàm như sau

gam(formula, family, data)

Với hồi quy chuẩn, dùng family = gaussian (mặc định).

Dữ liệu huấn luyện đậu nành, soybean_train, đã được nạp sẵn. Nó có hai cột: biến mục tiêu weight và biến Time. Để so sánh, mô hình tuyến tính model.lin, được fit với công thức weight ~ Time, cũng đã được nạp sẵn.

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

Học có giám sát với R: Hồi quy

Xem khóa học

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.

# soybean_train is available
summary(soybean_train)

# Plot weight vs Time (Time on x axis)
ggplot(soybean_train, aes(x = ___, y = ___)) + 
  geom_point() 
Chỉnh sửa và Chạy Mã