求教！如何将年度行业虚拟变量加入回归模型？

假如我有一个数据，是四个年度，21个行业的。我看别人的论文都是在模型最后加上年度和行业的控制变量。
这个要怎么实现啊？