Summarize the backpain{HSAUR3}
檢視資料型態
## 'data.frame': 434 obs. of 4 variables:
## $ ID : Factor w/ 217 levels "1","2","3","4",..: 1 1 2 2 3 3 4 4 5 5 ...
## $ status : Factor w/ 2 levels "case","control": 1 2 1 2 1 2 1 2 1 2 ...
## $ driver : Factor w/ 2 levels "no","yes": 2 2 2 2 2 2 1 1 2 2 ...
## $ suburban: Factor w/ 2 levels "no","yes": 2 1 2 2 1 2 1 1 1 2 ...
## ID status driver suburban
## 1 1 case yes yes
## 2 1 control yes no
## 3 2 case yes yes
## 4 2 control yes yes
## 5 3 case yes no
## 6 3 control yes yes
計算個數
dta.1 <- data.frame(table(dta$driver, dta$suburban, dta$status))
colnames(dta.1) <- c("driver", "suburban", "status", "frequency")
dta.1## driver suburban status frequency
## 1 no no case 26
## 2 yes no case 64
## 3 no yes case 6
## 4 yes yes case 121
## 5 no no control 47
## 6 yes no control 63
## 7 no yes control 7
## 8 yes yes control 100
繪製表格
driver <- c("no", "no", "yes", "yes")
suburban <- c("no", "yes", "no", "yes")
case <- c(26, 6, 64, 121)
control <- c(47, 7, 63, 100)
total <- c(73, 13, 127, 221)
data.frame(driver, suburban, case, control, total)## driver suburban case control total
## 1 no no 26 47 73
## 2 no yes 6 7 13
## 3 yes no 64 63 127
## 4 yes yes 121 100 221