r - R 调查包中如何使用 BRR 权重？

Question

如果您的数据集已经包含 BRR 权重，有人知道如何使用 Lumley 调查包中的 BRR 权重来估计方差吗？

我正在使用 PISA 数据，他们的数据集中已包含 80 个 BRR 重复项。我怎样才能让 as.svrepdesign 使用这些，而不是尝试创建自己的？我尝试了以下操作并收到了后续错误:

dstrat <- svydesign(id=~uniqueID,strata=~strataVar, weights=~studentWeight, 
                data=data, nest=TRUE)
dstrat <- as.svrepdesign(dstrat, type="BRR")

Error in brrweights(design$strata[, 1], design$cluster[, 1], ..., 
    fay.rho = fay.rho,  : Can't split with odd numbers of PSUs in a stratum

任何帮助将不胜感激，谢谢。

Answer 1

如果您已有包含重复权重的数据框，则无需使用 as.svrepdesign() :) 您可以直接从数据框创建重复加权设计。

假设您有一个名为 mainwgt 的主权重列和名为 repwgt1 到 repwgt80 的 80 个重复权重列的数据，您可以使用此 --

yoursurvey <-
    svrepdesign( 
    weights = ~mainwgt , 
    repweights = "repwgt[0-9]+" , 
    type = "BRR", 
    data = yourdata ,
    combined.weights = TRUE
)

——这样，您就不必识别确切的列号。然后你可以运行正常的调查命令，例如 --

svymean( ~variable , design = yoursurvey )

如果您想要另一个示例，这里是 some example code和 an explanatory blog post使用当前的人口调查。