※ 引述《lianhua (蓮華)》之銘言：
: [問題類型]:
: 程式諮詢
: [軟體熟悉度]:
: 使用者
: [問題敘述]:
: 因為我需要跑大量的貝氏模擬，大迴圈次數500次
: 大迴圈裡面的小迴圈20000次
: 每一次小迴圈都需要生成一個5XM (M大概100左右) 的矩陣
: 但是我的電腦大迴圈大概都跑個200次左右就會跳出
: 錯誤: 無法配置大小為 526.2 Mb 的向量
: 這種訊息 請問是我的電腦記憶體不足，還是語法太拙劣？
: [程式範例]:
: 程式很長，所以僅貼部分
: N = 500
: L0 = numeric()
: L1 = numeric()
: L2 = numeric()
: Alpha1 = numeric()
: Alpha2 = numeric()
關於numeric()這部分我可以給點小建議，
如果可以的話，還是盡量給定向量長度。
http://statlab.nchc.org.tw/rnotes/?page_id=233
這篇文章看到「盡量避免讓 R 自動增加向量長度」這段，
可以很明顯看到有給定向量長度及沒給定向量長度的差別。
: for(k in 1:N){
: set.seed(k)
: ．
: ．
: ．
: 這裡是我懷疑出錯的部分，因為如果我不是生成隨機矩陣的話
: 就不會有這個問題
: iter = 20000
: repeat{
: ### step 1 fill up Z ###
: g1 = function(x){
: ((l1+l0)*c1*(x^(c1-1)))/((l1+l0)*c1*(x^(c1-1))+l2*c2*(x^(c2-1)))
: }
: g2 = function(x){
: (l1*c1*(x^(c1-1)))/(l1*c1*(x^(c1-1))+(l2+l0)*c2*(x^(c2-1)))
: }
: Z1 = matrix(,ncol=2,nrow=length(t1))
: Z2 = matrix(,ncol=2,nrow=length(t2))
: for(i in 1:length(t1)){
: Z1[i,] = rmultinom(1,1,c(g1(t1[i]),(1-g1(t1[i]))))
: }
: for(i in 1:length(t2)){
: Z2[i,] = rmultinom(1,1,c(g2(t2[i]),(1-g2(t2[i]))))
: }
: }
: 想請各位先進為我解答，感謝～
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 114.34.27.68