svypredmeans() 的比例置信区间

Question

我试图获得比例的 95% CI，这些比例实际上是预测边际平均值，如使用 R 的survey package计算的那样。我包括这个可重复的例子，在内容方面没有意义，但希望能很好地说明我的目的：

library(survey)
data(api)
dstrat <- svydesign(id=~1,strata=~stype, weights=~pw, data=apistrat, fpc=~fpc)

# 1.- marginal means for groups according to a variable 
svyglm(I(sch.wide=="Yes") ~ awards+comp.imp, design=dstrat, family=quasipoisson()) %>% svypredmeans(., ~yr.rnd)

# 2.- the confidence intervals I'd like to compute by svyciprop(), in an obviously wrong way   
svyby(I(sch.wide=="Yes") ~ awards+comp.imp, ~yr.rnd, design=dstrat, svyciprop, vartype="ci", method="xlogit")

我无法弄清楚如何将正确的参数输入svyciprop() ，或者这是否可能。 函数svyciprop()采用单个公式，这似乎与计算预测边际均值的方式不兼容，也不与svypredmeans()的输出svypredmeans() 。

在此先感谢您的帮助！

编辑道歉，因为代码是正确的，但有一个错字。 但是，还有一个后续问题：第 1 步中预测边际均值的估计与第 2 步中的相同估计不匹配。有人可以解释为什么会出现这些差异吗？

Answer 1

您的代码不起作用的原因很简单：R 区分大小写。 将sch.wide=="yes"更改为sch.wide=="Yes" ，它应该可以工作。