How are copulas used in the real world? I have been reading about copula models. Essentially, copula models seem to be a creative method for creating a joint probability distribution from several variables, in which each individual variable can belong to a different probability distribution.
My Question: How are copula models actually used in the real world?
For example, suppose you have data corresponding to 5 variables $X_1, X_2, X_3, X_4$ and $X_5$. Let's say that you fit a copula model to these data.
Given values of $X_2, X_3, X_4$ and $X_5$, suppose you want to find the most probable value of $X_1$. At first glance, it seems like you could consider this as a conditional probability distribution $X_1 | X_2 = x_2, X_3 = x_3, X_4 = x_4, X_5 = x_5\dots$ therefore, you could use the Metropolis Hastings algorithm to draw many random samples from the above conditional probability distribution function, and then take the expected value of these samples as your answer?
Have I understood this correctly? Is this the main application of copula models?
Note: I have read about copula regression that relates more closely to GLM regression models – but it seems like the approach I outlined above might be closer to the intended use for copulas?
 A: Generalized estimating equations might be considered an example of this idea.  The joint likelihood of correlated data is specified through marginal distributions and a correlation structure.  If you were to verify the operating characteristics of confidence intervals for such a model you could simulate correlated normal observations, apply the the probability integral transformation to turn these into uniform observations, and apply the probability integral transformation once again to turn these into another multivariate distribution of interest.
If you are interested in a feature of the distribution of $X_1$ (e.g., the mean, a percentile, etc.) given values of $X_2$ through $X_5$, these other variables can be included as covariates in a univariate model for $X_1$.  When it comes to inference one can apply a link function to the estimator so that it is approximately normally distributed and use the delta method to estimate the standard error and invert a Wald test to produce a confidence interval.
A similar solution can be approximated through simulation, but it is computationally intensive.
