作者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