关于在stata的循环语句中引用mata环境

4907

收藏 2013-08-26

我想要解决这样一个问题：有一批矩阵，想要求得每个矩阵每行的最小值，输出成一批列向量。

因为stata里本来没有直接求矩阵行最小值的命令，所以引用mata环境；又因为有一批这样的矩阵，所以要用循环语句。

这里问题来了，要退出mata环境必须要用end，但是stata认为这个end表示我想要退出循环！因为循环也可以用end退出，这样的话我就没有办法在循环语句中引用mata环境。

在这里向各位高手求助应该如何解决这个问题？

我举了一个简单的例子，程序如下：

mat x1 = (1,8\3,4)
mat x2 = (6,2\3,9)
mat x3 = (4,6\9,11)

local p = 1
while `p' <= 3{
mat z = x`p'

mata
a=st_matrix("z")
a
y=rowmax(a)
y
st_matrix("m",y)
end

mat y`p' = m
matlist y`p'

local p = `p'+1
}

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

McvLJY

2013-8-26 23:54:30

运行后的结果如下：

. mat x1 = (1,8\3,4)

. mat x2 = (6,2\3,9)

. mat x3 = (4,6\9,11)

.
. local p = 1

. while `p' <= 3{
  2. mat z = x`p'
  3.
. mata
  4. a=st_matrix("z")
  5. a
  6. y=rowmax(a)
  7. y
  8. st_matrix("m",y)
  9. end
--Break--
r(1);

end of do-file

--Break--
r(1);

也就是说一旦读到end，stata就退出了循环语句。