https://github.com/cran/ape
Tip revision: 10898aebdf6661a0b81ba21bf24969336b544a60 authored by Emmanuel Paradis on 21 December 2021, 08:20:05 UTC
version 5.6
version 5.6
Tip revision: 10898ae
mantel.test.Rd
\name{mantel.test}
\alias{mantel.test}
\title{Mantel Test for Similarity of Two Matrices}
\description{
This function computes Mantel's permutation test for similarity of two
matrices. It permutes the rows and columns of the second matrix
randomly and calculates a \eqn{Z}-statistic.
}
\usage{
mantel.test(m1, m2, nperm = 999, graph = FALSE,
alternative = "two.sided", ...)
}
\arguments{
\item{m1}{a numeric matrix giving a measure of pairwise distances,
correlations, or similarities among observations.}
\item{m2}{a second numeric matrix giving another measure of pairwise
distances, correlations, or similarities among observations.}
\item{nperm}{the number of times to permute the data.}
\item{graph}{a logical indicating whether to produce a summary graph
(by default the graph is not plotted).}
\item{alternative}{a character string defining the alternative
hypothesis: \code{"two.sided"} (default), \code{"less"},
\code{"greater"}, or any unambiguous abbreviation of these.}
\item{\dots}{further arguments to be passed to \code{plot()} (to add a
title, change the axis labels, and so on).}
}
\details{
The function calculates a \eqn{Z}-statistic for the Mantel test, equal to
the sum of the pairwise product of the lower triangles of the
permuted matrices, for each permutation of rows and columns. It
compares the permuted distribution with the \eqn{Z}-statistic observed
for the actual data.
The present implementation can analyse symmetric as well as (since
version 5.1 of \pkg{ape}) asymmetric matrices (see Mantel 1967,
Sects. 4 and 5). The diagonals of both matrices are ignored.
If \code{graph = TRUE}, the functions plots the density estimate of
the permutation distribution along with the observed \eqn{Z}-statistic
as a vertical line.
The \code{\dots} argument allows the user to give further options to
the \code{plot} function: the title main be changed with \code{main=},
the axis labels with \code{xlab =}, and \code{ylab =}, and so on.
}
\value{
\item{z.stat}{the \eqn{Z}-statistic (sum of rows*columns of lower
triangle) of the data matrices.}
\item{p}{\eqn{P}-value (quantile of the observed \eqn{Z}-statistic in
the permutation distribution).}
\item{alternative}{the alternative hypothesis.}
}
\references{
Mantel, N. (1967) The detection of disease clustering and a
generalized regression approach. \emph{Cancer Research}, \bold{27},
209--220.
Manly, B. F. J. (1986) \emph{Multivariate statistical methods: a primer.}
London: Chapman & Hall.
}
\author{
Original code in S by Ben Bolker, ported to \R by Julien Claude
}
\examples{
q1 <- matrix(runif(36), nrow = 6)
q2 <- matrix(runif(36), nrow = 6)
diag(q1) <- diag(q2) <- 0
mantel.test(q1, q2, graph = TRUE,
main = "Mantel test: a random example with 6 X 6 matrices
representing asymmetric relationships",
xlab = "z-statistic", ylab = "Density",
sub = "The vertical line shows the observed z-statistic")
}
\keyword{multivariate}