我的模型中想控制行业虚拟变量和年度虚拟变量（使用固定效应模型），年度虚拟变量可以在estimate里的period中选择fixed，那行业虚拟变量呢？如何生成？

共16个行业，我试着手动输入了15个行业虚拟变量的数据，分别命名为D2,D3,D4......然后应该怎么做回归呢？请给出具体操作步骤（小白一个，，，若有朋友尽快回复，万分感谢）