def post_mean(self, chain=None, tune=None):
    """Calculate the mean of the posterior distribution
    """

    tune = tune or self.get_tune
    chain = chain or self.get_chain()[-tune:]

    return chain.reshape(-1, chain.shape[-1]).mean(axis=0)