作者wcchang (wcchang)
看板Statistics
标题Re: [问题] SAS关於PROC IML大型反矩阵运算
时间Sun Apr 27 10:12:14 2008
我要算的是Canonical correlation
有先试过5000*5000左右的矩阵
但因为singular的关系inverse无法找到
想用g-inverse 找但有点困难...
目前以变异数矩阵前後多乘 Diagonal matrix ##(-0.5)方式
切割可以算吗? GENE DATA...
※ 引述《granzi (乌木)》之铭言：
: ※ 引述《wcchang (wcchang)》之铭言：
: : 有没有人跑过SAS超过6000*6000的矩阵
: : 找CORR及INVERSE都有记忆体不够的问题耶!
: : "ERROR: (execution) Unable to allocate sufficient memory.
: : At least 437961640 more bytes required."
: : MATLAB 也似乎无法处理这麽大的矩阵
: : (CORR部分我用回圈一一处理元素 勉强可以找到
: : 但反矩阵就没法度了)
: 看样子你需要4.3G以上的记忆体，这在一般的PC上是不可能的。
: 为什麽要跑这麽大的Corr及inverse(我猜你要算的是多变量常态分布的covariance
: matrix 好算PDF的值）？
: 或许可以先切割一下
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.116.194.56
1F：推 tew:试试看SAS内建的package 04/27 13:00
2F：推 tew:proc cancorr; 04/27 13:00
3F：推 granzi:之前搞错了，是437MB的记忆体，换台电脑或许可以跑得动 04/27 14:56
4F：→ wcchang:谢谢两位!proc cancorr可以跑,可惜Var>>N,cancorr=1... 04/27 17:12
5F：→ wcchang:真的要用G-inverse了吗... 04/27 17:15