HW2 Q1 (generating a binomial distribution)
Hi All,
Here is the code I wrote for hw2_q1:
from random import randint
import seaborn as sns
import numpy as np
N = 100000
v_first = np.zeros(N)
v_min = np.zeros(N)
v_rand = np.zeros(N)
#Loop starts
for i in range(N):
n,p = 10, 0.5
flips = np.random.binomial(n, p, 1000)
first_coin = 0
min_coin = np.argmin(flips)
rand_coin = np.random.randint(0,1000)
v_first[i] = flips[first_coin]
v_min[i] = flips[min_coin]
v_rand[i] = flips[rand_coin]
#Loop ends
print np.average(v_min)
My answer turns out to be 0.367, while the correct answer is 0.01. Could someone please point out where I am going wrong?
|