Groups > comp.lang.java.programmer > #9332

RandomDirichlet?

Path	csiph.com!x330-a1.tempe.blueboxinc.net!usenet.pasdenom.info!gegeweb.org!de-l.enfer-du-nord.net!feeder1.enfer-du-nord.net!fu-berlin.de!uni-berlin.de!individual.net!not-for-mail
From	Giovanni Azua <bravegag@hotmail.com>
Newsgroups	comp.lang.java.programmer
Subject	RandomDirichlet?
Date	Tue, 01 Nov 2011 10:09:32 +0100
Lines	24
Message-ID	<CAD575DC.8551%bravegag@hotmail.com> (permalink)
Mime-Version	1.0
Content-Type	text/plain; charset="US-ASCII"
Content-Transfer-Encoding	7bit
X-Trace	individual.net IjTMsI9mnZiKD6Y40jukIw42KDM1eo8yA0DlGL0iZmoKhqlrXw
Cancel-Lock	sha1:TRgkeb3t78k4zQQ6aVexXrwHDss=
User-Agent	Microsoft-Entourage/12.31.0.110725
Thread-Topic	RandomDirichlet?
Thread-Index	AcyYdfegYCfNJx971kmIB/vYAIgt+g==
Xref	x330-a1.tempe.blueboxinc.net comp.lang.java.programmer:9332

Show key headers only | View raw

Hello,

I need to generate random TPCH workloads but I don't want a Uniform
distribution where Q1, RF1, RF2 would be equally probable i.e. using new
Random().nextInt(..)

I instead would like to draw samples from a Dirichlet distribution which
would give a very skewed PD depending on the input alphas (which specify how
skewed I want the PD to be). The Matlab code would be exactly:

function r = dirichlet_sample(a,n)
p = length(a);
r = gamrnd(repmat(a,n,1),1,n,p);
r = r ./ repmat(sum(r,2),1,p);
End

But there is no "gamrnd" RandomGamma in Java nor RandomDirichlet. Before I
spend time writing this, can anyone advice for a low-footprint library
implementing those functions?

TIA,
Best regards,
Giovanni

Back to comp.lang.java.programmer | Previous | Next — Next in thread | Find similar | Unroll thread

Thread

RandomDirichlet? Giovanni Azua <bravegag@hotmail.com> - 2011-11-01 10:09 +0100
  Re: RandomDirichlet? Eric Sosman <esosman@ieee-dot-org.invalid> - 2011-11-01 07:50 -0400
    Re: RandomDirichlet? Giovanni Azua <bravegag@hotmail.com> - 2011-11-01 19:50 +0100
  Re: RandomDirichlet? Wojtek <nowhere@a.com> - 2011-11-01 04:54 -0700
  Re: RandomDirichlet? Roedy Green <see_website@mindprod.com.invalid> - 2011-11-01 06:11 -0700
    Re: RandomDirichlet? Roedy Green <see_website@mindprod.com.invalid> - 2011-11-01 06:16 -0700
    Re: RandomDirichlet? Eric Sosman <esosman@ieee-dot-org.invalid> - 2011-11-01 20:23 -0400
      Re: RandomDirichlet? Roedy Green <see_website@mindprod.com.invalid> - 2011-11-02 01:39 -0700
        Re: RandomDirichlet? Eric Sosman <esosman@ieee-dot-org.invalid> - 2011-11-02 08:17 -0400
  Re: RandomDirichlet? Arne Vajhøj <arne@vajhoej.dk> - 2011-11-01 22:09 -0400

csiph-web