mahala: Mahalanobis Distance

In metan: Multi Environment Trials Analysis

Mahalanobis Distance

Description

Compute the Mahalanobis distance of all pairwise rows in .means. The result is a symmetric matrix containing the distances that may be used for hierarchical clustering.

Usage

mahala(.means, covar, inverted = FALSE)

Arguments

.means	A matrix of data with, say, p columns.
covar	The covariance matrix.
inverted	Logical argument. If TRUE, covar is supposed to contain the inverse of the covariance matrix.

Value

A symmetric matrix with the Mahalanobis' distance.

Author(s)

Tiago Olivoto [e-mail]

Examples

library(metan)
library(dplyr)
# Compute the mean for genotypes
means <- mean_by(data_ge, GEN) %>%
         column_to_rownames("GEN")

# Compute the covariance matrix
covmat <- cov(means)

# Compute the distance
dist <- mahala(means, covmat)

# Dendrogram
dend <- dist %>%
        as.dist() %>%
        hclust() %>%
        as.dendrogram()
plot(dend)

metan documentation built on March 7, 2023, 5:34 p.m.