WECCAhc: Weighted clustering of array CGH data.

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

Description

Weighted distance matrix calculation, followed by dendrogram construction for regioned hard called array CGH data.

Usage

1
WECCAhc(cghdata.regioned, dist.measure = "agree", linkage = "ward", weight.type = "all.equal")

Arguments

cghdata.regioned

A list-object as returned by the function regioning.

dist.measure

The distance measure to be used. This is either "agree" (the agreement similarity) or "conc" (the concordance similarity).

linkage

The linkage method to be used, like "single", "average", "complete", "ward", "total".

weight.type

Region weighting to be used in the calculation of the distance, either "all.equal" or "heterogeneity".

Details

The distance between the call profiles of all sample pairs are calculated using special distance measures suitable for this data type. The distance between the call signature of two regions is either defined via the agreement or concordance similarity. The distance between two call profiles is defined as a weighted average of the distances between individual regions. Region weighing currently implemented: no-weighing (all regions contribute equally to the overall distance), or heterogeneity weighing (weights are proportional to their Shannon's entropy). Once the distance matrix has been calculated, this is submitted to the hclust function, which clusters the samples hierarchically, returning a dendrogram.

Value

An object of class hclust which describes the tree produced by the clustering process. See hclust for a list of its components.

Author(s)

Wessel N. van Wieringen: w.vanwieringen@vumc.nl

References

Van Wieringen, W.N., Van de Wiel, M.A., Ylstra, B. (2008), "Weighted clustering of called aCGH data", Biostatistics, 9(3), 484-500.

See Also

dist, hclust, KLdiv, regioning

Examples

1
2
3
4
5
6
7
8
# generate object of cghCall-class
data(WiltingCalled)

# make region data (soft and hard calls)
WiltingRegioned <- regioning(WiltingCalled)

# clustering with soft.calls
dendrogram <- WECCAhc(WiltingRegioned)

tgac-vumc/WECCA documentation built on May 31, 2019, 9 a.m.