请问工具变量乘以时间的stata命令怎么写？以下代码报错
顺便问一下有没有相关使用工具变量*时间的文献推荐？

//生成时间虚拟变量
tab year,gen(dyear)
//河流密度*时间虚拟变量
gen rivery =riverm*(dyear)