在 R 上的 Gujarati 和 Porter 中從 logit model 示例中復制邊際效應

Question

我正在復制古吉拉特語和波特（西班牙語版）的計量經濟學書籍中的 logit model 示例。 我對 model 估計沒有任何問題，但我無法復制邊際效應。 在書中，回歸結果如下：

在以下幾行中，我展示了我的回歸結果：

> dat <- data.frame(
+   debito = c(0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 1, 
+              1, 0, 0, 0, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 
+              1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 1, 1, 1, 0, 1),
+   saldo = c(1756, 748, 1501, 1831, 1622, 1886,  740, 1593, 1169, 2125, 
+             1554, 1474, 1913, 1218, 1006, 2215,  137,  167, 343, 2557, 
+             2276, 1494, 2144, 1995, 1053, 1526, 1120, 1838, 1746, 1616,
+             1958, 634, 580, 1320, 1675, 789, 1735, 1784, 1326, 2051, 1044,
+             1885, 1790, 765, 1645, 32, 1266, 890, 2204, 2409, 1338, 2076,
+             1708, 2138, 2375, 1455, 1487, 1125, 1989, 2156),
+   cajero = c(13, 9, 10, 10, 14, 17, 6, 10, 6, 18, 12, 12, 6, 10, 12, 20, 
+              7, 5, 7, 20, 15, 11, 17, 10, 8, 8, 8, 7, 11, 10, 6, 2, 4, 4,
+              6, 8, 12, 11, 16, 14, 7, 10, 11, 4, 6, 2, 11, 7, 14, 16, 14,
+              12, 13, 18, 12, 9, 8, 6, 12, 14),
+   interes = c(1, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+               0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0,
+               0, 1, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0)
+ )
> mod <- glm(debito ~ saldo + cajero + interes,data = dat,
+            family = "binomial")
> summary(mod)

Call:
glm(formula = debito ~ saldo + cajero + interes, family = "binomial", 
    data = dat)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-1.6619  -0.9712  -0.6629   1.1440   1.7076  

Coefficients:
              Estimate Std. Error z value Pr(>|z|)  
(Intercept) -0.5749004  0.7857866  -0.732   0.4644  
saldo        0.0012481  0.0006973   1.790   0.0735 .
cajero      -0.1202251  0.0939842  -1.279   0.2008  
interes     -1.3520856  0.6809872  -1.985   0.0471 *
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 82.108  on 59  degrees of freedom
Residual deviance: 75.500  on 56  degrees of freedom
AIC: 83.5

Number of Fisher Scoring iterations: 4

如您所見，我得到了相同的結果。 我的問題是當我嘗試復制邊際效應時，書中的邊際效應如下：

我嘗試使用以下方法復制結果：

> coef(mod) * mean(dlogis(predict(mod, type = "link")))
  (Intercept)         saldo        cajero       interes 
-0.1262098468  0.0002740024 -0.0263934252 -0.2968279711 

如你所見，我失敗了。 我知道書中的結果是從 Stata 獲得的，所以我不知道計算的細節或者是否有可能在 R 中復制計算。 有辦法嗎？

Answer 1

使用 package“邊際”來獲得邊際效應。

install.packages("margins")
library("margins")

mod <- glm(debito ~ saldo + cajero + interes, data = dat, family = "binomial")

(m <- margins(mod))

以下是輸出：

Average marginal effects
glm(formula = debito ~ saldo + cajero + interes, family = "binomial",     data = dat)

    saldo   cajero interes
 0.000274 -0.02639 -0.2968

你的列名正確嗎？