July 13, 2014

100% accurate medical testing

The Wireless has a story about a fatal disease where there’s an essentially 100% accurate test available.

Alice Harbourne has a 50% chance of Huntington’s Disease. If she gets tested, she will have either a 0% or 100% chance, and despite some recent progress on the mechanism of the disease, there is no treatment.

July 11, 2014

Another prostate cancer study

Today’s prostate cancer risk factor, per the Herald, is vasectomy. The press release is here; the paper isn’t open-access.

This is a much more reliable study than the one earlier in the week about cycling, and there’s reasonable case that this one is worth a press release.

In 1986, the researchers recruited about 50000 men (health professionals: mostly dentists and vets), then followed them up to see how their health changed over time.  This research involves the 43000 who hadn’t had any sort of cancer at the start of the study. As the Herald says, about a quarter of the men had a vasectomy, and there have been 6000 prostate cancer diagnoses. So there’s a reasonable sample size, and there is a good chance you would have heard about this result if no difference had been found (though probably not via the Daily Mail)

The relative increase in risk is estimated as about 10% overall and about 20% for ‘high-grade’ tumours, which is much more plausible than the five-fold increase claimed for cycling.  The researchers had information about the number of prostate cancer tests the men had had, so they can say this isn’t explained by a difference in screening — the cycling study only had total number of doctor visits in the past year. Also, the 20% difference is seen in prostate cancer deaths, not just in diagnoses, though if you only consider deaths the evidence is borderline.  Despite all this, the researchers quite rightly don’t claim the result is conclusive.

There are two things the story doesn’t say. First, if you Google the name of the lead researcher and ‘prostate cancer’, one of the top hits is another paper on prostate cancer (and coffee, protective). That is, the Health Professionals Followup Study, like its sister cohort, the Nurses Health Study, is in the business of looking for correlations between a long list of interesting exposures and potential effects. Some of what it finds will be noise, even if it appears to pass sanity checks and statistical filters. They aren’t doing anything wrong, that’s just what life is like.

Second, there were 167 lethal prostate cancers in men with vasectomies. If the excess risk of 20% is really due to vasectomy, rather than something else, that would mean about 27 cancers caused by 12000 vasectomies. Combining lethal and advanced cases, the same approach gives an estimated 38 cases from 12000 vasectomies. So, if this is causation, the risk is 2 or 3 serious prostate cancers for every 1000 vasectomies. That’s not trivial, but I think it sounds smaller than “20% raised risk”.

July 10, 2014

Summaries of income

I don’t want to get into the general business of election fact-checking, but we have a Stat-of-the-Week  nomination for a statement that is (a) about a specifically statistical issue at the high-school level, and (b) unambiguously wrong.  From Richard Prebble’s “The Letter”:

 Cunliffe is basing Labour’s election campaign around the claim that inequality is growing. Fact check: inequality is falling and New Zealand remains a very equal country. The claim that around a quarter of a million children are in poverty is dubious, to say the very least. Cunliffe says households in poverty have less than 60 percent of the medium income after housing costs. If Bill Gates came to live in New Zealand, the medium income of the country would rise and, according to that logic, more children would be in poverty.

David Cunliffe, as you presumably know, talked about the median, not “medium”; the use of a fraction of median income as a relative poverty threshold is very common internationally. The reason for using the median is precisely that the median income of the country would not rise if a few billionaires were added to the population. The median, the income of the household in the middle of the income distribution, is very insensitive to changes in or additions of a few values. That’s what it’s for.

While I’m writing, I might as well mention the inequality statistics.  Mr Cunliffe isn’t making up his figures on children in poverty; they can be found in the 2014 Household Incomes Report from the Ministry of Social Development [update: that figure is 260000, which matches what The Letter reported was said, but the actual speech said 285000]. The report also gives trends in the Gini index of inequality and in the proportion of income spent on housing.  StatsNZ gives trends in the ratio of 80th to 20th percentile of income, before and after housing costs. The details of trends in inequality depend on how you measure it, but by these measures it is neither falling, nor notably low internationally.

July 9, 2014

Would I have heard if the results were different?

The story about cycling and prostate cancer in the Herald (or the Daily Mail) is a good opportunity to look at some of the rules of thumb for deciding which stories to read or believe:

Firstly, would you have heard if the results were the other way around? Almost certainly not: prostate cancer wasn’t the main point of this study, and there wasn’t a previously-suspected relationship.

Second, for cancer specifically, is this mortality or diagnosis data? That is, are we seeing an increase in detection or in cancer? This is diagnosis data; so it could be just an increase in detection. The researchers were confident it wasn’t, but we must remember the immortal words of Mandy Rice-Davies

Thirdly, what sort of study is it? Obviously it can’t be experimental, but a good study design would be to ask people about cycling (or even better, measure cycling) and then see whether it’s the bike fanatics who develop cancer. This study was a self-selected survey of cyclists, getting self-reported data about past cycling and past diagnosis of prostate cancer. It’s a fairly extreme sample, too: half of them cycle more than 5.75 hours per week.

Fourth, how strong is the evidence of association, and what sort of sample size are we looking at? The association is just barely statistically significant (p=0.046 in one model, p=0.025 in a second), and there are only 36 prostate cancer cases in the sample.  It’s pretty borderline.  The estimated relative risk is huge, because it has to be given the sample size, but the uncertainty range is also huge. The confidence interval on the relative risk of 5 reported by the Herald goes from 1.5 to 18.

Fifth, what does previous research say? This is in the story

‘To the best of our knowledge, this is the first study to demonstrate an association between prostate cancer and cycling, so there are no studies hypothesizing a pathophysiological mechanism for such a link.’

Sixth, what do other experts think? We don’t know. The closest thing to an independent comment is this in the press release

“Physicians should discuss the potential risks and health benefits of cycling with their patients, and how it may impact their overall health,” says Ajay Nehra, MD, Editor-in-Chief of Journal of Men’s Health and Chair, Department of Urology, Director, Men’s Health, Rush University Medical Center, Chicago, IL.

He could have said that without reading the paper.

In summary, there’s borderline evidence from a weak study design for a sensational finding that isn’t supported by any prior evidence. This is fine as research, but it shouldn’t be in the headlines.

You can read the research paper here for the next month, and the journal press release here.

Recycling

In March, I wrote

The Herald has a story about a potential blood test for dementia, which gives the opportunity to talk about an important statistical issue. The research seems to be good, and the results are plausible, though they need to be confirmed in a separate, larger sample before they can really be believed. …

 But it’s the description of the accuracy of the test that might be misleading.

There’s a Herald story today about new test; the same comments apply, except that the research paper is open-access

NRL Predictions for Round 18

Team Ratings for Round 18

The basic method is described on my Department home page. I have made some changes to the methodology this year, including shrinking the ratings between seasons.

Here are the team ratings prior to this week’s games, along with the ratings at the start of the season.

Current Rating Rating at Season Start Difference
Roosters 6.97 12.35 -5.40
Sea Eagles 6.56 9.10 -2.50
Rabbitohs 5.76 5.82 -0.10
Warriors 3.11 -0.72 3.80
Cowboys 2.96 6.01 -3.10
Panthers 2.96 -2.48 5.40
Bulldogs 2.50 2.46 0.00
Broncos 2.43 -4.69 7.10
Storm 0.40 7.64 -7.20
Knights -1.29 5.23 -6.50
Titans -3.51 1.45 -5.00
Dragons -4.01 -7.57 3.60
Eels -6.44 -18.45 12.00
Sharks -6.52 2.32 -8.80
Wests Tigers -6.64 -11.26 4.60
Raiders -7.04 -8.99 2.00

 

Performance So Far

So far there have been 122 matches played, 68 of which were correctly predicted, a success rate of 55.7%.

Here are the predictions for last week’s games.

Game Date Score Prediction Correct
1 Bulldogs vs. Sea Eagles Jul 04 23 – 16 -1.10 FALSE
2 Dragons vs. Cowboys Jul 05 27 – 24 -3.80 FALSE
3 Roosters vs. Sharks Jul 05 28 – 30 22.20 FALSE
4 Wests Tigers vs. Panthers Jul 06 10 – 26 -2.60 TRUE
5 Rabbitohs vs. Titans Jul 07 10 – 14 17.60 FALSE

 

Predictions for Round 18

Here are the predictions for Round 18. The prediction is my estimated expected points difference with a positive margin being a win to the home team, and a negative margin a win to the away team.

Game Date Winner Prediction
1 Sea Eagles vs. Wests Tigers Jul 11 Sea Eagles 17.70
2 Warriors vs. Eels Jul 12 Warriors 14.10
3 Storm vs. Bulldogs Jul 12 Storm 2.40
4 Titans vs. Raiders Jul 13 Titans 8.00
5 Sharks vs. Knights Jul 13 Knights -0.70
6 Panthers vs. Broncos Jul 14 Panthers 5.00

 

Super 15 Predictions for Round 19

Team Ratings for Round 19

The basic method is described on my Department home page. I have made some changes to the methodology this year, including shrinking the ratings between seasons.

Here are the team ratings prior to this week’s games, along with the ratings at the start of the season.

Current Rating Rating at Season Start Difference
Waratahs 8.40 1.67 6.70
Crusaders 7.96 8.80 -0.80
Sharks 4.75 4.57 0.20
Hurricanes 2.89 -1.44 4.30
Stormers 2.85 4.38 -1.50
Brumbies 2.13 4.12 -2.00
Blues 1.82 -1.92 3.70
Bulls 1.76 4.87 -3.10
Chiefs 1.63 4.38 -2.70
Highlanders -2.30 -4.48 2.20
Cheetahs -3.62 0.12 -3.70
Reds -3.72 0.58 -4.30
Force -3.84 -5.37 1.50
Lions -5.31 -6.93 1.60
Rebels -8.41 -6.36 -2.10

 

