**Journal of Mathematical Finance**

Vol.08 No.01(2018), Article ID:81873,13 pages

10.4236/jmf.2018.81001

Extended Model of Stock Price Behaviour

Nico Koning^{1*}, Daniel T. Cassidy^{2}, Rachid Ouyed^{1}^{ }

^{1}Department of Physics and Astronomy, University of Calgary, Calgary, Canada

^{2}Department of Engineering Physics, McMaster University, Hamilton, Canada

Copyright © 2018 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: October 30, 2017; Accepted: January 16, 2018; Published: January 19, 2018

ABSTRACT

We have developed an extended model for stock price behaviour that is able to accommodate fat-tailed distributions with support as large as $\left[-\infty \mathrm{,}\infty \right]$ . The “homogeneously saturated” (HS) model avoids exponential price changes for large fluctuations by means of a saturation parameter. In the limit where the saturation parameter is zero, the standard model of stock price behaviour (i.e., geometric Brownian motion) is recovered. We compare simulated stock price series generated for both the standard and HS model for the DJIA and five random stocks from the NYSE and NASDAQ exchanges. We find that in all cases, the HS model provides a better fit to the observed price series than the standard model. This has implications to many areas of finance including the Black-Scholes formula for option pricing.

**Keywords:**

Stock Price, Fat-Tail, Brownian Motion, Student’s T, Black-Scholes

1. Introduction

The standard model of stock price behaviour generates prices through geometric Brownian motion with a deterministic drift rate [1] [2] [3] [4] . This model is used throughout the financial world; most notably in deriving the Black-Scholes formula of option pricing. The standard model for generating stock prices is

$\frac{\text{d}S\left(t\right)}{\text{d}t}=\alpha S\left(t\right)+\sigma S\left(t\right)f\left(t\right)$ (1)

where $S\left(t\right)$ is the price of the stock at time t, $\alpha $ is the drift rate, $\sigma $ is the volatility of the stock price and $f\left(t\right)$ is a zero mean, normally distributed, stochastic, uncorrelated in time, noise driving term. Using $\eta \left(t\right)=\alpha +\sigma f\left(t\right)$ , we can simplify Equation (1) to

$\frac{\text{d}S\left(t\right)}{\text{d}t}=\eta \left(t\right)S\left(t\right)\mathrm{.}$ (2)

Integrating Equation (2) gives the predicted price of the stock at a later time t:

${S}_{s}\left(t\right)={S}_{0}{\text{e}}^{\omega \left(t\right)}\mathrm{,}$ (3)

where the subscript s denotes the standard model, ${S}_{0}$ is the price of the stock at $t=0$ , and for clarity, we have made the substitution

