NormalizeInfo: Normalize tree distances

In TreeDist: Calculate and Map Distances Between Phylogenetic Trees

Normalize tree distances

Description

NormalizeInfo() is an internal function used to normalize information against a reference, such as the total information present in a pair of trees.

Usage

NormalizeInfo(
  unnormalized,
  tree1,
  tree2,
  InfoInTree,
  infoInBoth = NULL,
  how = TRUE,
  Combine = "+",
  ...
)

Arguments

unnormalized	Numeric value, vector or matrix to be normalized.
tree1, tree2	Trees from which unnormalized was calculated.
InfoInTree	Function to calculate the information content of each tree.
infoInBoth	Optional numeric specifying information content of both trees independently. If unspecified (NULL), this will be calculated using the method specified by how.
how	Method for normalization, perhaps specified using the normalize argument to a tree distance function. See details for options.
...	Additional parameters to InfoInTree() or how().

Details

The unnormalized value(s) are normalized by dividing by a denominator calculated based on the how parameter. Valid options include:

FALSE

No normalization is performed; the unnormalized values are returned.

TRUE

Unless infoInBoth is specified, the information in each tree is computed using InfoInTree(), and the two values combined using Combine().

A numeric value, vector or matrix

how is used as the denominator; the returned value is unnormalized / how.

A function

Unless infoInBoth is specified, the information in each tree is computed using InfoInTree(), and the two values combined using how. NormalizeInfo(how = Func) is thus equivalent to NormalizeInfo(how = TRUE, Combine = Func).

Value

NormalizeInfo() returns an object corresponding to the normalized values of unnormalized.

Author(s)

Martin R. Smith ([e-mail])

Examples

library("TreeTools", quietly = TRUE)
pair1 <- c(BalancedTree(9), StarTree(9))
pair2 <- c(BalancedTree(9), PectinateTree(9))

# We'll let the number of nodes define the total information in a tree
Nnode(pair1)
Nnode(pair2)

# Let's normalize a unit distance
rawDist <- cbind(c(1, 1), c(1, 1))

# With `Combine = "+"`, the maximum distance is the sum of
# the information in each tree
denominator <- outer(Nnode(pair1), Nnode(pair2), "+")

NormalizeInfo(rawDist, pair1, pair2, InfoInTree = ape::Nnode, Combine = "+")
rawDist / denominator


# A denominator can be specified manually using `how`:
NormalizeInfo(rawDist, pair1, pair2, InfoInTree = ape::Nnode, how = 16)
rawDist / 16


# `how` also allows the denominator to be computed from trees:
outer(Nnode(pair1), Nnode(pair2), pmin)
NormalizeInfo(rawDist, pair1, pair2, InfoInTree = ape::Nnode, how = pmin)
rawDist / outer(Nnode(pair1), Nnode(pair2), pmin)

TreeDist documentation built on Sept. 11, 2024, 9:10 p.m.