Mercurial > hg > python

view distr_2.py @ 28:bb46148fab91

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
second half of old /c/Drive_F/merged_d/src/python/distr.py, modified to do coin flip experiment
author Henry S. Thompson <ht@inf.ed.ac.uk>
date Wed, 28 Jul 2021 17:54:52 +0100
parents
children
line wrap: on
line source

# experimenting with distributions
# Usage: distr.py N bucketWidth max
import sys, random, signif
from nltk import FreqDist

mu=float(sys.argv[1])
sigma=float(sys.argv[2])
n=int(sys.argv[3])
scale=float(sys.argv[4])
fd=FreqDist()
fd2=FreqDist()
for i in xrange(n):
  fd[round(random.normalvariate(mu,sigma)/scale)]+=1
  fd2[round(random.gauss(mu,sigma)/scale)]+=1
fd.bell(xtra=fd2)