小编pba*_*lis的帖子

有没有办法从标准误差中获得预测行为,lfe::felm如果使用投影方法扫除固定效果felm？这个问题与此处的问题非常相似,但该问题的答案都不能用于估计标准误差或置信度/预测间隔.我知道目前没有predict.felm,但我想知道是否有类似于上面链接的变通方法,可能也可用于估计预测间隔

library(DAAG)
library(lfe)

model1 <- lm(data = cps1, re74 ~ age + nodeg + marr)
predict(model1, newdata = data.frame(age=40, nodeg = 0, marr=1), se.fit = T, interval="prediction")$fit
# Result:        fit      lwr      upr
# 1 18436.18 2339.335 34533.03

model2 <- felm(data = cps1, re74 ~ age | nodeg + marr)
predict(model2, newdata = data.frame(age=40, nodeg = 0, marr=1), se.fit = T, interval="prediction")$fit
# Does not work

目标是估计yhat的预测间隔,我认为我需要计算完整的方差 - 协方差矩阵(包括固定效应).我一直无法弄清楚如何做到这一点,我想知道它是否在计算上是可行的.

我怀疑我在做一些非常愚蠢的事情，但我无法让该dst功能按lubridate预期工作。

library(lubridate)

x <- c("2016-01-01", "2016-06-01")

dst(x) # Returns c(FALSE, TRUE)
dst(as.Date(x)) # Returns c(FALSE, FALSE)

我在这两种情况下期望的结果都是c(FALSE, TRUE)。dst但是，只有当我传递字符向量而不是对象时，我才会得到预期的结果Date。我使用的是 OS X，我当前的时区是 PST ( America/Los_Angeles)。