使用survey()包对加权数据的子群体进行估算

时间:2016-12-16 04:03:07

标签: survey weighted

ftp://cran.r-project.org/pub/R/web/packages/survey/vignettes/domain.pdf

完整的数据集是tch2012。但是,我只对满足两个标准的tch2012子群感兴趣:age <= 5gender == "female".在该子群体中,我想比较那些疾病(disease == "1")并且没有疾病(disease == "0").

这是我写的代码:

library(survey)

tch2012.tsl.dsgn <- svydesign(id= ~HOSP_KID, strata= ~KID_STRATUM, weights = ~DISCWT, data = tch2012, nest = TRUE)

创建一个指向5岁及以下女性儿科年龄亚群的指针

tch2012_f_age5.tsl.dsgn <- subset(tch2012.tsl.dsgn, AGE <= 5 & gender == "female")

5岁及以下女性儿科患者有无疾病总数的加权数据

svyby(~count, ~disease, design=tch2012_f_age5.tsl.dsgn, svytotal)

但是,当我运行svyby()

时,我收到以下错误消息

sum(sapply(covmats, ncol))中的错误:无效'类型'(list)参数

由于我不太熟悉处理加权数据,我不知道如何解决问题。

提前感谢您的帮助!

1 个答案:

答案 0 :(得分:0)

此代码有效

var app = angular.module('app', []);

请添加一个可重复性最小的示例How to make a great R reproducible example?

来编辑您的问题