From: David Arnold <dwarnold45@suddenlink.net>
To: mailing list for ConTeXt users <ntg-context@ntg.nl>
Subject: Re: Random numbers in metapost
Date: Thu, 31 Dec 2009 10:17:02 -0800 [thread overview]
Message-ID: <D10D8E19-0D72-4B2B-BEE7-F46D4DC1B17B@suddenlink.net> (raw)
In-Reply-To: <4B3CA38B.7060602@elvenkind.com>
[-- Attachment #1.1: Type: text/plain, Size: 1582 bytes --]
Taco et al,
I think I've found a way. Anthony Phan has written a brilliant stats package for Metapost.
http://www-math.univ-poitiers.fr/~phan/statsmac.html
I think (haven't tried yet) I can select a number at random from the uniform distribution on [0,1], then use one of Anthony's functions such as:
poissonicdf(expr p, lambda). Inverse cumulative distribution function of P(lambda).
See Anthony's manual at:
http://www-math.univ-poitiers.fr/~phan/downloads/metapost/statsman.pdf
David.
On Dec 31, 2009, at 5:13 AM, Taco Hoekwater wrote:
> David Arnold wrote:
>> Hi,
>> Does anyone have any strategies for drawing random numbers from
>> various statistics distributions in Metapost? For example, say,
>> drawing 100 random numbers from a beta distribution?
>
> Assuming you are using mkiv: create a lua function and use its
> calculated results to create the metapost code. Metapost only knows
> about uniformdeviate and normaldeviate, and metapost is not a good
> language to program mathematics in.
>
> Best wishes,
> Taco
> ___________________________________________________________________________________
> If your question is of interest to others as well, please add an entry to the Wiki!
>
> maillist : ntg-context@ntg.nl / http://www.ntg.nl/mailman/listinfo/ntg-context
> webpage : http://www.pragma-ade.nl / http://tex.aanhet.net
> archive : http://foundry.supelec.fr/projects/contextrev/
> wiki : http://contextgarden.net
> ___________________________________________________________________________________
[-- Attachment #1.2: Type: text/html, Size: 2877 bytes --]
[-- Attachment #2: Type: text/plain, Size: 486 bytes --]
___________________________________________________________________________________
If your question is of interest to others as well, please add an entry to the Wiki!
maillist : ntg-context@ntg.nl / http://www.ntg.nl/mailman/listinfo/ntg-context
webpage : http://www.pragma-ade.nl / http://tex.aanhet.net
archive : http://foundry.supelec.fr/projects/contextrev/
wiki : http://contextgarden.net
___________________________________________________________________________________
prev parent reply other threads:[~2009-12-31 18:17 UTC|newest]
Thread overview: 3+ messages / expand[flat|nested] mbox.gz Atom feed top
2009-12-31 7:08 David Arnold
2009-12-31 13:13 ` Taco Hoekwater
2009-12-31 18:17 ` David Arnold [this message]
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to=D10D8E19-0D72-4B2B-BEE7-F46D4DC1B17B@suddenlink.net \
--to=dwarnold45@suddenlink.net \
--cc=ntg-context@ntg.nl \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for NNTP newsgroup(s).