stata的矩阵问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

stata中设立矩阵后，能否将该矩阵某行某列的数值赋予某个暂元吗？
如果能，具体怎么操作啊？？
还有stata中的矩阵，能够提取某一行或某一列，然后赋予到新的向量（矩阵）去吗？？就像matlab那样

请不要涉及mata，只在stata下能否完成上述任务
求教各位！！！！

扫码加我 拉你入群

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

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

分享1 收藏8 回帖

关键词：Stata tata MATLAB atlab matla matlab

唉，都问了n个问题了，没一个人理{:soso_e109:}

stata中设立矩阵后，能否将该矩阵某行某列的数值赋予某个暂元吗？

*把矩阵的第1行第2列元素赋予宏a
mat a=x[1,2]
loc a=det(a)

还有stata中的矩阵，能够提取某一行或某一列，然后赋予到新的向量（矩阵）去吗？？

*令矩阵x的第2列成为矩阵a
mat a=x[1...,2..2]

sungmoo 发表于 2012-5-17 02:42
*令矩阵x的第2列成为矩阵a
mat a=x[1...,2..2]

感谢斑竹！！
那个赋予宏的，我后来试出来了，原来以前一直是使用圆括号的，难怪老赋不上去。

但是那个提取行（列）向量的，中间的“..."是什么含义呢，没看太明白

中间的“..."是什么含义呢

x[1..2,3..4]即x的“第1行到第2行、第3列到第4列”（交集）的所有元素形成的矩阵。

若".."后面仍为"."，则其表示“（到）最后一行（列）”。

sungmoo 发表于 2012-5-17 15:02
x[1..2,3..4]即x的“第1行到第2行、第3列到第4列”（交集）的所有元素形成的矩阵。

若".."后面仍为" ...

真心感谢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'"
        }
}