$\omega \left(t\right)={\displaystyle {\int}_{0}^{t}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\eta \left({t}^{\prime}\right)\text{d}{t}^{\prime}.$ (4)

Equations (3) and (4) suggest that the price of the stock depends exponentially on the integral of the noise driving term, $f\left(t\right)$ . Since this term is assumed to be normally distributed, the probability of a large price is essentially zero and the predicted price remains bounded. If, however, the underlying noise is not normally distributed, Equation (3) might predict wild price swings that are unrealistic and not observed on the market.

Stock returns are generally assumed to follow a normal distribution, in part owing to mathematical simplicity. It has been known for some time, however, that this assumption is not supported by actual stock prices (e.g., [5] [6] [7] ). For example, daily returns of the DJIA and the S&P 500 indices are described by a fat-tailed distribution [3] [8] . Prices predicted by the standard model, using a normal distribution as the noise driving term, will therefore be inaccurate and simply substituting a fat-tailed distribution will permit the infinite prices mentioned above.

A simulation of stock prices using Equation (2) and drawing $\eta \left(t\right)$ from different distributions is given in Figure 1. The fat-tailed distribution used in

Figure 1. Simulated stock price over 1000 days using the standard model, Equation (2), with ${S}_{0}=1000$ , $t=1$ day, and $\eta \left(1\right)$ drawn from different distributions but for the same sequence of pseudo random numbers. The black line shows daily increments from a normal distribution with $\sigma =0.01,\mu =0.001$ ; the blue line shows draws from a Student’s T with $\nu =2.5$ ; the green line, a Student’s T with $\nu =2.0$ ; and the red line, a Student’s T with $\nu =1.5$ . Each Student’s T distribution has $\beta =0.01$ and $\mu =0.001$ . It is evident that when the noise used in the standard model is a fat-tailed distribution (e.g., Student’s T), the predicted price is subject to large fluctuations. The extent of these fluctuations increases with fatter tails (i.e., smaller $\nu $ ).

this comparison is the Student’s T distribution (see Appendix A), where the “fatness’’ of the tails is governed by the shape parameter $\nu $ . It is clear from Figure 1 that when a fat-tailed distribution is used in the standard model, large price fluctuations can result.

There exist several approaches to pricing stocks when the noise driving term is a fat-tailed distribution. One approach is to modify the tails of the distribution such that the contributions far into the tails are negligible while not affecting significantly the central portion of the distribution, which fits well the observed data (e.g. [9] [10] [11] [12] [13] ). Capping the value of the stock [8] [14] , or truncating the distribution [13] are alternative possibilities.

A different approach to pricing stocks when the underlying distribution is fat-tailed is to allow for saturation of the stock price by depletion of the resource that supports the price (i.e., by depletion of the reservoir of money that is available to purchase the stock). This is the approach we investigate in this paper. This “homogeneously saturated” (HS) model for the price of a stock is constructed and compared to the standard model to gain insight into the pricing of financial assets when the underlying distribution is fat-tailed.

2. Homogeneously Saturated (HS) Model

The aim of this section is to develop a pricing model which can handle fat-tailed distributions as the noise driving term. Let $M\left(t\right)$ be the amount of money available in a reservoir to buy the stock, N the rate at which money is added to the reservoir, $\gamma S\left(t\right)M\left(t\right)$ the rate at which money is removed from the reservoir due to the purchase of the stock, and $\zeta M\left(t\right)$ the rate at which money is removed from the reservoir due to the purchase of other goods. We then have

$\frac{\text{d}}{\text{d}t}M\left(t\right)=N-\gamma S\left(t\right)M\left(t\right)-\zeta M\left(t\right)+\delta f\left(t\right)$ (5)

where $\delta f\left(t\right)$ is a noise driving term. It is interesting to note that the noise in the homogeneously saturated model is ascribed to fluctuations in the amount of money available to invest in the stock. All parameters ( $\zeta $ , $N$ , $\gamma $ ) in Equation (5) have a time dependence, but it is assumed that these parameters vary slowly enough that they can be treated as constants. As well, we assume the system is in

steady state, such that $\frac{\text{d}}{\text{d}t}M\left(t\right)=0$ , leading to

$M\left(t\right)=\frac{\alpha +\rho f\left(t\right)}{1+\lambda S\left(t\right)}$ (6)

where we have made the substitutions: $\alpha =N/\zeta $ , $\rho =\delta /\zeta $ and $\lambda =\gamma /\zeta $ .

The next step in our derivation is inspired from laser physics, where coupled rate equations are used to describe the interaction between the laser output (analogous to $S\left(t\right)$ in our case) and the inversion (analogous to $M\left(t\right)$ in our case) [15] . The saturation of the inversion in a laser enforces conservation of energy, thereby forcing the output to be finite and to track the input; a trait we desire in our model. Therefore coupling $S\left(t\right)$ to $M\left(t\right)$ will allow for saturation of the stock price such that it cannot continually increase exponentially as allowed by the standard model. We therefore make the following assumption:

$\frac{\text{d}}{\text{d}t}S\left(t\right)=M\left(t\right)S\left(t\right).$ (7)

The validity of Equation (7) will be determined by how well it fits the available data. Using Equation (6) in (7) gives

$\frac{\text{d}}{\text{d}t}S\left(t\right)=\frac{\alpha S\left(t\right)+\rho f\left(t\right)S\left(t\right)}{1+\lambda S\left(t\right)}.$ (8)

Making the substitution $\eta \left(t\right)=\alpha +\rho f\left(t\right)$ we have

$\frac{\text{d}}{\text{d}t}S\left(t\right)=\frac{\eta \left(t\right)S\left(t\right)}{1+\lambda S\left(t\right)}.$ (9)

The price of the stock in our HS model is therefore

$S\left(t\right)=\frac{{S}_{0}{\text{e}}^{\omega \left(t\right)}}{{\text{e}}^{\lambda \left(S\left(t\right)-{S}_{0}\right)}}=\frac{{S}_{s}\left(t\right)}{{\text{e}}^{\lambda \left(S\left(t\right)-{S}_{0}\right)}}.$ (10)

The HS model is similar in form to the standard model, Equation (3), except for the ${\text{e}}^{\lambda \left(S\left(t\right)-{S}_{0}\right)}$ term in the denominator of the former. Indeed, when $\lambda =0$ the HS model reduces to the standard model. The $S\left(t\right)$ dependency in the denominator of Equation (10) effectively saturates the price of the stock; without it wealth would not be conserved and $S\left(t\right)$ could continually increase exponentially with time.

To demonstrate the behaviour of the saturation parameter, we generate simulated prices over a 1000 day period using the HS model, Equation (9), by adding different and independent one day solutions for
$S\left(1\right)$
for each t to obtain the price for the next day, i.e.,
$S\left(t+1\right)=S\left(t\right)+{S}_{t+1}\left(1\right)$
, for
$t=0,\cdots ,1000$
with
$S\left(0\right)=1000$
and
${S}_{t+1}\left(1\right)$
being the one day solution over the time interval t to
$t+1$
. These generated price series are shown in Figure 2 for various values of
$\lambda $
and for the same sequence of draws from a pseudo random number generator (PRNG) to determine
${\eta}_{t}\left(1\right)$
and hence
${S}_{t}\left(1\right)$
for each t.
${\eta}_{t}\left(1\right)$
is
$\eta \left(1\right)$
created using the t^{t}^{h} value from the PRNG sequence for a Student’s T distribution with shape parameter
$\nu =2.5$
, scale parameter
$\beta =0.01$
, and drift
$\mu =0.001$
. When
$\lambda $
is zero, the price of the stock increases exponentially with the noise driving term and leads to large price swings. The case
$\lambda =0$
is equivalent to the standard model with a fat-tailed Student’s T noise source. When
$\lambda >0$
, the price of the stock approximately linearly follows the noise; i.e., it becomes saturated. Variations in the price of the stock are increasingly damped with higher values of
$\lambda $
. See Appendix B for a justification of the approach
$S\left(t+1\right)=S\left(t\right)+{S}_{t+1}\left(1\right)$
.

To reiterate, $\lambda $ describes the rate at which money is removed from the money reservoir owing to purchase of the stock. $\lambda =0$ (as in the standard model) means that buying the stock has no effect on the money supply; essentially the standard model assumes an infinite reservoir of money is available to buy the stock. From Equation (9), it can be observed that the reciprocal of the saturation parameter, ${\lambda}^{-1}$ , has the same units as the stock price $S\left(t\right)$ . One can identify ${\lambda}^{-1}$ as the stock price at saturation. When $S\left(t\right)={\lambda}^{-1}$ , the instantaneous rate of change with time of $S\left(t\right)$ is one-half of what it would be with $\lambda =0$ ; see Equation (9). For $S\left(t\right)\ll {\lambda}^{-1}$ , the reservoir of money is not depleted (or saturated) by the rate of transactions, and the rate of change of $S\left(t\right)$ is similar in magnitude to the standard model. For $S\left(t\right)\gg {\lambda}^{-1}$ , the reservoir of money that is available to purchase the stock is saturated (or depleted) by the rate of transactions and the time rate of change of $S\left(t\right)$ is greatly diminished.

3. Results

The goal in developing the HS model is to provide an extension to the standard model that can accommodate assets whose returns are fat-tail distributed. In this section we compare the HS model to the standard model using real data in an effort to corroborate our claim.

Our metric for how close simulated prices match the observed will be the mean absolute percentage error (MAPE) ( [16] , e.g.), defined as

$M\left(A,O\right)=100\times \frac{1}{n}{\displaystyle \sum _{t=1}^{n}}\left|\frac{{A}_{t}-{O}_{t}}{{A}_{t}}\right|$ (11)

where n is the number of simulated days, ${A}_{t}$ is the observed price on day t, and ${O}_{t}$ is the simulated price on day t. Of course since we are drawing randomly from distributions, any given simulated price series will yield a different M. We therefore create a set of k = 100,000 simulated price series for each model and compare the average M from those trials. That is, we compare

${M}_{\text{avg},D}\left(\lambda \right)=\frac{1}{k}{\displaystyle \sum _{i=1}^{k}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}M\left({\text{SIM}}_{D}\left(\lambda \right),\text{OBS}\right)$ (12)

where D is the distribution (either N for normal or ST for Student’s T), ${\text{SIM}}_{D}{\left(\lambda \right)}_{i}$ is a simulated price series obtained using Equation (9) (for a given $\lambda $ and distribution), and OBS is the observed price series.

As an example, we use closing values of the Dow Jones Industrial Average (DJIA) for 1000 days starting on January 1 2010 as our observed price series. The returns along with the Normal and Student’s T best fit are shown in the first panel of Figure 3 and the best fit parameters to the 1000 one-day returns are given in the first row of Table 1. It is clear from both Figure 3 and the ${\chi}^{2}$ in the table that the Student’s T is the superior fit to the DJIA one-day returns.

We begin by determining the saturation parameter, ${\lambda}_{0}$ , that gives the best fit between the Student’s T simulated price series and the observed DJIA price series. We do this by using Equation (12) as a function of $\lambda $ in a standard optimization technique (i.e., Brent’s Method [17] ). We find that ${\lambda}_{0}=2.04\times {10}^{-4}$ , with ${M}_{\text{avg}}=7.5\%$ for the DJIA data. We then proceed to determine

Figure 2. Simulated stock prices $S\left(t\right)$ over a 1000 day period using the HS model for $t=1$ day, ${S}_{0}=1000$ , and for various values of $\lambda $ . The same sequence of draws from a Student’s T distribution with $\nu =2.5$ , $\beta =0.01$ , and $\mu =0.001$ were used to create the $S\left(t\right)$ for the different values of $\lambda $ . As $\lambda $ increases, the large price swings seen with the standard model ( $\lambda =0$ ) are suppressed.

Figure 3. Normal and Student’s T fits to daily returns from six different assets. The best fit parameters for each distribution are given in Table 1. In all cases the Student’s T provides a superior fit.

• ${M}_{\text{avg}\mathrm{,}N}\left(0\right)$ (i.e., standard model, Equation (9) with $\lambda =0$ and $\eta $ drawn from the normal distribution with parameters given in Table 1).

• ${M}_{\text{avg}\mathrm{,}ST}\left(0\right)$ (i.e., Student’s T distribution directly inserted in the standard

Table 1. Best fit parameters for Normal and Student’s T distributions to one-day returns with ±uncertainties (confidence level = 0.95) in parentheses. Each stock data set starts on January 1 2010 and runs for 1000 days. The reduced c^{2} for each fit is also given. b is the scale parameter for the Student’s T distribution; see Appendix A.

model).

• ${M}_{\text{avg}\mathrm{,}ST}\left({\lambda}_{0}\right)$ (i.e., HS model, Equation (9) with $\lambda ={\lambda}_{0}$ and $\eta $ drawn from the Student’s T distribution with parameters given in Table 1).

The results for the DJIA are given in the first row of Table 2. There are several things to notice. First, the standard model performs relatively poorly for the DJIA with ${M}_{\text{avg}}=33\%$ . This can be expected since the daily returns of the DJIA are fit much better by a fat-tailed distribution. Second, when a fat-tailed distribution is substituted into the standard model directly, ${M}_{\text{avg}}$ is actually larger than with a normal distribution despite the better fit to the one-day returns. This, as we mentioned earlier, is expected due to the large probability of large price swings when directly using a fat-tailed distribution in the standard model. These two points taken together justify the need for a better pricing model that can handle fat-tailed distributions. Third, the HS model with $\lambda =2.04\times {10}^{-4}$ , clearly gives superior results to the standard model with an ${M}_{\text{avg}}$ of only 7.5%.

To test these ideas further, we have randomly selected five stocks from the NYSE and NASDAQ exchanges; TRP, PAGG, VBF, PNW and HBCP. The results from these additional stocks are given in Figure 3, Table 1, and Table 2. Specifically, Figure 3 displays best fits of normal and Student’s T distributions to the daily returns for the five random stock selections and for the DJIA. The parameters determined for the best fit distributions are used in the simulations with the standard and HS models for the stock price.

In 4 out of the 6 cases, the standard model performs relatively poorly. The exceptions are TRP with ${M}_{\text{avg}}=19.3\%$ and VBF with ${M}_{\text{avg}}=13.5\%$ . For TRP this is expected since its returns are actually well fit by a normal distribution. VBF, on the other hand, is poorly fit by a normal distribution so its low ${M}_{\text{avg}}$ is surprising.

In all cases, except marginally PNW, simply substituting a fat-tailed distribution into the standard model results in a higher ${M}_{\text{avg}}$ than the standard model. This, as we mentioned previously, is due to the large price swings inherent

Table 2. Comparison between simulated and observed price series for different models. The average MAPE (Equation (12)) generated from 100,000 simulated price series and the observed price series (see Figure 5) for various assets are given. ${M}_{\text{avg}\mathrm{,}N}\left(0\right)$ represents the standard model (i.e., Equation (9) with $\lambda =0$ and $\eta $ drawn from the normal distribution), ${M}_{\text{avg}\mathrm{,}ST}\left(0\right)$ Equation (9) with $\lambda =0$ and $\eta $ drawn from the Student’s T distribution, ${M}_{\text{avg}\mathrm{,}ST}\left({\lambda}_{0}\right)$ the HS model with $\lambda ={\lambda}_{0}$ and $\eta $ drawn from the Student’s T distribution. The best fit parameters used for each distribution are given in Table 1. The minimum and maximum M are given in parenthesis.

in fat-tailed distributions, and is the motivation for this paper.

Finally in all 6 cases, the HS model performs better than the standard model at predicting future stock prices. The most surprising result here is that for VBF the HS model ( ${M}_{\text{avg}}=13.2\%$ ) outperforms the standard model ( ${M}_{\text{avg}}=13.5\%$ ) by only a slight margin. This is unexpected since VBF returns are fit much better by a fat ( $\nu =1.69$ ) Student’s T distribution. A clue may lie in the exceptionally high maximum M score for VBF (~600%) indicating that large price fluctuations are still present. This, however, is not nearly as large as the maximum M score of ~500,000% for a Student’s T distribution in the standard model without saturation. This may signal a limitation of the HS model in that although the addition of the saturation parameter, $\lambda $ , reduces the effects of large price swings for fat-tailed distributions, it does not eliminate them altogether. This limitation is most pronounced for exceptionally fat-tailed distributions where the probability of large returns is likely.

Figure 4 shows the distribution of 100,000 M values for each stock, with a bin width of $M=1\%$ for the standard and HS models. In all cases except HBCP, the HS model distribution is thinner and further left than the standard model. This implies that the HS model is more likely to match the observed stock price. VBF and HBCP have nearly overlapping distributions which is affirmed by their similar ${M}_{\text{avg}}$ scores between the standard and HS models. Figure 5 shows example price series for each stock taken from those that make up the peak bin of the distribution in Figure 4. This gives an indication of how well each model may be expected to perform.

Figure 4. Distribution of MAPE scores using both the standard and HS model for six different assets. In all cases, except HBCP, the HS model has a narrower distribution with a lower peak MAPE score.

Figure 5. Examples of price series using both the standard and HS model. The price series for each model was chosen randomly from those that make up the peak bin (MAPE score given in parentheses) of their respective distributions shown in Figure 4. The observed price series for each asset is shown in black.

4. Conclusions

The standard model, Equation (1), is the most popular model for stock price behaviour. It has wide reaching influence, most notably in the development of the Black-Scholes formula for option pricing ( [1] , e.g.). A primary assumption of the standard model is that daily stock returns follow a normal distribution; an assumption unfounded in the actual market. Fat-tailed distributions fit the returns much better, but when the support for the fat-tailed distributions is assumed to be $\left[-\infty \mathrm{,}+\infty \right]$ and these fat-tailed distributions are used in the standard model, the standard model allows for extreme (infinite) price swings that are not observed. In this paper we have developed a model for stock price behaviour that adequately allows for fat-tailed distributions with support of $\left[-\infty \mathrm{,}+\infty \right]$ . Using data from the DJIA and five random stocks from the NYSE and NASDAQ exchanges, we have shown/reaffirmed:

• Daily stock returns are indeed better fitted by a fat-tailed distribution (Student’s T in this case) than by a normal (Gaussian) distribution.

• The standard model, in most cases, performs poorly when predicting future stock prices.

• Fat-tailed distributions, when substituted directly into the standard model ( $\lambda =0$ ), perform poorly when predicting future stock prices.

• A non-zero saturation parameter, $\lambda $ , when added to the standard model is able to effectively suppress large price swings inherent in fat-tailed distributions for daily returns. $\lambda $ represents the fraction of money removed from the reservoir due to the purchase of the stock. In the standard model $\lambda =0$ which implies an infinite amount of money available for the purchase of the stock.

• The HS model ( $\lambda \ne 0$ ) developed in this paper is able to incorporate fat-tailed distributed returns, and consistently outperforms the standard model when predicting future stock prices.

We have shown that a model for stock price behaviour using fat-tailed distributions with support as large as $\left[-\infty \mathrm{,}+\infty \right]$ for daily returns is attainable. This model is realistic in that it allows for a finite supply of money to be saturated by the net rate of transactions and does not require truncation or capping [8] [9] [13] [14] , or modification of the distribution of returns [10] [11] [12] to avoid very large and unobserved stock prices. This homogeneously saturated model, which borrows from laser physics [15] , will have widely spread repercussions, most importantly for the Black-Scholes formula for option pricing, which will be investigated in future work. The homogeneously saturated (HS) model requires that a value for the saturation parameter $\lambda $ , which is defined in Equation (9), must be estimated from historical data or assumed. The HS model ascribes fluctuations in the price of a stock to fluctuations in the amount of money available to invest in the stock.

Acknowledgements

This work was funded in part by the Natural Sciences and Engineering Research Council of Canada.

Cite this paper

Koning, N., Cassidy, D.T. and Ouyed, R. (2018) Extended Model of Stock Price Behaviour. Journal of Mathematical Finance, 8, 1-13. https://doi.org/10.4236/jmf.2018.81001

References

- 1. Hull, J.C. (2006) Options, Futures, and Other Derivatives. 6th Edition, Pearson Education Inc., New Jersey, 269-272.
- 2. McDonald, R.L. (2006) Derivatives Markets. 2nd Edition, Pearson Education Inc., Boston, Ch. 20.
- 3. Bouchard, J.-P. and Potters, M. (2003) Theory of Financial Risk and Derivative Pricing. 2nd Edition, Cambridge University Press, Cambridge, Ch. 13. https://doi.org/10.1017/CBO9780511753893
- 4. Lax, M., Cai, W. and Xu, M. (2006) Random Processes in Physics and Finance. Oxford University Press, New York, Ch. 16.4. https://doi.org/10.1093/acprof:oso/9780198567769.001.0001
- 5. Mandelbrot, B. (1963) The Variation of Certain Speculative Prices. Journal of Business, 36, 394-419. https://doi.org/10.1086/294632
- 6. Fama, E. (1965) The Behavior of Stock Market Prices. Journal of Business, 38, 34-105. https://doi.org/10.1086/294743
- 7. Clark, P. (1973) A Subordinated Stochastic Process Model with Finite Variance for Speculative Prices. Econometrica, 41, 135-155. https://doi.org/10.2307/1913889
- 8. Cassidy, D.T., Hamp, M.J. and Ouyed, R. (2010) Pricing European Options with a log Student’s t-Distribution: A Gosset Formula. Physica A, 389, 5736-5748. https://doi.org/10.1016/j.physa.2010.08.037
- 9. Cassidy, D.T. (2011) Describing n-Day Returns with Student’s t-Distributions. Physica A, 390, 2794-2802. https://doi.org/10.1016/j.physa.2011.03.019
- 10. Moriconi, L. (2007) Delta Hedged Option Valuation with Underlying Non-Gaussian Returns. Physica A, 380, 343-350. https://doi.org/10.1016/j.physa.2007.01.018
- 11. Lim, G.C., Martin, G.M. and Martin, V.L. (2006) Pricing Currency Options in the Presence of Time-Varying Volatility and Non-Normalities. Journal of Multinational Financial Management, 16, 291-314. https://doi.org/10.1016/j.mulfin.2005.08.004
- 12. Lye, J.N. and Martin, V.L. (1993) Robust Estimation, Non-Normalities and Generalized Exponential Distributions. Journal of the American Statistical Association, 88, 261-267.
- 13. Cassidy, D.T. (2012) Effective Truncation of a Student’s t-Distribution by Truncation of the Chi Distribution in a Chi-Normal Mixture. Open Journal of Statistics, 2, 519-525. https://doi.org/10.4236/ojs.2012.25067
- 14. Cassidy, D.T., Hamp, M.J. and Ouyed, R. (2013) Log Student’s t-Distribution Based Option Sensitivities: Greeks for the Gosset Formulae. Quantitative Finance, 13, 1289-1302. https://doi.org/10.1080/14697688.2012.744087
- 15. Cassidy, D.T. (1984) Analytic Description of a Homogeneously Broadened Injection Laser. IEEE Journal of Quantum Electronics, QE-20, 913-918. https://doi.org/10.1109/JQE.1984.1072481
- 16. Hyndman, R.J. and Athanasopoulos, G. (2013) Forecasting: Principles and Practice. OTexts, Melbourne, Ch. 2.5.
- 17. Brent, R.P. (1973) Algorithms for Minimization without Derivatives. Prentice-Hall, Englewood Cliffs, NJ, Ch. 4.

Appendix A: Student’s T Distribution

A Student’s T distribution is a fat-tailed distribution with a probability density function

$f\left(t\right)=\frac{1}{\beta \sqrt{\text{\pi}\nu}}\frac{\Gamma \left(\frac{\nu +1}{2}\right)}{\Gamma \left(\frac{\nu}{2}\right)}{\left(1+\frac{{\left(t-\mu \right)}^{2}}{\nu {\beta}^{2}}\right)}^{-\frac{\nu +1}{2}}$ (A.1)

where $\Gamma $ is the gamma function, $\mu $ is the mean, $\beta $ is the scale parameter, and $\nu $ is the shape parameter. Figure A1 shows the effect of varying $\nu $ (left panel) and $\beta $ (right panel). The standard deviation of Equation (A.1) is

$\sigma =\beta \sqrt{\frac{\nu}{\nu -2}}$ for $\nu >2$ .

Appendix B: Additivity of Solutions

We wish to show that solutions $S(\cdot )$ to the differential Equation (9) have the property $S\left(t+1\right)=S\left(t\right)+{S}_{t+1}\left(1\right)$ .

Add a second argument to $S\left(t\right)$ and $w\left(t\right)$ to indicate the time origin:

$w\left(a\mathrm{,}b\right)={\displaystyle {\int}_{a}^{b}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\eta \left(\xi \right)\text{d}\xi $ (B.1)

and

$S\left(a\mathrm{,}b\right)=\frac{S\left(a\mathrm{,}a\right){\text{e}}^{w\left(a\mathrm{,}b\right)}}{exp\left(\lambda \left(S\left(a\mathrm{,}b\right)-S\left(a\mathrm{,}a\right)\right)\right)}\mathrm{.}$ (B.2)

$S\left(a\mathrm{,}a\right)$ equals $S\left(a\right)$ and is the value of S at a. $S\left(a\mathrm{,}b\right)$ is the value of ${S\left(t\right)|}_{t=b}$ given the initial value at $t=a$ was $S\left(a\right)={S}_{a}$ .

Figure A1. Effect of varying $\nu $ (left panel, $\beta =1$ , $\mu =0$ ) and $\beta $ (right panel, $\nu =1$ , $\mu =0$ ) on the Student’s T distribution (Equation (A.1)). $\nu =\infty $ corresponds to the normal distribution.

Start with $S\left(a\mathrm{,}b\right)+S\left(b\mathrm{,}c\right)$ and use Equation (B.2) to simplify:

$S\left(a\mathrm{,}b\right)+S\left(b\mathrm{,}c\right)$ (B.3)

$=\frac{{S}_{a}{\text{e}}^{w\left(a,b\right)}}{\mathrm{exp}\left(\lambda \left(S\left(a,b\right)-{S}_{a}\right)\right)}+\frac{{S}_{b}{\text{e}}^{w\left(b,c\right)}}{\mathrm{exp}\left(\lambda \left(S\left(b,c\right)-{S}_{b}\right)\right)}$ (B.4)

$=\frac{\mathrm{exp}\left(\lambda \left(S\left(b,c\right)-{S}_{b}\right)\right){S}_{a}{\text{e}}^{w\left(a,b\right)}+\mathrm{exp}\left(\lambda \left(S\left(a,b\right)-{S}_{a}\right)\right){S}_{b}{\text{e}}^{w\left(b,c\right)}}{\mathrm{exp}\left(\lambda \left(S\left(a,b\right)+S\left(b,c\right)-\left({S}_{a}+{S}_{b}\right)\right)\right)}$ (B.5)

$=\frac{\frac{{S}_{b}{\text{e}}^{w\left(b,c\right)}}{S\left(b,c\right)}{S}_{a}{\text{e}}^{w\left(a,b\right)}+\frac{{S}_{a}{\text{e}}^{w\left(a,b\right)}}{S\left(a,b\right)}{S}_{b}{\text{e}}^{w\left(b,c\right)}}{\mathrm{exp}\left(\lambda \left(S\left(a,b\right)+S\left(b,c\right)-\left({S}_{a}+{S}_{b}\right)\right)\right)}$ (B.6)

$=\frac{{S}_{a}{S}_{b}{\text{e}}^{w\left(a,c\right)}\left(\frac{1}{S\left(b,c\right)}+\frac{1}{S\left(a,b\right)}\right)}{\mathrm{exp}\left(\lambda \left(S\left(a,b\right)+S\left(b,c\right)-\left({S}_{a}+{S}_{b}\right)\right)\right)}$ (B.7)

$=\frac{\left(S\left(a,b\right)+S\left(b,c\right)\right){\text{e}}^{w\left(a,c\right)}\left(\frac{{S}_{a}{S}_{b}}{S\left(a,b\right)S\left(b,c\right)}\right)}{\mathrm{exp}\left(\lambda \left(S\left(a,b\right)+S\left(b,c\right)-\left({S}_{a}+{S}_{b}\right)\right)\right)}.$ (B.8)

But ${S}_{b}=S\left(a,b\right)$ , by definition, since ${S}_{b}$ is the value of $S\left(t\right)$ at time $t=b$ , $a<b<c$ , given that the initial value for the time series is ${S}_{a}$ . Then,

$S\left(a,b\right)+S\left(b,c\right)$ (B.9)

$=\frac{\left(S\left(a,b\right)+S\left(b,c\right)\right){\text{e}}^{w\left(a,c\right)}\left(\frac{{S}_{a}}{S\left(b,c\right)}\right)}{\mathrm{exp}\left(\lambda \left(S\left(b,c\right)-{S}_{a}\right)\right)}$ (B.10)

or

$S\left(b,c\right)=\frac{{S}_{a}{\text{e}}^{w\left(a,c\right)}}{\mathrm{exp}\left(\lambda \left(S\left(b,c\right)-{S}_{a}\right)\right)}.$ (B.11)

$S\left(b,c\right)$ is the value of ${S\left(t\right)|}_{t=c}$ given a value of ${S}_{b}$ at time $t=b$ . The

interim value b does not appear in Equation (B.11). If one replaces $S\left(b,c\right)$ in Equation (B.11) with $S\left(a,c\right)={S}_{c}$ , since the initial value for the time series is ${S}_{a}$ , then one has

$S\left(a,b\right)+S\left(b,c\right)=S\left(a,c\right)={S}_{c}=\frac{{S}_{a}{\text{e}}^{w\left(a,c\right)}}{\mathrm{exp}\left(\lambda \left({S}_{c}-{S}_{a}\right)\right)}.$ (B.12)

This completes the demonstration and illustrates a useful property of the homogeneously saturated solution. It is possible to simulate a time series by adding increments to the series or by solving for the value over the full time interval. If one chooses to add increments, then one must be careful to use the appropriate initial condition for the increment. This appropriate initial condition is the value just before the increment starts.

^{1}
${M}_{\text{avg}\mathrm{,}ST}\left({\lambda}_{0}\right)$
does not yield a minimum value for the VBF data. This is because the low
$\nu $
parameter for VBF ensures that extreme price swings are realized even when the saturation parameter is large (see text). Essentially
$\lambda =\infty $
, a straight line, gives the lowest average. The minimum M was used in the optimization to find
${\lambda}_{0}$
for this case only.