ftp://cran.r-project.org/pub/pub/r/web/packages/survey/vignettes/domain.pdf
完整的数据集为TCH2012。但是，我只对满足两个标准的TCH2H2H2H2H2H2H2012感兴趣：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) of参数
由于我对处理加权数据的处理不太熟悉，所以我不知道如何麻烦拍摄。
事先感谢您的帮助！
此代码工作
library(survey)
data(api)
dclus1<-svydesign(id=~dnum, weights=~pw, data=apiclus1, fpc=~fpc)
x <- subset( dclus1 , sch.wide == 'Yes' )
svyby(~api00, ~stype, design=x, svytotal)
请通过添加最小可重现的示例来编辑您的问题，如何制作一个出色的R可重现示例？