相关疑难解决方法(0)

任何有关这方面的帮助将非常感激.我正在使用Lumley调查包,我正在尝试简化我的代码,但是遇到了一些麻烦.

包中的svymean函数在我的代码中调用如下,其中第一个参数是指示我想要哪些变量的公式,第二个参数是该数据集:

svymean(~hq_ehla, FraSvy, na.rm=TRUE)

我正在尝试创建一个函数来提取分类变量的均值(比例)和标准误差,所以我做了以下函数:

stats <- function(repstat, num) {
    estmean <- as.numeric(round(100 * repstat[num], digits=0))
    estse <- round(100 * sqrt(attributes(repstat)$var[num,num]), digits=1)
    return(list(mean=estmean, se=estse))
}

这是有效的,所以当我拿出我的第一个类别的平均值和se时,我使用:

stats(svymean(~hq_ehla, FraSvy, na.rm=TRUE), 1)$mean
stats(svymean(~hq_ehla, FraSvy, na.rm=TRUE), 1)$se

我希望能够做到的是将其简化为更短的东西,也许我只需要写:

stats(FraSvy, "hq_ehla", 1)$mean

或类似的东西.问题是我无法弄清楚如何使用变量名将公式传递给函数.