815102137 发表于 2014-8-25 21:45
执行V=st_view(X, ., .)时错误提示是“X not found”。之前已经按步骤创建了X矩阵啊。看了help文件,仍然 ...
不好意思,刚才弄错了,应该用v=st_matrix("x")把stata中矩阵转换到mata中,所以改写为:
v=st_matrix("x")
p=.
l=.
eigensystem(v, p,l)
p
l
结果为
p
-.723241249 .129430781 + .330324585i .129430781 - .330324585i
2 | -.580525181 .792993021 .792993021
3 | -.284953867 -.007865144 - .383946807i -.007865144 + .383946807i
4 | -.220148611 -.29553534 + .019452345i -.29553534 - .019452345i
5 | -.10123484 .008825556 + .100104836i .008825556 - .100104836i
。。。
l
1 2 3
+-------------------------------------------------------------------------------------
1 | 5.70476568 -.068718011 - 1.70348216i -.068718011 + 1.70348216i
+-------------------------------------------------------------------------------------
4 5
---------------------------------------------------------+
1 -.28366483 - .67491643i -.28366483 + .67491643i |
---------------------------------------------------------+