2016-10-20 6 views
0

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) 

을하지만 :

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)

는 나는 그들이 사용하는 동일한 데이터와 코드를 사용하고 생각 샘플에 대한 진단 플롯을 수행합니다.

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). 최신 버전을 사용하여 코드가 작동합니다.