*To*: Debian Bug Tracking System <submit@bugs.debian.org>*Subject*: Bug#249546: ITP: libranlip -- generates random variates with multivariate Lipschitz density*From*: Anibal Monsalve Salazar <A.Monsalve.Salazar@IEEE.org>*Date*: Tue, 18 May 2004 14:03:20 +1000*Message-id*: <20040518040320.GF2577@tairona.apana.org.au>*Reply-to*: Anibal Monsalve Salazar <A.Monsalve.Salazar@IEEE.org>, 249546@bugs.debian.org

Package: wnpp Severity: wishlist * Package name : libranlip Version : 1.0 Upstream Author : Gleb Beliakov <gleb@deakin.edu.au> * URL : http://www.deakin.edu.au/~gleb/ranlip.html * License : GPL Description : generates random variates with multivariate Lipschitz density RanLip generates random variates with an arbitrary multivariate Lipschitz density. While generation of random numbers from a variety of distributions is implemented in many packages (like GSL library http://www.gnu.org/software/gsl/ and UNURAN library http://statistik.wu-wien.ac.at/unuran/), generation of random variate with an arbitrary distribution, especially in the multivariate case, is a very challenging task. RanLip is a method of generation of random variates with arbitrary Lipschitz-continuous densities, which works in the univariate and multivariate cases, if the dimension is not very large (say 3-10 variables). Lipschitz condition implies that the rate of change of the function (in this case, probability density p(x)) is bounded: |p(x)-p(y)|<M||x-y||. From this condition, we can build an overestimate of the density, so called hat function h(x)>=p(x), using a number of values of p(x) at some points. The more values we use, the better is the hat function. The method of acceptance/rejection then works as follows: generatea random variate X with density h(x); generate an independent uniform on (0,1) random number Z; if p(X)<=Z h(X), then return X, otherwise repeat all the above steps. RanLip constructs a piecewise constant hat function of the required density p(x) by subdividing the domain of p (an n-dimensional rectangle) into many smaller rectangles, and computes the upper bound on p(x) within each of these rectangles, and uses this upper bound as the value of the hat function. See http://www.deakin.edu.au/~gleb/ranlip.html for more information. Anibal Monsalve Salazar -- .''`. Debian GNU/Linux | Building 28C : :' : Free Operating System | Monash University VIC 3800, Australia `. `' http://debian.org/ | http://www-personal.monash.edu/~anibal/ `- |

**Attachment:
pgpI7I3XPH28H.pgp**

- Prev by Date:
**Re: @debian.org email forwarding and SPF** - Next by Date:
**Re: dh_testdir: I have no package to build.** - Previous by thread:
**O: imwheel -- Program to support the "wheel" on most new mice** - Next by thread:
**is there a procedure for reporting discouraging maintainer bugresponses?** - Index(es):