egen ye_coun = group(year country)        这是年份国家固定效应
egen g = group(product_code6)               这是产品固定效应
su g                                                           q指产品数量，p指产品价格， Lnp_nagative=-lnp
loc n=r(max)
g m=.
g e_p=.
forv i=1/`n'{
reg lnq Lnp_nagative pop i.ye_coun  if g==`i'
mat b=e(b)
replace m=b[1,1] if g==`i'
predict r if e(sample),r
replace e_p=r if e(sample)
drop r
}

上面是我原来的代码
但是现在为解决内生性问题，想拿工具变量（两个工具变量：lne，lncost）来替代lnp进行回归。
上面的语句该怎么修改？[smile][smile]