真心感谢stata版中sungmoo斑竹，蓝色斑竹，连老师等老师，为我们这些素不相识的人细心解答，让我们少走了不少弯路！！！！

最后再请教斑竹一事，
因为stata处理矩阵的能力太弱，但mata中的循环语句还不太熟悉：
所以，想在是在stata的循环语句或者if语句中调用mata，但不知能否实现，该如何操作？
比如（以下例子乱举的），
forvalue i=1/100{
mata
1+1
end
}
有end后，stata就会终止循环

mata是专门处理矩阵的，stata处理的主要对象并不是矩阵。

可以使用mata的循环功能，再利用mata与stata间的数据传递功能，把mata的信息转化成stata的信息。

另外，这里涉及的一个问题是，应该说，stata的“使用思维”主要不是“矩阵式”的（这与其他软件不太一样），应该尽可能利用stata的各种对象（数据存储形式）与命令来实现既定目标。

Stata有现成的命令，例子1如下：

capt prog drop seq3 //简短的程序seq3
program seq3
drop _all
rndbin 10 0.5 1
gen z=0
forvalues i=3/10{
replace z=1 if xb==xb[_n-1] & xb==xb[_n-2] in `i'
}
sum z
end

//模拟抽样，将均值存储到矩阵a中
matrix a=J(100,100,-1.000)
set more off
forvalues i=1/100{
        forvalues j=1/100{
                quietly simulate max=r(max), reps(1000) nodots:seq3
                quietly sum
                //di "r(mean): "`r(mean)'
                matrix a[`i',`j']=(`r(mean)')
                di %5.4g as text "matrix a: "a[`i',`j']
        }
}


例子2
//逐行逐列的显示矩阵元素
matrix a=(1,0,1\2,1,0\-3,2,-5)
forvalues i=1/`=rowsof(a)'{
        forvalues j=1/`=colsof(a)'{
                local d=el(a,`i',`j')
                di "`d'"
        }
}