Mv: Continuous Jump heights function

In konkam/BNPdensity: Ferguson-Klass Type Algorithm for Posterior Normalized Random Measures

Continuous Jump heights function

Description

This function evaluates the M(v) function that determines the jump heights in the "continuous" part of an increasing additive process.

Usage

Mv(u, alpha, beta, gama, low, upp, N)

Details

For internal use.

Examples

## The function is currently defined as
function(u = 0.5, alpha = 1, beta = 1, gama = 1 / 2, low = 1e-04,
         upp = 10, N = 5001) {
  x <- -log(seq(from = exp(-low), to = exp(-upp), length = N))
  f <- alpha / gamma(1 - gama) * x^(-(1 + gama)) * exp(-(u +
    beta) * x)
  dx <- diff(x)
  h <- (f[-1] + f[-N]) / 2
  Mv <- rep(0, N)
  for (i in seq(N - 1, 1)) Mv[i] <- Mv[i + 1] + dx[i] * h[i]
  return(list(v = x, Mv = Mv))
}

konkam/BNPdensity documentation built on March 14, 2024, 7:15 a.m.