HW2 Q1 (generating a binomial distribution)

[vahuja4]

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?

[willgao]

nu is the fraction, not the number of the head count.