Estoy trabajando con un modelo de regresión de efectos fijos usando plm.

El modelo se ve así:

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" es una fórmula que definí anteriormente. Tengo muchas variables independientes, por lo que esto lo hizo más eficiente.

Lo que quiero hacer es obtener mis valores ajustados (mis yhats) y unirlos a mi conjunto de datos base; data.reg2

Pude obtener los valores ajustados usando este código:

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

Sin embargo, esto solo me da un vector de una columna de valores ajustados únicamente: no tengo forma de unirlo a mi conjunto de datos base.

Alternativamente, he intentado algo como esto:

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

Sin embargo, me sale este error con eso:

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

¿Hay alguna otra forma de obtener mis valores ajustados en mi conjunto de datos base? ¿O alguien puede explicar el error que estoy recibiendo y tal vez una forma de solucionarlo?

Debo señalar que no quiero calcular manualmente los yhats basados en mis betas. Tengo demasiadas variables independientes para esa opción y mi fórmula definida (fml) puede cambiar para que esa opción no sea eficiente.

¡¡Muchas gracias!!