Ultrametrics.HAR {TML}

Hit-and-Run Sampler for the space of ultrametrics

Description

This sampler samples a point in the space of ultrametrics where each point represents an equidistant tree on n leaves

Usage

Ultrametrics.HAR(x0, n, I = 1, h = 1)

Arguments

x0

an equidistant tree defined as ultrametric

n

number of leaves for the equidistant tree

I

number of states in the Markov chain

h

height of phylogenetic tree

Value

point in the space of ultrametrics over n leaves

Author(s)

Ruriko Yoshida [e-mail]

References

Yoshida, Ruriko, Keiji Miura and David Barnhill (2022). Hit and Run Sampling from Tropically Convex Sets.

Examples

x0 <-Sim_Trees15[1,]
n<-10

Ultrametrics.HAR(x0, n, I = 50, h = 1)