[問題] 在R上使用R2OpenBUGS如何使程式跑得更快?

[phil5566 (5566)]

[問題類型]:效能諮詢


[軟體熟悉度]:新手(沒寫過程式)


[問題敘述]:

我想在R上使用R2OpenBUGS這個package,使得在R上操作openbugs進行模擬,

之前測試過跑完一次要5~6天,但這個耗時好像太久,小弟剩的時間不多

想請教版上高手和先進們,要如何讓程式運作的更快?

在下面附上的code裡,主要定義了兩個函數分別是boot.CI(用來計算信賴區間)

和mcmc.fun(用來操作openbugs),其中mcmc.fun裡

1.POW(a,b)的意思就是a的b次方,打a^b的話,openbugs會無法辨別

2.dat[dat==1/0]=0,就是把data裡inf換成0,否則inf,openbugs會無法辨別

大概是這樣,謝謝




[程式範例]:

http://pastie.org/10991607



[環境敘述]:

R x64 3.3.1,

"R2OpenBUGS"package



[關鍵字]:

R2OpenBUGS



--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 36.224.29.126
※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1485022460.A.801.html
→ celestialgod: 它是call openbug幫你跑，你要去問怎麼加速openbug
→ celestialgod: 才行
→ celestialgod: 我記得我學校老師最後是自己用fortran寫Bayesian的
→ celestialgod: 模擬
→ phil5566: 感謝C大回答,有問題我會再發問~
→ clickhere: 不用加速openbug. 這種iid情況,用多台跑即可.
推 psinqoo: fortran<-完全不想回憶~~天阿~~
→ HumuHumu: R20