0

関数を使用してパッケージ出力ggmcmcの要約pdfファイルを生成するためにパッケージを使用しています。ただし、次のエラー メッセージが表示されます。rjagsggmcmc()

> 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 は、以前のバージョンとは異なる方法でビン数の計算を解決します。以前のバージョンは ggplot2:::bin に依存していましたが、0.5.1 はビンとそのビン幅を単独で計算します。

一部のパラメータの範囲が非常に極端であるため、丸め誤差によって一部のパラメータのビンが 1 つ多いまたは 1 つ少なくなるため、このエラーが発生する可能性があります。

于 2013-10-04T09:15:28.367 に答える