stata中如何得到一个现有矩阵的特征值和特征向量

如题，在以前帖子中看到1.可将Excel中数据粘到stata的data editor 中，再执行命令“mkmat var*，mat(X)”即可得到矩阵X

2.计算手动输入的矩阵特征值和特征向量的命令为：
   mata
x=(1,2,3,4,5\1/3,1,5,3,4\ 1/2,1/5,1,3,2\1/3,1/3,1/3,1,5\1/5,1/4,1/2,1/5,1)   //手工输入5*5矩阵
p=.
l=.
eigensystem(x,p,l)
p   //显示特征向量
l    //显示特征值
end
我的是已经粘贴好的矩阵X，就没有执行手工输入矩阵的命令，直接执行下边的命令，结果出错了。请问如何解呀！求大神解救。

好像要先讲stata中的矩阵转换成mata中的矩阵，现在mata中设置一个只有缺失值的矩阵，然后用st_view把stata中X矩阵的数据转换到mata中你设置的新矩阵中，然后就可以用求特征根的命令了。
思路是这样的，不知对你有没有帮助，自己查查help文档。
mata
V=st_view(X, ., .)
p=.
l=.
eigensystem(V, p,l)

pkl 发表于 2014-8-25 10:19
好像要先讲stata中的矩阵转换成mata中的矩阵，现在mata中设置一个只有缺失值的矩阵，然后用st_view把stata中 ...

非常感谢你，我好好学习一下。

pkl 发表于 2014-8-25 10:19
好像要先讲stata中的矩阵转换成mata中的矩阵，现在mata中设置一个只有缺失值的矩阵，然后用st_view把stata中 ...

执行V=st_view(X, ., .)时错误提示是“X not found”。之前已经按步骤创建了X矩阵啊。看了help文件，仍然没弄好。能麻烦您再指点一下吗？多谢了。

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  |
     ---------------------------------------------------------+

pkl 发表于 2014-8-25 22:31
不好意思，刚才弄错了，应该用v=st_matrix("x")把stata中矩阵转换到mata中，所以改写为：
v=st_matrix(" ...

成功了！非常感谢你帮助我，谢谢！你太赞了！

pkl 发表于 2014-8-25 22:31
不好意思，刚才弄错了，应该用v=st_matrix("x")把stata中矩阵转换到mata中，所以改写为：
v=st_matrix(" ...

能再问您一个问题吗？得到的特征矩阵怎样才能使得在mata模式之外也可用。退出mata模式后，我执行mat list p，显示的是“p not found”。

815102137 发表于 2014-8-27 08:22
能再问您一个问题吗？得到的特征矩阵怎样才能使得在mata模式之外也可用。退出mata模式后，我执行mat list ...

Mata是stata的一个编译程序，其实是相当于两个编译器，所以，stata中的变量在Mata中不能直接引用，Mata中的变量在stata中也不能直接引用，所以要转换，你看看Mata手册吧。

pkl 发表于 2014-8-27 12:58
Mata是stata的一个编译程序，其实是相当于两个编译器，所以，stata中的变量在Mata中不能直接引用，Mata中 ...

好的，谢谢。

这是我找到的，供后来读贴的人们使用吧。