This article has a correction
- Chris Parker, medical statistician1,
- Carol Coupland, associate professor in medical statistics2,
- Julia Hippisley-Cox, professor of clinical epidemiology and general practice2
- 1Nottinghamshire County Teaching Primary Care Trust, Hucknall Health Centre, Hucknall, Nottingham NG15 7JE
- 2Division of Primary Care, University Park, Nottingham NG2 7RD
- Correspondence to: J Hippisley-Cox
- Accepted 22 June 2010
Objective To determine whether antipsychotic drugs are associated with an increased risk of venous thromboembolism, and to examine risks by type of antipsychotic, potency, and dose.
Design Population based nested case-control study.
Setting The UK QResearch primary care database.
Participants Patients (cases) with a first ever record of venous thromboembolism between 1 January 1996 and 1 July 2007; each was matched with up to four controls by age, calendar time, sex, and practice.
Main outcome measures Odds ratios for venous thromboembolism associated with antipsychotic drugs adjusted for comorbidity; concomitant drug exposure.
Results There were 25 532 eligible cases (15 975 with deep vein thrombosis and 9557 with pulmonary embolism) and 89 491 matched controls from a study population of 7 267 673. Individuals prescribed antipsychotic drugs in the previous 24 months had a 32% greater risk of venous thromboembolism than non-users, despite adjustment for potential risk factors (odds ratio 1.32, 95% confidence interval 1.23 to 1.42). Patients who had started a new drug in the previous three months had about twice the risk (1.97, 1.66 to 2.33). The risk was greater for individuals prescribed atypical rather than conventional drugs (adjusted odds ratio 1.73, 1.37 to 2.17, for atypical drugs; 1.28, 1.18 to 1.38, for conventional drugs). It also tended to be greater for patients prescribed low rather than high potency drugs (1.99, 1.52 to 2.62, for low potency; 1.28, 1.18 to 1.38, for high potency). The estimated number of extra cases of venous thromboembolism per 10 000 patients treated over one year was 4 (3 to 5) in patients of all ages and 10 (7 to 13) for patients aged 65 and over.
Conclusions There is an association between use of antipsychotic drugs and risk of venous thromboembolism in a large primary care population. The increased risk was more marked among new users and those prescribed atypical antipsychotic drugs.
Venous thromboembolism is an important and preventable cause of morbidity and mortality.1 Up to a quarter of affected patients die within a week,2 and almost a third of survivors experience long term effects.3 4 Some research has indicated that antipsychotic drugs, some of which are also widely prescribed for nausea, vomiting, and vertigo, might be associated with an increased risk of venous thromboembolism. Early case reports led to several studies, but the findings have been inconsistent and the possible association has received little attention.
One case-control study found a sevenfold risk of venous thromboembolism among current users of antipsychotic drugs,5 and another a 13-fold risk of death from pulmonary embolism.6 The studies were based on small numbers, excluded older people, and were conducted before the widespread use of atypical antipsychotic drugs. A hospital based case-control study found a 3.5-fold risk associated with antipsychotics, with a lower risk for atypical drugs than for conventional drugs.7 Two large cohort studies focused on people aged 65 or over who were taking antipsychotic drugs.8 9 One found a non-significant overall association,8 although the subgroup taking butyrophenone conventional drugs showed a 43% increased risk. The other, based on residents of nursing homes, found a doubled risk associated with atypical drugs but no effect with conventional drugs.9 A further cohort study focused on the atypical clozapine and found that current users had a fivefold risk of death from pulmonary embolism compared with past users.10
There are therefore grounds for concern, but considerable uncertainty remains on an association that, if proved, would have important implications. Several plausible biological mechanisms have been suggested, including enhanced aggregation of platelets, raised concentrations of anticardiolipin antibodies, and exacerbation of venous stasis.11
We investigated whether antipsychotic drugs are associated with an increased risk of venous thromboembolism and examined risks by type, potency, and dose, adjusting for comorbidity and concomitant drug exposure. We carried out a nested case-control study in a large population based primary care cohort in the United Kingdom.
The QResearch database version 16 (www.qresearch.org) holds the anonymised primary care clinical records of over 11 million people registered at any time in the past 16 years with 525 UK general practices. The database includes data on patients’ demographics, characteristics, consultations, symptoms, diagnoses, referrals, results of investigations, and prescribed drugs. It has been subject to detailed analyses including age-sex distributions, birth rates, death rates, and consultation rates and shows good correspondence with other sources.12 Recording of clinical diagnoses including psychotic illness13 and venous thromboembolism14 has good levels of completeness and accuracy in UK general practice databases.
The study population consisted of an open cohort of patients registered with participating practices between 1 January 1996 and 1 July 2007. Cases were all patients aged between 16 and 100 with a first ever record of venous thromboembolism (deep vein thrombosis or pulmonary embolism) during the study period, including postmortem diagnoses, identified from computer recorded Read codes (diagnostic codes) for venous thromboembolism recorded within the patient’s electronic record.
We used incidence density sampling to identify up to four controls for each case, matched by single year of age, calendar time, sex, and practice. All controls were alive and registered with the practice at the date of the first recorded diagnosis of venous thromboembolism in their matched case: this was the index date for each case and their controls. Controls had no diagnosis of venous thromboembolism in their record up to this date.
Potential participants (cases and controls) were excluded if they had less than 24 months of data before the index date or if area census data on socioeconomic status were missing from their record (usually temporary residents). Controls were excluded if they had had any prescriptions for warfarin before the index date, as this could be treatment for a previous, unrecorded venous thromboembolism. For cases, prescriptions for warfarin in the six weeks before diagnosis of venous thromboembolism did not trigger exclusion as they could plausibly be treatment for the index event itself. Cases in which there was any use of warfarin earlier than this six week period, however, were excluded.
We extracted demographic information including year of birth, sex, and fifths of Townsend score (a measure of socioeconomic status), together with the most recent body mass index (BMI) and smoking status before the index date. Individuals were coded as having schizophrenia, bipolar disorder, and/or dementia if these diagnoses were recorded before the index date. We extracted data recorded before the index date on comorbid conditions that could increase the risk of venous thromboembolism (coronary heart disease, cardiac failure, stroke, cancer, inflammatory bowel disease, liver disease, varicose veins, gastrointestinal bleed, Parkinson’s disease, renal disease, asthma, diabetes, hypertension, hyperlipidaemia).1 15 16 For several events associated with increased risk in the short term (hip surgery, hip or lower limb fractures, acute infections, pregnancy),1 15 16 we extracted data for the six months before the index date. We also extracted information on whether there was any computer recorded evidence of a hospital admission in the preceding 31-183 days (admissions in the preceding 30 days might have been for the thrombosis itself) and included this in the analysis as a binary variable.
We assessed exposure to drugs on the basis of prescriptions on or before the index date. For the exposures of primary interest (antipsychotics, BNF, chapters 4.2.1 and 4.2.2), full prescription data were collected, including drug name, formulation, dose instructions, and date. For some drugs that increase risk of venous thromboembolism (oral contraceptives and hormone replacement therapy,16 17 tamoxifen18) and other commonly used drugs (statins, non-steroidal anti-inflammatory drugs (NSAIDs), aspirin, antimanics), we recorded only whether any prescriptions were issued in the 24 months before the index date.
The analyses were conducted in Stata version 9.2. Each individual was classified for exposure to antipsychotics as a current user (one or more prescriptions within three months before index date); recent user (four and 12 months before); past user (13 and 24 months before); or not exposed within the previous 24 months. We subclassified current users as new users (prescription for antipsychotic during the three months before the index date, after at least 12 months with no prescription for that antipsychotic) or as continuing users. Exposure during the past 24 months was also classified according to type (conventional only, atypical only, both) and potency (high potency only, low potency only, or both, with high potency defined as equivalent dose of more than 100 mg chlorpromazine). When possible, we calculated the daily dose equivalent of each script and classified individuals as low dose only (≤25% of maximum recommended in BNF), medium/low dose (under maximum), or any high dose (maximum and over). The number of different antipsychotic drugs prescribed to each individual, and the number of scripts each received, were ascertained.
Patients with more than one mental health indication were categorised according to a hierarchy: schizophrenia, bipolar disorder without schizophrenia, and dementia without schizophrenia or bipolar disorder. BMI values outside the range 15-50 were treated as missing.
We undertook multiple conditional logistic regression to estimate the odds ratio (with 95% confidence interval) for risk of venous thromboembolism associated with any use of antipsychotics in the 24 months before the index date. The odds ratios are unbiased estimates of rate ratios as we used incidence density sampling to select controls. Odds ratios and 95% confidence intervals were estimated for subgroups based on timing, type, potency, duration, and dose of antipsychotics and for individual antipsychotics with a sufficiently high level of use. We compared the risk for current antipsychotic users (prescriptions in the previous three months) between new users and continuing users.
The analyses were adjusted for socioeconomic status, the comorbidity and drug variables listed above, and the number of complete months of data before the index date. Adjustment was not initially made for smoking or BMI, because of missing data, but after we imputed missing values we repeated the analyses with these variables added to the models. We tested for interactions between use of antipsychotics and age, sex, socioeconomic status, BMI, smoking status, and category of mental health.
We repeated the main analyses twice. The first analysis was restricted to idiopathic venous thromboembolism, and cases and controls with major risk factors for venous thromboembolism (previous cancer, coronary heart disease, stroke, congestive cardiac failure; hip surgery, hip or lower limb fracture, or pregnancy within the previous six months) were excluded. We then excluded individuals with diagnoses of manic depression or schizophrenia.
We calculated the numbers needed to harm per year by applying the incidence of first venous thromboembolism and using the adjusted odds ratio from current use of antipsychotic drugs by type, potency, timing, quantity, and dose. We calculated this separately for all patients aged 16 and over and for those aged 65 and over. We also estimated the number of additional cases of venous thromboembolism expected per 10 000 treated patients per year. We calculated approximate 95% confidence intervals for these estimates.
Multiple imputation was carried out with the Stata ICE programs,19 with all the main analysis variables entered in the prediction equation, together with the matching variables of age, sex, and index year. Five imputed datasets were created and combined for analysis.
A total of 453 QResearch practices met the inclusion criteria. The study population consisted of 7 267 673 patients with 26 702 594 person years of observation. We identified 31 612 incident cases of venous thromboembolism, an overall crude rate of 118 per 100 000 person years (table 1⇓, figure⇓). Incidence rates among women rose from 32 per 100 000 person years at age 16-24 to 403 at age 85 and over; for men the corresponding numbers were 13 and 380.
Cases and controls
The 31 612 cases were matched to 125 559 controls. We excluded 6016 cases (19%) and 17 595 controls (14%), principally because they had less than 24 months of data before the index date (table 2⇓). Cases without any matched controls, and controls without a matched case, were then excluded, leaving 25 532 eligible cases (15 975 with deep vein thrombosis and 9557 with pulmonary embolism) and 89 491 eligible controls.
Table 3⇓ shows that cases were more likely than controls to have a high BMI and marginally more likely to live in an area of socioeconomic deprivation. Overall, the prevalence of relevant mental health conditions was 0.4% for schizophrenia, 0.3% for bipolar disorder, and 1.0% for dementia; eight cases and 31 controls had more than one of these. All the suggested risk factors for venous thromboembolism were more prevalent among cases than controls. For many conditions the difference was small, but there was a greater disparity for cancer and for recent fractures, hip surgery, and acute infections. The use of drugs that could increase risk of venous thromboembolism (oral contraceptives, hormone replacement therapy, tamoxifen) was more common among cases than controls.
Risk associated with antipsychotics
New users—Table 4⇓ shows that 2126 cases (8.3%) and 4752 controls (5.3%) had received an antipsychotic drug in the previous 24 months. Overall, antipsychotic users had a 32% greater risk of venous thromboembolism than non-users (adjusted odds ratio 1.32, 95% confidence interval 1.23 to 1.42). Among antipsychotic users, 38% were classed as current users (a script within the previous three months): their increase in risk was 56% (1.56, 1.39 to 1.75) compared with 36% for recent users (1.36, 1.20 to 1.54); for past users the risk was not significantly increased (1.04, 0.91 to 1.18). Among current antipsychotic users, 15% had started a new drug within the three months before the index date, and this group of new users showed a greater increase in risk than continuing users (1.97, 1.66 to 2.33, for new users; 1.29, 1.11 to 1.51, for continuing users).
Dose and type of antipsychotic—Table 4⇑ also shows that 87% of controls (4156/4752) of those prescribed antipsychotics had received only conventional drugs, 90% (4293/4752) had received only high potency drugs, and 51% (2404/4752) had received only a single prescription in the past 24 months; 6% (256) of those taking antipsychotics had changed to a different antipsychotic drug, and 1% (44) had received the maximum daily dose (or greater) of any antipsychotic.
Atypical v conventional—The risk was greater for individuals prescribed atypical drugs rather than conventional drugs (adjusted odds ratios 1.73, 1.37 to 2.17, for atypical; 1.28, 1.18 to 1.38, for conventional). It was also greater for patients prescribed low rather than high potency drugs (1.99, 1.52 to 2.62, for low potency; 1.28, 1.18 to 1.38, for high potency).
Number of prescriptions—The risk was marginally higher for individuals receiving more than 12 scripts in the past 24 months (adjusted odds ratio 1.44, 1.18 to 1.75), but those receiving only one had a significantly higher risk than those receiving none (1.32, 1.17 to 1.49).
Multiple drugs—Individuals prescribed two or more different antipsychotics had a greater risk than those receiving one (adjusted odds ratio 1.99, 1.49 to 2.65, for two or more drugs; 1.29, 1.20 to 1.39, for one drug).
Dose response and route of delivery—There was an apparent dose-response gradient, but individuals receiving medium or high dose scripts were also more likely to have received two or more different drugs and to have received atypical drugs rather than conventional drugs. When we added the number of drugs or the types of antipsychotics to the models there was no significant relation between dose and risk. Individuals receiving any scripts for injections (short acting or depot) showed a substantially higher risk (adjusted odds ratio 3.24, 2.04 to 5.17) than individuals receiving only oral drugs or suppositories (1.30, 1.20 to 1.40), but the numbers receiving injections were small.
Individual drugs—The most commonly prescribed drug was prochlorperazine, a high potency phenothiazine, which, in addition to treating psychosis, is commonly prescribed for nausea, vomiting, and vertigo. Prochlorperazine was prescribed for 5133 individuals, 75% of the total number who had taken antipsychotics, and for 3243 of these only a single script had been given (not tabulated). The other most commonly prescribed drugs were (in order) risperidone, haloperidol, olanzapine, chlorpromazine, trifluoperazine, and quetiapine (table4⇑). Separate odds ratios were estimated for exposure to these drugs, and the highest risks were found among individuals prescribed quetiapine, chlorpromazine, and haloperidol (adjusted odds ratio 2.81 (1.75 to 4.50), 1.77 (1.27 to 2.48), and 2.17 (1.55 to 3.02), respectively).
Sensitivity analyses—All the results were similar when we carried out analyses on individuals having excluded those with either schizophrenia or manic depression and having excluded those with major risk factors for venous thromboembolism (see appendix on bmj.com for full results).
Missing data—Data on smoking or BMI, or both, were missing for 6534 cases (25.6%) and 25 720 controls (28.7%): they were marginally more likely to be men and to live in an area of social deprivation than those with complete data. Odds ratios based on the imputed datasets, with smoking and BMI added to the adjustment variables, were close to the ones reported in table 4. For example, the overall odds ratios after imputation were 1.32 (1.24 to 1.40) for any antipsychotic use in the previous 24 months, 1.57 (1.43 to 1.72) for current use, 1.31 (1.19 to 1.46) for recent use, and 1.06 (0.95 to 1.18) for past use.
Confounding variables—Table 5 shows associations between risk of venous thromboembolism and the variables treated as potential confounders⇓, adjusted for each other and for any use of antipsychotics in the previous 24 months. We used the imputed datasets so that smoking and BMI could be included, and the table shows significant associations with risk for both these variables and socioeconomic status. Individuals with a diagnosis of dementia were at higher risk than those with schizophrenia, bipolar disorder, or none of these conditions. Most of the comorbidities and health events considered as potential risk factors were associated with a relatively small increase in risk, but there was a more than threefold increase associated with cancer and about a 13-fold increase associated with recent surgery or fractures. Individuals prescribed statins or aspirin in the past 24 months had a lower risk of venous thromboembolism, and those prescribed NSAIDS, oral contraceptives, hormone replacement therapy, tamoxifen, or anti-manics had a higher risk.
Interactions—Inclusion of interaction terms showed that the risk associated with antipsychotic use did not differ significantly between subgroups based on age, sex, socioeconomic status, BMI, or mental health category. Smokers had a higher risk related to antipsychotics than non-smokers: the adjusted odds ratio for antipsychotic use was 1.68 (1.45 to 1.93) among smokers and 1.27 (1.18 to 1.36) among non-smokers.
Numbers needed to harm—Table 6⇓ shows the numbers needed to treat to harm (NNH) for each category of antipsychotic use and the number of excess cases per 10 000 patients treated over a year for all patients aged over 16 and for those aged 65 years and older. For example, the NNH for any antipsychotic use in the past 24 months for patients aged 65 year and older was 1044 (795 to 1452); for new users in the past three months it was 344 (251 to 506); and for continuing users it was 1152 (655 to 3037). The corresponding numbers of excess cases of venous thromboembolism per 10 000 treated patients were 10 (7 to 13), 29 (20 to 40), and 9 (3 to 15), respectively.
Summary of main findings
This nested case-control study, based on the primary care clinical records of 115 000 people in the UK, found a 32% increased risk of venous thromboembolism for individuals prescribed antipsychotic drugs in the past 24 months. The increase in risk was 56% for individuals with any antipsychotic use in the past three months, and 97% for those who had newly started on an antipsychotic in the past three months. The absolute risks, however, were low, with an excess of four extra cases of venous thromboembolism per 10 000 patients treated over one year in patients of all ages, and 10 for patients aged 65 and over.
Methodological considerations and strengths
The study was based on the primary care records of a large and representative population cohort, avoiding bias from selection, non-response, or poor recall. General practices in the UK have good levels of accuracy and completeness in recording clinical diagnoses and prescribed drugs,20 and practices that contribute to primary care databases such as QResearch are thought to perform even better.21 The recorded clinical diagnoses of venous thromboembolism were not independently verified for the study, but in other studies the inclusion of “possible” venous thromboembolism cases did not alter the estimates based on confirmed cases.14 22 The observed incidence rate in our population was close to the 117 per 100 000 reported elsewhere.23
We adjusted for many potential confounding factors, including comorbidities, health events, and drugs. We carried out an alternative analysis restricted to individuals with no major risk factors for venous thromboembolism because inclusion of non-idiopathic cases can bias an estimate towards the null.24 We found little difference between the estimates based on all individuals (with adjustment for risk factors) and those based on individuals without major risk factors. Clinicians prescribing antipsychotics might have no current perception of thrombotic risk and therefore do not take an individual’s other thrombotic risk factors into account in the decision. We also carried out an analysis excluding patients with a recorded diagnosis of schizophrenia or manic depression and again found little difference compared with the estimates based on all individuals.
Another strength of our study was that we had access to full prescription details for all antipsychotics, including drug name and formulation, dose instructions, and dates. This enabled us to look in detail at characteristics of the drug exposure in relation to risk. It also allowed for a comparison of risk between new users and those who had taken antipsychotic drugs for longer. Inclusion of all prevalent users can underestimate risks, which occur soon after the start of treatment,25 and we found that the risk was highest in the early months, in common with some previous studies.5 9
As with all observational studies, our study has limitations. About a quarter of the records had missing data on smoking status or BMI, or both. By imputing values for these variables we were able to ascertain that including them in the models had almost no effect on the estimates. There might have been some confounding by indication as in most patients we were unable to ascertain the reason for prescription of antipsychotics. This makes it difficult to distinguish the effects of the drugs from any effect of the condition for which the drug was prescribed. When we excluded people with schizophrenia or manic depression from the analyses, however, the results were largely unchanged. There could still be residual confounding and the true effect of this group of drugs is probably less than has been estimated in this study.
Our study was based on patients who have access to primary care and excluded any individuals not registered with a general practitioner—for example, those in prisons. In the overall database the most deprived socioeconomic group is marginally under-represented. Even in a dataset of this size, the number of cases and controls exposed to individual drugs was often too small to estimate separate odds ratios. The possibility of residual confounding by indication cannot be ruled out and is a further barrier to making definitive risk comparisons between different drugs. Furthermore, our study wasn’t designed to and is not able to offer hypotheses on individual mechanisms in the genesis of thromboembolism for each drug. In terms of the overall findings, however, the fact that past users of antipsychotics showed no significantly increased risk argues that the drugs, rather than any fixed characteristics of those prescribed them, are responsible for the increased risk in current users.
We have presented numbers needed to harm and number of estimated excess cases per 10 000 patients, which provides useful information at a population level. With this study design we were unable to derive individualised measures of absolute risk, which would take account of the patient’s characteristics.
Study in context of previous studies
Previous studies on antipsychotics and venous thromboembolism had substantial differences in sample and design and produced highly disparate findings. Our case-control study found smaller estimated effects than previous case-control studies,5 6 7 but our sample was much larger and more representative of the general population. Previous cohort studies were based on more restricted samples (such as older people,8 residents of nursing homes,9 and current or past users of clozapine10). The two studies based on older people produced conflicting findings, one finding a significant association only for a subgroup using conventional drugs,8 and the other an association for atypical drugs but not conventional drugs (though the latter were being prescribed at lower doses than previously found to increase risk).9 In our study, the risk was higher for individuals taking atypical drugs than for those taking conventional drugs. Our findings of a greater risk for individuals prescribed low rather than high potency drugs, and of no dose-response relation that could not be explained by other characteristics of the drug exposure, agree with the limited previous evidence.5 6
What this study adds
Our study adds to the accumulating evidence of adverse health events associated with antipsychotic drugs. Recent studies have confirmed the significantly increased risks of serious events and mortality among people with dementia who are treated for behavioural problems with these drugs.26 27 Nearly all use of antipsychotics in the primary care population we studied was for conditions such as nausea, vomiting, and vertigo, implying that the message on potential risk could have important implications for general practice. An increased risk of stroke associated with antipsychotics has recently been confirmed in primary care patients.28 Our findings would need to be replicated on another database before changes in clinical practice can be recommended, and larger numbers would be required to estimate the risks associated with individual antipsychotics. A study such as a self controlled case series analysis, which controls for all known and unknown confounders provided they are constant over time,29 would strengthen the evidence we found.29
If other studies replicate these findings, antipsychotic drugs should be used more cautiously for nausea and agitation, etc, especially among patients at high risk of thromboembolism. Patients need information on the balance of risks and benefits of these drugs before they start treatment. To achieve this, we need new algorithms designed to estimate a patient’s absolute risk of thromboembolism that take account of individual level factors such as age, sex, deprivation, smoking, comorbidities, and concurrent drugs. These can then be used to explain the likely risk-benefit ratio to patients considering antipsychotic treatment as well as to identify those at highest risk who are likely to need closer monitoring or interventions to reduce their level of risk.
What is already known on this topic
Antipsychotic drugs, some of which are also prescribed for nausea, vomiting, and vertigo, have been associated with an increased risk of venous thromboembolism
Previous studies have been small, restricted to certain population groups, or have not covered newer atypical antipsychotic drugs
What this study adds
In a large primary care population there was an association between use of antipsychotic drugs and risk of venous thromboembolism
The increased risk was more marked among new users and those prescribed atypical antipsychotic drugs
Though these findings add to the accumulating evidence of adverse health events associated with antipsychotic drugs, they should be confirmed with other data sources
Cite this as: BMJ 2010;341:c4245
We acknowledge the contribution of EMIS and EMIS practices contributing to the QResearch database.
Contributors: CP contributed to the design, reviewed the literature, undertook the primary analysis, contributed to the interpretation, and wrote the first and subsequent drafts of the paper. CC contributed to the development of the protocol, design and analysis, and interpretation and drafting of the paper. JH-C initiated the study, contributed to the design, obtained approvals, prepared the data, and checked and added to the analysis, interpretation, and drafting of the paper. All authors approved the final draft of the paper originally submitted to the BMJ. JH-C and CC approved the final draft of the paper revised after referees’ comments as CP had retired. JH-C is guarantor.
Funding: This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.
Competing interests: All authors have completed the Unified Competing Interest form at www.icmje.org/coi_disclosure.pdf (available on request from the corresponding author) and declare no support from any organisation for the submitted work. JHC is codirector of QResearch, a not-for-profit organisation that is a joint partnership between the University of Nottingham and EMIS. EMIS is the leading supplier of IT for 60% of UK general practices. This work and any views expressed within it are solely those of the co-authors and not of any affiliated bodies or organisations.
Ethical approval: The proposal was approved by the Trent multicentre research ethics committee.
Data sharing: The patient level data from the QResearch are specifically licensed according to its governance framework. See www.qresearch.org for further details. The Read codes groups used are available from the authors on request.
This is an open-access article distributed under the terms of the Creative Commons Attribution Non-commercial License, which permits use, distribution, and reproduction in any medium, provided the original work is properly cited, the use is non commercial and is otherwise in compliance with the license. See: http://creativecommons.org/licenses/by-nc/2.0/ and http://creativecommons.org/licenses/by-nc/2.0/legalcode.