0

我正在使用该ggmcmc程序包生成使用该函数的rjags程序包输出的摘要 pdf 文件。ggmcmc()但是,我收到以下错误消息:

> 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

谁能帮我确定错误是在哪里引起的以及如何纠正它?我错过了一些明显的东西吗?

4

1 回答 1

2

ggmcmc 0.5.1 以与以前版本不同的方式解决了 bin 数量的计算。以前的版本依赖于 ggplot2:::bin,而 0.5.1 自己计算 bin 和它们的 binwidth。

您的情况可能是某些参数的范围非常极端,以至于舍入误差会使它们中的一些具有更多或更少的 bin,因此会产生此错误。

于 2013-10-04T09:15:28.367 回答