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)

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

1 个答案: