all.jl

abstract type AbstractRejectionChains <: AbstractMCMC.AbstractChains end

include("rejection.jl")
include("tables.jl")
include("stats.jl")


function Base.show(io::IO, c::AbstractRejectionChains)
    print(io, "Chains (×$(nchains(c))) {", eltype(c), "} $(length(c)) elements. ")
end

function Base.display(c::AbstractRejectionChains)
    println(c)
    if is_multilevel(c)
        ls = maximum.(levels(c))
        if all( ls[2:end] .== ls[1] )
            println(" proxy level L=", ls[1])
        else
            println(" proxy levels ", ls)
        end
    end
    data = info(c)
    for k in keys(data)
        if nchains(c) == 1
            println(" ", k, " : ", data[k][1])
        else
            m, s = mean(data[k]), std(data[k])
            m = round.(m, sigdigits=5)
            s = round.(s, sigdigits=3)
            println(" ", k, " : ", m, " ± ", s)
        end
    end
end