Returns a set of indices of the original data set that maximizes the mean and minimizes the variance of the distances between pairs of plot records.

1 | ```
hcr(d, nout, nsampl=1000)
``` |

`d` |
An object of class |

`nout` |
The number of sites (plot records) to be chosen among those available in |

`nsampl` |
The number of resampling trials to be compared. |

Many subsets of the input data are selected randomly. These subsets are sorted by decreasing mean dissimilarity between pairs of plot records, and then sorted again by increasing variance of these dissimilarities. Ranks from both sortings are summed for each subset, and the subset with the lowest summed rank is considered as the most representative.

Returns a vector containing the indices of the selected sites (plot records) to be used for sub-setting the original table.

Miquel De Cáceres, Forest Science Center of Catalonia

Lengyel, A., Chytry, M., Tichy, L. (2011) Heterogeneity-constrained random resampling of phytosociological databases. Journal of Vegetation Science 22: 175-183.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 | ```
## Loads data (38 columns and 33 species)
data(wetland)
dim(wetland)
## Constructs the chord distance matrix
## (see also \code{\link{decostand}} in package vegan)
wetland.chord =dist(as.data.frame(sweep(as.matrix(wetland), 1,
sqrt(rowSums(as.matrix(wetland)^2)), "/")))
## Performs HCR resampling. Returns indices of objects
sel = hcr(wetland.chord, nout=20, nsampl=1000)
## Prints the names of the plot records
print(row.names(wetland)[sel])
## Subset the original distance matrix
sel.chord = as.dist(as.matrix(wetland.chord)[sel,sel])
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.