ｓこｐｐｐ - 経済統計練習帳

from scipy.stats import binom
 
p=[0.25,0.5,0.75]
pre=[0.3,0.4,0.3]
n = 8
z=[]
for i in range(3):
    rv =binom(n, p[i])
    z.append(pre[i]*rv.pmf(3))
 
z1=0
for i in range(3): 
   z1=z1+z[i]
 
for i in range(3):
    z[i]=z[i]/z1
 
print(z)