Running maximum for Geometric Brownian Motion

Can anyone provide the expression and source for the running maximum $M_t$ for geometric Brownian motion $X_t$ as a function of the initial value $X_0$, drift $\mu$ and diffusion $\sigma$? $X_t$ evolves as

$X_t = X_0*exp[(\mu-\sigma^2/2)t + \sigma W_t]$

where $W_t$ is a Wiener process

• What do you mean by "expression and source"? You can say a lot about this running maximum by noticing that an equivalent problem is to describe the maximum of $(\mu - \frac{\sigma^2}{2})t + \sigma W_t$. This is a Brownian Motion with drift and you can compute the CDF and Laplace transform of a BM with drift exactly by using Girsanov's theorem. – Chris Janjigian Jan 15 '14 at 3:48
• Chris, thanks for your response. By "expression and source", I am requesting the closed form solution for M(t) and a citation of the book or paper where this is published. I will explore the approach you described to find the BM maximum. My background in stochastic mathematics is not very strong, so if you can point me to the result, it would be greatly appreciated. – RK94 Jan 15 '14 at 5:11
• The maximum is highly path dependent, so you're not going to get a closed form expression. The best you can do is talk about its distribution. For a single time, the distribution of $M(t)$ is computable. This is an equivalent question to finding the the CDF of the maximum of a BM plus drift: math.stackexchange.com/questions/133628/… – Chris Janjigian Jan 15 '14 at 16:26
• Thanks Chris. I meant to say closed form expression for the density of M(t), so your response addresses my question. I am modeling a dataset where the event frequencies appear to evolve as the distribution of maximum integers of a GBM with drift, based on simulations. The simulations are slow and I am searching for a more efficient way to analyze the distributions for large datasets. Thanks again. – RK94 Jan 16 '14 at 1:30

1. First of all, it will be easier to find the distribution of the running maximum of $\log X$ which in your notation will be $\log M$.
2. Assume for a moment that $\log X$ has zero drift. Reflection principle then tells you that for any $a > \log X_0$ $$P(\log M_t > a) = 2 P(\log X_t > a ).\tag{1} \label{1}$$ Similarly for $a<\log X_0$, with reversed inequalities in \eqref{1}. Now since $\log X_t$ is Gaussian, one can evaluate the rhs of \eqref{1}. This gives you CDF and after differentiation a PDF of $\log M_t$. This should be enough for your purposes.