Univariate Multimodal Random Number Generator (Scripts) Publisher's description
from Avan Suinesiaputra
Sometimes we need samples that come from a mixture of pdf's, where more than 1 peaks are observed. E.g
Sometimes we need samples that come from a mixture of pdf's, where more than 1 peaks are observed. E.g. for analyzing goodness of a non-parametric kernel density estimation method.
By using the Rejection method, this function generates random numbers from a mixture (average sum) of N normal (Gaussian) distributions. It gives also the mixture probability density function (pdf) as an optional output.
An example of 500 samples generated from 3 Gaussian distributions are shown in the figure, i.e. by the following command:
>> X = umgrn([-4 0 5],[1 2 1.5],500);
>> figure; hist(X);
See 'help umgrn' for more details.
System Requirements:MATLAB 6.5.1 (R13SP1)
Program Release Status: New Release
Program Install Support: Install and Uninstall