R:如何使用调查包计算比例？

Question

这只是一个非常简单的问题，但我无法从网络和书籍中找到要使用的正确功能。

这是我从这里的一篇文章中得到的一个例子。

df <- data.frame(sex = c('F', 'M', 'F', 'M', 'M', 'M', 'F', 'F'),
             married = c(1,1,1,1,0,0,1,1),
             pens = c(0, 1, 1, 1, 1, 1, 0, 0),
             weight = c(1.12, 0.55, 1.1, 0.6, 0.23, 0.23, 0.66, 0.67))

d.s <- svydesign(ids=~1, data=df, weights=~weight)

我想计算已婚等百分比变量并计算标准误差？

我也想做已婚和钢笔的交叉表并获得结果比例的标准误差？

我怎么做？

我试过 svymean 但它会将数值视为整数而不是因子。

Answer 1

使用 svytable

summary(d.s)
svytable(~married+pens, d.s)
svytable(married~pens, d.s)
svytable(married~., d.s) #with all variable