mysum {mixstock}R Documentation

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)

Arguments

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

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)

[Package mixstock version 0.9 Index]