2016-10-20 3 views
0

Rのベータ回帰を実行するためにRパッケージ 'zoib'を使っています。41ページのサンプルR・ジャーナルに発表されたパッケージの作者:ベータ回帰のための 'zoib' Rパッケージを使用したときの 'ts'関数のエラー

Lui F and Kong Y. 2015. zoib: An R Package for Bayesian Inference for Beta Regression and Zero/One Inflated Beta Regression. The R Journal 7(2)

私は、彼らが使用してまったく同じデータとコードを使用しています信じて:

library(zoib) 
data("GasolineYield", package="zoib") 
GasolineYield$batch <- as.factor(GasolineYield$batch) 
d <- GasolineYield 
eg1.fixed <- zoib(yield ~ temp + as.factor(batch) | 1, data=GasolineYield, joint=FALSE, 
     random=0, EUID=1:nrow(d), zero.inflation=F, one.inflation=F, 
     n.iter=1050, n.thin=5, n.burn=50) 
sample1 <- eg1$coeff 
traceplot(sample1) 
autocorr.plot(sample1) 
gelman.diag(sample1) 

私がしようとしたときしかし、私はエラーを取得していますサンプルの診断プロットを行う。

Error in ts(seq(from = start(x), to = end(x), by = thin(x)), start = start(x), : 
    invalid time series parameters specified 

コードが機能していないか、私は問題を解決するために何ができるのか、なぜ私が理解することはできません:これはエラーメッセージです。私はzoibによって呼び出されるtime関数にエラーをトレースすることができます。サンプルオブジェクトにはtsp属性が含まれていないと思われるかもしれませんが、zoibパッケージの作成者はモデル出力がcodaと一緒に使用されるので、私は非常に混乱しています。私はMCMCや時系列オブジェクトでの作業経験があまりないので、明らかに何かが欠けているかもしれません。誰がパッケージの作者によって提供された例が失敗しているのか、そしてその解決策は何かを説明することはできますか?

答えて

1

私はパッケージの著者(Fang Liu)に電子メールを送って、実際に私が持っているパッケージのバージョンにバグがありましたが、最新のバージョンzoib(バージョ1.4.2)。最新のバージョンを使用して、コードが動作するようになりました。

関連する問題