Бивариат Probit / logit R: как найти ВСЕ коэффициенты и маргинальные эффекты с помощью пакета "zeligverse"

Question

Я использую двумерную модель logit в R с пакетом zeligverse. Я хочу рассчитать влияние моих независимых переменных на P (Y1 = 1), P (Y2 = 1), P (Y1 = 1, Y2 =0), P (Y1 = 1, Y2 = 1), P (Y1 = 0, Y2 = 1), P (Y1 = 0, Y2 = 0), P (Y1 = 1 | Y2 = 0) и все остальныеусловные вероятности (Y1 и Y2 - мои зависимые переменные. Они оба равны 0 или 1).Мне также нужны все предельные эффекты, связанные с этими вероятностями для каждой независимой переменной.

Знаете ли вы, как найти их в этом пакете (или в другом пакете, если он работает лучше)?

Answer 1

Не уверен, что это то, что вы ищете (не стесняйтесь, отметьте меня, если нет).Zelig пакеты действительно являются правильным выбором для вашего конкретного вопроса.

library(Zelig)

## Let X_i be independent variable
## Assume you are working with a univariate target variable Y where Y \in {0, 1} 

set.seed(123)
m <- 100
df <- data.frame(
  Y = rbinom(m, 1, 0.5),
  X1 = rbinom(m, 1, 0.95),
  X2 = rbinom(m, 1, 0.95) 
)

## Fit model once:
fit <- zelig(
  Y ~ .,
  model = "logit", 
  data = df, 
  cite = FALSE
)
summary(fit)

## Let's focus on the binomial predictor 2
x.out1 <- setx(fit, X2=1)

## Run estimation based on a posterior distribution:
postFit <- Zelig::sim(fit, x=x.out1)
summary(postFit)
# plot(postFit)