2016-01-15 58 views
0

我通過R使用BUGS軟件進行貝葉斯分析,並利用ggmcmc包進行貝葉斯推斷。修改ggmcmc圖中的家族參數

在我最近的例子中,我有一個監視器下的參數矩陣b,尺寸爲5x8。現在,如果我直接從ggmcmc包中使用一個圖,那麼參數非常多,以至於我不能在輸出後面的圖中看到一些東西。

例如ggs_histogram

enter image description here

現在ggmcmc積函數有一個參數調用的家人和你用這個來選擇的參數的子集的情節,包括。在官方軟件包頁面中,它表示你必須設置家族等於一個正則表達式,它與你想要的參數匹配,如果你有參數a,b並且你想繪製b(family ='b'), 。

現在我想從我提到的b矩陣中只繪製一列元素,例如b [1,1],b [2,1],b [3,1],...,b [ 8,1]

所以我試圖以通常的方式子集,像family ='b [,1]'。

錯誤seq.default(MN,MX,由= BW):

'從'

不能NA,NaN或無限

此外:警告信息:

1 :在min(x)中:沒有非缺少的參數min;返回Inf

2:在max(x)中:沒有非缺失參數爲max;返回-Inf

任何想法?也許正確的正則表達式或ggplot facet_grid運球?

+1

如果您包含可重複的示例,它將更容易幫助您。 – MLavoie

+1

添加最少的帖子,與您的代碼相同,並添加一些數據片段。我們無法讀懂頭腦。 –

+0

我故意跳過共享數據和可重複代碼的部分,因爲它會產生更多問題,並且會讓事情變得更加混亂。沒關係,我解決了它 –

回答

2

最後,ggmcmc包官方pdf文件有我想要的所有信息。我對正則表達式的需求是正確的,並且該包的教程對於正則表達式預期具有的形式非常有用。

所以,如果我想假設來推斷關於參數矩陣的第一列的元素,

family='b\\[.,1\\]' 

會做的工作就好了。這適用於ggmcmc軟件包的任何推理功能。