使用ggmcmc包時出錯

Question

我正在使用ggmcmc包生成使用ggmcmc()函數生成rjags包輸出的摘要pdf文件。不過，我得到了以下錯誤消息：

> ggmcmc(x, file = "Model0-output.pdf") 
Plotting histograms 
Error in data.frame(..., check.names = FALSE) : 
arguments imply differing number of rows: 160, 164 

當我檢查我與ggs()功能創建的輸入數據框的結構，一切看起來是正確的。

> str(x) 
'data.frame': 240000 obs. of 4 variables: 
$ Iteration: int 1 2 3 4 5 6 7 8 9 10 ... 
$ Chain : int 1 1 1 1 1 1 1 1 1 1 ... 
$ Parameter: Factor w/ 32 levels "N[1]","N[2]",..: 1 1 1 1 1 1 1 1 1 1 ... 
$ value : num 96 87 76 79 89 95 85 78 86 89 ... 
- attr(*, "nChains")= int 3 
- attr(*, "nParameters")= int 32 
- attr(*, "nIterations")= int 2500 
- attr(*, "nBurnin")= num 2000 
- attr(*, "nThin")= num 2 
- attr(*, "description")= chr "postout0" 
- attr(*, "parallel")= logi FALSE 

任何人都可以幫助我確定錯誤是在哪裏造成的，我如何糾正它？我錯過了明顯的東西嗎？

Answer 1

ggmcmc 0.5.1以與先前版本中不同的方式計算垃圾箱數量。以前的版本依賴於ggplot2 ::: bin，而0.5.1自己計算bin和binwidth。

很可能你的情況是，某些參數的範圍太過極端，舍入誤差會使其中一些有一個或多個分箱，因此產生這個誤差。