MATLAB Answers

sample from arbitrary continuous distribution

조회 수: 72(최근 30일)
Xin
Xin 8 Jan 2017
댓글: Xin 20 Feb 2017
Say I have a very complicated probability distribution function: x^(-2)*y^(-3)*exp(-z/t), where x y z t are all continuous positive numbers. The integration in x,y,z,t space is not 1, so there has to be some scaling factor. How can I devise a function, or is there a convenient function with some symbolic that allows me to sample a random, continuous set of numbers from this distribution? What if I want to sample a full conditional distribution, say x, with y,z,t as known?
Many thanks.

  댓글 수: 1

the cyclist
the cyclist 8 Jan 2017
Are your x, y, and z correlated with each other, or independent?

댓글을 달려면 로그인하십시오.

채택된 답변

the cyclist
the cyclist 8 Jan 2017
편집: the cyclist 9 Jan 2017
There are two key elements to your problem:
  • Sampling from an arbitrary distribution
  • Sampling from a multivariate distribution (and are those variables correlated?)
To do the first part -- for example just getting samples from x^(-2) -- there is a method using the inverse of the cdf. I suggest you take a look at this thread in the old MATLAB forum to get started.
If you have independent variables, then you can just multiply the distributions of each variable. If you have dependent variables, then your best bet is to use the copula approach.
There are examples in that documentation that will help you get started. I suggest that you try some baby steps for each of these two parts, so that you build your understanding, rather than trying to solve your whole problem at once.

  댓글 수: 1

Xin
Xin 20 Feb 2017
Thank you for your answer, I have got it working.

댓글을 달려면 로그인하십시오.

추가 답변(0개)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by