Displaying 1 result from an estimated 1 matches for "dul3".
Did you mean:
dul
2004 Feb 11
0
gelman.diag question
...doing wrong?
> dim(samp)
[1] 10000 86
> dul<-mcmc(b1<-samp[,1],start=1, end=10000, thin=10)
Error in data[1:nobs, , drop = FALSE] : incorrect number of dimensions
Then, as described below,I tried other way to create the mcmc object but I
get the message "Error in gelman.diag(dul3) : You need at least two chains".
I think I am not understanding how the Gelman-Rubin test works in R. Could
you give me a small example.
> dul2<-mcmc(b1<-as.matrix(samp[,1]),start=1, end=10000, thin=10)
> dul3<-matrix(1000,2,c(dul,dul2))
> is.mcmc(dul3)
[1] FALSE
> dul3...