我写了一个包装ftable
因为我需要计算许多变量的频率和百分比的平面表:
mytable <- function(...) {
tab <- ftable(...,
exclude = NULL)
prop <- prop.table(x = tab,
margin = 2) * 100
bind <- cbind(as.matrix(x = tab),
as.matrix(x = prop))
margin <- addmargins(A = bind,
margin = 1)
round(x = margin,
digits = 1)
}
mytable(formula = wool + tension ~ breaks,
data = warpbreaks)
A_L A_M A_H B_L B_M B_H A_L A_M A_H B_L B_M B_H
10 0 0 1 0 0 0 0.0 0.0 11.1 0.0 0.0 0.0
12 0 1 0 0 0 0 0.0 11.1 0.0 0.0 0.0 0.0
13 0 0 0 0 0 1 0.0 0.0 0.0 0.0 0.0 11.1
14 0 0 0 1 0 0 0.0 0.0 0.0 11.1 0.0 0.0
15 0 0 1 0 0 2 0.0 0.0 11.1 0.0 0.0 22.2
...
Sum 9 9 9 9 9 9 100.0 100.0 100.0 100.0 100.0 100.0
但是,我不能使用子集参数ftable
与我的功能,也不与最小mytable <- function(...) ftable(...)
:
mytable(formula = wool + tension ~ breaks,
data = warpbreaks,
subset = breaks < 20)
Error in eval(substitute(subset), data, env) :
..3 used in an incorrect context, no ... to look in
我知道我可以在数据参数中使用子集data = warpbreaks[warpbreaks$breaks < 20, ]
作为一种解决方法,但我希望提高我对 R 的了解。”高级R https://adv-r.had.co.nz/Computing-on-the-language.htm”帮助我理解错误是由于非标准评估造成的,但我没有设法纠正我的代码。
所以我的问题是:
- 我怎样才能告诉R去寻找
breaks
in warpbreaks
?
- 更一般地说,是否有更明显的基本 R 方法来计算单个变量和多个变量的垂直布局中的频率和百分比的平面表? (我可以使用以下命令获得单个变量的垂直布局
mytable(x = warpbreaks$tension, row.vars = 1)
.)