Journal of Applied Mathematics and Stochastic Analysis

Volume 2009 (2009), Article ID 215817, 16 pages

http://dx.doi.org/10.1155/2009/215817

## On Modelling Long Term Stock Returns with Ergodic Diffusion Processes: Arbitrage and Arbitrage-Free Specifications

School of Actuarial Studies, Australian School of Business, University of New South Wales, NSW 2052, Australia

Received 13 March 2009; Accepted 29 July 2009

Academic Editor: Vo Anh

Copyright © 2009 Bernard Wong. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We investigate the arbitrage-free property of stock price models where the local martingale component is based on an ergodic diffusion with a specified stationary distribution. These models are particularly useful for long horizon asset-liability management as they allow the modelling of long term stock returns with heavy tail ergodic diffusions, with tractable, time homogeneous dynamics, and which moreover admit a complete financial market, leading to unique pricing and hedging strategies. Unfortunately the standard specifications of these models in literature admit arbitrage opportunities. We investigate in detail the features of the existing model specifications which create these arbitrage opportunities and consequently construct a modification that is arbitrage free.

#### 1. Introduction

Ever since the fundamental work of Black and Scholes, there has been extensive work in the literature on alternative stock price models. In a continuous time setting, these include, for example, jump diffusions, Levy processes, stochastic volatility models, and regime switching models. Extensive references can be found in Cont and Tankov [1] and Fouque et al. [2]. These models and tools have proven invaluable for long term asset liability management, in particular with applications to insurance and pensions. These include the modelling and pricing of long term embedded guarantees (see, e.g., Sherris [3], Bauer et al. [4], Milevsky and Salisbury [5], Zaglauer and Bauer [6]) and also with the study of optimal asset allocation problems (see, e.g., Cairns [7], Gerber and Shiu [8], Stamos [9]). Extensive references to the vast literature can also be found in the monographs of Hardy [10], Schmidli [11], Milevsky [12], and Møller and Steffensen [13].

In this paper we investigate the arbitrage-free property of the class of Brownian based stock price models where the local martingale component of the (log) stock returns is assumed to be an ergodic diffusion. This class of models was first investigated by Bibby and Sørensen [14] and Rydberg [15, 16] (henceforth “BSR”) who reported good fit of their models to financial data. The use of ergodic diffusions imply that the marginal distributions over a long horizon will be approximately equal to that of the specified stationary distribution, such as the Student- (Bibby and Sørensen [14] provide additional discussion of this property). The dynamics of these models are time homogeneous and, as it is based on diffusions, tractable. The financial market under these models will be complete, and hence the valuation of options and guarantees can be performed without requiring extra assumptions regarding the market price of risk. In contrast, most alternative stock price models admit incomplete markets, with no unique pricing of options and guarantees available in general.

A significant drawback of the ergodic diffusion approach however was also noted by BSR (in particular, Bibby and Sørensen [14] and Rydberg [15]) who showed that no standard equivalent local martingale measure can exist for the ergodic diffusion model with a (generalized) hyperbolic ergodic distribution model they considered, and so the model they considered is only arbitrage free up to a stopping time. By the fundamental theorem of asset pricing (Harrison and Kreps [17], Harrison and Pliska [18], Delbaen and Schachermayer [19]) it follows that these models permit arbitrage. Furthermore the discussions of this issue in BSR for the (generalized) hyperbolic class of models further suggest that this feature may perhaps be present in many ergodic diffusion-based models in general.

Following on from the previous discussion in literature, in this paper we provide a detailed proof that any ergodic diffusion process used as a stock return model, and as specified in literature, will admit arbitrage. We further analyze in detail the cause for these arbitrage opportunities and consequently propose a modification that is arbitrage-free. This modification once again opens up the application of ergodic diffusion models to problems in insurance and finance.

