计算加权logit模型的边际效应

Question

我使用＆＃34;调查计算了加权logit模型＆＃34;封装

 library(survey)

 Mod_design <- svrepdesign(variables = Data,
                    weights = weights,
                    repweights = REP_WGT,
                    type = "JKn",
                    scale = 1,       
                    rscales = 1)

 Mod <- svyglm(Education ~ 
                 Sex + Age + Edu_Parents, 
                 family = quasibinomial, design = Mod_design)

现在我想计算这个模型的边际效应。没有权重，我通常会使用mfx包的logitmfx函数。遗憾的是，使用此软件包计算加权模型的边际效应是不可能的，到目前为止，我无法找到解决此问题的方法。

有没有办法计算加权模型的边际效应？

Answer 1

您可以使用以下代码从 R 包 svyglm 中通过 survey 估计的模型获得边际效应：

summary(margins(MODEL_NAME, variables = "VARIABLE_OF_INTEREST", design=YOUR_DESIGN))