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

Description Usage Arguments Details Value Author(s) References See Also Examples

Description

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

Usage

1
2
3
4
## 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, sara.wade@eng.cam.ac.uk

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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
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.