spacepaste

Paste details

reply | raw

  1.  
  2. data = """100m 14 Aug 2016
  3. Usain Bolt 21 AUG 1986
  4. Justin Gatlin 10 FEB 1982
  5. Andrew De Grasse 10 NOV 1994
  6. Yohan Blake 26 DEC 1989
  7. Akani Simbine 21 SEP 1993
  8. Ben Youssef Meite 11 NOV 1986
  9. Jimmy Vicaut 27 FEB 1992
  10. Trayvon Bromell 10 JUL 1995
  12. 200m 18 Aug 2016
  13. Usain Bolt 21 AUG 1986
  14. Andre De Grasse 10 NOV 1994
  15. Christophe Lemaitre 11 JUN 1990
  16. Adam Gemili 06 OCT 1993
  17. Churandy Martina 03 JUL 1984
  18. Lashawn Merritt 27 JUN 1986
  19. Alonso Edward 08 DEC 1989
  20. Ramil Guliyev 29 MAY 1990
  22. 400m 14 Aug 2016
  23. Wayde Van Niekerk 15 Jul 1992
  24. Kirani James 1 Sep 1992
  25. Lashawn Merritt 27 Jun 1986
  26. Machel Cedenio 6 Sep 1995
  27. Karabo Sibanda 2 Jul 1998
  28. Ali Khamis Khamis 30 Jun 1995
  29. Bralon Taplin 8 May 1992
  30. Matthew Hudson-Smith 26 Oct 1994
  32. 800m 15 Aug 2016
  33. David Lekuta Rudisha 17 Dec 1988
  34. TAOUFIK MAKHLOUFI 29 APR 1988
  35. CLAYTON MURPHY 26 FEB 1995
  36. PIERRE-AMBROISE BOSSE 11 MAY 1992
  37. FERGUSON CHERUIYOT ROTICH 30 NOV 1989
  38. MARCIN LEWANDOWSKI 13 JUN 1987
  39. ALFRED KIPKETER 28 DEC 1996
  40. BORIS BERIAN 19 DEC 1992
  42. 1500m 20 Aug 2016
  43. MATTHEW CENTROWITZ 18 OCT 1989
  44. TAOUFIK MAKHLOUFI 29 APR 1988
  45. NICHOLAS WILLIS 25 APR 1983
  46. AYANLEH SOULEIMAN 03 DEC 1992
  47. ABDALAATI IGUIDER 25 MAR 1987
  48. DAVID BUSTOS 25 AUG 1990
  49. BEN BLANKENSHIP 15 DEC 1988
  50. RYAN GREGSON 26 APR 1990
  52. 5000m 20 Aug 2016
  53. MOHAMED FARAH 23 MAR 1983
  54. PAUL KIPKEMOI CHELIMO 27 OCT 1990
  55. HAGOS GEBRHIWET 11 MAY 1994
  56. MOHAMMED AHMED 05 JAN 1991
  57. BERNARD LAGAT 12 DEC 1974
  58. ANDREW BUTCHART 14 OCT 1991
  59. ALBERT KIBICHII ROP 17 JUL 1992
  60. JOSHUA KIPRUI CHEPTEGEI 12 SEP 1996
  62. 10000m 13 Aug 2016
  63. MOHAMED FARAH 23 MAR 1983
  64. PAUL KIPNGETICH TANUI 22 DEC 1990
  65. TAMIRAT TOLA 11 AUG 1991
  66. YIGREM DEMELASH 26 JAN 1994
  67. GALEN RUPP 08 MAY 1986
  68. JOSHUA KIPRUI CHEPTEGEI 12 SEP 1996
  69. BEDAN KAROKI MUCHIRI 21 AUG 1990
  70. ZERSENAY TADESE 08 FEB 1982
  72. 42195m 21 Aug 2016
  73. Eliud Kipchoge 5 NOV 1984
  74. FEYISA LILESA 01 FEB 1990
  75. GALEN RUPP 08 MAY 1986
  76. GHIRMAY GHEBRESLASSIE 14 NOV 1995
  77. ALPHONCE FELIX SIMBU 14 FEB 1992
  78. JARED WARD 09 SEP 1988
  79. TADESSE ABRAHAM 12 AUG 1982
  80. MUNYO SOLOMON MUTAI 22 OCT 1992
  81. """
  83. data2 = """
  84. 100m 27 JUL 1996
  85. 30 DEC 1973
  86. 21 JUN 1972
  87. 12 DEC 1974
  88. 02 NOV 1975
  89. 31 AUG 1968
  90. 28 DEC 1973
  91. 11 AUG 1963
  92. 30 DEC 1970
  94. 200m 1 AUG 1996
  95. 02 OCT 1967
  96. 31 DEC 1965
  97. 30 MAR 1976
  98. 23 JAN 1970
  99. 24 FEB 1975
  100. 30 NOV 1962
  101. 16 APR 1967
  102. 01 NOV 1970
  104. 400m 29 JUL 1996
  105. 31 MAR 1966
  106. 15 NOV 1969
  107. 23 NOV 1972
  108. 16 NOV 1973
  109. 11 MAY 1970
  110. 03 MAR 1968
  111. 04 APR 1975
  112. 18 NOV 1971
  114. 800m 31 JUL 1996
  115. 16 SEP 1972
  116. 07 JUN 1972
  117. 27 MAY 1961
  118. 09 APR 1974
  119. 27 JUN 1969
  120. 06 OCT 1972
  121. 20 NOV 1968
  123. 1500m 3 AUG 1996
  124. 28 FEB 1970
  125. 16 FEB 1969
  126. 24 OCT 1970
  127. 20 JAN 1969
  128. 22 FEB 1964
  129. 28 DEC 1962
  130. 03 NOV 1972
  131. 24 APR 1976
  133. 5000m 3 AUG 1996
  134. 9 DEC 1973
  135. 26 JUN 1970
  136. 7 AUG 1969
  137. 9 FEB 1965
  138. 27 SEP 1971
  139. 18 AUG 1970
  140. 25 FEB 1968
  141. 15 APR 1968
  143. 10000m 29 JUL 1996
  144. 18 APR 1973
  145. 17 JUN 1969
  146. 16 JAN 1972
  147. 18 JUN 1966
  148. 12 DEC 1973
  149. 25 JUN 1969
  150. 29 JAN 1967
  151. 14 JUL 1964
  153. 42195m 4 AUG 1996
  154. 15 APR 1971
  155. 11 OCT 1970
  156. 19 DEC 1973
  157. 3 MAR 1963
  158. 6 JAN 1964
  159. 9 JAN 1968
  160. 26 SEP 1962
  161. 7 AUG 1962
  162. """
  163. # Sources: https://www.olympic.org/ and Wikipedia
  164. # (the official results site is some missing information on some finals and the race dates)
  166. from datetime import datetime
  167. from matplotlib import pyplot as plt
  168. import numpy as np
  170. def conv(x):
  171. return datetime.strptime(x, '%d %b %Y')
  173. def extract(data):
  174. lines = [x for x in data.split("\n") if x]
  175. dates = []
  176. dist = []
  177. runners = []
  179. for line in lines:
  180. if "\t" not in line: line = " \t" + line
  181. a, b = line.split("\t")
  182. if line[0].isdigit():
  183. dist.append(int(a[:-1]))
  184. dates.append(conv(b))
  185. runners.append([])
  186. else:
  187. runners[-1].append(conv(b))
  189. assert(all(len(x) in [7, 8] for x in runners))
  190. avg = [sum((dates[i] - y).days for y in runners[i]) / 365.2422 / len(runners[i]) for i in xrange(len(dist))]
  191. for i in xrange(len(dist)):
  192. print "%6sm: %.1f" % (dist[i], avg[i])
  194. even = sum(sum([x.year % 2 == 0 for x in y]) for y in runners)
  195. total = sum(len(x) for x in runners)
  196. print "%.1f%% (%d out of %d) born in even year" % (100.*even/total, even, total)
  197. return dates, dist, runners, avg
  199. dates, dist, runners, avg = extract(data)
  200. plt.semilogx(dist, avg, 'r.', mew=5, ms=10, label='Rio 2016')
  201. plt.plot(dist, np.poly1d(np.polyfit(np.log(dist), avg, 1))(np.log(dist)), 'r--')
  202. plt.tick_params(axis='x', which='minor', bottom='off', top='off')
  204. dates, dist, runners, avg = extract(data2)
  205. plt.semilogx(dist, avg, 'b.', mew=5, ms=10, label='Atlanta 1996')
  206. plt.plot(dist, np.poly1d(np.polyfit(np.log(dist), avg, 1))(np.log(dist)), 'b--')
  208. plt.xticks(dist, map(str, dist))
  209. plt.xlabel('Race distance (m)', fontsize=18)
  210. plt.xlim([plt.xlim()[0] - 20, plt.xlim()[1]])
  211. plt.ylim([plt.ylim()[0] - 0.2, plt.ylim()[1] + 0.2])
  212. plt.ylabel('Mean age of finalists (years)', fontsize=18)
  213. plt.title('Average age of Olympic finalists with race distance', fontsize=20)
  214. plt.legend(loc='upper left', numpoints=1)
  215. plt.tight_layout()
  216. plt.show()
  217.  
Powered by spacepaste. Paste removal.