如何使用 R 为 dhs 家庭暴力数据加权数据（科特迪瓦 2011-2012）答案

【问题标题】：How to weight data using R for dhs domestic violence data (Cote d'Ivoire 2011-2012)如何使用 R 为 dhs 家庭暴力数据加权数据（科特迪瓦 2011-2012）
【发布时间】：2020-10-30 17:22:17
【问题描述】：

我想使用 R (2011-2012) 为科特迪瓦 DHS 设置调查设计。我专注于家庭暴力。

我使用了这个代码：

wt <- ipv_studyf$d005/1000000 # creating the variable weight

ipvdesign <- svydesign(ids = ipv_studyf$v021,      # clustering- psu
                        strata = ipv_studyf$v023,   # strata variable
                        weights = wt,               # the weighting variable
                        data = ipv_studyf,          # the dataset
                        variance="HT" )             # Horvitz-Thompson estimator

但是，我不断收到此错误消息：

**Error in na.weight(data.frame(weights)) : missing values in `weights'**

如果有人对如何在 R 中设置数据集有任何见解并帮助我修复此错误消息，我将不胜感激！

【问题讨论】：

试试i1 <- !is.na(wt)，然后做svydesign(ids = ipv_studyf$v021[i1], strata = ipv_studyf$v023[i1], weights = wt[i1], data = ipv_studyf, variance = "HT")

标签： r survey weighted

【解决方案1】：

我们可以为NA创建一个逻辑向量，对数据进行子集化并应用函数

 i1 <- !is.na(wt)
 svydesign(ids = ipv_studyf$v021[i1], 
           strata = ipv_studyf$v023[i1], 
           weights = wt[i1], 
           data = ipv_studyf, variance = "HT")

【讨论】：

非常感谢，@akrun！它运行良好。