Прогноз вне выборки по модели с фиксированными эффектами

Рассмотрим модель:

library(plm)
data("Produc", package = "plm")
model <-  plm(pcap ~ hwy + water, data = Produc, model = 'within')

Чтобы рассчитать подогнанное значение модели, нам просто нужно использовать:

predict(model)

Однако при попытке сделать это из примера:

predict(model, newdata = data.frame('hwy' = 1, 'water' = 1))

Получит ошибку:

Error in crossprod(beta, t(X)) : non-conformable arguments

Что для меня довольно странно, потому что этот код будет работать для любой модели, кроме 'within'. Я ищу, что есть функция fixef, которая делает прогнозы по модели с фиксированным эффектом, но, к сожалению, только в образце. Итак: есть ли какое-либо решение, как мы можем предсказать выборку по модели с фиксированным эффектом?


r predict plm
person John    schedule 13.01.2021    source источник
comment
stackoverflow .com/questions/7123060/ может помочь   -  person user63230    schedule 13.01.2021

Ответы (1)


arrow_upward
0
arrow_downward

Просто удалите перехват для модели:

model <-  plm(pcap ~ 0 + hwy + water, data = Produc, model = 'within')
predict(model, newdata = data.frame('hwy' = 1, 'water' = 1))
3.980911
person John    schedule 14.01.2021