Performs clustering of a matrix of 0s and 1s, ie. the censoring indicator columns for multiple variables. If there are multiple detection limits within a column first convert the 0/1 designated to be Below vs Above the highest detection limit in the column. Detection limits may differ among columns.

binaryClust(dat.frame, method = "ward.D2", group = NULL, ncluster = NULL)

Arguments

dat.frame

A data frame containing only the 0/1 columns, one column per chemical parameter.

method

Method of forming clusters. The default is "ward.D2", which is appropriate for a variety of types of data. Another appropriate option is “average” – average distances between cluster centers. See the vegan package for other possible clustering methods.

group

Optional grouping variable. If used, sites being clustered will be represented by their group name, rather than by the row number.

ncluster

Optional number of clusters to be differentiated on the graph. Clusters are fenced off with rectangles.

Value

Prints a cluster dendrogram based on clustering method

References

Helsel, D.R., 2011. Statistics for Censored Environmental Data using Minitab and R, 2nd ed. John Wiley & Sons, USA, N.J.

Examples

data(PbHeron)

# without group specified
binaryClust(PbHeron[,4:15])


# With Group argument
binaryClust(PbHeron[,4:15],group=PbHeron$DosageGroup)