MathGroup Archive 2007

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: RandomArray from user defined distribution?


rob wrote:
> I'd like to use the RandomArray to produce some data from 
> what I think is a Poisson distribution in t
> P[t] = a^2 t Exp[-a*t]  where a is mean, sigma.
> 
> I see one can use RandomArray to produce sample data from a 
> lot of continuous distributions but the Poisson isn't among 
> them (it's only available in the discrete form).
> 
> I've made a bunch of crippled attempts to force my P[t] to 
> put out examples but have failed. Any suggestions? Thanks.
> 

Your density is a GammaDistribution[2,3], se e.g.

http://mathworld.wolfram.com/GammaDistribution.html

If you want to implement a Random from your own Distribution, the 
easiest way is to generate a 1000 point Table of (P(x),x) of paired 
values of the inverted cumulated Distribution function P^-1, construct 
an interpolating  function iP using Fit and use the Random generator 
choose:=ip[Random[Real,{0,1}]]

-- 

Roland Franzius


  • Prev by Date: Re: numerical_solution
  • Next by Date: Re: conditional plotstyles in ListPlot
  • Previous by thread: Re: RandomArray from user defined distribution?
  • Next by thread: Re: RandomArray from user defined distribution?