S4 defintiions of mcmc and mcmc.list

gistfile1.rebol

##' Redefine mcmc as S4 class
setClass("mcmc", contains="matrix",
         representation(start="numeric",
                        end="numeric",
                        thin="numeric"))
setOldClass("mcmc", S4Class="mcmc")

setClass("mcmc.list", contains="list",
         validity = function(object) {
             if (length([email protected]) == 0) {
                 TRUE
             } else {
                 allmcmc <- all(lapply([email protected], is, class2="mcmc"))
                 if (allmcmc) {
                     TRUE
                 } else {
                     "Not all elements in the list are mcmc objects."
                 }
             }
         })
setOldClass("mcmc.list", S4Class="mcmc.list")