预测忽略 feglm 估计中的 newdata 参数

Question

估计我的模型后

train_data <- df_est_hdfe[idx,]
test_data <- df_est_hdfe[-idx,]

est <- c("NATURAL", "GDPsim", "cap_lab_sim", "cap_lab_sim_sq", "cont")
form <- formula(paste0("rta ~ ",paste(est, collapse = " + "), "| exp_ind + imp_ind"))

train_data <- as.data.frame(train_data)

probithdfe <- feglm(form, data = train_data[,c("rta", est, "exp_ind", "imp_ind")], family = binomial(link = "probit"))

我想用它来预测给定一些新数据的输出

pred_LM <- predict(probithdfe, newdata = test_data[,c(est, "exp_ind", "imp_ind")], type = "response")

但是，我无法将我的新数据输入到预测函数中，因为

dim(test_data)
[1] 5910   13

而

length(pred_LM)
[1] 11546

这是训练数据的大小。如何将我的估计模型应用于新数据？我认为这可能特定于 alpaca 包的 feglm 规范。

Answer 1

我找到了一个使用 fixest 包的解决方法，它产生几乎相同的估计值，但允许我将新数据传递到 predict 函数中。该代码与上面显示的代码相同，因为两个包（fixest 和 alpaca）使用相同的 feglm 命令进行估算。

预测忽略 feglm 估计中的 newdata 参数

Predict ignores newdata argument in feglm estimation

r

predict

glm