Code covered by the BSD License  

Highlights from
Generation of Random Variates

image thumbnail

Generation of Random Variates



generates random variates from over 870 univariate distributions

transgam_pdf(x, alpha, mu, c)
% transgam_pdf.m - evaluates a Transformed Gamma Probability Density.
%   See "On Recursive evaluation of Mixed Poisson Probabilities and Related Quantities", 
%   G.E.Willmot, Scand. Actuarial J., vol.2, p.114, 1993. 
%  Created by Jim Huntley,  05/06/09
%   Vector Form of the PDF!!!

function [pdf] = transgam_pdf(x, alpha, mu, c)

%persistent coef

    %coef = mu^alpha * c / gamma(alpha);
    coef = exp(alpha*log(mu) +log(c) - gammaln(alpha));

pdf = coef .* x.^(c*alpha-1) .* exp(-mu.*x.^c);


Contact us