没问题！我来帮你搞定这个不用手动敲16个虚拟变量的麻烦事~ 你需要在回归中加入年份的二阶多项式，同时只纳入45-60年的虚拟变量，这里有两种高效的方法：

方法一：用因子变量直接指定范围（最简洁） #

不用提前生成任何虚拟变量，直接借助Stata的因子变量语法，精准提取你需要的年份区间：

首先生成年份的二阶项：

gen year2 = year^2

然后构造一个临时变量，仅在45-60年有值、其余年份为空，这样Stata只会为目标年份生成虚拟变量，同时保留1-70年的全部样本：

gen year_dummy = year if inrange(year,45,60)
regress y year year2 i.year_dummy if inrange(year,1,70)

这样一来，Stata会自动为year=45到60的每个年份生成对应虚拟变量，完全不用手动逐个创建。

方法二：用循环批量生成虚拟变量（适合需要单独调用变量的场景） #

如果你需要把这些虚拟变量保留在数据集中后续使用，可以用forvalues循环批量生成：

// 先生成年份二阶项
gen year2 = year^2

// 批量生成45到60年的虚拟变量，每个变量命名为year_xxx
forvalues y = 45(1)60 {
    gen year_`y' = (year == `y')
}

// 回归时直接用变量范围引用所有生成的虚拟变量
regress y year year2 year_45-year_60 if inrange(year,1,70)

小提醒 #

• 两种方法中，Stata都会自动省略一个虚拟变量作为基准组（默认是最小的年份45），这是回归中避免多重共线性的常规操作。如果你确实需要纳入所有45-60的虚拟变量（不省略），可以在回归命令末尾加上noconstant，但这种情况很少见，需要谨慎使用。
• 方法一不需要额外生成变量，代码更简洁；方法二生成了明确的虚拟变量，方便后续查看或单独调用。

内容的提问来源于stack exchange，提问作者xgolt01