random numbers from custom mixture distribution

2 次查看(过去 30 天)
I would like to generate random numbers from a mixture distribution defined as,
pdf_mix = p*gampdf(x,a,b) + (1-p)*lognpdf(x,mu,sigma);
I am able to obtain the best fit model using mle, but I also want random numbers generated from this mixture distribution. Is this possible? Can I get some help on this?
R Venkatessh

采纳的回答

Roger Stafford
Roger Stafford 2013-10-6
Here's another method you might consider that would allow you to directly use the 'gamrnd' and 'lognrnd' functions in the Statistics Toolbox. Suppose you want an m by n matrix, R, of random values in accordance with your mixed distribution.
U = rand(m,n)<=p; % U is true with probability p and false with probability 1-p
R = U.*gamrnd(a,b,m,n)+(1-U).*lognrnd(mu,sigma,m,n);
R should have the desired distribution.

更多回答(1 个)

Image Analyst
Image Analyst 2013-10-6
Calculate the CDF and then get a random number from the uniform distribution and plug it in. See the method described here: http://en.wikipedia.org/wiki/Inverse_transform_sampling I've also attached my program for doing it with a rayleigh distribution. Click on the filename below (not the picture, but the filename).

Community Treasure Hunt

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

Start Hunting!

Translated by