urplanck function

UNU.RAN Planck random variate generator

UNU.RAN Planck random variate generator

UNU.RAN random variate generator for the Planck distribution with shape parameter a. It also allows sampling from the truncated distribution.

[Special Generator] -- Sampling Function: Planck.

urplanck(n, a, lb = 1.e-12, ub = Inf)

Arguments

  • n: size of required sample.
  • a: (strictly positive) shape parameter.
  • lb: lower bound of (truncated) distribution.
  • ub: upper bound of (truncated) distribution.

Details

The Planck distribution with parameter a

has density proportional to

f(x)xaexp(x)1f(x) (xa)/(exp(x)1) f(x) \sim \frac{x^a}{\exp(x)-1}f(x) ~ (x^a)/(exp(x)-1)

for x0x \ge 0 and a>=1a >= 1.

The generation algorithm uses transformed density rejection TDR . The parameters lb and ub can be used to generate variates from the Planck distribution truncated to the interval (lb,ub).

See Also

runif and .Random.seed about random number generation and unuran for the UNU.RAN class.

References

W. H"ormann, J. Leydold, and G. Derflinger (2004): Automatic Nonuniform Random Variate Generation. Springer-Verlag, Berlin Heidelberg

Author(s)

Josef Leydold and Wolfgang H"ormann unuran@statmath.wu.ac.at .

Note

This function is wrapper for the UNU.RAN class in .

Examples

## Create a sample of size 1000 x <- urplanck(n=1000,a=2)
  • Maintainer: Josef Leydold
  • License: GPL (>= 2)
  • Last published: 2025-04-07