Modelng preference tme n mddle dstance trathlons Iztok Fster, 1 Andres Iglesas, 2 Suash Deb, 3, 4 Dušan Fster, 5 and Iztok Fster Jr. 6 1 Unversty of Marbor, Faculty of Electrcal Engneerng and Computer Scence, Smetanova 17, 2000 Marbor, Slovena; Emal: ztok.fster1@um.s 2 Unversty of Cantabra, Avda. de los Castros, s/n, 39005 Santander, Span 3 IT & Educatonal Consultant, Ranch, Jharkhand, Inda 4 Dstngushed Professoral Assocate, Decson Scences and Modellng Program, arxv:1707.00718v1 [cs.ne] 3 Jul 2017 Vctora Unversty, Melbourne, Australa 5 Unversty of Marbor, Faculty of Mechancal Engneerng, Smetanova 17, 2000 Marbor, Slovena 6 Unversty of Marbor, Faculty of Electrcal Engneerng and Computer Scence, Smetanova 17, 2000 Marbor, Slovena (Dated: 1 July 2017) Modelng preference tme n trathlons means predctng the ntermedate tmes of partcular sports dscplnes by a gven overall fnsh tme n a specfc trathlon course for the athlete wth the known personal best result. Ths s a hard task for athletes and sport traners due to a lot of dfferent factors that need to be taken nto account, e.g., athlete s abltes, health, mental preparatons and even ther current sports form. So far, ths process was calculated manually wthout any specfc software tools or usng the artfcal ntellgence. Ths paper presents the new soluton for modelng preference tme n mddle dstance trathlons based on partcle swarm optmzaton algorthm and archve of exstng sports results. Intal results are presented, whch suggest the usefulness of proposed approach, whle remarks for future mprovements and use are also emphaszed.
2 I. INTRODUCTION Mddle dstance trathlon (also called: Half-IRONMAN or IRONMAN 70.3) seems to be one of the most pleasant trathlon dstance. Typcally, t helps athletes as a preparaton for full dstance IRONMAN 7. Its overall dstance s half shorther than IRONMAN and conssts of 1.9 km of swmmng, 90 km of cyclng and 21.1 km of runnng. Usually, tranng for ths race s not as hard as tranng for the Full-IRONMAN. Nowadays, a lot of people are nvolved n these races, because dfferent groups of athletes can compete on the same course at the same tme, although the acheved results are managed separately for professonal and amateur groups of athletes as well as accordng to ther gender. Furthermore, the results of amateur athletes are also classfed nto dfferent age groups. Every year, the number of amateur athletes s ncreasng n these compettons 5. In the past years, tranng for mddle dstance trathlons has progressed sgnfcantly. On the one hand, varous equpments have arsen for mprovng the sports tranng, whle a lot of new IT technologes have also been appled to ths sport. At the moment, almost every athlete s preparng for these races by usng sports watches or smart phones for trackng the sports actvtes. Producers of these wearable moble devces 1 (e.g. Garmn) offer the applcatons for montorng the performed sports actvtes on the Internet that have became more and more advanced every year. Sports traners developed some technologes capable to compare the results of ther tranees wth the opponent s results. Almost all results of sport compettons are now avalable onlne. Therefore, an expert can see from these results, how ther tranees are prepared on the race, whle analyss of some detaled results can also enable them to see pace on partcular segment of race. By the same token, the profle of a race course s also a very nterestng nformaton for the sports traners. Tradtonally, trathlons are organzed consecutvely each year n the same place. Mostly race courses are not changed drastcally, but stays more or less the same. Sometmes only slght changes are needed on the course due to a bad route, reparng etc. As a matter of fact, sports traners and athletes gan a lot of nformaton about race course features from the earler years. In lne wth ths, sports traners are even capable of modelng (or plannng) a partcular fnsh tme for ther athletes n specfc course. Obvously, they can predct ths tme very precsely only when ther tranee adhere ther nstructons. Ths modelng has been done manually n the past. In fact, we dd not know any specal automatc tool for modelng the preference tme n trathlon untl recently. Hence, ths paper proposes a novel soluton based on Partcle Swarm Optmzaton 6 (PSO) that s capable to explore the data from earler edtons of the partcular trathlon competton.
3 In the next sectons, we present the modelng the preference tme n trathlons (secton II). Secton III provdes readers the background nformaton wth emphass on the nature-nspred populaton-based algorthms. In lne wth ths, the PSO algorthm s taken nto a deeper look. The PSO algorthm for modelng preference tme n mddle dstance trathlons s dscussed n Secton IV, whle experments and results are presented n Secton V. The paper s concluded wth Secton VI that summarzes the performed work and outlnes the future research. II. PREFERENCE TIME IN TRIATHLONS Modelng the preference tme n trathlons becomes an mportant task for athletes and ther sports traners. Usng the modelng, the sports traners can predct the fnsh tme on a specfc race course for ther athletes. Mostly, the modelng s performed some weeks or months before the compettons, because the athletes need to be famlar wth the strategy of the racng on a specfc course as proposed by ther sports traner. Modelng of preference tme n trathlon s much harder than n sngle sport dscplnes such as runnng half-marathon or cyclng tme-tral. In fact, trathlon s a mult-sport dscplne that conssts of three sports (.e., swmmng, cyclng, runnng), whle there are also two transtons. In the frst transton, athletes prepare themselves for cyclng after swmmng, whle n the second transton, they fnsh wth cyclng and prepare themselves for runnng. Usually, the transton areas are located at the same place, whle sometmes they are located on the dfferent places. A very mportant aspect of modelng preference tme s not only psycho-physologcal preparaton of the athlete on a race, but also studyng the strategy of racng accordng to course characterstcs. Especally, cyclng course has a huge nfluence on athlete s fnsh tme. Roughly speakng, most of the cyclng courses conssts of uphlls and downhlls. Typcally, some athletes are good n uphll, whle others n downhlls. Addtonally, wnd plays also a very mportant role n plannng the proper strategy for overcomng the cyclng course. As a result, dfferences n cyclng can have the hgher nfluence on the fnsh results than n other two sports. However, these dfferences are much more recognzed by amateur that by the professonal athletes. There s one unwrtten rule suggestng that all sports dscplnes should be fnshed n a lnear correlaton of acheved tmes between dfferent athletes. Ths means that the rato between the tmes acheved by swmmng and cyclng are correlated wth the tmes acheved by cyclng and runnng wthn the same categores of athletes. A motvaton behnd ths s that the racng course should be overcame wth the smlar relatve speed per sports dscplne regardng the other mem-
4 bers of the same group. Consequently, the rank of the specfc athlete remans constantly durng the race wthn the same group of compettors. In ths way, the athletes must acheve the good results n all three dscplnes. Obvously, these assumptons hold perfectly for the short-dstance trathlons, where draftng s permtted. As a result, athlete s fnal result depends strongly on an outcome of the frst dscplne (.e., swmmng). Normally, the lnear correlaton r between two varables x and y s calculated usng the Pearson s correlaton coeffcent that s expressed, as follows: n =1 r = (x x)(y ȳ) n =1 (x x) 2 n =1 (y ȳ), (1) 2 where r [ 1, 1] and r = 1 denotes the overall postve lnear correlaton, r = 1 s the overall negatve lnear correlaton and r = 0 determnes that no correlaton exsts between the observed varables. Table I presents results of top 5 professonal male athletes on IRONMAN 70.3 Tawan 2015. From the frst sght, we can see that the acheved results of swmmng and cyclng are lnearly TABLE I. IRONMAN 70.3 Tawan 2015 - TOP 5 PRO - Male Athlete SWIM [mn] BIKE [mn] RUN [mn] Guy Crawford 24.00 102.63 81.40 Chrstan Kramer 24.02 102.52 83.30 Fredrk Croneborg 24.95 107.38 80.83 Cameron Brown 25.73 110.47 82.77 Paul Ambrose 24.87 107.33 88.57 correlated, because ranks of the compettors after swmmng reman the same also after the cyclng. On the other hand, the results of cyclng and runnng are not too lnear correlated. Ths fndng s proven also by calculatng the Pearson s correlaton coeffcent for SWIM-BIKE r SWIM BIKE = 0.9938 (Fg. 1) and for BIKE-RUN only r BIKE RUN = 0.1804 (Fg. 2). As can be seen from Fgs. 1-2, acheved tmes for SWIM-BIKE correlate perfectly, whle the tmes for BIKE-RUN expose the low correlaton. The man reason for ths behavor s caused by the Fredrc Croneborg that acheved the best results n runnng. However, ths achevement ensured hm the thrd place n the general classfcaton, whle n contrary, the ffth classfed compettor Paul Ambrose lost hs thrd poston after two dscplnes due to the worse runnng. In summary, the man dea behnd the modelng preference tme n mddle dstance trathlons that was extracted from the already mentoned example can be formulated as follows: predctng
5 FIG. 1. IRONMAN 70.3 Tawan 2015-TOP 5 PRO-Male - SWIM-BIKE FIG. 2. IRONMAN 70.3 Tawan 2015-TOP 5 PRO-Male - BIKE-RUN the ntermedate tmes of dfferent sports dscplnes for compettors n mddle dstance trathlons should be based on the premse of keepng ther acheved rankng n the correspondng group constantly durng the whole competton. Ths can be mplemented by usng the mathematcal defnton of correlaton. III. BACKGROUND INFORMATION In ths secton, the background nformaton needed for understandng the subjects that follow are dscussed. Indeed, the followng subjects are dealt wth: nature-nspred algorthms the partcle swarm optmzaton In the remander of ths secton, the mentoned subjects are presented n detals.
6 A. Nature-nspred optmzaton algorthms Nature-nspred optmzaton algorthms are a knd of algorthms that are based on behavour of natural or bologcal systems for solvng NP-hard optmzaton problems. In a nutshell, they are dvded nto two categores 4,8 : Evolutonary Algorthms (EAs), Swarm Intellgence (SI) based algorthms. The former mmcs a Darwnan struggle for survval, where the ftter ndvduals have more chances to survve n the thoughtless crcumstances prevalng n the envronment. On the other hand, the latter, so-called SI-based algorthms, are nspred by socal lvng nsects and anmals, where each ndvdual represents only a smple creature able to perform smple actons. However, they are capable to perform complex tasks, when workng together. The partcle swarm optmzaton belongs to a famly of SI-based algorthms and s based on the socal behavor of brds. In the study, ths algorthm was appled for modelng preference tme n mddle dstance trathlons. B. Partcle swarm optmzaton Partcle swarm optmzaton (PSO) 6,9 belongs to a famly of stochastc populaton-based algorthms wth a populaton consstng of real-valued vectors (also partcles) x (t) = (x (t),1,..., x(t),d )T for = 1,..., Np and j = 1,..., D, where Np denotes the populaton sze and D s a dmenson of the problem. Thus, the algorthm explores new solutons by movng the partcles throughout a search space n drecton of the current best soluton. Indeed, two populatons of partcles are managed by the algorthm,.e., besde the populaton of partcles x (t) best solutons p (t) generaton 2. The new partcle poston s generated as presented n Eq. (2). also the populaton of local. Addtonally, the best soluton n the populaton g (t) s mantaned n each v (t+1) x (t+1) = v (t) + C 1 U(0, 1)(p (t) = x (t) + v (t+1), x (t) ) + C 2 U(0, 1)(g (t) x (t) ), (2) where U(0, 1) denotes a random value drawn from the unform dstrbuton n nterval [0, 1], and C 1 and C 2 are learnng factors. The pseudo-code of the orgnal PSO algorthm s llustrated n Algorthm 1.
7 Algorthm 1 Pseudocode of the basc partcle swarm optmzaton algorthm Input: PSO populaton of partcles x = (x 1,..., x D ) T for = 1... Np, MAX F E. Output: The best soluton x best and ts correspondng value f mn = mn(f(x)). 1: nt partcles; 2: eval = 0; 3: whle termnaton condton not met do 4: for = 1 to Np do 5: f = evaluate the new soluton(x ); 6: eval = eval + 1; 7: f f pbest then 8: p = x ; pbest = f ; // save the local best soluton 9: end f 10: f f f mn then 11: x best = x ; f mn = f ; // save the global best soluton 12: end f 13: x = generate new soluton(x ); 14: end for 15: end whle The PSO algorthm has also been used for generaton of eatng plans for athletes 2. IV. PSO FOR MODELING PREFERENCE TIME The basc PSO algorthm can not be used drectly for modelng preference tme n trathlons. In order to ft t to ths problem, the basc PSO needs to be modfed. There are three mportant modfcatons for mplementaton, as follows: representaton of ndvduals, ntalzaton of solutons, and ftness functon defnton. Indvduals n the PSO algorthm for modelng preference tme n trathlons are represented as a real-valued vectors x wth dmenson of D = 5 denotng the predcted tmes ntermedate of the specfc dscplnes for an observed athlete. Thus, each partcle conssts of fve elements representng tmes of swmmng, transton 1, cyclng, transton 2 and runnng, respectvely.
8 TABLE II. Constrant s handlng Dscplnes Constrants Swm t SWIM [25.0, 50.0] Transton T1 t T1 [2.0, 5.0] Cyclng t BIKE [140.0, 180.0] Transton T2 t T2 [2.0, 5.0] Runnng t RUN [85.0, 120.0] Intalzaton of the partcle elements s performed randomly. However, each element must satsfy the followng boundary constrants as presented n Table II. These constrants refer to a regular amateur athletes and can capture the slghtly hgher ntervals by the professonals. In common, these settngs must be talored regardng the observed athlete. Addtonally, constrants must not be volated also by movng the partcles through the search space accordng to Eq. (2). Ftness functon that evaluates the qualty of solutons s expressed as follows. At frst, the overall tme of the trathlon s calculated accordng to the followng equaton: D h(x ) = x,j, (3) where the elements of the partcles are smple summed. Then, the maxmum mprovement tme s j=1 expressed as follows: K max h(x ), K max < h(x ), g(x ) = MAX TIME, otherwse, (4) where K max denotes for the 5 % of mproved personal best fnsh tme acheved by the observed athlete. The motvaton behnd the functon g(x ) s to search for solutons that are the best ftted to the maxmum mprovement tme K max. Thus, solutons wth the overall fnsh tme better than the K max are gnored by assgnng t the constant tme MAX TIME denotng the worst overall fnsh tme. Fnally, the ftness functon f(x ) s defned as follows: h(x ), r 1 < r 2, f(x ) = MAX TIME, otherwse, (5) where r 1 and r 2 denote the Pearson correlaton coeffcents of an archve of the results for the professonal/age group of athletes Ar (n) and Ar (n+1) expressed as: r 1 = r SWIM BIKE (Ar (n) ) + r BIKE RUN (Ar (n) ), (6)
9 TABLE III. Parameter settng for the PSO algorthm Parameter Value NP 50 maxf es 10000 D 5 C 1, C 2 2.0 current best max 300.0 and r 2 the r 2 = r SWIM BIKE (Ar (n+1) ) + r BIKE RUN (Ar (n+1) ). (7) Varable Ar (n) n Eq. (6) denotes an archve of the best results as acheved by the best n-compettors n an observed professonal/age group n some competton. On the other hand, varable Ar (n+1) n Eq. (7) determnes the archve Ar (n), to whch the predcted results x are added. The dea behnd the calculaton of the ftness functon s searchng for a soluton that s the most ftted to correlaton of the result archve. In lne wth ths, the observed athletes would keep ther predcted ntermedate results wth the results acheved by the group of athletes, belongng to the specfc professonal/age group. V. EXPERIMENTS The purpose of experments was to valdate the PSO algorthm for modelng preference tme. In lne wth ths, experments were conducted on two dfferent archves of the exstng results, where the predcted ntermedate achevements must correlate wth the results acheved by the athletes n the correspondng professonal/age groups. The PSO algorthm for modelng preference tme was mplemented n Python programmng language. The task of the algorthm was to predct the ntermedate results for a specfc athlete such that overall fnal tme of 300 mn (or 5 hours) s acheved for someone, whch has the personal best result around 3 mn above the predcted tme (t best > 300 mn) and would lke to acheve the new personal best fnsh tme slghtly under 5 hours (t new best 300 mn). Parameter settng of the PSO as presented n Table III were used durng the expermental work. Thus, the 5 ndependent runs were conducted n order to evaluate the results of the stochastc SI-based algorthm. In ths study, result archves have been employed that were acheved n the prevous IRONMAN compettons worlwde. Precsely, the results of athletes tracked durng IRONMAN 70.3 races n
10 2015 from ronman.com webstes were used that are publcly avalable as datasets 3. In these datasets, data are saved n CSV and JSON fle formats and can be used n the proposed PSO algorthm wthout any specal preprocessng except all tmes represented as character strngs must be converted to mnutes (floatng-pont values). Indeed, these data consst of: the name of an athlete, the naton, the category, the fnsh place, splt tmes of swmmng, bkng and runnng, splt tmes of both transtons, and the overall tme. Although the mentoned datasets nclude the dfferent sets of results, here, we were focused especally on two IRONMAN compettons,.e., IRONMAN 70.3 Taupo 2015 and Wesbaden 2015. Indeed, the results of the frst 30 compettors n age group 25-29 years represent a bass, from whch the archves Ar (n) and Ar (n+1) were bult. A. Results The results of predcted ntermedate tmes for the observed 26 year old man athlete were based on archves of results acheved on two dfferent mddle dstance trathlon courses,.e., IRONMAN 70.3 Wesbaden 2015 and Taupo 2015. Each archve conssts of 30 best results from the age group 25-29 years. Interestngly, the calculated Pearson s correlaton coeffcent for the frst competton (W esbaden) (T aupo) was r 1 = 0.7256, whle for the second r 1 = 0.2051. The results of experments are depcted n Tables IV and V, where columns represent the predcted ntermedate tmes for three sports dscplnes as well as two transton areas. Furthermore, the predcted overall fnsh tme column s added to both tables. Thus, the frst fve rows llustrate the results for each ndependent runs, whle the last two rows denote the mean and standard devaton of the results n aforementoned rows, respectvely. Mostly, the dfference between the Pearson s correlaton coeffcents are reflected n the correspondng standard devaton values depcted n the tables. Ths means that predcted ntermedate
11 TABLE IV. Predcted ntermedate tmes based on results acheved n IRONMAN 70.3 Wesbaden Run Swmmng T1 Cyclng T2 Runnng Total 1 32:45.65 3:25.88 2:45:11.42 4:35.33 1:34:01.53 4:59:59.82 2 34:11.84 2:00.00 2:49:04.71 4:41.07 1:30:02.36 4:59:59.99 3 35:47.80 3:13.38 2:50:35.38 3:22.77 1:27:00.64 4:59:59.99 4 31:23.28 2:14.76 2:44:14.28 3:46.04 1:38:21.47 4:59:59.86 5 35:07.18 3:10.33 2:48:43.17 2:36.88 1:30:22.41 5:00:00.00 Mean 33:51.15 2:48.87 2:47:33.79 3:48.42 1:31:57.68 4:59:59.93 Stdev 1:47.20 0:38.67 2:42.88 0:51.84 4:21.36 0:00.08 TABLE V. Predcted ntermedate tmes based on results acheved n IRONMAN 70.3 Taupo Run Swmmng T1 Cyclng T2 Runnng Total 1 41:02.92 2:41.81 2:47:35.78 3:39.45 1:25:00.00 4:59:59.98 2 33:59.64 4:02.92 2:37:53.48 2:46.39 1:41:17.54 4:59:59.99 3 38:11.37 2:00.00 2:43:28.88 5:00.00 1:31:19.73 4:59:59.98 4 37:59.19 3:31.27 2:43:33.20 2:40.58 1:32:15.73 4:59:59.99 5 39:37.81 3:04.22 2:43:11.92 4:35.71 1:29:30.30 4:59:59.98 Mean 38:10.19 3:04.04 2:43:08.65 3:44.43 1:31:52.66 4:59:59.99 Stdev 2:38.43 0:46.91 3:27.19 1:02.86 5:57.51 0:00.00 tmes n dfferent runs of the PSO algorthm dstngush lesser n Table IV than n Table V, because the Pearson s coeffcent underlyng the results n the frst table s much hgher than the same coeffcent underlyng the results n the second table. For nstance, predcted values for IRONMAN 70.3 Wesbaden are for swmmng wthn 31-35 mnutes (4 mnutes), cyclng wthn 2:44-2:50 hours (6 mnutes) and runnng wthn 1:27-1:38 hours (11 mnutes). In contrary, the values for the IRONMAN 70.3 Taupo are more dspersed, because these are stuated for swmmng wthn 33-41 mnutes (8 mnutes), for cyclng wthn 2:37-2:47 hours (10 mnutes) and runnng wthn 1:25-1:41 hour (16 mnutes). In general, the Pearson s correlaton coeffcent has a bg nfluence on modelng the preference tme n trathlons. Analyss of the coeffcent s especally nterestng for the professonal athlete s groups, where these values are much hgher. Consequently, the predcton of the overall fnsh tme n mddle trathlons s much easer n these groups, because the rght champons must be excellent n all dscplnes.
12 B. Lessons learned Before the experments were fnshed, we dd not have any expectatons only a feelng about how postve effects should have an usage of the PSO algorthm wth Pearson s correlaton coeffcent appled to the real data. However, experments revealed an nterestng story. They confrmed, on one hand that t s possble to model the preference tme n mddle trathlons, whle on the other hand that ths predcton can be very relable. From the results of experments, we can nfer that the more than the data are correlated, the more relable predcton can be expected. For nstance, the frst 30 compettors n age group 25-29 years at IRONMAN 70.3 Wesbaden 2015 have ther tmes much more correlated between dscplnes than the frst 30 compettors at IRONMAN 70.3 Taupo 2015. Hence, predcted results for the race n Wesbaden s much more relable than the second one. Addtonally, the observed athlete can also select, whch of the offered preference tmes wll be hs target tme on the same race course. However, the relablty of the preference tmes can addtonally be mproved, when the personal best tme of the observed athletes wll be taken nto account. In lne wth ths, the frst 15 results better than ther personal best and the frst 15 results worse than ther personal best can be selected from the archve n place to take the best 30 results. VI. CONCLUSION In ths paper, we have proposed a new soluton for modelng preference tme n mddle dstance trathlon usng the PSO algorthm. In lne wth ths, we have systematcally outlned a problem and proposed a theoretcal model to solve t. The theoretcal model was later valdated by some experments, where real data from IRONMAN 70.3 compettons were taken nto consderaton. There s also a bunch of possbltes for future extenson of the proposed soluton. On the other hand, the nfluence of data sze should be studed, whlst the desgn of other ftness functons would also be nterestng. 1 Andrew Best and Barry Braun. Usng a novel data resource to explore heart rate durng mountan and road runnng. Physologcal Reports, 5(8):e13256, 2017. 2 Dušan Fster, Samo Rauter, Iztok Fster, and Iztok Fster Jr. Generatng eatng plans for athletes usng the partcle swarm optmzaton. In 17th Internatonal Symposum on Computatonal Intellgence and Informatcs (CINTI), pages 193 198, 2016.
13 3 Iztok Fster Jr. and Dušan Fster. A collecton of ronman, ronman 70.3 and ultra-trathlon race results, verson 0.1, techncal report 0110, 2016. 4 Iztok Fster Jr., Xn-She Yang, Iztok Fster, Janez Brest, and Dušan Fster. A bref revew of naturenspred algorthms for optmzaton. Elektrotehnšk vestnk, 80(3):116 122, 2013. 5 Anna M Kaleta, Ewa Lewcka, Alcja Dabrowska-Kugacka, Zuzanna Lewcka-Potocka, Elżbeta Wabch, Wojcech Potock, and Grzegorz Raczak. Intensve exercse and ts effect on the heart: Is more always better? Cardology Journal, 24(2):111 116, 2017. 6 J Kennedy and R Eberhart. Partcle swarm optmzaton. In Neural Networks, 1995. Proceedngs., IEEE Internatonal Conference on, volume 4, pages 1942 1948. IEEE, 1995. 7 Beat Knechtle, Pantels T Nkolads, Thomas Rosemann, and Chrstoph A Rüst. Der ronman-trathlon. Praxs (16618157), 105(13), 2016. 8 S Salcedo-Sanz. Modern meta-heurstcs based on nonlnear physcs processes: A revew of models and desgn procedures. Physcs Reports, 655:1 70, 2016. 9 Yudong Zhang, Shuhua Wang, and Genln J. A comprehensve survey on partcle swarm optmzaton algorthm and ts applcatons. Mathematcal Problems n Engneerng, 2015, 2015.