hk0330_inclass_1

把HSAUR3叫出來

library(HSAUR3)

## Loading required package: tools

dta <- HSAUR3::backpain

看一下資料型態

str(dta)

## '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 ...

head(dta)

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