Re: [問題] SAS關於PROC IML大型反矩陣運算

看板Statistics作者 (wcchang)時間17年前 (2008/04/27 10:12), 編輯推噓3(302)
留言5則, 3人參與, 最新討論串3/3 (看更多)
我要算的是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

04/27 13:00, , 1F
試試看SAS內建的package
04/27 13:00, 1F

04/27 13:00, , 2F
proc cancorr;
04/27 13:00, 2F

04/27 14:56, , 3F
之前搞錯了,是437MB的記憶體,換台電腦或許可以跑得動
04/27 14:56, 3F

04/27 17:12, , 4F
謝謝兩位!proc cancorr可以跑,可惜Var>>N,cancorr=1...
04/27 17:12, 4F

04/27 17:15, , 5F
真的要用G-inverse了嗎...
04/27 17:15, 5F
文章代碼(AID): #184-3_hc (Statistics)
文章代碼(AID): #184-3_hc (Statistics)