I want to create a function which will perform panel regression with 3-level dummies included.
Let\'s consider within model with time effects :
library