我想为复杂的调查设计提供数据。我有重复重量(约50)和采样重量。虽然我找到了一个如何在R中进行插补的例子(见下文),但是它有一个取样权重的地方(weight=~WTPFQX6
,见下面有“设计”的行“)但是没有重复权重的地方。你能帮我解决一下如何用重复权重来估算数据吗?
library(mitools)
library(RSQLite)
impdata <- imputationList(c("set1","set2","set3","set4","set5"),
dbtype="SQLite", dbname="~/nhanes/imp.db")
designs <- svydesign(id=~SDPPSU6, strat=~SDPSTRA6,
weight=~WTPFQX6, data=impdata, nest=TRUE)
答案 0 :(得分:2)
您需要将svrepdesign
而不是svydesign
用于支持复制的设计,并且您不再需要指定id=
或strat=
参数,因为该信息应该是包含在复制权重中。在asdfree.com上最接近乘以nhanes的两个例子应该是
和
但svrepdesign
有一堆其他参数,您需要通过匹配某些cdc发布的统计信息来确认您正确使用这些参数。祝你好运!