Estou trabalhando com um modelo de regressão de efeitos fixos usando plm.

O modelo fica assim:

FE.model <-plm(fml, data = data.reg2,
           index=c('Site.ID','date.hour'), # cross section ID and time series ID
           model='within', #coefficients are fixed
           effect='individual')
summary(FE.model)

"fml" é uma fórmula definida anteriormente. Eu tenho muitas variáveis independentes, então isso tornou mais eficiente.

O que eu quero fazer é obter meus valores ajustados (meus yhats) e juntá-los ao meu conjunto de dados base; data.reg2

Consegui obter os valores ajustados usando este código:

 Fe.model.fitted <- FE.model$model[[1]] - FE.model$residuals

No entanto, isso me fornece apenas um vetor de uma coluna de valores ajustados - não tenho como associá-lo ao meu conjunto de dados base.

Como alternativa, tentei algo assim:

 Fe.model.fitted <- cbind(data.reg2, resid=resid(FE.model), fitted=fitted(FE.model))

No entanto, eu recebo este erro com isso:

 Error in as.data.frame.default(x[[i]], optional = TRUE) : cannot coerce class ""pseries"" to a data.frame

Existem outras maneiras de obter meus valores ajustados no meu conjunto de dados base? Ou alguém pode explicar o erro que estou recebendo e talvez uma maneira de corrigi-lo?

Devo observar que não quero calcular manualmente os yhats com base nos meus betas. Eu tenho muitas variáveis independentes para essa opção e minha fórmula definida (fml) pode mudar para que a opção não seja eficiente.

Muito Obrigado!!