*这边简单利用Stata的mata来处理
*因为用到putmata与getmata所以限用Stata12版以上
inp ID1 X1
1 10
2 20
3 30
end
putmata ID1X1=(ID1 X1)
mata: ones=J(3,1,1)
mata: a=ID1X1#ones
clear
inp ID2 X2
1 10
2 20
3 30
end
putmata ID2X2=(ID2 X2)
mata: b=ones#ID2X2
clear
mata: c=(a,b)
getmata (ID1 X1 ID2 X2)=c