Предельные эффекты после панели пробит в R

Question

Я использую R.

Я оценил коэффициенты панельной пробитной модели следующим образом:

Model_1 <- pglm(dependent variable ~ var_1 +
                       var_2 +
                       var_3+
                       var_4,
                     model=("random"), 
                     effect=("twoways"), 
                     index=c("indiv", "time"), 
                     family=binomial(link="probit"), data=data)

Что я хотел бы сделать, это вычислить предельные эффекты в среднем для каждой из этих переменных (как в stata после xtprobit). Но я не нашел никакого эквивалента в R.

Заранее спасибо, Марк