Merge plm fitted values to dataset
问题 I'm working with a fixed effects regression model using plm. The model looks like this: 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" is a formula I defined previously. I have many independent variables, so this made it more efficient. What I want to do is get my fitted values (my yhats) and join them to my base dataset; data.reg2 I was able to get