Я пытаюсь получить подходящие значения из линейной модели со многими факторами, которые я хотел бы оценить, используя функцию felm
из пакета R lfe
. Если я неправильно интерпретирую значение fitted.values
, возвращаемого функцией, похоже, что эти значения не соответствуют выходным данным, которые я получаю, когда создаю их вручную. Вот пример, адаптированный из документации пакета:
library(lfe)
set.seed(42)
nn = 10
n1 = 3
x <- rnorm(nn)
f1 <- sample(n1, length(x), replace=TRUE)
y <- 2.13*x + cos(f1) + rnorm(length(x), sd=0.5)
est <- felm(y ~ x | f1)
estb <- lm(y~x+factor(f1)-1)
# we have exactly the same coefficients
getfe(est)['effect']/estb$coefficients[2:(n1+1)]
est$coefficients/estb$coefficients[1]
# but different fitted values -- in fact all having the same group offset
estb$fitted.values-est$fitted.values
Что это за смещения? Намерен ли felm
возвращать подходящее значение другого типа? Спасибо за поиск