r - 双变量 Probit/logit R：如何使用“zeligverse”包找到所有系数和边际效应

Question

我正在使用 zeligverse 包在 R 中运行双变量 logit 模型。我想计算我的独立变量对 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)