如何将mata 的数据传回stata

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

小弟求教，如果在mata中有一个n*k的矩阵，我想把这个矩阵的数据传回到stata中，我知道使用st_addvar() 和st_store() 但是，如果我想一次生成k个变量，即矩阵的每一列都成为一个在stata中的变量而且列数如果很多， 这样应该如何来做呢？ 谢谢。。

扫码加我 拉你入群

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

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

分享0 收藏4 回帖

关键词：Stata tata MATA Mat Store 如何 而且

http://www.pinggu.org/bbs/thread-547419-1-1.html

2# sungmoo
谢谢楼上，这个贴我之前看过了，不过我的问题是我需要同时从矩阵中产生多个变量储存在stata中，有点类似于一个循环语句，我现在这样写，假设我在mata中已经存储了一个n*k的矩阵x，我想生成k个变量，即x的每一列
i=1
n=k
while (i<=k) {
y=.
st_subview(y,x,.,i)
st_addobs(rows(y))
st_addvar("float", "t") /*关键就是这步和下一步，如何用st_addvar循环生成多个变量，我现在的做法只能生成一个变量t,如果我想通过循环生成t1 t2....到tk,应该怎样做
st_store(., "t",y)
i++
}

superzhousi 发表于 2010-8-30 22:47 小弟求教，如果在mata中有一个n*k的矩阵，我想把这个矩阵的数据传回到stata中，我知道使用st_addvar() 和st_store() 但是，如果我想一次生成k个变量，即矩阵的每一列都成为一个在stata中的变量而且列数如果很多， 这样应该如何来做呢？ 谢谢。。

*一个比较粗浅的作法：

*（1）在stata中
. set obs n

*（2）用generate生成k个变量（变量名自己随意定），不妨令它们都只有缺失值
*（3）进入mata，生成n*k矩阵m
*（4）使用st_store()
: st_store(.,.,m)
: end

*mata中的m的值便把stata中的k个变量的值代换了。

4# sungmoo
果然厉害，问题解决，没想到st_store(.,.,m)这样用，非常感谢

*一个比上面稍微复杂的例子（注意st_store()中各值对结果的作用）：

clear all
inp a b c d e
. . . . .
. . . . .
. . . . .
. . . . .
end

mata
z=(1,2\3,4)
st_store((2,3),("a","b"),z)
end

list

st_store(a,b,c)的关键是：
（1）a、b、c都是mata中的向量；
（2）以a与b在stata当前数据库中定义出一个“观测值区域”，该区域须与c同型；
（3）以c中的值替换该区域的值。

6# sungmoo
也在看这个东西，谨受教

set obs 30
mata
r1=30
c1=5
a=uniform(r1,c1)
d1=J(r1,1,0)
for (i = 1; i <= c1; i++) {
d1=a[.,i]
name = "x" + strofreal(i)
st_store(.,st_addvar("float",name),d1)
}
end

set obs 30
mata
r1=30
c1=5
a=uniform(r1,c1)
d1=J(r1,1,0)
for (i = 1; i <= c1; i++) {
d1=a[.,i]
name = "x" + strofreal(i)
st_store(.,st_addvar("float",name),d1)
}
end

*不用循环
clear all
set ob 30
mata
a=uniform(30,5)
st_store(.,st_addvar("float",("x1","x2","x3","x4","x5")),a)
end

*用循环
clear all
set ob 30
mata
a=uniform(30,5)
for(i=1;i<=5;i++) st_store(.,st_addvar("float","x"+strofreal(i)),a[.,i])
end