calcIterativeDifference: Evaluate Bernstein Function differences

In hsloot/rmo: A package for the Marshall-Olkin distribution

Evaluate Bernstein Function differences

Description

Calculate the values for a Bernstein function and its higher-order, alternating iterated forward differences, possibly scaled by a binomial coefficient, i.e.,

\binom{n}{k} {(-1)}^{j-1} \Delta^{j}{ \psi(c x) } , \quad x > 0 .

The evaluation of Bernstein functions using this formula is usually not numerically stable. Consequently, the various alternative approaches are used dependent on the class of the Bernstein function.

Usage

calcIterativeDifference(
  object,
  x,
  difference_order = 0L,
  n = 1L,
  k = 0L,
  cscale = 1,
  ...
)

Arguments

object	An object deriving from the class BernsteinFunction.
x	A nonnegative numeric vector at which the iterated difference of the Bernstein function is evaluated.
difference_order	A nonnegative integer with the order of the alternating iterated forward differences taken on the Bernstein function.
n, k	Nonnegative numbers for the binomial factor.
cscale	A positive number for the composite scaling factor.
...	Pass-through parameter.

See Also

Other Bernstein function generics: calcExShockArrivalIntensities(), calcExShockSizeArrivalIntensities(), calcMDCMGeneratorMatrix(), calcShockArrivalIntensities(), calcValue(), getDefaultMethodString(), getLevyDensity(), getStieltjesDensity()

Examples

bf <- AlphaStableBernsteinFunction(alpha = 0.7)
calcIterativeDifference(bf, 1:5)

hsloot/rmo documentation built on May 1, 2024, 4:28 a.m.