march and june there were more anorexics born than averageand more in june itself let' look at that worrisome statistic for those women born in june the team studied women who had been diagnosed as anorexicso the mean number of births per month was slightly more than this suggests that the number born in june was ( * let' write short program to see if we can reject the null hypothesis that this occurred purely by chance def juneprob(numtrials)june for trial in range(numtrials)june for in range( )if random randint( , = june + if june > june + jprob june /float(numtrialsprint 'probability of at least births in june ='jprob figure probability of anorexics being born in june when we ran juneprob( it printed probability of at least births in june it looks as if the probability of at least babies being born in june purely by chance is around so perhaps those researchers in aberdeen are on to something wellthey might have been on to something had they started with the hypothesis that more babies who will become anorexic are born in juneand then run study designed to check that hypothesis but that is not what they did insteadthey looked at the data and thenimitating the texas sharpshooterdrew circle around june the right statistical question to have asked is what is the probability that there was at least one month (out of in which at least babies were born the program in figure answers that question def anyprob(numtrials)anymonth for trial in range(numtrials)months [ ]* for in range( )months[random randint( , )+ if max(months> anymonth + aprob anymonth /float(numtrialsprint 'probability of at least births in some month ='aprob figure probability of anorexics being born in some month