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

我试图用一个模型来预测一个二元结果，该模型使用调查数据包含一个随机效应。我在下面介绍了抽样设计，所以请随意评论我的调查加权方法。我的主要问题是如何在调查加权模型中包含随机效应。以下是到目前为止的代码：

# 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)

#库
图书馆（调查）
#生成dataframe对象，其中d是工作数据帧
dfobj这在调查包中不可用，我相当肯定它在R中不可用

当随机效应的聚类与抽样单位相同时，有一些技术可用于对调查数据进行广义线性混合模型拟合。在Stata（-gllamm-）和MLwin中有实现，可能还有其他实现

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