我尝试遍历各列,如果满足条件,则打印比例表,如果不满足条件,则打印均值表。
数据
library(survey)
dat <- data.frame(id=c(1,2,3), weight=c(0,2,0.1), var1=c(2,3,4), var2=c(2,6,7))
design <- svydesign(id=~1, weights=~weight, data=dat)
cols <- c("var1", "var2")
type <- c("prop", "mean")
我尝试过:
for(i in seq_along(cols)){
ifelse(type[i]=="prop",
print(prop.table(svytable(bquote(~.(as.name(cols[i]))), design))),
print(svymean(bquote(~.(as.name(cols[i]))), design))) }
我收到此错误消息:
Error in array(x, c(length(x), 1L), if (!is.null(names(x))) list(names(x), :
'data' must be of a vector type, was 'language'
我不太确定为什么bquote会给您这么多问题,但是如果我使用公式,它可以正常工作:
for(i in seq_along(cols)){
FORMULA= as.formula(paste("~",cols[i]))
ifelse(type[i]=="prop",
print(prop.table(svytable(FORMULA, design))),
print(svymean(FORMULA, design))) }
var1
2 3 4
0.00000000 0.95238095 0.04761905
mean SE
var2 6.0476 0.0786