liao083681 发表于 2015-1-15 18:47 
你解决了没??
我解决了,谢谢关心。
我给出我的算法,比如说做一个3*1向量的自回归,就应该算出每一个OLS回归式子的残差,然后再求这三个残差的变异数共变数矩阵,最后把它呈现出来。
predict resid1, resid eq(#1)
predict resid2, resid eq(#2)
predict resid3, resid eq(#3) //上面三个式子分别求残差
cor resid1 resid2 resid3, c //这个式子是计算变异数,最后那个c是表示求选择求变异数
mat a = r(C) //刚才求得变异数共变数矩阵被保存在r(C)矩阵里面,将其赋值给矩阵a