使用的指令是:
program modelselect1
1. scalar aic=ln(e(rss)/e(N))+2*e(rank)/e(N)
2. scalar list aic
3. end
forvalues q=1/24{
2. regress d.lninr l.lninr L(1/`q').d.lninr trend
3. display "q=`q'"
4. modelselect1
5. }
结果总是在q=3就停止了,显示no observations(r2000).
然后我利用另外两个变量(lncny,lnwti,lninr)做出来的结果也是这样,滞后为4就显示no observations了。很郁闷
这里附上数据请大牛帮我看下,刚用stata不久,请多指教