Tuesday, 20 January 2015

A bunch of papers

The beginning of the new year has been particularly busy, as I'm working on several interesting projects. On the bright side, some of these are starting to give their fruits and, coincidentally, in the last few days we've had a few papers finalised (ie published, accepted for publication or submitted to the arxiv in an advanced status).

The first one has been published in Cost Effectiveness and Resource Allocation (the open access version is here). I've been involved in this paper with colleagues at UCL. The paper is an economic evaluation of an interesting and rather complex community trial conducted in Malawi, a country with particularly low life expectancy and high rates of HIV. In the paper, we did most of the economic analysis using BCEA.

The second one has also just been published in Pharmacoeconomics and is an "educational" piece that I co-wrote with several colleagues at UCL. I think this too was an interesting piece of work, in that we tried to focus on several statistical issues that are of concern in many economic evaluations $-$ the idea that health economics is in many ways inextricably related to statistics is of course one of my pieces de resistance (I guess I've showed off enough complicated words for a post...).

The third one is the RDD paper, which we had submitted ages ago to Statistics in Medicine. I had a brilliant experience with the Structural Zeros paper $-$ I submitted the first version in August and the paper was online by November. This time around, we had to struggle a lot more (apparently they couldn't find suitable reviewers, then the reviews arrived but took some time, then we responded to the comments $-$ long story short, it's been almost one year). Anyway, finally, they seem to have accepted the paper (which we previously arxived a similar version); we need a couple more changes and we should be good to go (I hope I'm not jinxing it!).

Finally, the last one is part of one of my PhD students (technically, I'm only second-supervising him). In fact, the paper develops a nice Bayesian non-parametric model to perform clustering and model selection simultaneously. We developed the model to handle a real clinical dataset, which records data on patients with lower urinary tract infection. I only knew little about Bayesian NP before working on this, so it was a nice opportunity. William has done a very good job in sorting this out and we have also submitted the full paper to Statistics in Medicine (hopefully, we'll get a quick turnaround!).

Monday, 19 January 2015

DIA Joint Adaptive Design and Bayesian Statistics Conference

This is my first real contribution to the ISBA Section on Biostatistics and Pharmaceutical Statistics, in my new role of secretary. Our section has formally endorsed this very interesting conference $-$ the timeline is very short, as the conference will start on February 11th.

The conference has two very interesting tutorials (the first one on Bayesian Methods for Drug Safety Evaluation and Signal Detection, given by David Olhsen and Amy Xia; and the second one on Use of Historical Data in Clinical Trials, taught by Heinz Schmidli). There'll also be other interesting talks and discussion.

In the Biostats/Pharma section, we'll try to organise (either directly, or through endorsement, or some other form) a few similar meetings (Bayes Pharma is of course another interesting one). 

Monday, 12 January 2015

No-go zones?

I think this is really serious: that's the story of an "expert" on terrorism, who's commented on Fox News about the situation in Europe. If you really want to laugh about it, you just take this at face value, as, erm, let's just say, "somebody" who just said a bunch of nonsense on TV.

But I think there's so much more to it, particularly after hearing an interview with this guy (on the BBC Radio 4 programme PM $-$ I'm sure there'll be a link to today's broadcast shortly). The point is that this guy is kind of defending himself by just saying that he's sorry for not having done his homework properly, and just taking as true oral information coming from some sources. He said he had used these sources in the past, so effectively didn't even bother to stop and think whether he was making a fool of himself.

It's astounding that you can get away, on TV, by just making a claim such as "in London, Muslim religious police beat anyone who doesn't dress according to Muslim, religious Muslim attire" (which is just not true!) without pointing out to any hard evidence to substantiate your claim. 

In fact, I think there's a wider problem: when I was reading the book on Enron scandal, I had noticed that quite often the authors stated something like "a source says...", or "one of the people involved in this said...", without going in to the details (NB: I should say that this is not an attack to the authors of that book $-$ just the way it is, I guess). Of course, in our line of work, you could never do anything like that. I suppose it may be time the media start following suit...

Monday, 5 January 2015


Today is the first day of my sabbatical term, which I have asked (and obtained) so I could work on writing the books $-$ the BCEA book had been commissioned for a while and it's not far from a nearly finalised version, while the 2nd edition of Bayesian Approaches to Clinical Trials and Health Care Evaluation (which I am co-authoring with relatively minor and unknown researchers, as some may say...) is slightly behind.

Although it wasn't a formal sabbatical, on another occasion I spent a relatively large chunk of time working from home $-$ that's back when I was working in Siena and was also collaborating to the write up of another book. Back then, it was a strange experience $-$ after a while, I wasn't really enjoying very much not going to the office every day and was actually, actively missing it. Part of it is probably that the office was bang in the middle of Siena (and this was the view)...

This time it will be probably quite different, as I will still go to the office once/twice a week (mostly for meetings and stuff), but at least I got to dump (in the friendliest way possible, I hope!) most of my teaching to Aidan, who's been so kind to cover for me, this term. Anyway, I keep telling people that, because I'm on sabbatical, this term will be good as I'll be able to do lots of work $-$ I just hope I've not already overbooked (pun intended!) myself...

Sunday, 4 January 2015

Advance notice

As a few people have already inquired (more or less formally) about the next edition of the short course on Bayesian Methods in Health Economics, I thought it would be nice to start the new year off by posting some details on that...

While we're exploring avenues for world domination and as yet unjustified expansion that will possibly, eventually lead us to be the Enron of statistical education (I'll post about this when I have a less confused idea about it), we know that the next edition will be held at the MRC Biostatistics Unit in Cambridge on 24-26 November, later this year (that's Tuesday to Thursday).

I think we'll very much keep the same format, which I think has worked fairly well in the two previous editions $-$ with only very minor tweaks, particularly to the first and, possibly last two talks. We may revise some of the practicals, but again I think that they kind-of work well, so I'm not sure about that... We'll re-advertise more widely closer to the event, but you can put the dates in your diary (and yes: we did have to book the room so much in advance!).

Wednesday, 31 December 2014

New Year's fireworks

I have to thank Xi'an for this: as he rightly guesses in his comments to my original post (earlier this year $-$ well, for a few more hours still, at least!), my spam filter has worked a treat and if it weren't for him, I would have missed this incredible spectacle of pre-new Year's eve fireworks by Vincent Granville

In his blog, Granville takes on evil statisticians who constantly talk down on data scientists. He seems to have a particular grudge with Andrew Gelman, who, Granville says, is our "leader and only influencer". 

What I found particularly amusing (and I needed this, as I'm still fighting off a nasty flu, so thank you very much, Xi'an!) is the bit when VC boasts about not even reading Gelman's publications, as they "are disseminated to a very small audience in obscure journals that pretty much no mainstream people read".

Now: I know that Google Scholar cannot be taken at face value and there are better ways of measuring how influential one's research or publications are. But you can see Gelman's numbers here. 14,000 citations on one book are sort-of suggesting a little less obscurity than that, I'd say. Sadly, Google Scholar can shed no light on how "mainstream" the people citing him are...

Happy new year!!

Tuesday, 16 December 2014


It's nearly the Christmas break and as I was writing the previous post (on our Workshop on cost-effectiveness thresholds), I just noticed the post-counter in the blog archive. While the decrease in the number of posts from 2012 to 2013 was minimal (98 to 90), this year I seem to have been fairly lazy (at posting, that is...) $-$ unless I frantically manage to write about 20 posts from now to the end of the year, of course, which I don't quite see happening...

I guess most of the "missing posts" (ie those I haven't written) are probably those about actually doing some stats $-$ which are fun to do, but also take time. For instance, the other day I found some data on Italian MPs attendance which I thought would be interesting to analyse, but then didn't have the time to actually do it. So, new year's resolution (slightly in advance): try and make some time for this.

By the way: while the one in the picture is not my cat, that's exactly what he's been doing pretty much for the whole day...

NICE and the cost-effectiveness thresholds: Can good intentions compensate for bad practice?

That's the title of the workshop we've held at UCL yesterday (I'd mentioned it in a previous post). I think it went remarkably well (OK $-$ as I'm the organiser, I may be over-enthusiastic, but I really think it was a very good day!). Despite the fact that we've purposely limited advertisement, as we wanted to keep it simple to start with, we had a good turnout (about 30 people!).

All the speakers have agreed to make their talks available (which are here) and we've agreed that I'd try and summarise some of the main discussion points coming out of the presentations, so that we can perhaps move this forward, somehow. Here's my attempt (any inconsistency is obviously due to my poor recollection of the events!):
  1. Matt gave a very interesting talk focussing on the international comparison (of 5 main countries); he started by discussing the different approaches at defining the cost-effectiveness threshold (which was also reprised by James in the second talk). Then, he moved to give some evidence of how in general agencies (with the notable exception of NICE) tend to resist a clear definition of what the threshold is. He presented some interesting work that has tried to elicit the underlying value of the threshold by analysing a set of past decisions in the Australian setting. What I found interesting of the Dutch case is the idea that the threshold can be set as a function of the disease severity, which made me think of potential links to my talk on risk aversion (see later). One of the interesting points of his discussion of NICE's case is the idea of a range of thresholds (rather than an absolute value). We had a little less time to discuss France and Japan, which are all interesting cases anyway (France because the formal use of cost-effectiveness methods is newly established; Japan because of the mix of private and public funding as well as the recent move towards formal inclusion of cost-effectiveness considerations, mainly driven by the need of capping health expenditure). Probably, in the end, the main message is to wonder whether strict adherence to a very specific guidance (such as NICE's) is a good idea. As I was listening to the talk, I also thought that in a more expanded version of this, it would be perhaps interesting to look at North America.
  2. James's talk concentrated on different ways of determining the threshold, effectively (and quite amusingly) pointing out at the many flaws in virtually all of these (which in itself is evidence of how complex this problem is!). I found the discussion of the "historical precedent" method fascinating and the fact that this is actually used in some countries is kind of strange. Then he gave a very clear account of the "budget exhaustion" method, linking back to some applied examples (eg Netherlands and Ireland), showing the incredible range of variability present in many cases. I found his comments on why it may not be ideal to have a range of thresholds very sensible and while I was listening to the talk I kind of changed my mind about that 2 or 3 times... James also linked to the Claxton et al paper on value based assessment and its potential impact in setting some plausible values. I think this is a very important piece of work that we ought to consider very careful (if we are to discuss the problem in general terms). Finally, I think another very good point that James made is about the idea that in fact, all this procedure is dynamic --- setting a threshold based on current budget and then ranking interventions accordingly may need continuous revisions, which in turn may also lead to modifications of the threshold itself, to add to the complexity of the problem.
  3. I gave arguably the most confused talk of the lot. The first part discussed the meaning of PSA and what are the quantities that should be involved (and why so). This, I feel, is something that I have thought about quite a lot and I think I have a clear idea of what the problem is and how it should be tackled. I was trying to lead this back to the main issue of selecting the threshold but I think Chris Jackson's comment that, as statisticians we kind of take that for granted, sums up nicely the fact that it is difficult to embed this aspect in the statistical model. In the second part of my talk, I have tried to (still confusingly) address the problem of risk aversion and its implications to PSA in general and the threshold in particular. I have presented a couple of potential utility functions that could be used to include a risk aversion parameter, together with a brief description of their characteristics. I liked James' comment that while difficult to think about it, it would be possible to elicit reasonable values of the risk aversion parameter, which in turn would make the analysis more applicable. While I was talking (and relating to the previous presentations) I also thought that may be risk-propensity is just as important (eg in cases where a disease is perceived to be so important as to be granted some sort of special treatment). Also, Matt asked the very fundamental question of whether decision-makers should be risk averse, in the first place $-$ I think nobody knew the answer to that one...
  4. Finally, Mike presented to us his model addressing the impact of the several assumptions underlying the standard definition of the threshold. The model attempts at transforming this into some more complex object, embedding the ideas of "value" and "information" available to the (different) decision-makers. The model is quite complex, but Mike explained it brilliantly (one of my PhD students came to me after the talk still shaking with the excitement!). The idea is to try and include formally several aspects that characterise the fixed assumptions underlying the definition of the threshold (eg the divisibility of the technologies involved, the level of information available to the decision-makers and, crucially, the fact that the process is driven by sequential decisions made, in general, by a set of actors, rather than the vaguely defined "decision-maker"). Model complexity aside (which I quite liked), the results are also very interesting and show that in fact, the straight line in the cost-effectiveness plane can be thought of as a special case of a more general threshold which can assume very different shapes (including kinks and steps). The very interesting point made by Mike is that if you start accounting for value and information, then you (read: NICE) may find yourself in the situation where it may be cost-effective to replace an existing technology with something that, while non cost-effective in the conventional sense is still a better option. This is kind of intuitive, but the model actually formalises the intuition and produces thresholds that have suitable shapes to accommodate this situation.
As I said, I don't really know where we go from here. But it was a very interesting start to this discussion!