mysum: summaries of MCMC output

In bbolker/mixstock: Functions for Mixed Stock Analysis

summaries of MCMC output

Description

calculates summary statistics of one or more columns of MCMC output

Usage

mysum(x, names = NULL)
mysumvec(x,names=NULL,levels=c(0.95,0.9))

Arguments

x	a numeric matrix representing MCMC output
names	vector of column names
levels	numeric vector of confidence levels

Details

a slight variant on the standard numeric summary: calculates mean, median, standard deviation, and lower and upper 90 and 95 percent quantiles. mysum gives a table (original data as rows, various stats as columns), while mysumvec collapses the answer to a vector

Value

a numeric vector of results, appropriately named

Examples

  x = matrix(runif(1000),ncol=2,dimnames=list(NULL,c("A","B")))
  mysum(x)
  mysumvec(x)

bbolker/mixstock documentation built on July 23, 2024, 12:18 p.m.