Rのベータ回帰を実行するためにRパッケージ 'zoib'を使っています。41ページのサンプルR・ジャーナルに発表されたパッケージの作者:ベータ回帰のための 'zoib' Rパッケージを使用したときの 'ts'関数のエラー
私は、彼らが使用してまったく同じデータとコードを使用しています信じて:
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や時系列オブジェクトでの作業経験があまりないので、明らかに何かが欠けているかもしれません。誰がパッケージの作者によって提供された例が失敗しているのか、そしてその解決策は何かを説明することはできますか?