Составьте некоторые данные:
set.seed(42)
Diuretic <- sample(c("yes", "no"), 100, replace=TRUE)
Year <- sample(c(2014, 2015, 2016), 100, replace=TRUE)
Beneficiaries <- round(rnorm(100, 35, 5))
dta <- data.frame(Diuretic, Year, Beneficiaries)
Теперь используйте xtabs
:
(dta.tbl <- xtabs(Beneficiaries~Diuretic+Year, dta)
# Year
# Diuretic 2014 2015 2016
# no 741 888 295
# yes 448 649 429
Добавьте итоги:
addmargins(dta.tbl)
# Year
# Diuretic 2014 2015 2016 Sum
# no 741 888 295 1924
# yes 448 649 429 1526
# Sum 1189 1537 724 3450