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

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

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

因为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
}

扫码加我 拉你入群

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

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

分享0 收藏3 回帖

关键词：Stata 循环语句 tata MATA Mat 程序 如何

运行后的结果如下：

. 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就退出了循环语句。

需要编写一个Mata函数，然后在程序中调用。

McvLJY 发表于 2013-8-26 23:54
运行后的结果如下：

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

mata计算部分写在一个do文件
循环里面do即可

kerrydu 发表于 2014-1-4 19:17
mata计算部分写在一个do文件
循环里面do即可

谢谢，这个方法可行！

使用mata mosave先将mata保存成.mo file，然后直接调用