Panel data model estimation with dummy variables In "R - project" I am trying to estimate the panel data lm model with plm function. When I include 3 dummy variables into the regression it doesn't appear in the summary of the model, but when I estimate a simple lm model it appears.
Why is it so? What should I do to estimate the statistics for those dummy variables?
 A: A possible reason might be that your dummies do not vary over time. In this case, the fixed effects estimator and first differencing will remove such variables from the model. The reason is that these estimators cannot identify variables that do not vary over time because those variables will be eliminated together with the unobserved fixed effects.
As an alternative you can try to use random effects and use a Hausman test to see whether random effects gives significantly different results from fixed effects. If not, you can use random effects. One problem with this approach is that the models are only comparable if they include the same variables. So for the Hausman test you should first exclude the dummies from the random effects model, too.
A: Another possibility is that the dummy variables are (very close to) co-linear. In this case, plm will automatically exclude the co-linear variables from its output. One way to check this is to run the model, and check the aliased object:
model.FE<- plm(income~area+weight+dproduct+dummy1+dummy2+dummy3,data=pa‌​nel, model = "within")
model.FE$aliased

If model.FE$aliased reports TRUE for any of the dummies, then they are linearly dependent (or very close being so). In that case, go back and check that dummy1,dummy2, and dummy3 are sufficiently different.
