关于DID分析，设计了虚拟变量DID（0，1），政策存在时间t,进行stata命令时报错

数据是面板格式的，想探究政策实施对地区产业结构的影响，是非平衡面板，各个政策的实施时间也不同。
设计模型：air(i,t)=a0+a1did*t+bx(i,t)+Ct+Ut+E
air是产业结构，did是政策实施的虚拟变量，实施前是0，之后是1，t是政策实施后到2015年的存在年份，如2005年实施，t=2015-2005=10,x(i,t)是各个控制变量（gov pergdp phone stud indus invest medical），Ct是时间固定效应，Ut是个体固定效应，E是遗漏。didt=did*t
开始 命令是： quietly xtreg ais didt gov pergdp phone stud indus invest medical i.year,fe r
must specify panelvar; use xtset
改完变成
. quietly xtset ais didt gov pergdp phone stud indus invest medical i.year,fe r
报错是：factor variables and time-series operators not allowed

请问这个分析应该怎么做？
如果模型改成air(i,t)=a0+a1did+bx(i,t)+Ct+Ut+E，又怎么做？

本人刚刚学STATA,错误如果比较低级，希望大神指导