ftp://cran.r-project.org/pub/R/web/packages/survey/vignettes/domain.pdf
完整的数据集是tch2012。但是,我只对满足两个标准的tch2012子群感兴趣:age <= 5
和gender == "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)
参数
由于我不太熟悉处理加权数据,我不知道如何解决问题。
提前感谢您的帮助!
答案 0 :(得分:0)
此代码有效
var app = angular.module('app', []);
请添加一个可重复性最小的示例How to make a great R reproducible example?
来编辑您的问题