Performance So Far

So far there have been 113 matches played, 73 of which were correctly predicted, a success rate of 64.6%.

Here are the predictions for last week’s games.


Game Date Score Prediction Correct
1 Chiefs vs. Hurricanes Jul 04 24 – 16 0.20 TRUE
2 Lions vs. Rebels Jul 04 34 – 17 5.70 TRUE
3 Crusaders vs. Blues Jul 05 21 – 13 8.80 TRUE
4 Force vs. Reds Jul 05 30 – 20 1.30 TRUE
5 Stormers vs. Bulls Jul 05 16 – 0 1.90 TRUE
6 Cheetahs vs. Sharks Jul 05 27 – 20 -7.70 FALSE
7 Waratahs vs. Highlanders Jul 06 44 – 16 12.90 TRUE

 

Predictions for Round 19

Here are the predictions for Round 19. The prediction is my estimated expected points difference with a positive margin being a win to the home team, and a negative margin a win to the away team.

Game Date Winner Prediction
1 Blues vs. Chiefs Jul 11 Blues 2.70
2 Brumbies vs. Force Jul 11 Brumbies 8.50
3 Bulls vs. Rebels Jul 11 Bulls 14.20
4 Crusaders vs. Highlanders Jul 12 Crusaders 12.80
5 Reds vs. Waratahs Jul 12 Waratahs -9.60
6 Lions vs. Cheetahs Jul 12 Lions 0.80
7 Stormers vs. Sharks Jul 12 Stormers 0.60

 

July 7, 2014

Stat of the Week Competition: July 5 – 11 2014

Each week, we would like to invite readers of Stats Chat to submit nominations for our Stat of the Week competition and be in with the chance to win an iTunes voucher.

Here’s how it works:

  • Anyone may add a comment on this post to nominate their Stat of the Week candidate before midday Friday July 11 2014.
  • Statistics can be bad, exemplary or fascinating.
  • The statistic must be in the NZ media during the period of July 5 – 11 2014 inclusive.
  • Quote the statistic, when and where it was published and tell us why it should be our Stat of the Week.

Next Monday at midday we’ll announce the winner of this week’s Stat of the Week competition, and start a new one.

(more…)

July 5, 2014

Once is accident, twice is coincidence

Back in 2010, a piece in Slate pointed out that a country’s success in the 2010 and 2006 World Cup knock-out rounds was strongly correlated with the proportion of the population infected by Toxoplasma gondii.  In 2010, Toxoplasma seroprevalence predicted all eight knockout-round wins; in 2006 it predicted seven of eight.

Toxoplasma, in case you weren’t introduced when you met it, is a single-celled organism that can live, and reproduce asexually, in pretty much any warm-blooded animal, but can only reproduce sexually in the guts of cats. That’s not the interesting part. The interesting part is that in rodents the parasite has effects on the brain, making the animal less cautious and more likely to end up in the gut of a cat. There’s some evidence Toxoplasma also has effects on human behaviour, though that’s still controversial.

Now, in 2014, I see a Tweet from Australian biologist Michael Whitehead

So, three times is enemy action?

There are good reasons to be sceptical: the football rankings haven’t changed all that much since 2006, so this isn’t really three independent tests. Also, the seroprevalence data is for the countries as a whole, not for the team members.  Still, in contrast to the predictions using Octopus vulgaris in the last World Cup, it’s not completely out of the question that there could be a real effect.

July 4, 2014

Briefly

  • “There are advantages and disadvantages for reporting the median, but over time it has become common practice worldwide to report market benchmark prices as the median.” from a good explanation at interest.co.nz
  • The maximum, on the other hand, is rarely a good summary on its own. This Herald story on slow licence suspensions is not an exception. A median or 90th percentile, or proportion taking longer than a reasonable duration would have been good additions. Also, how many licence suspensions are there for accumulated demerits?
  • TheWireless is having a ‘theme’ on Risk. It’s pretty much non-quantitative, which I think misses something, but they aren’t trying to draw unwarranted generalisations from qualitative data. I liked the story on yellow-stickered earthquake-risk buildings; an interesting counterpoint is Eric Crampton’s post on house-hunting in Wellington.
  • “Measure twice; cut once” is the old saying. It’s good that the government’s welfare reform program is being evaluated. Not so good that the evaluation plan is secret even under OIA.
  • Interestingly, the retracted and recently republished paper on GMOs and Roundup  (previous StatsChat coverage)wasn’t peer-reviewed. Or, rather, it wasn’t peer-reviewed again — the journal decided that the initial review before the retraction was enough. This was not made very clear in the paper or press material.