Расчет предельных эффектов для взвешенной логит-модели

Я рассчитал взвешенную логит-модель с помощью пакета «Survey».

 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)

Теперь я хотел бы рассчитать предельные эффекты для этой модели. Без весов я обычно использовал бы функцию logitmfx пакета mfx. К сожалению, с помощью этого пакета невозможно рассчитать предельные эффекты для взвешенных моделей, и до сих пор я не мог найти способ решить эту проблему.

Есть ли способ рассчитать предельные эффекты для взвешенных моделей?


r logistic-regression glm
person Joachim Schork    schedule 04.03.2015    source источник
comment
Это может быть вопрос CrossValidated, потому что вам нужно знать, как их вычислить, учитывая, что в R., вероятно, нет готового решения.   -  person Thomas    schedule 04.03.2015

Ответы (1)


arrow_upward
1
arrow_downward

Вы можете получить предельные эффекты от моделей, оцененных через svyglm из пакета R survey, используя следующий код:

summary(margins(MODEL_NAME, variables = "VARIABLE_OF_INTEREST", design=YOUR_DESIGN))
person Caserio    schedule 08.03.2021
comment
Спасибо, что повторно активировали эту ветку, и за ответ! - person Joachim Schork; 08.03.2021