rflexscan (version 0.1)

plot.rflexscan: Graph plotting of flexscan results

Description

Display detected clusters by a graph representation.

Usage

# S3 method for rflexscan
plot(x, rank = 1:length(x$cluster), pval = 1,
  vertexsize = max(x$coordinates[, 1]) - min(x$coordinates[, 1]),
  xlab = colnames(x$coordinates)[1], ylab = colnames(x$coordinates)[2],
  xlim = c(min(x$coordinates[, 1]), max(x$coordinates[, 1])),
  ylim = c(min(x$coordinates[, 2]), max(x$coordinates[, 2])), ...)

Arguments

x

An rflexscan object.

rank

An integer vector which specifies ranks of clusters to be displayed.

pval

A threshold of P-value. Clusters with P-values of <pval will be displayed.

vertexsize

Size of vertex of the graph.

xlab

A label of the x axis.

ylab

A label of the y axis.

xlim

The x limits of the plot.

ylim

The y limits of the plot.

...

Other parameters to be passed to plot.igraph function.

Details

Clusters are colored using the current palette. Please use palette function to specify colors of each cluster. Note that clusters with ranks larger than the number of colors in the palette are not highlighted.

See Also

flexscan

Examples

# NOT RUN {
# load sample data (North Carolina SIDS data)
library(spdep)
data("nc.sids")

# calculate the expected numbers of cases
expected <- nc.sids$BIR74 * sum(nc.sids$SID74) / sum(nc.sids$BIR74)

# run FleXScan
fls <- flexscan(x = nc.sids$x, y = nc.sids$y,
                observed = nc.sids$SID74,
                expected = expected,
                name = rownames(nc.sids),
                clustersize = 10,
                nb = ncCR85.nb)

# display all clusters
plot(fls)

# display clusters with rank 1, 2 and 3
plot(fls, rank = c(1, 2, 3))

# display clusters of P-value <= 0.05
plot(fls, pval = 0.05)

# }