The outline of this paper is as follows. In Section 2 we briefly review the construction of local martingales based on ergodic diffusion processes. Section 3 sets out the financial market we consider and defines the economic notions of portfolios and arbitrage. The standard model specification considered in the existing literature is investigated, and associated arbitrage opportunities identified, in Section 4. This analysis is consequently used to construct an alternative, arbitrage-free specification in Section 5. An extension of the arbitrage-free specification to stock markets with a stochastic term structure of interest rates can be found in Section 6. Section 7 concludes.

#### 2. Local Martingales Based on Ergodic Diffusions

Let be a Brownian motion. Consider a local martingale of the form where is a continuous, strictly positive deterministic function. By Engelbert and Schmidt [20] it follows that a nonexploding solution to the stochastic differential equation (2.1) will always exist.

Consider an interval with . Let be a continuous, bounded probability density function which is strictly positive on and zero outside . Typically in applications to stock returns we are interested in the case when is , and we will make this assumption for the analysis in Sections 4–7.

Using standard diffusion theory (cf. Karlin and Taylor [21]), it was noted by Bibby and Sørensen [14] and Rydberg [15] that can be constructed to be an ergodic diffusion with stationary density by selecting for some arbitrary constant . In particular they considered from the (generalized) hyperbolic class, which includes as special cases the Student- and Normal Inverse Gaussian distributions. As an example, a local martingale with coefficient will possess a Student- (with degrees of freedom) stationary distribution, while a local martingale with will possess a Cauchy stationary distribution. Other distributions can also be used, dependent on the characteristics of the process being modelled.

It is worth noting that the existence of a stationary distribution for (2.1) is closely related to the concept of volatility-induced stationarity in the interest rate modelling literature (Conley et al. [22], Nicolau [23]).

Note that the above is not the only method of constructing a Brownian-based model with a specified distribution. There are two alternative approaches. The first alternative (Bibby et al. [24], Borkovec and Klüppelberg [25]) constructs semimartingales with stationary density by considering diffusion processes with nonzero drift. We do not pursue this approach as our stock price construction considers local martingales, which necessarily have zero drift. The second alternative (Dupire [26], Madan and Yor [27]) uses the Fokker-Plank equation of a diffusion, with an additional time scaling assumption, to construct a diffusion with a specified marginal distribution. A significant drawback of this second alternative approach however is that the dynamics of the resulting local martingales will be time inhomogeneous in general. In contrast, the dynamics of the local martingale (2.1) is time homogeneous.

#### 3. Financial Market

For the financial market we consider a probability space and the time interval , the filtration being generated by -dimensional Brownian Motion, augmented to satisfy the usual conditions. In cases where we consider two nonequivalent measures we will augment with respect to the null sets of both measures. Denote by the equivalence class of all progressively measurable processes satisfying For clarity of presentation, in Sections 3–5 we will assume that the only source of uncertainty in the financial market arises from the stock. Specifically, assume that there are two primary securities traded in the time interval . The first is a savings account with for some constant short rate . The second is a strictly positive stock price , with . We consider alternative specifications for in Sections 4 and 5.

An extension of our framework to include stochastic interest rates can be found in the Section 6.

Portfolios are formed by holding an amount of of the savings account and of the stock, with being progressively measurable. The value of a self-financing portfolio should satisfy for all , To ensure that there are no pathological doubling strategies we also require that the discounted value of a portfolio is bounded below by . Portfolios satisfying this constraint are referred to as tame (Dybvig and Huang [28]). Finally, an arbitrage opportunity is defined as a tame portfolio with value process satisfying

#### 4. Arbitrage Opportunities under the Standard Specification

Ergodic diffusion-based stock price models were first considered by Bibby and Sørensen [14] and Rydberg [15, 16]. Their models for the stock price process are of the form with , constant , and being a driftless ergodic diffusion with stationary density and associated state space . In the following we will call ergodic diffusion based models of the form (4.1) as being of the standard specification.

As is nonexplosive by construction, it follows that The model can equivalently be represented by the following stochastic differential equation by an application of Ito's formula: Representation (4.3) is sometimes more convenient for calculation purposes.

