summary.c.estimate: Summarize and plot the estimate of the partition

In muschellij2/mcclust.ext: Point estimation and credible balls for Bayesian cluster analysis

Description

Summarizes and plots the estimate of the partition in a Bayesian cluster analysis model.

Usage

## S3 method for class 'c.estimate'
summary(object, ...)
## S3 method for class 'c.estimate'
plot(x,data=NULL,dx=NULL,xgrid=NULL,dxgrid=NULL,...)

Arguments

object	an object of class "c.estimate", i.e. a clustering estimate .
x	an object of class "c.estimate", i.e. a clustering estimate .
data	the dataset contained in a data.frame with ncol(x$cl) rows of data points.
dx	for ncol(data)=1, the estimated density at the observed data points.
xgrid	for ncol(data)=1, a grid of data points for density estimation.
dxgrid	for ncol(data)=1, the estimated density at the grid of data points.
...	other inputs to summary or plot.

Details

Summarizes and plots the clustering estimates returned by the functions minVI and minbinder.ext. In plots, data points are colored according to cluster membership. For nrow(x)=1, the data points are plotted against the density (which is estimated via a call to density if not provided). For nrow(x)=2 the data points are plotted, and for nrow(x)>2, the data points are plotted in the space spanned by the first two principal components.

Value

method	the optimization method used to obtain the point estimate.
k	(a vector of) the number of clusters in the point estimate. Returns a vector if n.c>0.
n.c	the number of point estimates in the object.
t	a list of length n.c of the table(s) with cluster sizes.
value	(a vector of) the posterior expected loss. Returns a vector if n.c>0.

Author(s)

Sara Wade, [e-mail]

References

Wade, S. and Ghahramani, Z. (2015) Bayesian cluster analysis: Point estimation and credible balls. Submitted. arXiv:1505.03339.

See Also

minVI and minbinder.ext

Examples

data(galaxy.draw)
data(galaxy.fit)
data(galaxy.pred)
x=data.frame(x=galaxy.fit$x)

# Find representative partition of posterior
psm=comp.psm(galaxy.draw)
galaxy.VI=minVI(psm,galaxy.draw,method=("all"),include.greedy=TRUE)
summary(galaxy.VI)
plot(galaxy.VI,data=x,dx=galaxy.fit$fx,xgrid=galaxy.pred$x,dxgrid=galaxy.pred$fx)
galaxy.B=minbinder.ext(psm,galaxy.draw,method=("all"),include.greedy=TRUE)
summary(galaxy.B)
plot(galaxy.B,data=x,dx=galaxy.fit$fx,xgrid=galaxy.pred$x,dxgrid=galaxy.pred$fx)

muschellij2/mcclust.ext documentation built on May 26, 2019, 9:36 a.m.