Ошибка в функции ts при использовании пакета «zoib» R для бета-регрессии

Question

Я работаю с пакетом «zoib» R для выполнения бета-регрессии в R. Я пытаюсь воспроизвести пример, приведенный на странице 41 в статье авторы пакета, опубликованные в газете R Journal:

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 

Я не могу понять, почему код не работает или то, что я могу сделать, чтобы исправить эту проблему. Я могу отследить ошибку до функции time, которая вызывается по zoib, и кажется, что, возможно, проблема заключается в том, что объект-образец не имеет атрибута tsp, но авторы пакета zoib дают понять, что их вывод модели предназначен для можно использовать с coda, поэтому я очень смущен. У меня нет большого опыта работы с MCMC или объектами временных рядов, поэтому, возможно, я просто пропустил что-то очевидное. Может ли кто-нибудь объяснить, почему неудачный пример, предоставленный авторами пакета, и что такое решение?

Answer 1

Я отправил электронное письмо автору пакета (Fang Liu), и она сообщила мне, что на самом деле была ошибка в версии пакета, который у меня есть, но что ошибка исправлена ​​в последней версии zoib (версия 1.4.2). Используя самую последнюю версию, код теперь работает.