Bibby and Sørensen [14] and Rydberg [15, 16] report a good fit of the generalized hyperbolic ergodic diffusion to financial data. Unfortunately they also noted that no standard equivalent martingale measure can exist for the models they considered. By the fundamental theorem of asset pricing (Harrison and Kreps [17], Harrison and Pliska [18], Delbaen and Schachermayer [19]) it follows that arbitrage opportunities exist in these models. Furthermore, their discussion and proof suggests that the same problem may also apply to any ergodic diffusion-based model of the form (4.1). In the following we provide a detailed proof of the arbitrage opportunity for general ergodic diffusion models and perform further analysis on the technical features of this model specification which created these opportunities. Consequently in Section 5 we construct an alternative model specification that is arbitrage free.

Theorem 4.1. *The financial market with stock price modelled by (4.1) admits arbitrage opportunities.*

*Proof. *By Ito's lemma the process satisfies the stochastic differential equation
In a Brownian setting, the fundamental theorem of asset pricing states that there is no arbitrage if and only if there exists an equivalent measure such that is a local martingale under the equivalent measure. As there is only one source of uncertainty, the market is complete. It follows that an equivalent local martingale measure exists if and only if there exists a strictly positive martingale with , and
where is the market price of risk process, with
Notice that as is a supermartingale as it is a local martingale that is bounded below. Consequently will be absorbing if reached. As the processes and are continuous and non-explosive by construction, , and hence
(cf. Kazamaki [29], Liptser and Shiryaev [30, 6.1.1]).

The martingale property of can be investigated by considering a candidate measure (cf. Kadota and Shepp [31], Delbaen and Shirakawa [32], Wong and Heyde [33], Rogers and Veraart [34] for an application of this technique in different settings), which is not assumed to be equivalent to a priori, with Brownian motion , and a process satisfying
with . Under we assume that will be stopped if explosion (to or ) occurs. Denote this stopping time as . If and are equivalent measures, then and will be equivalent in law under .

Note that can also be represented in terms of a process , with
with also being stopped at the explosion time . Note also that if
then the stochastic integral in (4.10) should be interpreted in the Liptser-Shiryaev [30, 4.2.9] sense.

For calculation purposes it is convenient to consider, under the measure , the related process defined by
which, by the definition of , satisfies
Correspondingly, consider under the measure a process , defined by
By Ito's lemma satisfies
Define the local martingale by
Observe that, by the boundedness of and Novikov's condition, is a strictly positive, true martingale. Hence we can define a measure equivalent to by
for some measurable event , and where represents the indicator function for . Girsanov's theorem shows that under there exists a Brownian motion and where satisfies
Consider the probability of explosion to in finite time. This can be checked with a criterion of Delbaen and Shirakawa [32]. Denote as the cumulative distribution function corresponding to . By the assumptions on in Section 2 we have, for all , with ,
as required by Delbaen and Shirakawa [32, Theorem 1.4]. Consequently the probability of hitting is decided by the convergence or divergence of the integral
with divergence being equivalent to a probability of the process reaching . As it follows that (4.20) is finite, and hence can hit in finite time. As is a local martingale bounded below by , it is also a supermartingale. Hence is absorbing, and we have
and by equivalence
which by comparison to (4.13) implies that and cannot be equivalent measures.

As this is a complete market we can identify a specific arbitrage portfolio and strategy using the techniques of Levental and Skorohod [35]. Firstly note that as is not a true martingale, we have for some constant . Define a -martingale by which is a strictly positive martingale. It follows by Ito's lemma that there exist such that we can represent in stochastic exponential form Consider the process , with Applying Ito's lemma shows that is the value process of a self financing portfolio, with As , we have implying that is an arbitrage opportunity.

#### 5. Arbitrage-Free Specification

In the previous section we have shown that the model specification (4.1) will always admit arbitrage opportunities. In this section we further investigate the technical features of this specification that create the arbitrage and consequently construct a modification that is arbitrage free.

As is not an equivalent measure to , the -local martingale , where , and is not strictly positive, with and defined by (4.10). By Kazamaki [29] and Liptser and Shiryaev [30, 6.1.1] we have where all processes are stopped at the explosion time of , which we will denote as . Note that, by comparing (4.9) and (4.10) with (5.2), and as is bounded, we have by construction.

