[統計][R]
>>> import math
>>> cmb = lambda m,n:math.factorial(m) / (math.factorial(n) * math.factorial(m-n))
>>> cmb(10,2)
45.0
>choose(10, 2)
[1] 45
>>> dbinom = lambda k,n,p:cmb(n,k) * math.pow(p,k) * math.pow((1-p),(n-k))
>>> dbinom(3,10,0.5)
0.1171875
> dbinom(3,10,0.5)
[1] 0.1171875
>>> pbinom = lambda q,n,p:sum([dbinom(k,n,p) for k in range(0,q+1)])
>>> pbinom(3,10,0.5)
0.171875
> pbinom(3,10,0.5)
[1] 0.171875
import math
cmb = lambda m,n:math.factorial(m) / (math.factorial(n) * math.factorial(m-n))
dbinom = lambda k,n,p:cmb(n,k) * math.pow(p,k) * math.pow((1-p),(n-k))
pbinom = lambda q,n,p:sum([dbinom(k,n,p) for k in range(0,q+1)])
YAGI Hiroto (piroto@a-net.email.ne.jp)
twitter http://twitter.com/pppiroto
Copyright© 矢木 浩人 All Rights Reserved.