mahala: Mahalanobis Distance
In TiagoOlivoto/WAASB: Multi Environment Trials Analysis

mahala

R Documentation

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 tiagoolivoto@gmail.com

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)

TiagoOlivoto/WAASB documentation built on March 20, 2024, 4:18 p.m.