假设你现在是2000年的一份截面数据，以下面一个简单例子给出解答，按照你的实际变量按需修改。如果没有理解错误，可以尝试一下。

试试：
clear

local N = 31  //研究对象的个体数
local T = 108  //研究的期数
local NT = `N'*`T'
set obs `NT'
set seed 12579

egen id =seq(),from(1) to(`N') block(`T') //设定T,N的值
sort id
* 生成地区
egen region = ntos(id), from(1/31) to("北京" "天津" "河北" "山西" "内蒙古" "辽宁" "吉林" "黑龙江" "上海" "江苏" "浙江" "安徽" "福建" "江西" "山东" "河南" "湖北" "湖南" "广东" "广西" "海南" "重庆" "四川" "贵州" "云南" "西藏" "陕西" "甘肃" "青海" "宁夏" "新疆")
* ntos命令通过 help egenmore 安装

* 生成日期
egen year = seq(), from(2010) to(2018) block(12)        // 先分块重复 12 次，再循环加 100，直到所有观察值数量
gen month=year*100
sort id year
by id year: gen months = _n + month

不是用xtreg year province就可以了吗？