请教各位大牛用stata得到动态回归系数的问题。

我想运行一个时间序列的线性回归，但希望系数是动态的。也就是：

y(t)=a + d1(t) + [ b+ d2(t)] *x1(t) + [c + d3(t)] *x2(t) + epsilon(t)

其中，系数中的随机变化部分d1(t), d2(t), d3(t)假定都服从AR（1）模型。也就是：
d1(t+1)=k* d1(t)+ e(t)
d2(t+1)=k* d2(t)+ e(t)
d3(t+1)=k* d3(t)+ e(t)

在stata里有dfactor命令，但不知道该怎样写出完整的指令？或者有其他方法可以解决这个问题吗？非常感谢大家的帮助！！