Tk Library Source Code

View Ticket
Login
Ticket UUID: 1859693
Title: Submission: Poisson distribution
Type: RFE Version: None
Submitter: erickb Created on: 2007-12-28 13:50:47
Subsystem: math Assigned To: arjenmarkus
Priority: 5 Medium Severity:
Status: Closed Last Modified: 2008-01-11 18:20:02
Resolution: Closed By: arjenmarkus
    Closed on: 2008-01-11 11:20:02
Description:
Hi,

I'm submitting a file that implements the standard 3 functions for a Poisson distribution:

pdf-poisson (*)
cdf-poisson
random-poisson

(*) For discrete distributions, this is properly a probability mass function (pmf). Would it be helpful or annoying to call this pmf-poisson?

The random-poisson function implements a pair of methods for generating Poisson-distributed random variates so it's reasonably fast at both small and large mu.

The pdf-poisson function has a very direct implementation.

The cdf-poisson function uses the incomplete gamma function that I submitted in a separate submission. I can re-send it if you like.

Thanks!
Eric
User Comments: arjenmarkus added on 2008-01-11 18:20:02:
Logged In: YES 
user_id=400048
Originator: NO

Incorporated this code

erickb added on 2008-01-01 06:08:48:

File Deleted - 259984: 



File Added - 260327: poisson.zip

Logged In: YES 
user_id=816411
Originator: YES

I'm submitting a collection of files as a zip archive. These include:
 * poisson.tcl
    - pdf-poisson
    - cdf-poisson
    - random-poisson
    - Definitions of tests & helper procs
 * possion_test.tcl
    - The code that runs the tests, including a visual comparison of a histogram produced by random-poisson and the pdf
 * simpletest.tcl
    - The code to run the non-visual tests
 * incgamma.tcl
    - Implementation of the incomplete gamma function

This code is more explicit about the tests.

Also, it includes one bug fix: the cdf should have been evaluated at floor(x), rather than at x. This is now corrected.

Eric

File Added: poisson.zip

erickb added on 2007-12-28 20:50:47:

File Added - 259984: random-poisson.tcl

Attachments: