plot.DiscrFact: plot Method for DiscrFact Objects

Description Usage Arguments Details Value Author(s) References Examples

Description

The plot method for class DiscrFact: Next to a plot of the tclust object which has been used for creating the DiscrFact object, a silhouette plot indicates the presence of groups with a large amount of doubtfully assigned observations. A third plot similar to the standard tclust plot serves to highlight the identified doubtful observations.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
## S3 method for class 'DiscrFact'
 plot(x, enum.plots = FALSE, ...)
plot_DiscrFact_p2 (x, xlab = "Discriminant Factor", 
                   ylab = "Clusters", main, xlim, 
                   print.Discr = TRUE, main.pre, ...)
                  
                  
                  
plot_DiscrFact_p3 (x, main = "Doubtful Assignments", col, pch, 
                   col.nodoubt = grey (0.8), by.cluster = FALSE, 
                   ...)
                   

Arguments

x

An object of class "DiscrFact" as from DiscrFact ().

enum.plots

A logical value indicating whether the plots shall be enumerated in their title ("(a)", "(b)", "(c)").

xlab, ylab, xlim

Arguments passed to funcion plot.tclust.

main

Argument passed to funcion plot.

print.Discr

A logical value indicating whether each clusters mean discriminant factor shall be plotted

main.pre

An optional string which is appended to the plot's caption.

pch, col

Arguments passed to function plot.

col.nodoubt

Color of all observations not considered as to be assigned doubtfully.

by.cluster

Logical value indicating whether parameters pch and col refer to observations (FALSE) or clusters (TRUE).

...

Arguments to be passed to or from other methods.

Details

plot.DiscrFact.p2 displays a silhouette plot based on the discriminant factors of the observations. A solution with many large discriminant factors is not reliable. Such clusters can be identified with this silhouette plot. Thus plot.DiscrFact.p3 displays the dataset, highlighting observations with discriminant factors greater than the given threshold. Function plot.DiscrFact combines the standard plot of a tclust object, and the two plots introduced here.

Value

No return value is provided.

Author(s)

Agustin Mayo Iscar, Luis Angel García Escudero, Heinrich Fritz

References

García-Escudero, L.A.; Gordaliza, A.; Matrán, C. and Mayo-Iscar, A. (2010), "Exploring the number of groups in robust model-based clustering." Statistics and Computing, (Forthcoming).
Preprint available at www.eio.uva.es/infor/personas/langel.html.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
sig <- diag (2)
cen <- rep (1, 2)
x <- rbind(mvtnorm::rmvnorm(360, cen * 0,   sig),
	       mvtnorm::rmvnorm(540, cen * 5,   sig * 6 - 2),
	       mvtnorm::rmvnorm(100, cen * 2.5, sig * 50)
)

clus.1 <- tclust (x, k = 2, alpha=0.1, restr.fact=12)
clus.2 <- tclust (x, k = 3, alpha=0.1, restr.fact=1)

dsc.1 <- DiscrFact (clus.1)
plot(dsc.1)

dsc.2 <- DiscrFact (clus.2)
plot (dsc.2)

dev.off ()
plot_DiscrFact_p2 (dsc.1)
plot_DiscrFact_p3 (dsc.2)

Example output

Loading required package: mvtnorm
Loading required package: sn
Loading required package: stats4

Attaching package: 'sn'

The following object is masked from 'package:stats':

    sd

Loading required package: mclust
Package 'mclust' version 5.4.1
Type 'citation("mclust")' for citing this R package in publications.
Loading required package: cluster
Warning message:
In .tclust.warn(O, O$ret) :
  The result is artificially constrained due to restr.fact = 1.
null device 
          1 

tclust documentation built on May 2, 2019, 3:20 a.m.