R中的概率函数和差函数的标准偏差

Question

我在R中很新。

我目前正在分析三个或三个以上孩子的单身女性的劳动力参与（LFP），以应对三个或三个以上孩子的单身女性的税收抵免时间表的扩大。

（治疗组：有三个或三个以上孩子的单身女性）

（对照组：有两个孩子的单身女性）

由于LFP是二进制变量，因此我使用了一个probit函数，并且为了隔离税收抵免扩展的影响，我使用了差异差异方法。

为了用代码说明这一点，

probit_model <- glm(lfp ~ three_kids + post_expansion + I(three_kids * post_expansion), family=binomial(link = "probit"), data=Mcps_s)
summary(probit_2)

其中three_kids是一个二进制变量，如果我有三个或三个以上的孩子，则该变量等于1，否则为零； post_expansion是一个二进制变量，如果是税收抵免扩展后的时间段，则等于1。

现在我感兴趣的是，我想获得差分差异估计及其标准偏差。

所以我使用predict()函数和diff()函数：

treat_prediction <- predict(probit_model,
                            newdata = data.frame("three_kids" = 1,
                                                 "post_expansion" = c(0, 1),
                                                 type = "response"))

diff(treat_prediction)

它的值为0.0216。

但是如何找到它的标准偏差？

如果您发现我犯的任何错误，请告诉我。

感谢您阅读我的文章。