4
$\begingroup$

Let $X\sim \mathcal{N}(\mu,\Sigma)$, where $\mu\in\mathbb{R}^n$ and $\Sigma\in\mathbb{R}^{n\times n}$.

How can I efficiently sample from $X | {\min{X}\le 0}$? (I.e. from the distribution of $X$ conditional on $\min{\{X_1,\dots,X_n\}}\le 0$.)

This seems harder than the usual problem of sampling from a truncated Normal, where Gibbs sampling works OK.

To see that Gibbs sampling is unlikely to perform well in my case, suppose that $\mu=[10,\dots,10]'$ and $\Sigma=I$, and suppose the sampler is initialized with $[0,10,\dots,10]'$. At the next step of the Gibbs sampler, the first coordinate must remain negative (as all others are positive), but with one coordinate negative, the conditional distributions of the other coordinates are just $\mathcal{N}(10,1)$, which is very unlikely to be negative. Thus the sampler will remain stuck with only the first coordinate negative.

$\endgroup$
0
$\begingroup$

One option (though I expect there's a better alternative):

Perform importance sampling, with the following proposal distribution:

  1. Draw $i$ from $\{1,\dots,n\}$ with probability proportional to $\Pr(X_i\le 0)$.
  2. Draw $x_i$ from the distribution of $X_i|X_i\le 0$.
  3. Draw $x_1,\dots,x_{i-1},x_{i+1},\dots,x_n$ from the distribution of $X_1,\dots,X_{i-1},X_{i+1},\dots,X_n|X_i=x_i$.

The advantage of this scheme is that there's no serial correlation in the draws, and all of the conditioning can be performed exactly. The disadvantage is that the weights used in step 1 may give a poor approximation to the true distribution in the presence of correlations, so the importance sampling may be inefficient.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.