import random
# flip a coin
random.choice([0, 1])
0
N = 100
flips = []
for i in range(N):
flip = random.choice([0, 1])
flips.append(flip)
print flips
[1, 1, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 1, 0, 1, 1, 1, 1, 1, 0, 0, 1, 0]
print sum(flips) / float(len(flips))
0.47
hist(flips)
(array([53, 0, 0, 0, 0, 0, 0, 0, 0, 47]), array([ 0. , 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1. ]), <a list of 10 Patch objects>)
# write a function to perform many flips
def flip_N(N):
flips = []
for i in range(N):
flip = random.choice([0, 1])
flips.append(flip)
return sum(flips) / float(N)
avg = flip_N(100)
print avg
0.46
# now perform many trials, and show the distribution
n_trials = 500
trials = []
N = 100
for i in range(n_trials):
avg = flip_N(N)
trials.append(avg)
hist(trials)
(array([ 8, 19, 54, 112, 134, 94, 54, 15, 6, 4]), array([ 0.37 , 0.399, 0.428, 0.457, 0.486, 0.515, 0.544, 0.573, 0.602, 0.631, 0.66 ]), <a list of 10 Patch objects>)
ls
coin-flips-evokerr.ipynb coin-flips-hljordt.ipynb coin-flips.ipynb coin-flips-lnj204.ipynb coin-flips-mherrmannsfeldt.ipynb coin-flips-qchrist2.ipynb coin-flips-Richard.ipynb filter-blast-csv-hljordt.ipynb filter-blast-csv.ipynb filter-blast-csv_qchrist2.ipynb graph-blast-bitscore-Copy0.ipynb graph-blast-bitscore-Copy6-Copy0.ipynb graph-blast-bitscore-Copy6.ipynb graph-blast-bitscore-evokerr.ipynb graph-blast-bitscore.ipynb graph-blast-bitscore-leonard.ipynb graph-blast-bitscore-lnj204.ipynb graph-blast-bitscore-mherrmannsfeldt.ipynb graph-blast-bitscore-qchrist2.ipynb graph-blast-bitscore-Richard.ipynb graph-blast-bitscore-sr320.ipynb hljordt.ipynb kmer-abundance-hljordt.ipynb kmer-abundance.ipynb kmer-abundance-lnj204.ipynb kmer-abundance-mherrmannsfeldt.ipynb kmer-abundance-qchrist2.ipynb kmer-abundance-Richard.ipynb kmer-abundance-titus.ipynb monty-hall-evokerr.ipynb monty-hall-hljordt.ipynb monty-hall.ipynb monty-hall-lnj204.ipynb monty-hall-mherrmannsfeldt.ipynb monty-hall-qchrist2.ipynb monty-hall-Richard.ipynb monty-hall-titus.ipynb README.txt