1 | initial version |
I usually use scipy for this:
sage: import scipy.stats
sage: scipy.stats.[TAB]
Display all 213 possibilities? (y or n)
scipy.stats.Tester scipy.stats.dweibull scipy.stats.gmean scipy.stats.linregress scipy.stats.pointbiserialr scipy.stats.t
scipy.stats.alpha scipy.stats.entropy scipy.stats.gompertz scipy.stats.loggamma scipy.stats.poisson scipy.stats.test
scipy.stats.anderson scipy.stats.erlang scipy.stats.gumbel_l scipy.stats.logistic scipy.stats.powerlaw scipy.stats.threshold
scipy.stats.anglit scipy.stats.expon scipy.stats.gumbel_r scipy.stats.loglaplace
[.. etc.]
sage: B = scipy.stats.binom(10, 0.5)
sage: B
<scipy.stats.distributions.rv_frozen object at 0x6aba7d0>
sage: B.[TAB]
B.args B.cdf B.dist B.entropy B.interval B.isf B.kwds
B.mean B.median B.moment B.pdf B.pmf B.ppf B.rvs
B.sf B.stats B.std B.var
sage: B.mean()
5.0
sage: B.rvs(15)
array([6, 5, 5, 5, 4, 7, 4, 7, 6, 4, 6, 7, 4, 7, 5])