data(chess)
aylmer.test(chess)
a <- matrix(1,9,9) # See Sloane's A110058
plot(randomprobs(a,1000),type="b",main="Importance of burn-in")
set.seed(0)
b <- diag(rep(6,6))
plot(randomprobs(b,B=1000,n=1000), type="b",main="Importance of burn-in, part II")
data(purum)
randomboards(purum,10)
Run the code above in your browser using DataLab