Predict 忽略 feglm 估计中的 newdata 参数答案

【问题标题】：Predict ignores newdata argument in feglm estimationPredict 忽略 feglm 估计中的 newdata 参数
【发布时间】：2021-02-07 06:22:13
【问题描述】：

估计我的模型后

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 规范。

【问题讨论】：

标签： r glm predict

【解决方案1】：

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

【讨论】：