调查R中的加权随机效应Logit模型

Question

我试图用包含使用调查数据的随机效应的模型来预测二元结果。我在下面介绍了抽样设计，请随时对我的调查加权方法发表评论。我的主要问题是如何在调查加权模型中包含随机效应。这里是代码达到这一点：

# Libraries 
library(survey) 
# Make dataframe object where d is the working dataframe 
dfobj <- svydesign((id = ~cluster+household, strata = ~interaction(region, urban) weights = ~chweight, strata = ~strata, data = d) 

# Run a logit model 
formula1 <- stunting ~ modern_toilet + diarrhoea + fever + insurance + sex + age + region_code 
model1 <- svyglm(formula=formula1,design=dfobj,family = quasibinomial) 

我想随机效应是在该地区。谢谢，

取样说明：

的MICS 2006年采用两阶段分层抽样设计。在抽样的第一阶段，选择了300 人口普查查点区（124个城市和176个农村EA）。这些是为GLSS 5选择的660个EAs（281个城市和379个农村）的子样本 。使用与其大小成比例的概率的系统抽样选择每个 区域中的聚类。

Answer 1

这不是在调查封装，我相当肯定它不是提供R.

有技术广义线性混合模型与调查数据拟合当集群的随机效应与抽样单位相同。 Stata（-gllamm-）和MLwin以及其他可能的实现。

混合模型的基于设计的推理是当前研究的一个领域。