2014-12-04 71 views
2

我想爲複雜勘測設計插補數據。我有重複的重量(約50)和取樣重量。雖然我發現了一個如何在R中進行插值的例子(見下文),但它有一個抽樣權重的地方(weight=~WTPFQX6,參見下面的「設計」一行),但沒有重複權重的地方。你能幫我弄清楚如何使用複製權重來推算數據嗎?複雜勘測設計的插補(nhanesIII)

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) 

回答

2

你需要使用svrepdesignsvydesign複製支持的設計,你不再需要指定id=strat=參數,因爲這些信息應包含在複製權。兩個最接近的例子來乘估算上asdfree.com NHANES應該

https://github.com/ajdamico/usgsd/blob/3ff408b6acc0971129d0435388ff13749d3f5129/Programme%20for%20the%20International%20Assessment%20of%20Adult%20Competencies/download%20import%20and%20design.R#L188

https://github.com/ajdamico/usgsd/blob/d300884bd63dd05c61e8a6fa76ed7293adae55c2/Consumer%20Expenditure%20Survey/2011%20fmly%20intrvw%20-%20analysis%20examples.R#L497

svrepdesign有一堆,你需要確認你的其它參數的使用正確匹配一些cdc發佈的統計信息。祝你好運!