Lặp lại và phân khối

Thiết kế thí nghiệm với R

Joanne Xiong

Data Scientist

Lặp lại

Cần lặp lại thí nghiệm để đánh giá đầy đủ biến thiên
Nếu chỉ thử hiệu quả thuốc trên một người, làm sao khái quát kết quả? (Không thể!)

library(dplyr)
mtcars %>%
    count(cyl)

    cyl  n
1   4     11
2   6      7
3   8     14

Phân khối

Giúp kiểm soát biến thiên bằng cách làm các nhóm điều trị giống nhau hơn
Trong nhóm, khác biệt sẽ nhỏ; giữa nhóm, khác biệt sẽ lớn
Ví dụ: phân khối nhóm điều trị theo giới tính

Biểu đồ hộp

# Boxplot of MPG by Car Cylinders
ggplot(mtcars, aes(x=as.factor(cyl), 
                   y=mpg)) + 
  geom_boxplot(fill="slateblue", 
               alpha=0.2) + 
  xlab("cyl")

văn bản thay thế

Hàm dùng để mô hình hóa

Mô hình tuyến tính
```
lm(formula, data, na.action,...)
```
Mô hình ANOVA một nhân tố
```
aov(formula, data = NULL, ...)
```
Mô hình ANOVA lồng nhau
```
anova(object,...)
```

Hãy luyện tập!

Thiết kế thí nghiệm với R