source('read_superbowl_data.R')
sb = read_superbowl_data()
##print(sb)

## Did one team begin closer to the goal:
##
print(t.test(sb[sb$Team =='N', ]$Starting_Distance, sb[sb$Team =='C', ]$Starting_Distance) )


## Build a large model to anwser many questions:
##
mp = glm(Score ~ Drive_Number + Starting_Distance, data=sb, family=binomial)
print(summary(mp))


## An interaction between Team and Drive_Number:
##
mp = glm(Score ~ Team * Drive_Number, data=sb, family=binomial)
print(summary(mp))


## An interaction between Team, Drive_Number, and Starting_Distance:
##
mp = glm(Score ~ Team * ( Drive_Number + Starting_Distance ), data=sb, family=binomial)
print(summary(mp))