tropical.centroid {TML}

Tropical centroid-based sampling about a center of mass

Description

This function is a centroid-based HAR sampler about a center of mass denoted by a location parameter with scale parameter in terms of the tropical distance

Usage

VE.HAR.centroid(D_s, x0, I = 1, m, s, tadd = max)

Arguments

D_s

matrix of vertices of a tropical simplex; each row is a vertex

x0

initial point for sampler, numeric vector

I

number of states in Markov chain

m

location parameter; numeric vector indicating centroid

s

scale parameter; in terms of tropical distance

tadd

function; max indicates max-plus addition, min indicates min-plus addition. Defaults to max

Value

next sampled point from the tropical polytope

Author(s)

David Barnhill [e-mail]

Examples

D_s <-matrix(c(0,10,10,0,10,0,0,0,10),3,3,TRUE)
x0 <- c(0,0,0)
m <- c(0,5,5)
s <- 1
VE.HAR.centroid(D_s, x0, I = 50,m,s)
VE.HAR.centroid(D_s, x0, I = 50,m,s,tadd=min)