我一直在研究如何计算两个分类变量中比例的CI。我已经看到了这个answer,它接近我正在寻找的东西,但我希望在两个变量的所有可能组合之间计算比例。我使用svymeans
和confint
让它工作,这是我正在寻找的输出,但是使用svyciprop
。此示例来自使用apiclus1数据框的调查包中的api数据。
library(survey)
data(api)
dclus1<-svydesign(id=~dnum, fpc=~fpc, data=apiclus1)
cbind(svymean(~interaction(sch.wide, stype), design = dclus1),
confint(svymean(~interaction(sch.wide, stype), design = dclus1)))
2.5 % 97.5 %
interaction(sch.wide, stype)No.E 0.06557377 0.031454436 0.09969311
interaction(sch.wide, stype)Yes.E 0.72131148 0.635732262 0.80689069
interaction(sch.wide, stype)No.H 0.01639344 -0.002458860 0.03524575
interaction(sch.wide, stype)Yes.H 0.06010929 0.018347771 0.10187081
interaction(sch.wide, stype)No.M 0.04371585 0.005196326 0.08223537
interaction(sch.wide, stype)Yes.M 0.09289617 0.050121072 0.13567128