
Compute the regression slope for diet 1, get s1.
Compute the regression slope for diet 2, get s2.
Consider the difference d = s1 - s2

For each point we have diet1 weights and diet 2 weights.

10,000 times
Shuffle the label of whether a weight is diet1 or diet2 at each time point.
Then recalculate the difference in the slopes and see how often it's
as big as s1-s2. 


# How do we select a team of 6 students for a quidditch team?
# mixing Gryffindors and Slytherins??



import random

Gryffindor = ["Harry", "Ron", "Hermione", "Fred", "George"]
Slytherin  = ["Draco", "Crabbe", "Goyle", "Snape", "Zabini"]


# SOLUTION 3
QuidditchTeam = []
pool = Gryffindor + Slytherin
random.shuffle(pool)
for i in range(6):
    QuidditchTeam.append(pool.pop())

print ("solution 3 quidditch team: ", QuidditchTeam)

