如何使用R评估dhs家庭暴力数据的权重科特迪瓦2011-2012

问题描述

我想使用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中设置数据集有任何见解并帮助我修复此错误消息,我将不胜感激!

解决方法

我们可以为NA创建一个逻辑向量,将数据子集并应用函数

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