BIBD - chức năng thận của mèo
Để thật sự hiểu BIBD trông như thế nào, hãy tự xây dựng một bộ dữ liệu chứa một BIBD từ đầu.
Giả sử bạn muốn kiểm tra sự khác biệt giữa bốn loại thức ăn ướt trong khẩu phần của mèo đối với chức năng thận. Tuy nhiên, thức ăn cho mèo khá đắt, nên mỗi khối ta chỉ thử 3 loại để tiết kiệm chi phí. Yếu tố chặn (blocking factor) là color của mèo, vì chúng ta không quan tâm yếu tố này trong thí nghiệm. Kết quả đầu ra sẽ là nồng độ creatinine trong máu được đo lường, một chỉ dấu về chức năng và rối loạn chức năng thận ở cả mèo và người.
Bài tập này là một phần của khóa học
Thiết kế thí nghiệm với R
Hướng dẫn bài tập
- Hàm tùy chỉnh
lambda()đã được nạp sẵn. Hãy tínhlambdavớit= 4,k= 3, vàr= 3 để đảm bảo BIBD là khả thi. - Chạy mã để lắp ráp bộ dữ liệu. Bạn có thể xem thứ tự các nghiệm thức (loại thức ăn) được dùng trong từng khối.
- Tạo
cat_modelvớiaov()theo mô tả thí nghiệm ở trên và xem kết quả bằngsummary(). Liệu loại thức ăn ướt có tạo ra khác biệt về mức creatinine không?
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.
# Calculate lambda
___
# Build the data.frame
creatinine <- c(1.98, 1.97, 2.35, 2.09, 1.87, 1.95, 2.08, 2.01, 1.84, 2.06, 1.97, 2.22)
food <- as.factor(c("A", "C", "D", "A", "B", "C", "B", "C", "D", "A", "B", "D"))
color <- as.factor(rep(c("Black", "White", "Orange", "Spotted"), each = 3))
cat_experiment <- as.data.frame(cbind(creatinine, food, color))
# Create cat_model and examine with summary()
___ <- ___ (___ ~ ___ + color, data = ___)
___