summary.swap: Summarize Gibbs samples for a k-gene model
Description
Calculate marginal posteriors for location, posterior means for
coefficients, and the Bayes Factor for k vs k-1 genesUsage
## S3 method for class 'swap':
summary(object, method=NULL, ncoef=length(sc.obj$alt.coef), nloc=sc.obj$nloc)
Arguments
ncoef
Optional. The number of coefficients in the class of
models. Typically, 2*nloc for method=="F2" and
nloc for all other methods
nloc
Optional. The number of loci in the sample space.
Value
- A list with components:
- loc.posteriorA vector of (marginal) posterior odds for each
locus compared to a no gene model
- coefsPosterior means of coefficients.
- ratioA list with components
mean, an estimate of the
Bayes Factor for k versus k-1 gene models, and var, an estimate
of its variance.
synopsis
summary.swap(object, method, ncoef, nloc, ...)