Expectation, Variance and Correlation of a bivariate Lognormal distribution

Question

If $Y \sim N(\mu,\sigma^2)$ is normally distributed, then $X=\mathrm{e}^Y$ is lognormally distributed. To get the log-$\mu$ and log-$\sigma$ of this lognormal distribution you calculate $$\sigma^2 = \ln\left( \frac{\mathit{Var}}{E^2} + 1 \right)$$ and $$\mu = \ln(E)-\frac{\sigma^2}{2}$$.

But how is this in the case of a bivariate Normaldistribution, $Y \sim N((\mu_1, \mu_2),\Sigma)$, $\Sigma=((\sigma_1,\rho),(\rho,\sigma_2))$ How do I get $\mu, \sigma, \rho$? And more general in the case of a multivariate distribution?

Answer 1

$\text{Cov}(X_1,X_2)=E(X_1X_2)-E(X_1)E(X_2)$

$E(X_1X_2)=E(e^{Y_1+Y_2})$

Now the distribution of $Y_1+Y_2$ is normal (and straightforward), so $E(e^{Y_1+Y_2})$ is just the expectation of a univariate lognormal.

The $E(X_1)E(X_2)$ term you can already do.

As a result, it's straightforward to write $\text{Cov}(X_1,X_2)$ in terms of $\mu,\sigma$ and $\rho$ and thereby to solve for $\rho$.

---

$Y_1+Y_2\sim N(\mu_1+\mu_2,\sigma_1^2+\sigma_2^2+2\rho\sigma_1\sigma_2)$, so

$e^{Y_1+Y_2}\sim logN(\mu_1+\mu_2,\sigma_1^2+\sigma_2^2+2\rho\sigma_1\sigma_2)$,

which has expectation $\exp[\mu_1+\mu_2+\frac{1}{2}(\sigma_1^2+\sigma_2^2+2\rho\sigma_1\sigma_2)]$.

$E(X_i)=\exp(\mu_i+\frac{1}{2}\sigma_i^2)$

So $\text{Cov}(X_1,X_2)=E(X_1)E(X_2)[\exp(\rho\sigma_1\sigma_2)-1]$

And hence:

$\exp(\rho\sigma_1\sigma_2)-1=\frac{\text{Cov}(X_1,X_2)}{E(X_1)E(X_2)}$

$\rho=\log(\frac{\text{Cov}(X_1,X_2)}{E(X_1)E(X_2)}+1)\cdot\frac{1}{\sigma_1\sigma_2}$

You can extend this approach to calculating $\rho_{ij}$ from $\text{Cov}(X_i,X_j)$ and the other quantities.

However, if you're trying to do this to estimate parameters from a sample, using sample moments of a lognormal to do parameter estimation (i.e. method-of-moments) doesn't always perform all that well. (You might consider MLE if you can.)