0

我使用anesrakeandweights包在我的调查中创建了加权数据。但是，除了包中的功能之外，我不确定之后如何使用权重wpct。我如何计算描述性统计数据并将加权数据与其他功能/包集成？

包中可重现的数据anesrake：

data("anes04")
anes04$caseid <- 1:length(anes04$age)
anes04$agecats <- cut(anes04$age, c(0, 25,35,45,55,65,99))
levels(anes04$agecats) <- c("age1824", "age2534", "age3544",
"age4554", "age5564", "age6599")
marriedtarget <- c(.4, .6)
agetarg <- c(.10, .15, .17, .23, .22, .13)
names(agetarg) <- c("age1824", "age2534", "age3544",
"age4554", "age5564", "age6599")
targets <- list(marriedtarget, agetarg)
names(targets) <- c("married", "agecats")
outsave <- anesrake(targets, anes04, caseid=anes04$caseid,
verbose=TRUE)
caseweights <- data.frame(cases=outsave$caseid, weights=outsave$weightvec)

这将为我提供一个带有数据框权重的新向量。所以，我的问题是，我怎样才能知道分析数据？如何将这些权重与汇总统计数据结合起来？

4

1 回答 1

1

您可以将权重作为weights=参数提供给survey::svydesign. 理想情况下，您会在调查包中进行 raking，以便您可以考虑 raking 的方差减少，但是将 raked 权重分析为只是抽样权重是非常标准的（至少在公共使用数据中）。

或者，如果您最终得到的 raking 规范足够简单，可以在survey::rakeor中重现survey::calibrate，您可以在调查包中重做 raking。

使用调查包的原因是它允许进行非常广泛的其他分析（甚至更多svyVGAM）。

于 2020-11-03T06:20:16.763 回答