作者lianhua (蓮華)
看板R_Language
標題[問題] 記憶體不足
時間Sat Jun 22 08:00:34 2013
[問題類型]:
程式諮詢
[軟體熟悉度]:
使用者
[問題敘述]:
因為我需要跑大量的貝氏模擬，大迴圈次數500次
大迴圈裡面的小迴圈20000次
每一次小迴圈都需要生成一個5XM (M大概100左右) 的矩陣
但是我的電腦大迴圈大概都跑個200次左右就會跳出
錯誤: 無法配置大小為 526.2 Mb 的向量
這種訊息 請問是我的電腦記憶體不足，還是語法太拙劣？
[程式範例]:
程式很長，所以僅貼部分
N = 500
L0 = numeric()
L1 = numeric()
L2 = numeric()
Alpha1 = numeric()
Alpha2 = numeric()
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: 140.115.45.220
※ 編輯: lianhua 來自: 140.115.45.220 (06/22 08:01)
1F：→ celestialgod:會跳出那個訊息 基本上是記憶體不足 06/22 20:02