Skip to content
All library documents

Calculating a Stock Price Threshold Probability under Geometric Brownian Motion

Article Quant Q&A · Author: Gus Montano

Summary

The document corrects a probability calculation for a stock modeled by geometric Brownian motion. The attempted solution approximates the future price itself as normally distributed, but under GBM the future price is lognormally distributed and its logarithm is normal. The accepted derivation writes the terminal price as its current value multiplied by an exponential involving drift, volatility, time, and a standard normal variable.

Taking logarithms of the price threshold converts the event into a standard normal tail probability. The document gives both the direct tail expression and an equivalent expression using the opposite sign, and another answer illustrates the calculation with the stated inputs. This method relies on the GBM assumptions and specified drift and volatility; it does not assess whether those assumptions or inputs fit a real stock. It also concerns a physical probability under the stated drift, rather than an option-pricing probability under a risk-neutral measure.

Key ideas

  • Under geometric Brownian motion, the terminal stock price is lognormally distributed.
  • The logarithm of the terminal price is normally distributed with drift adjusted by half the variance rate.
  • A price threshold event can be converted to a standard normal tail probability by taking logarithms.
  • Treating the stock price itself as normally distributed gives a different and unsuitable calculation under GBM.
  • The result depends on the assumed drift, volatility, horizon, and probability measure.

Tags

Full text
# Confidence Intervals of Stock Following a Geometric Brownian Motion


# Confidence Intervals of Stock Following a Geometric Brownian Motion












In preparation for my Options, Future's and Risk Management examination next week, I have been presented with a series of questions and their answers. Unfortunately, my lecturer, one of the less organised, does not respond to emails and attempts for consultation. I have resorted to these forums to relieve some stress.

My question is presented as follows:

> The share price of company XYC Inc. exhibits an instantaneous drift of 7% per year with return volatility of 45%. What is the probability that XYZ shares exceed \$95 after 10 months when they cost $55 today

Of course, I will display my attempted solution.

First, I assume that the change in stock price follows a geometric brownian motion (GBM). That is,

$$\frac{\Delta S}{S_{0}}=\mu \Delta t+\sigma\sqrt{\Delta t}\cdot \varepsilon.$$

Following some algebra,

$$ \begin{align*} \frac{\Delta S}{S_{0}} &=\mu \Delta t+\sigma\sqrt{\Delta t} \cdot \varepsilon \\ \frac{S-S_{0}}{S_{0}} &= \mu \Delta t+\sigma\sqrt{\Delta t} \cdot \varepsilon \\S &= \left(S_{0} + \mu S_{0} \Delta t\right) + \sigma S_{0} \sqrt{\Delta t} \cdot \varepsilon \end{align*} $$

Therefore the distribution of future stock price is given by

$$S \sim \phi\left(S_{0} + \mu S_{0} \Delta t,\left(\sigma S_{0} \sqrt{\Delta t}\right)^{2}\right).$$

Substituting appropriate figures,

$$S \sim \left(58.21, \left(22.59\right)^2\right).$$

For probabilistic problems regarding normal distributions, I relate to standardised scores. I calculate that

$$z_{95} = 1.63.$$

Using Microsoft Excel, the probability that the z-score is greater than 1.63, and therefore, the price of the stock is greater than 95 is given by

$$1- \mathrm{NORMDIST(95,58.21,22.59,TRUE)}.$$

The answer I get is 5.17%. The answer states it is 8.23%.

I would be beyond thankful for any help and advice on how to properly solve this problem.

Thank you in advanced,

Gustavo.

## Answer by Gordon (score 2, accepted)

https://quant.stackexchange.com/a/18439

As the stock price process $S$ follows a geometric Brownian motion, we have that \begin{align*} S_T &= S_0 e^{(\mu-\frac{1}{2}\sigma^2)\, T + \sigma\, W_T}\\ &= S_0 e^{(\mu-\frac{1}{2}\sigma^2)\, T + \sigma\, \sqrt{T}\, \xi}, \end{align*} where $\xi$ is a standard normal random variable. Then, we have the probability \begin{align*} P(S_T > 95) &= P\Big( S_0 e^{(\mu-\frac{1}{2}\sigma^2)\, T + \sigma\, \sqrt{T}\, \xi} > 95\Big)\\ &= P\bigg(\xi > \frac{\ln \frac{95}{S_0} - (\mu-\frac{1}{2}\sigma^2)\, T}{\sigma\, \sqrt{T}} \bigg)\\ &= 1- NORMSDIST\left(\frac{\ln \frac{95}{S_0} - (\mu-\frac{1}{2}\sigma^2)\, T}{\sigma\, \sqrt{T}} \right)\\ &=NORMSDIST\left(\frac{\ln \frac{S_0}{95} + (\mu-\frac{1}{2}\sigma^2)\, T}{\sigma\, \sqrt{T}} \right). \end{align*}

## Answer by Egodym (score 1)

https://quant.stackexchange.com/a/18440

Given its price today, the stock price at time T is lognormally distributed, whereas $lnS_T$ is normally distributed, that is

$lnS_T$ ~ $N \Bigr(lnS_0 + (\mu- \frac{\sigma^2}{2}T),\sigma^2T \Bigl)$

see for example Hull - Options, Futures, and other Derivatives.

Plugging in the numbers you get

$lnS_T$ ~ $N(3.981291519,0.16875)$

Then the probability you want is $P[lnS_T>ln95]=1-P[lnS_T<ln95]=1-P\Bigr[Z<z=\frac{ln95-3.981291519}{\sqrt{0.16875}}\Bigl] \\=1-P[Z<z=1.391]=1-0.9177=0.0823$

Shown in full with attribution under the source's licence. Licence: CC BY-SA 4.0 (Stack Exchange)

This summary was written by Stratmill's research agent from the original; it is not a copy of the source.