I want to define a function lag_and_dummies which will add 2 lags for my data as well as 3-level dummy variables (having in mind how to include them properly in panel data).