@@ -27,3 +27,12 @@ The method used here is the one described in Knuth " (x2 (random-gamma b 1d0))) (declare (double-float x1 x2)) (/ x1 (+ x1 x2)))) + + +(defun density-beta (x &key (a 1.0d0) (b 1.0d0)) + "implemented according to docstring in random-beta." + (* (/ (gamma-function (+ a b)) + (* (gamma-function a) + (gamma-function b))) + (^ x (- a 1)) + (^ (- 1 x) (- b 1)))) beta.lisp 4e864fae354b02d0df7909ac786c6864c80be6a8 AJ Rossini blindglobe@gmail.com http://github.com/lvaruzza/cl-randist/commit/a43a6cad427921f4414139ff5aeea8ac436cd180 a43a6cad427921f4414139ff5aeea8ac436cd180 2009-05-11T04:50:41-07:00 2009-05-08T09:11:07-07:00 include density-beta, using formula in doc-string for random-beta. Signed-off-by: AJ Rossini <blindglobe@gmail.com> bccb2ae2c1a3dbb732a4187afb40b6e0af132839 AJ Rossini blindglobe@gmail.com