

The distribution determines an homogeneous Levy process, and this process is representable through subordination of Brownian motion by the inverse Gaussian process. The normal inverse Gaussian distribution is defined as a variance-mean mixture of a normal distribution with the inverse Gaussian as the mixing distribution. A numerical scheme, based on the collocation method for boundary value problems, is further illustrated, in order to get precise approximations of the free-boundary problem solution, which seems to be very hard to derive analytically, because of the particular structure of the Lévy measure of an inverse Gaussian process.

We show that the initial optimal stopping problem for the posterior probability of one of the hypotheses can be reduced to a free-boundary problem, whose unknown boundary points are characterized by the principles of the continuous or smooth fit and whose unknown value function solves a linear integro-differential equation over the continuation set. This problem arises when we are willing to test the positive drift of an unobservable Brownian motion, for which only the first passage times over positive thresholds can be recorded. We analyze the Bayesian formulation of the sequential testing of two simple hypotheses for the distributional characteristics of an inverse Gaussian process.