As is not strictly positive, we have and as is bounded it follows that the term in (5.2) causing (5.4) is of the form . In comparison to the standard model specification (4.1) it is apparent that this term arose from a drift adjustment for the quadratic variation of . This suggests that an alternative stock price model of the form or, equivalently, where is the quadratic variation process of , defined by will be arbitrage free.

Theorem 5.1. *The financial market with stock model (5.6)-(5.7) does not admit arbitrage opportunities.*

*Proof. *Note firstly that we have, by the assumptions on in Section 2,
Hence by Kazamaki [29] and Liptser and Shiryaev [30, 6.1.1] it follows that is strictly positive over .

From Ito's lemma the discounted stock price process satisfies
An equivalent local martingale measure exists if and only if there exists a strictly positive martingale with , and
where is the market price of risk process, with
By construction the process is continuous and bounded. Hence by Novikov's condition it follows that is a strictly positive martingale and an (unique) equivalent local martingale measure exists. By the fundamental theorem of asset pricing it follows that the model is arbitrage free.

Under the unique equivalent local martingale measure the discounted stock price satisfies with

#### 6. Stock Markets with Stochastic Term Structure of Interest Rates

In Sections 3 to 5 we considered a financial market with a constant interest rate for clarity of presentation. In practical applications however the long term nature of many problems (e.g., in insurance and pensions) imply that such an assumption may be inappropriate. In this section we extend Theorem 5.1 to stock markets with a stochastic term structure of interest rates.

To allow for imperfect correlation between stock and interest rates we now consider a probability space and the time interval , the filtration being generated by 2 dimensional Brownian motion , augmented to satisfy the usual conditions. Following Section 5, we consider a stock price model of the form where the local martingale is assumed to be for constants , and stationary density .

Interest rate variability will be introduced via the second Brownian motion and its augmented filtration , with . Following Heath et al. [36] (cf. Musiela and Rutkowski [37, Chapter 11]), assume that for every , the instantaneous forward rate follows for a Borel measurable function , and progressively measurable processes and satisfying Assume , , and are differentiable with respect to , with bounded first derivatives , , and . It is known (cf. Musiela and Rutkowski [37, Proposition 11.1.1]) that the short rate process is a continuous semimartingale, with where Finally, assume that there exists a progressively measurable process such that, for any , we have and such that the process , with is a strictly positive martingale. This assumption is standard in literature and can intuitively be interpreted as assuming that the interest rate market is internally arbitrage free.

Under the above setup the savings account satisfies and the price of a -maturity Zero Coupon Bond process, with initial value , satisfies where

Theorem 6.1. *The financial market with stock model (6.1)-(6.2) and interest rate term structure modelled by (6.3)–(6.8) does not admit arbitrage opportunities.*

*Proof. *The discounted stock price process satisfies
under the measure . As is a strictly positive martingale by assumption, we can define a measure equivalent to by
for some measurable event , and where represents the indicator function for . By Girsanov's theorem, under we have a -dimensional Brownian motion , with
with defined by (6.7). Under the short rate process satisfies
which is continuous and is independent of by the assumptions on the coefficients of (6.3).

Consider now the nonnegative local martingale , with
and where is the market price of risk process corresponding to , with
Notice that by the continuity of . As we also have , and hence by equivalence
Consequently, by the boundedness of , we have
and, by noting that and are independent of , we have (cf. Liptser and Shiryaev, [30, Example 6.2.4])
It follows by Novikov's condition that
and in particular,
Equations (6.19) and (6.22) imply that is a strictly positive martingale, and hence we can define a measure equivalent to by
for some measurable event , and where represents the indicator function for . By Girsanov's Theorem it follows that discounted asset prices are local martingales under . As is equivalent to which is in turn equivalent to , is an equivalent local martingale measure. By the fundamental theorem of asset pricing it follows that the model is arbitrage free.

#### 7. Conclusions

