Mercurial > hg > python
view distr_2.py @ 40:f13017bb4502
works on markup with -v ecclerig...
author | Henry S. Thompson <ht@inf.ed.ac.uk> |
---|---|
date | Mon, 03 Jan 2022 19:19:20 +0000 |
parents | bb46148fab91 |
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)