klastorin <- function(phy)
    .Defunct(msg = '\'klastorin\' has been removed from ape,
    see help("ape-defunct") for details.')

heterozygosity <- function (x, variance = FALSE)
    .Defunct(msg = '\'heterozygosity\' has been moved from ape to pegas,
    see help("ape-defunct") for details.')

H <- function(x, variance = FALSE)
    heterozygosity (x, variance = FALSE)

nuc.div <- function(x, variance = FALSE, pairwise.deletion = FALSE)
    .Defunct(msg = '\'nuc.div\' has been moved from ape to pegas,
    see help("ape-defunct") for details.')

theta.h <- function(x, standard.error = FALSE)
    .Defunct(msg = '\'theta.h\' has been moved from ape to pegas,
    see help("ape-defunct") for details.')

theta.k <- function(x, n = NULL, k = NULL)
    .Defunct(msg = '\'theta.k\' has been moved from ape to pegas,
    see help("ape-defunct") for details.')

theta.s <- function(s, n, variance = FALSE)
    .Defunct(msg = '\'theta.s\' has been moved from ape to pegas,
    see help("ape-defunct") for details.')