tag:blogger.com,1999:blog-34204451214858230342024-03-13T23:38:48.279-07:00You Guys are Tweakers.A cycling blog that was taken over by math.Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.comBlogger27125tag:blogger.com,1999:blog-3420445121485823034.post-8757470455085393592009-07-30T16:21:00.000-07:002009-07-31T15:24:34.568-07:00Domestiques: Who Needs Them?In previous posts <a href="http://tweakersonwheels.blogspot.com/2009/07/best-domestiques-podium-edition.html?pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print">here</a> and <a href="http://tweakersonwheels.blogspot.com/2009/07/andera-tonti-best-domestique-ever.html?pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print">here</a> I looked at some ways to measure the relative value of domestiques. We can also turn this question around and consider which leading riders are most dependent on having certain teammates around to win. Over his career, for instance, has Alessandro Petacchi needed certain lead out men around to win? Have grand tour winners relied on help from specific teammates?<br /><br />The simplest way I can think to do this is to do a <a href="http://en.wikipedia.org/wiki/Linear_regression">linear regression</a> on the team leader’s results. Linear regression finds the model that best fits all of the leader’s results in terms of a linear combination of parameters from each teammate. These parameters are fit from results data. As usual, I will take the logarithm of all results to put more value on higher placings. Essentially, each result is represented by “adding up” all of the contributions from the teammates that were in that race. This is formally defined for each race as:<br /><br />Log(<i>R</i>) = <i>β<sub>0</sub> + β<sub>1</sub> x<sub>1</sub> + β<sub>2</sub>x<sub>2</sub> + β<sub>3</sub>x<sub>3</sub> …<br /><br />R </i>is the leader’s result in the race and the<i> x<sub>i</sub></i> correspond to each teammate the leader has ever raced with. If teammate 1 was present in the race, <i>x<sub>1</sub></i> is 1. If not, <i>x<sub>1</sub></i> is 0. Writing this equation for every individual race, we get a big series of algebraic equations in which we know all the results and all the <i>x<sub>i</sub></i>. We then find the best fit for each of the regression coefficients <i>β<sub>i</sub></i>, which correspond to how much each teammate contributes. Helpful teammates will have negative <i>β<sub>i</sub></i> since they reduce the result. Teammates with a positive <i>β<sub>i</sub></i> tend to make the leader’s results worse.<i><br /></i><br />The coefficient <i>β<sub>0</sub></i>, called the intercept, can be thought of as the leader’s base result before teammates get factored in. It is the same for every race, as chosen to best fit all races. A rider with a large intercept relies on specific teammates to bring his result down to a top placing, whereas a rider with an intercept near zero generally does well regardless of which teammates are present. Note that these calculations depend on having a lot of results with a variety of teammates in order to tease out the contributions from each domestique.<br /><br />I identified 33 riders with at least 15 podiums from <a href="http://www.cqranking.com/men/asp/gen/start.asp">Cycling Quotient</a> and performed regression on their career results (excluding ITT and when riding for national teams). Here are the riders, number of races I used in calculations, their intercept, and their most valuable domestique (minimum 50 races together):<br /><br /><a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhMWalSDd7lj-Se4s8zJddEBOKXH-Y7nCcDbYfz8LvpwGmrxLAKOuHwch72dVlGXkxEekIkEIb1na7xpbpRf4iwCHIuVRLK4vbzq1v4JRqJFt5IHVYPKeqvnoswI_fgPL-Eok6ncHcayItO/s1600-h/topstab.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 348px; height: 400px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhMWalSDd7lj-Se4s8zJddEBOKXH-Y7nCcDbYfz8LvpwGmrxLAKOuHwch72dVlGXkxEekIkEIb1na7xpbpRf4iwCHIuVRLK4vbzq1v4JRqJFt5IHVYPKeqvnoswI_fgPL-Eok6ncHcayItO/s400/topstab.png" alt="" id="BLOGGER_PHOTO_ID_5364400825224114578" border="0" /></a>The precise value of the intercept isn't very meaningful in itself since I fit log-transformed results, but for reference a value of zero would mean the rider always wins independent of their teammates. So the riders at the top of this list have needed less help to get their results, in the sense that they do well regardless of who they’re racing with. Riders at the bottom are those whose frequency of a good result is dependent on having certain teammates present. Some interesting points:<br /><ul><li>Greipel and McEwen are the sprinters who are high on the list. Greipel has been successful on a team that primarily supports other sprinters – both Cavendish and Henderson have large and positive regression coefficients, meaning they systematically harm Greipel’s results when they’re around. McEwen has made his career winning grand tour stages for teams busy supporting a GC contender. Regression correctly identifies these guys as riders who don’t rely on team support.</li><br /><li>Similarly, Kirchen and Pellizotti are GC riders who have never had the benefit of a dedicated support team in grand tours. They race relatively independently, as the analysis shows.</li><br /><li>Recent Astana drama aside, I interpret Contador’s place high on this list to mean he is strong enough to win regardless of who happens to be in the same kit. So we shouldn’t doubt Contador’s grand tour chances in the future, no matter where he ends up in 2009.</li><br /><li>The riders in the middle (Boonen, Bettini, Menchov, etc.) are all good candidates for guys whose support has depended on age and the specific race. As these riders became more experienced and targeted their races their team support increased, but many of their early results were achieved without a team built around them.</li><br /><li>Intercepts of 0.6-1.5 are most common, suggesting that it is standard for both GC riders and stage hunters to rely quite a bit on their teammates. This is not surprising.<br /></li><br /><li>It appears that O’Grady needs specific teammates present in order to do well. I suspect his intercept is so extreme because the majority of his podiums are from grand tours that he has ridden with a core set of teammates, so the regression associates them with his success. This could be coincidence, but we can’t say for sure.</li><br /><li>Many of the leader-domestique pairings are very sensible. Contador with Paulinho, Petacchi with Velo, and Armstrong with Rubiera are just a few of the well-known combinations that appear on the list. </li><br /></ul>There are a lot of other ways to judge how much a team leader relies on teammates for success, so I will probably try other methods in the future. Regression, however, is a relatively simple and common way to address problems like this so I decided to start with it.<br /><br /><span style="font-size:85%;"><br />Technical notes: Data source is Cycling Quotient and includes about 900 races from 2002 to the present. To avoid partial result listings, I considered approximately 900 races in which more than 100 riders are listed in the results. Roughly half of the races were stages from grand tours, and the remaining results are mostly the major one-day races and lesser stage races. Individual time trials and national team events were excluded from the analysis. Regressions were fit using ordinary least squares.</span><i><span style="font-style: italic;"><br /></span></i><i><br /></i>Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-75237586905374945512009-07-28T21:28:00.000-07:002010-07-23T13:28:22.900-07:00Best Domestiques (Podium Edition)I recently posted <a href="http://tweakersonwheels.blogspot.com/2009/07/andera-tonti-best-domestique-ever.html?pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print">a statistical analysis</a> that identified domestiques who are associated with better team results. For example, I found that when Quick Step started Kevin Hulsmans in a race last year, their best finish was an average of 10 places higher than when they did not. So you might say Hulsmans was worth 10 places to his best Quick Step teammate. I also calculated how significant the effects were in terms of the statistical likelihood that such an effect might be a random fluctuation. In doing this, I used log-transformed results to put more weight on better placings. This basically made the difference between 1st and 10th as important as the difference between 10th and 100th. Although this approach is fine for some purposes, I think it still underestimates the importance of a top finish.<br /><br />This post will propose an alternate method that focuses on podium placings. In a bike race, top 10 results are satisfying only in that they suggest the potential for a 1st, 2nd, or 3rd place finish down the line. So here I will ask if certain riders increase the frequency of their team achieving a podium position. As with the previous analysis I will do this on both a year-by-year and career basis, including races between 2002 and 2009.<br /><br />As an example, consider Marco Velo. From 2002 to 2008 Velo was a leadout man for Alessandro Petacchi, one of the era's dominant field sprinters, and now performs similar duties on Quick Step. Over that time, Velo has appeared in the results of 281 races, 69 of which have a teammate on the podium (not Velo). His teams also raced 461 times without him, with 57 podiums. So Velo's team has achieved more podiums in far fewer races Velo contested: 69/281 versus 57/461. This corresponds to an odds ratio of 2.3, meaning that it was 2.3 times more likely that Velo's team made the podium when he was in the race. That sounds pretty good, right? But, of course, you also want to know if this a significant difference given these sample sizes. We can use a statistical test to determine that the likelihood of this effect in random data is <span style="font-style: italic;">P</span> = 2e-5, or 0.0002%. Quite significant, suggesting that Marco Velo is an excellent domestique. Good for him.<br /><br />Using <a href="http://www.cqranking.com/men/asp/gen/start.asp">CQ data</a> for all riders (see the riveting technical notes below and on previous posts for more details), I went searching for other extraordinarily valuable domestiques. I identified every rider/year combination with a <span style="font-style: italic;">P</span> less than 0.01. Here are the rider, year, team, odds ratio, <span style="font-style: italic;">P</span> value, and most common teammate on the podium for each significant finding:<br /><br /><a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjv_ruKdcVouO6s2vThkDLZ-AG5wnVxEPUAKm8Sm-lNb6CHoCiaZeWV-HARh9pUuEZ4GTsDaNSiHENTZpdeDr9jezy5m5a4Z8H591wqO_FccmGNN3C1M5qAbDgaf01PD-J7EJbZttZM-eHj/s1600-h/byyear.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 360px; height: 400px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjv_ruKdcVouO6s2vThkDLZ-AG5wnVxEPUAKm8Sm-lNb6CHoCiaZeWV-HARh9pUuEZ4GTsDaNSiHENTZpdeDr9jezy5m5a4Z8H591wqO_FccmGNN3C1M5qAbDgaf01PD-J7EJbZttZM-eHj/s400/byyear.png" alt="" id="BLOGGER_PHOTO_ID_5368844971029588226" border="0" /></a>The odds ratio is how many times more likely it is that a teammate gets on the podium when the listed rider is racing (larger is better). Infinite results (INF) occur when the team never placed on the podium without the rider present. The <span style="font-style: italic;">P</span> value is the chances that this result might have arisen from random noise (lower is better). I also did the same calculation for each rider's career -- at least using the results I have from 2002-2009:<br /><br /><a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgCp_wa9V2gjkXEt6_ujheftkOjygDPR-g_bbd54ay6-tnKJXHJ_3yTdD10p8vvZO7NjIv4pVszSz_a_Nt3etZlHwyZt3tjpiG6dHD7N-1QMv-O0dQufsk6McJzBNS7wOfG-K0ZMo5-RYJq/s1600-h/allyears.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 400px; height: 163px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgCp_wa9V2gjkXEt6_ujheftkOjygDPR-g_bbd54ay6-tnKJXHJ_3yTdD10p8vvZO7NjIv4pVszSz_a_Nt3etZlHwyZt3tjpiG6dHD7N-1QMv-O0dQufsk6McJzBNS7wOfG-K0ZMo5-RYJq/s400/allyears.png" alt="" id="BLOGGER_PHOTO_ID_5368851453738140082" border="0" /></a>We can compare these two tables with the previous results and see that there is a fair amount of overlap. For instance, the 2008 season for Kevin Hulsmans is still significant, but now instead of saying he's worth 10 placings we can credit him with a three-fold increase in podium spots. Notably missing is the 2003 incarnation of Andrea Tonti, whom I previously declared to be the best domestique ever. Although his presence corresponded to an astounding gain of 33 placings, he wasn't around for enough teammates' podiums to make this list. So he might be an example of moving teammates into the top 10, but not all the way to the big money.<br /><br />As before, I'm not implying that a domestique whose specific presence doesn't yield enhanced podium returns isn't doing his job well. He might be on a team that is always putting riders on the podium, or a team with second-rate team leaders who rarely crack the top three. Basically all I'm doing here is identifying domestiques who have shown a pattern of association with good team results. Determining whether the domestique is actually <span style="font-style: italic;">causing </span>the better results is a judgment call that the statistics cannot make.<br /><br />I prefer this method to my previous one, primarily because it's easier to understand and focuses better on top results. However, this it's bedeviled by some of the same issues. A couple major ones are:<br /><ul><li>False positives. The significance levels appear to be quite low, but since I've done thousands of tests there may be many false positives here. However, I'm not sure how independent these tests are so I can't easily compute a correction. I would have to do a large number of permutation tests to get an empirical idea of the precise false positive rate. </li><br /><li>Disregarded cofactors. As we all know, correlation does not necessarily imply causation. An analysis like this may be fraught with causal variables that have been ignored in the analysis. For example, it is difficult to separate the contribution of one domestique from another, and from that of the team leader. Many of the riders on the list are Alessandro Petacchi's leadout train (Velo, Ongarato, Tosatto). Were these guys extraordinarily suited to leading out their man, did one of them carry the weight for them all, or were they just lucky to be working for the fastest guy around? It might be impossible to separate the contributions of Petacchi and his leadout train with the results I have, but it's worth thinking about. This analysis doesn't really try. Another cofactor is the nature of specific event. Since pack finishes are so common, domestiques that aid in sprints will have more significant results due to the greater sample size of sprints.</li></ul><br /><span style="font-size:85%;">Technical Notes: Data source is Cycling Quotient. To avoid partial result listings, I considered approximately 900 races in which more than 100 riders are listed in the results. Roughly half of the races were stages from grand tours, and the remaining results are mostly the major one-day races and lesser stage races. Individual time trials and national team events were excluded from the analysis. To avoid small sample sizes, odds ratios and <span style="font-style: italic;">P</span> values were only computed if there were five results in every test set. <a href="http://en.wikipedia.org/wiki/Odds_ratio">The odds ratio</a> is defined as <span style="font-style: italic;">p_r</span>(1-<span style="font-style: italic;">p_r</span>)/<span style="font-style: italic;">p_nr</span>(1-<span style="font-style: italic;">p_nr</span>), where <span style="font-style: italic;">p_r</span> and <span style="font-style: italic;">p_nr</span> are the frequencies of a team podium place when the rider is and is not in the race, respectively. <span style="font-style: italic;">P</span> values are calculated using <a href="http://en.wikipedia.org/wiki/Fisher%27s_exact_test">Fisher's exact test</a>, which assumes a hypergeometric distribution for the null hypothesis.<br /><br /><br /></span>Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-26720624299339954852009-07-27T17:23:00.000-07:002009-07-27T17:34:42.314-07:00UPDATED: Should Unprecedented Tour Success Warrant Suspicion?<a href="http://tweakersonwheels.blogspot.com/2009/07/should-unprecedented-tour-success.html?pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print">Last week</a> I considered whether a big improvement in a rider's results is reasonable grounds for suspicion of doping. Bradley Wiggins is the current poster boy for this sort of skepticism, having just finished fourth in the Tour de France. After last year's Tour, the news that Bernhard Kohl and Stefan Schumacher had been caught using CERA was easy to believe due to the sense that their performances had improved to an extent that wasn't natural.<br /><br />However, before jumping to conclusions about Wiggins we should ask whether dopers really show significantly improved results and, if so, whether such improvements also occur for non-dopers. This can be answered with some fairly straightforward statistical testing to compare a rider's current results with their previous results. I previously defined two parameters to quantify the improvement and significance of a given rider's results during a given year:<br /><ul><li><span style="font-style: italic;">R</span>: The difference between a rider’s mean placings in a given year and mean placings in previous years. Larger numbers mean greater improvement.</li><br /><li><span style="font-style: italic;">P</span>: The likelihood that this difference is real and not simply a result of random fluctuations (statistical significance). Smaller numbers mean greater significance.</li></ul>When computing <span style="font-style: italic;">P</span>, I first take the logarithm of all results in order to enhance the value of top placings and minimize differences between mid-pack finishes (e.g. the difference between 2nd and 12th is much more important than the difference between 102th and 112nd). See below for the enthralling technical details.<br /><br />I computed <span style="font-style: italic;">R</span> and<span style="font-style: italic;"> P</span> for 934 riders over the years 2003-2009. For each rider, I only considered years in which 10 or more results were in the Cycling Quotient database. Forty-five rider/year pairings showed statistically significant improvements (see technical note for the definition of "significant"). Here they are, with 2009 cases in red:<br /><br /><a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhIXY4Fb9ZTCZIJ8MMcuPdc8-wAeSOrLrkmXEmsRoX0fYMuAfdcN8Mk5fxdGNhsq9lQgt-9m3-iOMnHUi-mgwOzxcRTS_3fYgUIdCiHJNamXtocyzIBliQ0XenYxKViSMKNjj8AeM9UWc7q/s1600-h/wiggins_tabs_2a.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 154px; height: 320px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhIXY4Fb9ZTCZIJ8MMcuPdc8-wAeSOrLrkmXEmsRoX0fYMuAfdcN8Mk5fxdGNhsq9lQgt-9m3-iOMnHUi-mgwOzxcRTS_3fYgUIdCiHJNamXtocyzIBliQ0XenYxKViSMKNjj8AeM9UWc7q/s320/wiggins_tabs_2a.png" alt="" id="BLOGGER_PHOTO_ID_5363301496295537938" border="0" /></a>Having improved an average of 50 places per race, Wiggins's 2009 is on this list. Columbia's Tony Martin is here as well, and has actually gained more than Wiggins this year. But there are very few convicted dopers here; where are our naughty friends? Adding Di Luca's 2009 to the list I showed on my previous post, the results for recent doping positives look like this:<br /><br /><a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjSacNjh9_xesLskWwMtOikTTiNFjCrsRm9NKK1deF7MV-3rCaHLAOnx7cAt0PDdQzYzHxkFYhXtSIBeFclGuH-YG12rVnnU5_9KkkrwJRxgAnUb9xnm1xeFYlSdvu8LJ5K9aqHtv3VtJhj/s1600-h/wiggins_tabs_2c.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 320px; height: 232px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjSacNjh9_xesLskWwMtOikTTiNFjCrsRm9NKK1deF7MV-3rCaHLAOnx7cAt0PDdQzYzHxkFYhXtSIBeFclGuH-YG12rVnnU5_9KkkrwJRxgAnUb9xnm1xeFYlSdvu8LJ5K9aqHtv3VtJhj/s320/wiggins_tabs_2c.png" alt="" id="BLOGGER_PHOTO_ID_5363301621853256882" border="0" /></a>Although a few of these riders show large gains in average results and fairly low <span style="font-style: italic;">P</span> values, none of these riders appear on the above list of significant cases. So I'd have to say there isn't much support to the idea that a big improvement in results is a sign of doping.<br /><br /><br /><span style="font-size:85%;">Technical Notes: I defined significance as having a </span><span style="font-style: italic;font-size:85%;" >P</span><span style="font-size:85%;"> less than 2e-5. This might sound overly conservative -- it means the chances that the rider's current and previous results are the same is only 0.002%. The problem is that I've done 2400 tests, so a cutoff of 5% would give me over 100 false positives. Dividing 0.05 by the number of tests, I get a </span><span style="font-style: italic;font-size:85%;" >P</span><span style="font-size:85%;"> cutoff of 2e-5 and don't need to worry about false positives. Additional technical notes of possible relevance <a href="http://tweakersonwheels.blogspot.com/2009/07/should-unprecedented-tour-success.html?pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print">here</a> and <a href="http://tweakersonwheels.blogspot.com/2009/07/andera-tonti-best-domestique-ever.html?pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print">here </a>(scroll down to the fine print).</span>Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-34796042024812703302009-07-25T19:45:00.000-07:002009-07-25T19:46:51.881-07:00TdF Stage 20I think it's fair to say that Armstrong kept his word on what he would do to last year's top 5. Will he now demand that Sastre apologize for making him apologize during the first week? Oh the drama.Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-48729288994888125812009-07-24T17:11:00.000-07:002009-07-27T10:09:26.996-07:00Andrea Tonti: Best Domestique Ever?Now that Mr Lance Armstrong is embracing his role as a domestique (well, sort of), I expect domestiequery to become the hottest trend in cycling chatter. Internet forum people will now get into heated arguments over whether Armstrong is clearly the most awesome domestique ever to ride, or obviously the embodiment of all that is wrong with cycling teamwork.<br /><br />However, I’m afraid this argument will be even more fruitless than those surrounding the achievements of team leaders. Leaders, after all, win races, which is relatively easy to remember. How often can one recall, much less judge, the efforts of domestiques? Sure, most of us immediately conjure up images of Jens Voigt hammering for his CSC/Saxo leaders, Yaroslav Popovych riding at the limit for Armstrong (but not Cadel Evans), and Johan Van Summeren spending tens of kilometers chasing down breaks for Robbie McEwen. But how essential were those efforts in the end? And, when it comes to grinding out the kilometers, how interchangeable are these guys?<br /><br />Ultimately I think a good domestique is a rider who makes their teammates better riders. And by better riders, I mean they achieve better results. It doesn’t matter what the domestique does – whether he chuffs out 60 km to chase down the break or takes a fall in the final kilometer to let his man escape, the only criterion is his teammates’ results.<br /><br />Best of all, a results-based approach enables a quantitative method. First, I take all races that a domestique’s team contested in a given year. At this point, any rider is a potential domestique so I do this for everyone. I then separate this set of races into two groups. The first group is races the domestique finished, and the second are races he did not finish. Ideally this would be based on races the domestique did or didn’t start, but I don’t have that data. I then took the result from each of these races achieved by the best-placed teammate. This gives me two sets of race results, corresponding to the team’s top finisher (other than the rider) in each event the domestique did or did not finish. If the rider’s teammates have significantly better results when he is present, that rider is a particularly valuable domestique.<br /><br />I calculate two quantities from these two sets of results:<br /><br /><ul><li><span style="font-style: italic;">D</span>: The difference between the team’s average best placing when the rider is present and the team’s average best placing when the rider is absent. Positive values mean the rider is a good domestique, negative numbers suggest he is not. <span style="font-style: italic;">D</span> is for the domestique value (it’s even the same in French!).</li><br /><li><span style="font-style: italic;">P</span>: The likelihood that this difference is real and not simply a result of random fluctuations (statistical significance). Smaller numbers mean greater significance.</li></ul><br />Data, as usual, are from the fantastic <a href="http://www.cqranking.com/men/asp/gen/start.asp">Cycling Quotient</a>. When computing <span style="font-style: italic;">P</span>, I first took the logarithm of all results in order to enhance the value of top placings and minimize differences between mid-pack finishes (e.g. the difference between 2nd and 12th is much more important than the difference between 102th and 112nd). See below for more enthralling technical notes.<br /><br />I calculated <span style="font-style: italic;">D</span> and <span style="font-style: italic;">P</span> for riders in the pro peloton over about 900 races from 2001-2009 (see those technical notes again). I computed this on both a year-by-year basis, since certain riders may be better domestiques on certain teams (looking at you, Popo), as well over an entire career (to the extent I have their results). Here are the top domestiques for 2002-2009, using a significance cutoff of 1E-4:<br /><br /><a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjSFNi2FkuKM4q0Bu9J-pNkYOdH9hXTtSb8Gl8NY34M26M22_J2-vS5BUeVTYW6dzzAsoOHZ6X-AjHxWMmdQV4LofIua99n9eIhxOb1674b3YgDe2teaPBsfUfQCOuKCgRhfranwtoGVGlg/s1600-h/bestdomestiques.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 320px; height: 295px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjSFNi2FkuKM4q0Bu9J-pNkYOdH9hXTtSb8Gl8NY34M26M22_J2-vS5BUeVTYW6dzzAsoOHZ6X-AjHxWMmdQV4LofIua99n9eIhxOb1674b3YgDe2teaPBsfUfQCOuKCgRhfranwtoGVGlg/s320/bestdomestiques.png" alt="" id="BLOGGER_PHOTO_ID_5362184063456693106" border="0" /></a>And the prize goes to Italian Andrea Tonti, whose 2003 season was the best performance as a domestique in the data I have. Tonti worked for Gilberto Simoni in his Giro win that season, and presumably when Tonti was not around Saeco did not enjoy the same success. Career performances were led by Sergio Barbero, whose presence in a race was typically worth 24 placings for his team's top finisher. Like Tonti, Barbero had a long career without many wins for himself -- the model domestique. The notables from 2009 thus far are the Milram riders Johannes Frohlinger, Peter Velits, and Fabian Wegmann.<br /><br />Overall, I was a little surprised at how short this list is. Very few riders, it appears, actually produce a significant improvement in team results. Of course, this is not to say that domestiques don’t earn their pay. Instead I interpret it to mean that domestiques are generally interchangeable, and there are very few riders who have an extraordinary ability to help a team leader.<br /><br />We can also determine the worst domestiques, riders whose teammates have better results when they’re not around. This isn’t necessarily a bad thing. I interpret these guys as the riders who shoulder the responsibility for winning, and when they’re not around someone else has to step up. Or maybe they’re just bad domestiques. Either way, there are many more significant anti-domestiques than domestiques when the same significance criterion is used:<br /><br /><a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjcliDVuEgvVakwmVDBkiFYLZnEV4fJfa6ZTDEWa6qoABp8aMKM_9j_sdI7ow4c1UrZ_kpX4f0boowK_mCjWlr212YtATFqKFRIjinu51-MTFW4v59FSQoVVTGd6Lh9hVkgkAUWBvsF3mVY/s1600-h/worstdomestiques.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 110px; height: 320px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjcliDVuEgvVakwmVDBkiFYLZnEV4fJfa6ZTDEWa6qoABp8aMKM_9j_sdI7ow4c1UrZ_kpX4f0boowK_mCjWlr212YtATFqKFRIjinu51-MTFW4v59FSQoVVTGd6Lh9hVkgkAUWBvsF3mVY/s320/worstdomestiques.png" alt="" id="BLOGGER_PHOTO_ID_5362184126921839458" border="0" /></a>This riders on this list make a lot of sense to me. Guys like Bettini and Zabel were true team leaders in that when they were in the race, no one else on the team needed worry about getting a result. Incidentally, Bettini was close to making the best domestiques list above for his 2002 season with Mapei, early in his career. GC riders tend not to appear here, possibly because their teammates often finish fairly high on mountain stages and hence the results are always pretty good.<br /><br />There are some obvious caveats with this analysis. For instance, there might be confounding factors like a domestique that is always paired in races with a very successful teammate (or teammates). This would associate the domestique with the results, even if the leader was independently a great rider. But if that leader rarely rode without this domestique, how would we know his greatness was independent of the domestique? This question is impossible to answer by looking at results alone. As a result of this, domestiques tend to appear in groups (e.g. Lampre in 2003). Furthermore, I only consider a team’s top placing in every race. There may well be domestiques that raise all of their teammates’ results, perhaps through their deft delivery of bottles from the team car.<br /><br /><span style="font-size:85%;">Technical Details: The data source is <a href="http://www.cqranking.com/men/asp/gen/start.asp">Cycling Quotient</a>. To avoid partial result listings, I considered approximately 900 races in which more than 100 riders are listed in the results. Significance was assessed using the Student’s t-test, a way of computing the likelihood that the two sets of race results are different. When performing a t-test for two sets of results, I require that each set has at least five results. Choosing a threshold for significance is a notoriously difficult problem. Naively one would take something like all P less than 0.05, but this implies a false positive every twenty t-tests. Over hundreds of riders this adds up to a lot of spurious positives. The conservative way to deal with this is to divide 0.05 by the number of t-tests (Bonferroni correction), but this assumes all tests are independent and that isn’t the case here. I chose an intermediate cutoff of P less than 0.0001.</span>Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-23877867579730731632009-07-23T21:58:00.000-07:002009-07-23T22:05:20.599-07:00SVD Analysis of the 2006 Tour<a href="http://tweakersonwheels.blogspot.com/2009/07/svd-analysis-of-2007-tour.html?pagewanted=print&pagewanted=print&pagewanted=print">As promised</a>, here is more of this <a href="http://tweakersonwheels.blogspot.com/2009/07/tour-de-france-has-34-stages.html?pagewanted=print&pagewanted=print&pagewanted=print">dorky, generally pointless stuff</a> – Part 3 of a series on using singular value decomposition to study grand tour results. In this post, I will look at the 2006 Tour de France, a year (in)famous for the Floyd Landis doping spectacle and Oscar Pereiro being declared winner in part due to gaining 30 minutes in a flat-stage breakaway that was not chased. Despite these odd circumstances, the overall complexity of this Tour was less than most grand tours, containing only 3.1 effective stages.<br /><br />As I discussed in the <a href="http://tweakersonwheels.blogspot.com/2009/07/svd-analysis-of-2007-tour.html?pagewanted=print">post on the 2007 tour</a>, SVD basically rearranges all of the results from all of the stages into a series of composite stages, the modes, that appear in the data with decreasing weights. Here are the modes and weights for the 2006 Tour:<br /><br /><a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEi8hWI4puVHfN4bcjoHVaK09yNWp2E2MX49ACyWkJMI6sWHyBRZeoA5BHq3BjiCcINEODfLIAL8aAbnRbnvSt6aPEXUFvJrKaLa-6j3k0d0jNg9olooGOEpijbkL8cDzO8LFu7kQ0bKalC-/s1600-h/tour06.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 400px; height: 170px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEi8hWI4puVHfN4bcjoHVaK09yNWp2E2MX49ACyWkJMI6sWHyBRZeoA5BHq3BjiCcINEODfLIAL8aAbnRbnvSt6aPEXUFvJrKaLa-6j3k0d0jNg9olooGOEpijbkL8cDzO8LFu7kQ0bKalC-/s400/tour06.png" alt="" id="BLOGGER_PHOTO_ID_5361887412764410722" border="0" /></a>The SVD modes are read as columns in the raster plot, with red and green corresponding to greater and lesser times for each stage. Only riders who finished the race are included in the results, and among those excluded is Floyd Landis since his results have been removed from the CQ record.<br /><br />The first striking feature of the mode patterns is how uneventful the first nine stages were. None of the large SVD modes has a large signal for these stages, meaning very few time gaps occurred. The major climbing stages were Stages 11, 15, 16, and 17, and sure enough these are the primary contributors to Mode 1. Stages 7 and 19 were individual time trials, the first of which had minor effects on the GC. Mode 2 primarily encodes large time losses on Stage 13 – the stage in which Voigt, Pereiro, Chavanel, Quinziato, and Grivko gained enormous time in the peloton. Mode 3 mostly encodes time gains on the Stage 17 to Morzine, which saw the peloton shatter in the wake of Landis’s shady escapade. Modes 4 and 5 are more corrections to mountain stage placings, and the latter modes are minor time gaps on flat stages and time trials.<br /><br />Each riders’ individual results can be recomputed by summing up these patterns, with each pattern separately weighted according to the individual rider. Looking at the weights for the first two modes for each rider, we see something quite unique in this Tour. This plot shows the extent to which each rider’s results (dot) exhibited Mode 1 (x-axis) and Mode 2 (y-axis), with final GC placing running from red to blue:<br /><br /><a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhG1BRlKKsfE9rIKkpylUomhl3Dt4XrnhYVxLEGREQqiwZbja6CdIm_YS3jUr69GCZlYJEzgm9fTVJcA3JfLpld3dEXEnHHr6cmybvKb8lCqKaPey3ZvLvRLurcWwmc3hy5Zdh20Ucx9edh/s1600-h/pca06.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 320px; height: 312px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhG1BRlKKsfE9rIKkpylUomhl3Dt4XrnhYVxLEGREQqiwZbja6CdIm_YS3jUr69GCZlYJEzgm9fTVJcA3JfLpld3dEXEnHHr6cmybvKb8lCqKaPey3ZvLvRLurcWwmc3hy5Zdh20Ucx9edh/s320/pca06.png" alt="" id="BLOGGER_PHOTO_ID_5361887544163316962" border="0" /></a>Comparing this with the results from the 2007 Tour we see a major rotation such that the results tend to spread from upper left to lower right rather than simply left to right. This, of course, is due to the Stage 13 breakaway. The four points isolate below are Pereiro, Voigt, Chavanel, and Quinziato (Grivko is not included since he did not finish the Tour). Their special status as breakaway survivors is clearly shown in their isolation on the plot. Note that Pererio is to the upper left of his breakaway companions, meaning that his Tour did have the makings for GC success even without the break. Looking at the top 10, we see how he compares to the others:<br /><br /><a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjlEnIp3T6_9AgNudsVYzL3L-uNXdzjYNLZQnMS_sFkwoFY4unLpOQ9_XRSg2igO2vGpFAx_y8S5b-7e2BUsiz0uaEdoHrUDuRMXfQz-2ek2ETpSmI_l0dGB2klwVCVnS38LYHIWh9XJUI5/s1600-h/pca_gc06.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 320px; height: 297px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjlEnIp3T6_9AgNudsVYzL3L-uNXdzjYNLZQnMS_sFkwoFY4unLpOQ9_XRSg2igO2vGpFAx_y8S5b-7e2BUsiz0uaEdoHrUDuRMXfQz-2ek2ETpSmI_l0dGB2klwVCVnS38LYHIWh9XJUI5/s320/pca_gc06.png" alt="" id="BLOGGER_PHOTO_ID_5361887671286808850" border="0" /></a>What are you doing all the way down there, Oscar? Winning the Tour de France, apparently.Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-53007747003118161292009-07-23T21:44:00.001-07:002009-07-23T21:52:19.658-07:00TdF Stage 18: The Problem with Greg LemondAlberto Contador scored a mighty impressive win in today's ITT. His post-race press conference was clouded by <a href="http://www.cyclingnews.com/news/contadors-climbing-credibility-questioned">questions about doping</a>, specifically generated by a <a href="http://www.lemonde.fr/sports/article/2009/07/23/alberto-prouve-moi-qu-on-peut-croire-en-toi-par-greg-lemond_1221871_3242.html?pagewanted=print">newspaper column</a> written by Greg Lemond (link is in French). In short, Lemond demands that Contador prove he is clean, otherwise we can only assume he is a doper. The very idea of computing a rider's physiological parameters from TV images is a bit silly, yet Lemond claims a specific VO2 max value for Contador that he considers impossible. When making arguments like this Lemond implicitly assumes arbitrary limits on the capabilities of cyclists, and those limits tend to align with his own performances.<br /><br />Personally I think Lemond is fine when he talks about the rise of EPO and its effects in the early 90’s. Clearly the EPO users got fast and there was indeed a peloton at two speeds, at least until everyone started using. When Lemond and others bemoan this I think they’re justified, because what they’re talking about is a change in ability across the entire peloton. The average rider got a lot faster, and while technology, training, and team support also improved I don’t think there’s any question that pharmaceuticals were the primary cause. <p>However, this is much different than looking at the performance of a single rider and calling them dirty, especially if that rider is the best in the peloton. Like everything else, rider ability is spread across a distribution and you can rarely say anything about the extremes of a distribution even if you can characterize the average. By definition, extreme cases like Contador are unlike anyone else in the population. They are literally in a class by themselves, so even if everyone they’re beating is doping it is impossible to know if they are as well (e.g. Armstrong in his prime). Unless they fail a test, of course. One would think that Greg Lemond, himself one of these extreme cases, would recognize this.<br /></p> <p>If we really did understand physiology to the extent Lemond thinks he does, then we could move past these statistical arguments. But the fact that statistics-based methods (epidemiology, statistical genetics) are still the source of most biomedical knowledge suggests otherwise.</p> <p>So when Lemond rants about the general use of drugs in cycling he’s got a point, but when he demands that a specific rider justify their results he’s just being a nut.</p>Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-57447650133950079852009-07-22T20:45:00.000-07:002009-07-24T16:10:30.583-07:00SVD Analysis of the 2007 TourI <a href="http://tweakersonwheels.blogspot.com/2009/07/tour-de-france-has-34-stages.html?pagewanted=print&pagewanted=print&pagewanted=print&pagewanted=print">previously posted</a> about using singular value decomposition (SVD) to analyze stage race results and estimated that based on the patterns in the results, grand tours could be reduced to about 3 or 4 stages without losing much information. The character of these “composite stages” depends on the results of each race on a case-by-case basis. Here I consider the 2007 Tour de France, which was the “simplest” tour, effectively having only 2.8 stages.<br /><br />SVD basically rearranges all of the results from all of the stages into a series of composite stages, the modes, that appear in the data with decreasing weights. The first mode is the most prominent pattern in the data, the second mode the second most prominent, etc, and the weight (aka the singular value) of each mode is that composite stage’s importance in the global data set. The rearrangement can be represented as a raster plot. Here are the modes and weights for the 2007 Tour:<br /><br /><a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhBmaG6k9NKPXlj6asuYI2jKabR8cDa9McIAIA2Umj8H0tKFJYlc_W10dHabY1PW5vdhC7D5pNvfbGhHsdQdKsO7J_yE_ovJ57gg59bsGjPxEvc7c-lcVrAhXhAZR_CjEfXaJRkCtr65naz/s1600-h/tour07.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 400px; height: 170px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhBmaG6k9NKPXlj6asuYI2jKabR8cDa9McIAIA2Umj8H0tKFJYlc_W10dHabY1PW5vdhC7D5pNvfbGhHsdQdKsO7J_yE_ovJ57gg59bsGjPxEvc7c-lcVrAhXhAZR_CjEfXaJRkCtr65naz/s400/tour07.png" alt="" id="BLOGGER_PHOTO_ID_5361497481427914242" border="0" /></a>The SVD modes are read as columns in the raster plot, with red and green corresponding to greater and lesser times for each stage.<br /><br />Looking up the first mode, it is not surprising that Stages 8, 14, 15, and 16 were mountain stages in the 2007 tour. In most grand tours, the mountain stages dominate the major modes because the time gaps in the mountains are so much greater than those in other stages. Mode 1 in particular has an much larger weight than any other mode, and it describes 93% of the results in the Tour. It can be though of as the base pattern for any rider’s results. This first mode alone correctly determines the three riders on the final podium (Contador, Evans, Leipheimer), although it interchanges Evans and Leipheimer’s final placings. It primarily encodes the climbing stages, but also includes some information on the ITTs in Stages 13 and 19. The second mode encodes a correction to the first, which is a time lost in the Alps relative to the Pyrenees, presumably representing riders who grow stronger in the last week. Mode 3 is mostly gains on Stage 9 paired with losses on Stage 15. Modes 4-8 contain further corrections to the mountain stages. Modes 9-14 generally account for time gaps in stages where breaks played a role and details of time trials. Modes 15-21 are slight differences in sprint stage and prologue performances – these weights are extremely small since most of the field finished sprint stages with the same time.<br /><br />Each riders’ individual results can be recomputed by summing up these patterns, with each pattern separately weighted according to the individual rider. A GC rider, therefore, will have a small weight for Mode 1 since they didn’t lose much time on the decisive stages, whereas a consistent member of the autobus will have a large weight for that mode corresponding to their large time losses in the mountains. We can, for instance, look at the first two modes for each rider. This plot shows the extent to which each rider’s results (dot) exhibited Mode 1 (x-axis) and Mode 2 (y-axis), with final GC placing running from red to blue:<br /><br /><a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEj2PU9NSj-tJ5LqMDm4g7dzBeK4GWRlwjeDglUtyNxjXk5uup8h82IeDjd0k65jHyp7tmmSR6dmkDBXASwUSX9pDyM_Gy5N0kXROt1-YEv70U326hp9yYaD6IObjLpfiWH1ko2u5bfBfpLg/s1600-h/pca07.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 320px; height: 312px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEj2PU9NSj-tJ5LqMDm4g7dzBeK4GWRlwjeDglUtyNxjXk5uup8h82IeDjd0k65jHyp7tmmSR6dmkDBXASwUSX9pDyM_Gy5N0kXROt1-YEv70U326hp9yYaD6IObjLpfiWH1ko2u5bfBfpLg/s320/pca07.png" alt="" id="BLOGGER_PHOTO_ID_5361497707457655554" border="0" /></a>GC riders have lower values for Mode 1 since they lost the least time on important stages. The blob of blue on the far right is the autobus riders who consistently lost a log of time. Zooming in on the top 10 we get an idea of how the GC is scattered:<br /><br /><a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjbf6n0ZIGKxkoejCARJEykJ_sEEZ94xgjuUAzuP01jWG7qPQxIFoVJ5vhOPkNz1Lw45tSt5EYc7Yv9qsQO1ufLfVTcub9Rk0QZ84pENYGOrlPHWfmUkZFuKU06cbgcA2HNC2U_0aljKzlf/s1600-h/pca_gc07.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 320px; height: 312px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjbf6n0ZIGKxkoejCARJEykJ_sEEZ94xgjuUAzuP01jWG7qPQxIFoVJ5vhOPkNz1Lw45tSt5EYc7Yv9qsQO1ufLfVTcub9Rk0QZ84pENYGOrlPHWfmUkZFuKU06cbgcA2HNC2U_0aljKzlf/s320/pca_gc07.png" alt="" id="BLOGGER_PHOTO_ID_5361497629088282146" border="0" /></a>Contador has the smallest Mode 1 component – he lost the least time in the big stages. Along with his teammates Leipheimer and Popovych, he also did relatively well later in the race, as signified by his large Mode 2 component. Although Leipheimer did well overall (Mode 1), his early losses in the Alps (Mode 2) were not great enough to overcome Evans. We could continue to add modes and represent the data in higher dimensions, but I think two is enough for a blog post.<br /><br />The 2007 Tour is a fairly typical case for SVD. Most Tours show a similar dominant mode that combines climbing and time trial stages, encompassing most of the GC. It usually takes the addition of one or two more modes to work out the precise order of the podium and a couple more to fill in the details of the top 10, but the majority of activity encoded in the smaller modes are mid-GC reassortments due to breaks.<br /><br />I think this generally means Tour results can be summarized by a few patterns in the data and therefore major changes in riders’ performances in the course of a grand tour are very rare. This is, however, a global and quantitative view of the data and should be taken as such; the difference between fourth and seventh place for an individual rider is still be quite important to them and their supporters! Vinokourov’s miniscule time gains on the Champs Élysées in 2005 do not appear until Mode 19 in a 2005 SVD analysis, but it certainly mattered to him and Levi Leipheimer.<br /><br />Next I’ll consider the 2006 Tour de France. How does Oscar Pereiro’s unorthodox victory look through SVD glasses?Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-33010743156944057652009-07-22T20:38:00.000-07:002009-07-22T20:44:34.228-07:00TdF Stage 17The GC activity today was all very exciting, but the ride of the day has to go to Thor Hushovd. <a href="http://www.cyclingnews.com/races/96th-tour-de-france-gt/stages/stage-19/results">Said Thor</a>:<br /><br /><blockquote>"Today everything just felt perfect. I attacked over the first climb, did a good descent and then had an amazing day in the front," Hushovd said. "I think this is the best day I've ever had on the bike."<br /></blockquote><br />This is being played up as a smackdown to Cavendish (as if this bike race was lacking in manufactured histrionics), but I think that distracts from what an awesome ride this was. It is great to see my geeky blogging disrupted by news of manly riding!Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-55916624505833460532009-07-22T16:05:00.000-07:002009-07-22T16:13:26.661-07:00The Tour de France has 3.4 StagesWith a time trial and mountain stage remaining in this year’s Tour de France, there are no lack of assurances that <a href="http://tour-de-france.velonews.com/article/95662/third-week-a-mystery-for-contador">the Tour is not over</a>. But once the GC has been sorted out by the early stages, how much do the results actually change? Is there any major variance in the results between each type of stage, or can the Tour be reduced to one representative mountain stage, a representative time trial, and a typical sprint stage?<br /><br />To determine how complex the Tour is, I calculated the <a href="http://en.wikipedia.org/wiki/Singular_value_decomposition">singular value decomposition</a> (SVD) of stage finish times. SVD is a linear algebra technique that rearranges data into a series of features, called modes, weighted from the most important to least important. For example, it is used in image compression to reduce a matrix of pixels to an approximate matrix that is reconstructed from the a handful of modes. Furthermore, by looking at how quickly the mode weights decrease, one can estimate the effective number of independent components in the data. A famous example was a <a href="http://www.pnas.org/content/100/13/7432.abstract">study of correlated voting records</a> in the US Supreme Court, which concluded that the nine justices could be approximated well by 4.7 independent justices with uncorrelated voting patterns. Another way of stating this is that the information content of US Supreme Court decisions is the same as a court of 4.7 justices rather than nine.<br /><br />Looking at Tour results, we can ask the same question. The Tour is usually 21 stages, but we can perform SVD analysis and determine the effective number of stages. I did this for a number of recent grand tours using results for time lost on each stage from <a href="http://www.cqranking.com/men/asp/gen/start.asp">Cycling Quotient</a>. Here are the results:<br /><br /><a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhHXxugqGDybDEqaIjylkzwwdFM8vsC2e7BbdWx3Rc0mQE8e_9yg-ge72frXpDa1KaUuZcxGTlHsQll4rhleSzVcB9MhLyYJFea3AVq4lxnuBQLlWzvmhiHiAqBlIQ5QDDRKKcbrK3sj7z9/s1600-h/effective_stages.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 320px; height: 259px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhHXxugqGDybDEqaIjylkzwwdFM8vsC2e7BbdWx3Rc0mQE8e_9yg-ge72frXpDa1KaUuZcxGTlHsQll4rhleSzVcB9MhLyYJFea3AVq4lxnuBQLlWzvmhiHiAqBlIQ5QDDRKKcbrK3sj7z9/s320/effective_stages.png" alt="" id="BLOGGER_PHOTO_ID_5361425287787076450" border="0" /></a>The average Tour effectively has 3.4 stages, the average Giro 3.3, and the average Vuelta 4.3. The most complex tour here is the 2003 edition, which involved multiple long breaks and a closely contested GC. The Vuelta is consistently more complex than the Tour and Giro, perhaps because breaks are allowed freer reign or riders in the lead are more likely to crack after a long racing season.<br /><br />So does the SVD analysis reduce the race to a single climbing, time trial, and sprinting stage, with perhaps a “half stage” for hilly transition stages? To fully answer this, we need to look at each stage race on a case-by-case basis. I will look at a few races in forthcoming posts. Briefly, since the consistently good climbers also tend to be the better time trialists, time trials and climbing stages are often combined into a single representative “GC mode”. Additional modes encode variations around this main trend, such as a difference for some riders between the Alps and Pyrenees, and breakaway results. It should be kept in mind that this analysis is entirely data driven – the modes correspond to how the results panned out rather than any preconceived notion of which stages are important in a grand tour. Since I do not have a handy source for complete results during the Indurain years I could not do the analysis, however I imagine the balance of time trialing and climbing might have been different then.<br /><br /><span style="font-size:85%;"><br />Technical notes: SVD is a linear algebra routine that produces a unique solution for a given data matrix. The results matrix was composed of stages x riders, such that each matrix element is the time lost in seconds for a given rider on a given stage. Time bonuses were not factored in because of personal laziness. Riders who did not finish all stages were excluded since SVD chokes on missing data. The effective number of stages is computed as the Shannon entropy for the fractional singular values squared. Missing grand tours are due to incomplete results in the CQ database.</span>Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-14721064463969620102009-07-21T12:21:00.000-07:002009-07-24T22:43:41.074-07:00Should Unprecedented Tour Success Warrant Suspicion?In recent years at the Tour de France we’ve seen breakthrough results from riders like Floyd Landis and Bernhard Kohl nullified by doping convictions. In this light, it is tempting to view the current success of Bradley Wiggins with skepticism (see <a href="http://www.podiumcafe.com/2009/7/19/954602/when-is-it-appropriate-to-be">here</a> for a generally civil discussion). Wiggins, now third on GC, is a multiple gold medalist and world champion on the track, but has never had much success in road or stage racing. He credits his new climbing ability to substantial weight loss paired with a new focus on road racing, but a hardcore skeptic who believes <a href="http://www.cyclingnews.com/news/kohl-tells-all-about-doping-1">Bernhard Kohl</a> is not going to be so easily convinced.<br /><br />Of course, at the moment there is no way we can know if Wiggins is on the good juice or not. But we can look at his pattern of improvement and ask how extraordinary it actually is. Similarly, we can look at dopers and determine if their results were entirely unprecedented. To this end, I will consider two quantities:<br /><ul><li><span style="font-style: italic;">R</span>: The difference between a rider’s mean placings in a given year and mean placings in previous years. Larger numbers mean greater improvement.</li><br /><li><span style="font-style: italic;">P</span>: The likelihood that this difference is real and not simply a result of random fluctuations (statistical significance). Smaller numbers mean greater significance. </li></ul>Data, as usual, are from the indispensable <a href="http://www.cqranking.com/men/asp/gen/start.asp">Cycling Quotient</a>. When computing <span style="font-style: italic;">P</span>, I first take the logarithm of all results in order to enhance the value of top placings and minimize differences between mid-pack finishes (e.g. the difference between 2nd and 12th is much more important than the difference between 102th and 112nd). See below for more exciting technical notes.<br /><br />Here are the career results for Bradley Wiggins:<br /><br /><a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiUkFhMZjs3uXF20CbW4e_dodUuFXWttL3iBGjyStZXXbDcncIkvWrtvuTPUM9T29ddbsDmLBzH2s15QZ7Zncdsb3Dfe273JCJuSDtCz0vAo2GSFsdUzmqR18nNSRDr-qZZ7rM6QoOzoBav/s1600-h/wiggins.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 360px; height: 255px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiUkFhMZjs3uXF20CbW4e_dodUuFXWttL3iBGjyStZXXbDcncIkvWrtvuTPUM9T29ddbsDmLBzH2s15QZ7Zncdsb3Dfe273JCJuSDtCz0vAo2GSFsdUzmqR18nNSRDr-qZZ7rM6QoOzoBav/s400/wiggins.png" alt="" id="BLOGGER_PHOTO_ID_5360996481712728690" border="0" /></a>Wiggins results so far in 2009 are definitely the best of his road career, with a mean placing of 66 compared to 108 over previous years. This corresponds to <span style="font-style: italic;">R</span> = 42.5. The likelihood of this result in random fluctuations is <span style="font-style: italic;">P</span> = 0.00012, suggesting that this is indeed statistically significance. Here is how Wiggins compares to a few other riders:<br /><br /><a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhc6aCVGhk5ijvU9K_WhHyjxoJq00I3-sQOUNf34Jx5C2tvI2ctoG5y3yH835E_GX71xcgdHHDdk-Gs5N87Iy2-QPawPSgjy85Xt8eri8_iexLYXS7D1Nq8t-mS593QW_EzLeaB___ji9n1/s1600-h/wiggins+tab2.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 200px; height: 54px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhc6aCVGhk5ijvU9K_WhHyjxoJq00I3-sQOUNf34Jx5C2tvI2ctoG5y3yH835E_GX71xcgdHHDdk-Gs5N87Iy2-QPawPSgjy85Xt8eri8_iexLYXS7D1Nq8t-mS593QW_EzLeaB___ji9n1/s200/wiggins+tab2.png" alt="" id="BLOGGER_PHOTO_ID_5361019295224298770" border="0" /></a>Sastre’s results appear to have improved slightly, but the <span style="font-style: italic;">P</span>-value tells us this is not significant. Haussler is having a very good season but not to the extent that Wiggins is. Vande Velde was impressive in 2008, but Wiggins has shown an even greater improvement in 2009.<br /><br />So how rare is this magic season of Bradley Wiggins? Looking back across the CQ data, I identified other riders who have had a season with an improvement greater than or better than Wiggins, defined as a <span style="font-style: italic;">R</span> greater than 42.5 and <span style="font-style: italic;">P</span> less than 0.00012:<br /><br /><a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiFp68AWje0houPbJwAawmBVLazI7pXfivu3PR1qgdUEZG7Wt0w9R-pU0pPFrHTEUQcWXdfuswDs7Gm7EI7ZD6Sy94qnywnScHmZVgCia-S-Ywcd7_LsyisfksgqCBLhxJRHHm_jBS583ER/s1600-h/wiggins+tab1.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 203px; height: 320px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiFp68AWje0houPbJwAawmBVLazI7pXfivu3PR1qgdUEZG7Wt0w9R-pU0pPFrHTEUQcWXdfuswDs7Gm7EI7ZD6Sy94qnywnScHmZVgCia-S-Ywcd7_LsyisfksgqCBLhxJRHHm_jBS583ER/s320/wiggins+tab1.png" alt="" id="BLOGGER_PHOTO_ID_5361017375858570018" border="0" /></a>Should we be skeptical of Wiggins based on this company? There are a few suspicious characters there, but overall this list is notably light on convicted dopers. We can compare these riders with some of the recent bad boys:<br /><br /><a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjg53uIsP-HagB3rIxri_2qdLZ_t_v32-a-87gG0Iodm_WmzRMrRI9b3wqNhIvECPld3jW3-9J9vfsT0_sT2WIf3MabN-MNwEGE3OCSeIMh-_ot_iPoKZCoQ9tB5zg1k6YxsKxzTsoKU0Se/s1600-h/wiggins+tab3.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 200px; height: 140px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjg53uIsP-HagB3rIxri_2qdLZ_t_v32-a-87gG0Iodm_WmzRMrRI9b3wqNhIvECPld3jW3-9J9vfsT0_sT2WIf3MabN-MNwEGE3OCSeIMh-_ot_iPoKZCoQ9tB5zg1k6YxsKxzTsoKU0Se/s200/wiggins+tab3.png" alt="" id="BLOGGER_PHOTO_ID_5361017536791572898" border="0" /></a>These improvements are generally more modest. Kashechkin, Landis, and Basso are the only riders in this lot whose gains were in the neighborhood of Wiggins's 2009. Three of these riders actually had a <span style="font-style: italic;">drop </span>in results when doping (negative <span style="font-style: italic;">R</span>), although these were not significant.<br /><br />The data fairly clearly contradict the idea that significant performance improvements are necessarily, or even likely, the result of doping. As always, we should apply the standard caveat that the absence of a positive test does not always imply clean riding. However, without evidence to the contrary I think we should conclude that it is common for riders to significantly improve their performances without the aid of doping. Likewise, it is common for dopers not to have significantly improved results. <span style="font-size:85%;"><br /><br /><br />Technical notes: Most race results on Cycling Quotient are partial, often listing only the top 10 or 20 riders. For calculation of <span style="font-style: italic;">R</span> and <span style="font-style: italic;">P</span>, only races with 100 or more listed riders were used. Calculations were therefore based on about 900 races from 2000 to 2009, with the later years having more races listed. Roughly half of the races were stages from grand tours, and the remaining results are mostly the major one-day races and lesser stage races. To avoid small sample sizes, <span style="font-style: italic;">R</span> and <span style="font-style: italic;">P</span> were only computed if the rider at hand had more than 10 race results for the year in question and more than 10 results prior to that year. The likelihood <span style="font-style: italic;">P</span> corresponds to the probability that there is no difference between a rider’s mean results in a given year and their mean results in previous years, and was computed using the Student’s t-test without assuming known or equal variances between the two samples.<br /></span>Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-68984518141670955772009-07-21T12:18:00.000-07:002009-07-21T12:21:17.194-07:00TdF Stage 16Well, wasn't that a blast from the past. Nice to see the old man throw down and bridge on a climb like that. He might just hold that second place yet.Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-23004937845686927372009-07-19T18:22:00.000-07:002009-07-19T22:16:44.031-07:00TdF Stage 15Perhaps it's a little early to call it, but the Tour GC tends to be decided by mountaintop finishes and we've probably seen enough. I think the major remaining questions are if Armstrong and Wiggins can maintain their podium spots on Ventoux.<br /><br />Armstrong got <a href="http://www.cyclingnews.com/news/sastre-respects-his-rivals-armstrong-should-do-the-same">a lot of grief</a> for a comment regarding the quality of the top finishers in last year's Tour, but after today it looks like he might have had a point. Running down the top 5 from '08:<br /><ol><li>Sastre is 11th at 3:52,</li><li>Evans is 14th at 4:27,</li><li>Kohl is flipping burgers in Vienna,</li><li>Menchov is 29th at 11:23,<br /></li><li>Vande Velde is 12th at 3:59.</li></ol>I am now waiting for Armstrong to "kick their asses". Team time trials do not count.Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-14877043405817237952009-07-18T13:34:00.000-07:002009-07-18T13:35:50.151-07:00TdF Stage 14My goodness, who would have ever thought a win by Sergei Ivanov would get everyone so hot and bothered?Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-33847628240613845022009-07-17T15:50:00.000-07:002009-07-20T17:09:53.761-07:00TdF Stage 13<span style="font-family:arial;">I have read that the other racers call Heinrich Haussler "Barbie" for his attention to fashion and grooming. Crying as he won today's stage probably won't help that image, although a 195 km break ahead of the cold, whimpering peloton should earn him some respect. After he took a painful-to-watch second place at Milano-San Remo and another second at De Ronde, today's result isn't a surprise. Looking at Haussler's pro results prior to the start of the Tour shows how he has been steadily marching down the top ten over the past few years:</span><br /><br /><br /><a style="font-family: arial;" onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiINJTW7BWhmgJPx3I_VNaGeUJX7vZUrwTwZkbUq9_3_CbSfGfXOx9HS0ZwC3chNTnqIf00xL5-4KpHmN8atDwWZv7SvYu2_m4-tc_V4yzh9I-snA46OLIKkdMqmrO1HnSTZbcKgbeZZ4wK/s1600-h/haussler.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 360px; height: 248px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiINJTW7BWhmgJPx3I_VNaGeUJX7vZUrwTwZkbUq9_3_CbSfGfXOx9HS0ZwC3chNTnqIf00xL5-4KpHmN8atDwWZv7SvYu2_m4-tc_V4yzh9I-snA46OLIKkdMqmrO1HnSTZbcKgbeZZ4wK/s400/haussler.png" alt="" id="BLOGGER_PHOTO_ID_5359565543574195234" border="0" /></a><span style="font-family:arial;">[Data from </span><a style="font-family: arial;" href="http://www.cqranking.com/men/asp/gen/start.asp">CyclingQuotient</a><span style="font-family:arial;">]<br /><br /></span>Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-91796041690777209422009-07-16T14:46:00.001-07:002009-07-17T15:54:27.347-07:00TdF Stage 12<span style="font-family:arial;">I've often heard that recovery becomes more difficult with age. If that's true, Lance Armstrong must be pleased with the way the hard stages are spread out in this year's Tour.</span><span style="font-family:arial;"><br /><br /><br />I think Nocentini is going to lose his fancy shirt tomorrow.</span>Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-83330809802674916022009-07-15T12:38:00.000-07:002009-07-15T15:52:33.342-07:00TdF Stage 11<span style="font-family:arial;">After years of working to chase down breaks for Robbie McEwen, gangly Belgian Johan Van Summeren got himself off the front today. Great to see him on the other side of that situation for once. Unfortunately for him, Columbia-HTC had plenty of help in bringing the race back together for a bunch sprint. There was a rumor that the finishing incline would be difficult for Cavendish, and I guess it was -- he could only beat Farrar by a wheel. </span><span style="font-family:arial;"><br /><br /><br />It must be really fun to be Eisel, Grabsch, or any of the other Columbia-HTC riders right now. There are few things more satisfying in cycling than working your ass off for a teammate that has the talent, smarts, and drive to deliver the win.<br /><br /><br />In </span><a style="font-family: arial;" href="http://www.cyclingnews.com/news/cavendish-defends-against-charge-of-racism">very important news</a><span style="font-family:arial;">, apparently the UCI leaked to L'Equipe that some unnamed rider thinks Cavendish is an anti-French racist! Cavendish dealt with this matter appropriately:</span><span style="font-family:arial;"><br /></span><blockquote><span style="font-family:arial;">"For sure I’m going to get arsy at some riders, because, you know, I’m an asshole," said Cavendish. "But it’s irrelevant their nationality, and irrelevant what they look like, or where they come from. Because, like I said, I’m an asshole."</span><span style="font-family:arial;"></span></blockquote><span style="font-family:arial;">Now Cavendish could have escaped this controversy by pointing out that as a sort-of-British person he is expected to dislike the French, but clearly he has come to the Tour in better form than that. Bravo!</span>Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-4268729548570301032009-07-14T12:43:00.000-07:002009-07-20T17:09:06.454-07:00Genetics of Athletic Performance (Dog Version)<span style="font-family:arial;"><a href="http://www.genome.gov/12513335">Dr. Elaine Ostrander</a> is a geneticist at the NIH. She studies genetic variation within and between populations, and tries to identify genetic mutations that cause disease. Her lab studies dog genetics as a model for understanding the genetic structure of populations. Dogs are fairly interesting in this regard because although they are all members of the same species, different breeds show very diverse phenotypes. The size difference between a great dane and a chihuahua, for instance, is far greater than the size difference between any two people. The cause of this size difference is genetic, and the Ostrander lab traced the variation to mutations in the gene IGF-1.</span><span style="font-family:arial;"><br /><br />As part of this research program, Ostrander’s lab published a </span><a style="font-family: arial;" href="http://www.plosgenetics.org/article/info%3Adoi%2F10.1371%2Fjournal.pgen.0030079">paper in 2007</a><span style="font-family:arial;"> on how a genetic mutation in whippets increased their athletic performance. Whippets are lean racing dogs that were bred from greyhounds in the late 1800s and, like amateur bike racers, they compete in categories based on racing ability. Dogs move up and down categories, ranked A to D from fast to slow, as they proceed through their careers.</span><span style="font-family:arial;"><br /><br />Not surprisingly, an individual whippet’s physiology has a large effect on their racing ability. In particular, there is a subpopulation called “bully whippets” that show larger than normal musculature. This condition is known as “double muscling” in some breeds of cattle and, in one known case, a human. These individuals have a substantially greater number of skeletal muscle fibers. Furthermore, in bully whippets this enhanced musculature tends to appear at two different degrees – bigger and biggest. To a geneticist, this strongly suggests that the cause is a single genetic mutation that is partially recessive. Dogs with one copy of the mutant gene and one normal gene (heterozygotes) are somewhat more muscular than normal, and dogs with two copies of the mutant gene (homozygotes) are extremely muscular. Following up on this, the Ostrander lab identified a mutation in the myostatin gene (MSTN) as the cause of the bully whippet phenotype. Sure enough, dogs with two mutant copies of MSTN were double muscled, and dogs with a single mutant MSTN had enhanced musculature.</span><span style="font-family:arial;"><br /><br />The MSTN mutation had a positive effect on racing ability, judging by the fact that mutant dogs were more likely to be Grade A racers. The MSTN gene was genotyped in 85 individuals and 13 were found to have the mutation. Nine of these 13 were Grade A racers, three were Grade B, and one was Grade D. If this mutation didn’t have an effect on race performance we would expect about 3 mutant dogs in each category; the likelihood we would get 9 or more in Grade A is 0.03%. Although it would be nicer to have more individuals and hence better statistics, this is pretty solid evidence that the MSTN mutation confers an advantage in racing.</span><span style="font-family:arial;"><br /><br />This is interesting but not very surprising. Clearly, certain genetic variants are going to improve athletic performance, especially ones that increase muscle mass. One might ask if there are other genetic determinants of racing ability that might explain the 13 Grade A dogs that had a normal MSTN gene. Are they also genetically superior to the other grades, in a way that doesn’t directly involve MSTN?</span><span style="font-family:arial;"><br /><br />Indeed, they are. To study this, the researchers genotyped each of the 85 dogs at 32 sites scattered across the genome and then compared the results. Based on this fairly sparse view of the genome (a dog has 39 chromosome pairs, so this is less than one site per chromosome), they computed a representative number they call pc1 for each individual (technically, the first principle component of the genotype). While one individual pc1 isn’t meaningful, we can say that two or more dogs with similar pc1 values are genetically similar. The results looked like this:</span><span style="font-family:arial;"><br /><br /></span><a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiTDT1O0AtMAsCQ3P07sdKT3TMcWtoV2CZJGQDctFptxjSvqP0Au0sWVK28bvud0-lIvqhslxyrbgU8hL7HqaykJA9uGgvjsXLs8LJQUUR9MIPL1wiCfU4AUGFheUUFhFoFneAXiXJ2942W/s1600-h/whippets1.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 400px; height: 286px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiTDT1O0AtMAsCQ3P07sdKT3TMcWtoV2CZJGQDctFptxjSvqP0Au0sWVK28bvud0-lIvqhslxyrbgU8hL7HqaykJA9uGgvjsXLs8LJQUUR9MIPL1wiCfU4AUGFheUUFhFoFneAXiXJ2942W/s400/whippets1.png" alt="" id="BLOGGER_PHOTO_ID_5358404887656011890" border="0" /></a><span style="font-family:arial;">The x-axis, pc1, represents the genotype and the y-axis shows the race grade (plus a small random number to scatter the points a little). It is clear that Grade A dogs tend to have low pc1 and Grade D dogs generally have high pc1. Grade B dogs cluster in the middle, and Grade C dogs are spread into a middle and low group. This population structure is partially due to the greater inbreeding within a race grade – fast dogs are bred with other fast dogs, so fast dogs will tend to be genetically similar. It is difficult to tell how great of an effect this has on the current study since the dog’s ancestries were not reported. That caveat aside, the cluster of points in the upper left suggests that the genotypes of the Grade A dogs provide a genetic advantage over the other race grades.</span><span style="font-family:arial;"><br /><br />This effect becomes even more clear when we consider only the dogs without the MSTN mutation:</span><span style="font-family:arial;"><br /><br /></span><a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgpyDENnW6h5TbSJwkNjniaSmC9CEK73UFl2fL1Fi3FpF-MEEB0UmAhkvbo-X9HeR-Ef0t_cJDrVVrOZA5QFTm_fVfhWIslnfBYCmh7zS7ZhAgA3GIsUzePbh7ZzMvld_Wa84qHCJin6aeE/s1600-h/whippets2.png"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 400px; height: 286px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgpyDENnW6h5TbSJwkNjniaSmC9CEK73UFl2fL1Fi3FpF-MEEB0UmAhkvbo-X9HeR-Ef0t_cJDrVVrOZA5QFTm_fVfhWIslnfBYCmh7zS7ZhAgA3GIsUzePbh7ZzMvld_Wa84qHCJin6aeE/s400/whippets2.png" alt="" id="BLOGGER_PHOTO_ID_5358405016447515762" border="0" /></a><span style="font-family:arial;">In this second plot points have been faded for dogs with at least one MSTN mutation (bully whippets). We see that all but one of the Grade A dogs with the highest pc1 values carry the mutation, suggesting that this single genetic abnormality has turned these five slower dogs into Grade A racers. Looking only at the dark points, we see strong evidence that MSTN is not the only genetic contributor to racing ability; there are 13 Grade A dogs without the mutation, and all except for a single outlier have a genetic makeup (negative pc1 values) that sets them apart from most dogs in the other grades. The authors suggest that about 70 other genes are contributing to this genetic advantage. </span><span style="font-family:arial;"><br /><br />So a major result of this study is a genetic signature for a fast whippet, comprised of 32 markers across the genome. In principle, if one were to give these researchers the DNA from a rookie whippet, they could probably predict to reasonable accuracy whether that dog will be Grade A based on its genotype. The test won’t be perfect, but one might significantly improve the test by including more than 32 genetic markers. Current technology can read millions of genetic markers for a few hundred dollars, so it certainly could be done.</span><span style="font-family:arial;"><br /><br />The paper concludes by noting that the MSTN mutation also confers greater musculature in humans, and therefore is a candidate for gene doping. With only one known human case of this mutation it is impossible to predict side effects, rendering this proposition quite dangerous even if it were technically feasible. However, it is probably a matter of time before such ideas will make their way into athletics. Personally I propose an alterative use for genetics in cycling. I would like to collect DNA samples from a few hundred of my fellow amateur racers, genotype them all, and determine the genetic profile of a Cat 1/2, Cat 3, and Cat 4 racer. New racers will be forced to do their ten Cat 5 races, after which they will be genotyped and assigned to their appropriate category. No more unnecessary thrashings of the Cat 3 and 4 fields by clearly superior racers who are amassing upgrade points – just put everyone where they belong and let their training, race smarts, teams, and luck determine the outcome. It will be a new age of bike racing, brought to us by science! Are you reading this, Shawn Farrell?</span>Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-43222965316368476932009-07-14T09:15:00.000-07:002009-07-15T17:14:30.692-07:00TdF Stage 10<span style="font-family:arial;">ASO's plans to liven up le Tour by taking everyone's radios away seems to have backfired. In fact, the race was so boring I suspect we saw a subtle protest at the finish -- Hushovd and Farrar agreed not to try coming around Cavendish so that even the sprint would be tedious. Or maybe they were waiting for their DS to point out the finish line in the last 200 meters. In any case the stage made the rest day press conferences look action packed.</span><span style="font-family:arial;"> Maybe next time assign all riders and DS random radios?<br /><br /><br />In other, equally exciting news it seems that Giro has their </span><a style="font-family: arial;" href="http://www.cyclingnews.com/races/96th-tour-de-france-gt/stages/stage-7/photos/78666">new lid</a><span style="font-family:arial;"> on a number of teams (Astana, Garmin, etc). Apparently it's very light and weighs not so many grams, but it looks to me like the helmet Postal was wearing in 1999 (maybe now with a Greek name?). I think this will be a marketing challenge for Giro; will the weight savings overcome the dated styling in the eyes of consumers? If my riding cohort is representative of the market, the answer is no. Lightness is fine, but cycling is all about looking good these days. In white, like a bride.</span>Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-69867735509907141092009-07-12T17:24:00.000-07:002009-07-14T12:25:13.036-07:00TdF Stage 9<span style="font-family:arial;">Egoi Martinez took the mountains jersey today, with Franco Pellizotti 23 points behind. I expected David Moncoutie to make a stronger run for the jersey, but he did enough too keep himself in contention by getting in today's early break and fighting for a few points. Yes, this is a minor matter, but I don't have much more to say about the fairly boring stage today. Whatever happened to that plucky little mountain scamp, Riccardo Ricco? I miss him.</span><span style="font-family:arial;"><br /><br /><br />I also missed most of the TV coverage since I was "enjoying" a 40km time trial today. Unfortunately I did just well enough to encourage myself to put some more work into TTs. My current muscle aches provide evidence that Fabian <a href="http://www.marlyspearson.com/cancellara%20time%20trial.jpg">Cancellara's massive glutes</a> have contributed greatly to his success.</span>Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-66806780982684997812009-07-11T12:49:00.000-07:002009-07-14T12:25:04.737-07:00TdF Stage 8<span style="font-family:arial;">I do love a small group sprint, and today Luis Leon Sanchez made it look easy. He waited form Casar to do most of the work to bring Efimkin back, closed the gap himself with a few pedal strokes, then paused and let Casar come around to lead him out. He glanced over his right shoulder as if he was giving Casar permission to go, then came around him easily for the win.</span><span style="font-family:arial;"><br /><br /><br />I hung around to watch the Versus interviews this morning. They really need to give George Hincapie a while to recover before doing post-stage interviews. Let him have something to drink, maybe some food to get some glucose back in his brain. Ambushing him while non-functional just isn't fair. Then Thor Hushovd had trouble getting his massive Viking hands through the tiny bike-racer sleeves on that green jersey.</span><span style="font-family:arial;"><br /><br /><br /><a href="http://www.cyclingnews.com/news/evans-attacks-but-fails-to-gain-time">Johan Bruyneel claims confusion</a> at the early attack from Cadel Evans. Don't play coy with us, you Belgian minx! Everyone knows you are the evil genius of le Tour, the puppetmaster who controls all.<br /><br /><br />Oscar Pereiro apparently gave up today. Is his career Beloking?</span>Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-23722569715671380812009-07-10T12:55:00.000-07:002009-07-10T13:13:24.531-07:00TdF Stage 7<span style="font-family:arial;">Fantastic win by Brice Feillu as France takes first and second on a mountaintop finish. Cue the "Is this what a clean Tour looks like?" speculation. It's another great result for France and a big day for Feillu's title sponsor, <a href="http://tour-de-france.velonews.com/photo/94860">le ribcage</a>.</span><span style="font-family:arial;"><br /><br /><br />In GC news, so concludes the Great Astana Leadership Drama of 2009. Contador now leads Armstrong and I doubt that is going to change. Armstrong should race for second, donate millions to the UCI anti-doping program, and maybe have a longshot chance for an eight Tour win as we head into fall. Leipheimer looked more comfortable than I expected, and Evans continues to loosen up his style nicely. Andy Schleck was disappointing; I expected him to follow Contador but he couldn't even drop his fat brother. And how about track star Bradley Wiggins climbing with the GC contenders and attacking in the final kilometer? Suspicalicious! </span>Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-68462093696395241442009-07-09T14:05:00.000-07:002009-07-09T14:10:09.327-07:00TdF Stage 6<span style="font-family:arial;">Over 20 riders crashed in today's stage. That looked far more dangerous than the TTT; I wonder who will be the first rider to call for canceling races when it rains. Angry about being left out, René Haselbacher crashed at the </span><a style="font-family: arial;" href="http://www.cyclingnews.com/races/61st-tour-of-austria-2-hc/stages/stage-5/results">Tour of Austria</a><span style="font-family:arial;">. </span><span style="font-family:arial;"><br /><br /><br />The UCI </span><a style="font-family: arial;" href="http://www.cyclingnews.com/news/uci-and-afld-running-targeted-testing-at-tour">sorta promises</a><span style="font-family:arial;"> next-day turnover for doping results. This could be a problem, since it doesn't allow much time to leak the positives to L'Equipe.</span><span style="font-family:arial;"><br /><br /><br />Tomorrow ends with a bit of switchbacky goodness:</span><br /><br /><a style="font-family: arial;" onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiK6jsyhj_QNhkZlm4S5HmlQ9ysYJb1pJ5hNsZdHQaI1cXvk4xxrmHhm1TC6QkiEus6oSiVT3f2ccw7OqhoGEIM4T0n1F3qU7g8DuG4fZ_bHw1ZE2V3OC4TScUn4l6-Rb-6WD3cRhzrtWOO/s1600-h/arcalis.jpg"><img style="margin: 0pt 10px 10px 0pt; float: left; cursor: pointer; width: 320px; height: 256px;" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiK6jsyhj_QNhkZlm4S5HmlQ9ysYJb1pJ5hNsZdHQaI1cXvk4xxrmHhm1TC6QkiEus6oSiVT3f2ccw7OqhoGEIM4T0n1F3qU7g8DuG4fZ_bHw1ZE2V3OC4TScUn4l6-Rb-6WD3cRhzrtWOO/s320/arcalis.jpg" alt="" id="BLOGGER_PHOTO_ID_5356570058602386914" border="0" /></a><span style="font-family:arial;"><br /><br /><br /><br /><br /><br /><br /><br /><br /><br /><br /><br /><br /><br /><br />Hard for me to see anyone other than Contador winning.</span>Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-48260201353144060172009-07-08T18:22:00.000-07:002009-07-10T13:13:35.747-07:00TdF Stage 5<span style="font-family:arial;">Well done, Mr Voeckler! I wonder how many other people didn't know he had never won a TdF stage. Maybe he was given an honorary win at some point.<br /><br /><br />I agree with <a href="http://www.cyclingnews.com/news/aldag-frustrated-by-pelotons-tactics">whiny Rolf Aldag</a>. If you're sending a team to the Tour, it only makes sense to try winning a stage.<br /><br /><br />And finally, <a href="http://www.cyclingnews.com/news/wiggins-in-the-shape-of-his-life">Bradley Wiggins is saying</a>:</span><br /><blockquote><span style="font-family:arial;"></span>[T]here are so many different levels at this Tour: there are 100 blokes who are good and 100 blokes who are not so good, and it's all mixed up.</blockquote><span style="font-family:arial;">Is he saying it's a <span style="font-style: italic;">peloton at two speeds</span>, and he's in the fast group? My goodness, that team is suspicious. </span>Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0tag:blogger.com,1999:blog-3420445121485823034.post-62587205267731194952009-07-08T17:20:00.000-07:002009-07-10T13:12:54.282-07:00TdF Stage 4<span style="font-family:arial;">In general, team time trials are only slightly less boring to watch than individual time trials (unless Rasmussen is involved). But this was pretty interesting. Two points:<br /></span><span style="font-family:arial;"><br /></span><ul><li><span style="font-family:arial;">I think bike handling is a respectable skill and should be rewarded. I was happy to see Paris-Roubaix decided by <span style="font-family:arial;">bike handling and I don't mind it being part of the</span> Tour.</span></li><br /><li><span style="font-family:arial;">Garmin had four guys working for the majority of the race, compared to almost twice that many for most teams. And yet they beat almost everyone, finishing second by only 18 seconds. That's an amazing effort. One might even say </span><span style="font-style: italic;font-family:arial;" >superhuman!</span><span style="font-family:arial;"> </span><span style="font-family:arial;"> </span><span style="font-family:arial;"> </span></li></ul>Pendletonhttp://www.blogger.com/profile/11007470830120510626noreply@blogger.com0