In this paper we investigated the arbitrage-free property of the class of stock price models where the local martingale component is based on an ergodic diffusion with a specified stationary distribution. The dynamics of these models are time homogeneous and, as it is based on Brownian motion, tractable. The financial market under these models will be complete, and hence the valuation of options and guarantees can be performed without requiring extra assumptions regarding the market price of risk. In this paper we provided a detailed proof that any ergodic diffusion process used as a stock return model, and as specified in the existing literature, will admit arbitrage in general. We further analyzed the technical cause for these arbitrage opportunities and consequently constructed a modification that is arbitrage-free. This arbitrage free property is shown to be true in financial markets both with and without stochastic interest rates. Our modification once again opens up the application of ergodic diffusion models to problems in insurance and finance.

#### References

- R. Cont and P. Tankov,
*Financial Modelling with Jump Processes*, Chapman & Hall/CRC Financial Mathematics Series, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2004. View at MathSciNet - J.-P. Fouque, G. Papanicolaou, and K. R. Sircar,
*Derivatives in Financial Markets with Stochastic Volatility*, Cambridge University Press, Cambridge, UK, 2000. View at MathSciNet - M. Sherris, “The valuation of option features in retirement benefits,”
*The Journal of Risk and Insurance*, vol. 62, no. 3, pp. 509–534, 1995. View at Google Scholar - D. Bauer, R. Kiesel, A. Kling, and J. Ruß, “Risk-neutral valuation of participating life insurance contracts,”
*Insurance: Mathematics & Economics*, vol. 39, no. 2, pp. 171–183, 2006. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. A. Milevsky and T. S. Salisbury, “Financial valuation of guaranteed minimum withdrawal benefits,”
*Insurance: Mathematics & Economics*, vol. 38, no. 1, pp. 21–38, 2006. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. Zaglauer and D. Bauer, “Risk-neutral valuation of participating life insurance contracts in a stochastic interest rate environment,”
*Insurance: Mathematics & Economics*, vol. 43, no. 1, pp. 29–40, 2008. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Cairns, “Some notes on the dynamics and optimal control of stochastic pension fund models in continuous time,”
*ASTIN Bulletin*, vol. 30, no. 1, pp. 19–55, 2000. View at Google Scholar - H. U. Gerber and E. S. W. Shiu, “Geometric Brownian motion models for assets and liabilities: from pension funding to optimal dividends,”
*North American Actuarial Journal*, vol. 7, no. 3, pp. 37–56, 2003. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Z. Stamos, “Optimal consumption and portfolio choice for pooled annuity funds,”
*Insurance: Mathematics & Economics*, vol. 43, no. 1, pp. 56–68, 2008. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Hardy,
*Investment Guarantees: Modeling and Risk Management for Equity-Linked Life Insurance*, John Wiley & Sons, River Edge, NJ, USA, 2003. - H. Schmidli,
*Stochastic Control in Insurance*, Probability and Its Applications, Springer, London, UK, 2008. View at MathSciNet - M. A. Milevsky,
*The Calculus of Retirement Income: Financial Models for Pension Annuities and Life Insurance*, Cambridge University Press, Cambridge, UK, 2006. View at MathSciNet - T. Møller and M. Steffensen,
*Market-Valuation Methods in Life and Pension Insurance*, Cambridge University Press, Cambridge, UK, 2007. - B. Bibby and M. Sørensen, “A hyperbolic diffusion model for stock prices,”
*Finance and Stochastics*, vol. 1, pp. 25–41, 1997. View at Google Scholar - T. Rydberg, “A note on the existence of equivalent martingale measures in a Markovian setting,”
*Finance and Stochastics*, vol. 1, pp. 251–257, 1997. View at Google Scholar - T. Rydberg, “Generalized hyperbolic diffusion processes with applications in finance,”
*Mathematical Finance*, vol. 9, no. 2, pp. 183–201, 1999. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. M. Harrison and D. M. Kreps, “Martingales and arbitrage in multiperiod securities markets,”
*Journal of Economic Theory*, vol. 20, no. 3, pp. 381–408, 1979. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. M. Harrison and S. R. Pliska, “Martingales and stochastic integrals in the theory of continuous trading,”
*Stochastic Processes and Their Applications*, vol. 11, no. 3, pp. 215–260, 1981. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - F. Delbaen and W. Schachermayer, “A general version of the fundamental theorem of asset pricing,”
*Mathematische Annalen*, vol. 300, no. 3, pp. 463–520, 1994. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. J. Engelbert and W. Schmidt, “On solutions of one-dimensional stochastic differential equations without drift,”
*Zeitschrift für Wahrscheinlichkeitstheorie und Verwandte Gebiete*, vol. 68, no. 3, pp. 287–314, 1985. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Karlin and H. M. Taylor,
*A Second Course in Stochastic Processes*, Academic Press, New York, NY, USA, 1981. View at MathSciNet - T. Conley, L. Hansen, E. Luttmer, and J. Scheinkman, “Short-term interest rates as subordinated diffusion,”
*Review of Financial Studies*, vol. 10, pp. 525–577, 1997. View at Google Scholar - J. Nicolau, “Processes with volatility-induced stationarity: an application for interest rates,”
*Statistica Neerlandica*, vol. 59, no. 4, pp. 376–396, 2005. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. Bibby, I. M. Skovgaard, and M. Sørensen, “Diffusion-type models with given marginal distribution and autocorrelation function,”
*Bernoulli*, vol. 11, no. 2, pp. 191–220, 2005. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Borkovec and C. Klüppelberg, “Extremal behavior of diffusion models in finance,”
*Extremes*, vol. 1, no. 1, pp. 47–80, 1998. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. Dupire, “Pricing with a smile,”
*Risk*, vol. 7, pp. 18–20, 1994. View at Google Scholar - D. B. Madan and M. Yor, “Making Markov martingales meet marginals: with explicit constructions,”
*Bernoulli*, vol. 8, no. 4, pp. 509–536, 2002. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - P. Dybvig and C. Huang, “Nonnegative wealth, absence of arbitrage, and feasible consumption plans,”
*Review of Financial Studies*, vol. 1, pp. 377–401, 1988. View at Google Scholar - N. Kazamaki, “Continuous exponential martingales and BMO,” in
*Lecture Notes in Mathematics*, vol. 1579, Springer, Berlin, Germany, 1994. View at Google Scholar - R. Liptser and A. Shiryaev,
*Statistics of Random Processes*, vol. 1, Springer, Berlin, Germany, 1977. - T. T. Kadota and L. A. Shepp, “Conditions for absolute continuity between a certain pair of probability measures,”
*Zeitschrift für Wahrscheinlichkeitstheorie und Verwandte Gebiete*, vol. 16, pp. 250–260, 1970. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - F. Delbaen and H. Shirakawa, “No arbitrage condition for positive diffusion price processes,”
*Asia-Pacific Financial Markets*, vol. 9, pp. 159–168, 2002. View at Google Scholar - B. Wong and C. C. Heyde, “On the martingale property of stochastic exponentials,”
*Journal of Applied Probability*, vol. 41, no. 3, pp. 654–664, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. C. G. Rogers and L. A. M. Veraart, “A stochastic volatility alternative to SABR,”
*Journal of Applied Probability*, vol. 45, no. 4, pp. 1071–1085, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Levental and A. V. Skorohod, “A necessary and sufficient condition for absence of arbitrage with tame portfolios,”
*The Annals of Applied Probability*, vol. 5, no. 4, pp. 906–925, 1995. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. Heath, R. Jarrow, and A. Morton, “Bond pricing and the term structure of interest rates: a new methodology for contingent claim valuation,”
*Econometrica*, vol. 60, pp. 77–105, 1992. View at Google Scholar - M. Musiela and M. Rutkowski,
*Martingale Methods in Financial Modelling*, vol. 36 of*Stochastic Modelling and Applied Probability*, Springer, Berlin, Germany, 2nd edition, 2005. View at MathSciNet