你只需要在我那段code之前做额外的一步,create a new dataset with i-j pairs:
use olddata
rename id id1
sort ind qu id1
save newdata1
rename id1 id2
rename earning1 earning2
sort ind qu id2
joinby ind qu using newdata1
drop if id1==id2
这样你的data就是i-j成对的了,然后用我上面回帖的代码就可以了,根本不用loop。
一般而言,能不用loop就不要用loop,stata的loop比较耗时。
其实用SAS的proc sql更有效率,基本一步就出来了