Article Text

Download PDFPDF

Determinants of community pharmacists’ quality of care: a population-based cohort study using pharmacy administrative claims data
  1. Nancy Winslade1,2,
  2. Robyn Tamblyn3
  1. 1 Department of Medicine, McGill University, Montreal, Quebec, Canada
  2. 2 Winslade Consultants, Ottawa, Ontario, Canada
  3. 3 Department of Epidemiology and Biostatistics, McGill University, Montreal, Quebec, Canada
  1. Correspondence to Dr Nancy Winslade; nancy.winslade{at}


Objective To determine if a prototype pharmacists’ services evaluation programme that uses linked community pharmacy claims and health administrative data to measure pharmacists’ performance can be used to identify characteristics of pharmacies providing higher quality of care.

Design Population-based cohort study using community pharmacy claims from 1 November 2009 to 30 June 2010.

Setting All community pharmacies in Quebec, Canada.

Participants 1742 pharmacies dispensing 8 655 348 antihypertensive prescriptions to 760 700 patients.

Primary outcome measure Patient adherence to antihypertensive medications.

Predictors Pharmacy level: dispensing workload, volume of pharmacist-provided professional services (eg, refusals to dispense, pharmacotherapy recommendations), pharmacy location, banner/chain, pharmacist overlap and within-pharmacy continuity of care. Patient level: sex, age, income, patient prescription cost, new/chronic therapy, single/multiple antihypertensive medications, single/multiple prescribers and single/multiple dispensing pharmacies. Dispensing level: prescription duration, time of day dispensed and antihypertensive class. Multivariate alternating logistic regression estimated predictors of the primary outcome, accounting for patient and pharmacy clustering.

Results 9.2% of dispensings of antihypertensive medications were provided to non-adherent patients. Male sex, decreasing age, new treatment, multiple prescribers and multiple dispensing pharmacies were risk factors for increased non-adherence. Pharmacies that provided more professional services were less likely to dispense to non-adherent hypertensive patients (OR: 0.60; 95% CI: 0.57 to 0.62) as were those with better scores on the Within-Pharmacy Continuity of Care Index. Neither increased pharmacists’ services for improving antihypertensive adherence per se nor increased pharmacist overlap impacted the odds of non-adherence. However, pharmacist overlap was strongly correlated with dispensing workload. There was significant unexplained variability among pharmacies belonging to different banners and chains.

Conclusions Pharmacy administrative claims data can be used to calculate pharmacy-level characteristics associated with improved quality of care. This study supports the importance of pharmacist’s professional services and continuity of pharmacist’s care.

  • quality in health care
  • hypertension
  • health policy

This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See:

View Full Text

Statistics from

Strengths and limitations of this study

  • The trial directly measured community pharmacy characteristics using pharmacy claims and health administrative data.

  • The primary quality of care outcome used a standardised method for measuring patient adherence to medications.

  • The trial was population based and included a large sample of patients from community pharmacies in Quebec.

  • Performance on only one quality of medication-use indicator was evaluated and results may not apply to additional measures of pharmacists’ quality of care.

  • Administrative data are limited in the extent to which they can measure services provided by pharmacists that were not billed.



Misuse of prescription medications, ranging from inappropriate prescribing to patient non-adherence, remains a significant and costly challenge to health systems.1 The medication-related expertise and accessibility of community pharmacists have led policy-makers to re-evaluate the role community pharmacists play in managing medication misuse.2 3 Emphasis has been placed on the care provided by pharmacists both as part of medication dispensing and via expanded professional services that target specific medication-misuse problems.4 5 Although such care can improve patient’s medication use, community pharmacists struggle to incorporate expanded professional services into their daily practice.6 7 As a result, payers continue to seek evidence of the real-world impact of community pharmacists’ services on medication misuse,4 8 9 and quality indicators of unsafe or interacting medications and management of non-adherent patients have been established as standardised outcome measures of pharmacists’ quality of care.10–12 The services that pharmacists provide to achieve high performance on these quality indicators can vary across jurisdictions.5 Developments in the use of community pharmacy administrative claims data have enabled the measurement of both pharmacy-level performance on these standardised quality indicators and the impact of pharmacists’ professional services on patient outcomes.13 14

To date there have been no precise methods of determining pharmacy-level characteristics that consistently support high levels of pharmacists’ performance and that could inform directions for pharmacy policy. Pharmacy characteristics such as workload, continuity of care, culture, workflow and overlap of pharmacists have been evaluated through self-report and with varying definitions of quality performance.15–17 The few studies that used standardised quality measures employed a potentially biased ecological approach to estimate pharmacy characteristics by determining a population-based quality metric in the geographical area and then assigning these population-based results to all pharmacies within that area.18–21 More robust methodologies are needed to measure the characteristics of the patient, pharmacy and workload situation when the patient receives the medication.22

One potentially powerful option is to use pharmacy administrative claims data to measure salient pharmacy characteristics. Until now use of such data has been limited to identifying whether the pharmacy is a chain or independent, and the volume of dispensing.20 23 This is primarily due to challenges in using the large volume of pharmacy administrative data to create accurate measures, as well as challenges linking pharmacy claims data to other health administrative databases to obtain information on patient and pharmacy characteristics. Increasingly these linkages have been enabled through interest by payers in monitoring performance and researchers in conducting population-based studies.24 25 We developed a framework for pharmacists’ services evaluation that uses linked pharmacy administrative claims and health administrative data to measure and feed back pharmacy-level performance on quality indicators, followed by diagnostic on-site assessments of lower performing pharmacies.26 The objective of this study was to determine if the linked administrative health data used within this prototype pharmacists’ services evaluation programme could be used to identify characteristics of pharmacies providing higher quality of care.



This study was conducted in Quebec, with a population of 8 million patients of whom approximately 3.5 million receive government support for payment of their medications via the Régie de l’Assurance Maladie du Québec (RAMQ). Since the late 1970s, Quebec pharmacists have been authorised to bill RAMQ for professional services such as refusals to dispense medications and written pharmaceutical opinions for the management of specific medication-use problems.27 28 RAMQ requires the date, hour, drug identification number, therapeutic drug class, dosage form, strength, quantity, duration of treatment, specific type and reason for the pharmacist service (eg, previous adverse effect or management of underuse of antihypertensive medications) and costs to RAMQ, the patient and for the overall prescription. All data are coded and can be linked to other health administrative data using unique encrypted identifiers for patients, prescribers, pharmacists and pharmacies. For patients, age, sex, postal code and average household income are recorded. For pharmacies, the location (eg, shopping centre) and type of pharmacy (independent or not) are maintained, along with the specific chain or banner to which the pharmacy belongs.

Study design

A population-based prospective cohort of patients was assembled for whom Quebec pharmacists billed for dispensings of antihypertensive medications between 1 November 2009 to 30 June 2010. A dispensing was defined as the preparation and provision of medications to a patient pursuant to a prescription, regardless of quantity of medication dispensed. Each time there was a dispensing for an antihypertensive medication we determined whether the dispensing was to a patient who was adherent or not over the 90 days prior to the dispensing. Characteristics of each dispensing, the patient and the pharmacy were measured and a multilevel model used to identify predictors of dispensing to a non-adherent patient.


All 1891 pharmacies in Quebec were included unless they had opted out of participating in a previously reported randomised controlled trial, were open <61 days, or had dispensed >165 317 prescriptions over the 8-month study period, which represented outliers with Z-scores >2.5.14 29 Pharmacies with shorter open days did not have sufficient data for reliable calculation of characteristics, and very high dispensing volumes were not representative of traditional community pharmacy practice in Quebec. We had sufficient sample size to have 90% power to detect a difference in antihypertensive adherence of 5% for most potential predictors.

Primary outcome

The primary outcome was whether a dispensing of an antihypertensive medication was provided to an adherent or non-adherent patient. Antihypertensive adherence was selected for this initial evaluation as antihypertensive medications are widely used and non-adherence is common.30 Our previous research had also documented that almost all community pharmacies in Quebec (99.7%) dispense antihypertensive medications, thereby allowing a population-based cohort for the current study.12

For each antihypertensive dispensing, we created a record of all dispensings of the same antihypertensive medication to the same patient from all pharmacies in Quebec over the previous 180 days. ‘Same medication’ was defined as the same drug in the same dosage format, regardless of strength. Switches to a new medication in the same therapeutic class were treated as new therapies. Dispensings of antihypertensive medications were excluded if the patient had not been treated with the same medication for at least 90 days or had not had continuous insurance coverage over the previous 180 days. As dispensing pharmacists are responsible for obtaining information on medications supplied from other pharmacies when determining adherence, each eligible dispensing was attributed to the dispensing pharmacy. We calculated the proportion of previous 90 days covered (PDC) for the same medication using the previous dispensing dates and number of days of supply provided at each dispensing and adjusting for early refills. If the PDC over the 90 days prior to the dispensing was less than 72 days (80%), then the dispensing was to a non-adherent patient.31

Potential predictors

Dispensing-level characteristics included the type of antihypertensive medication dispensed, the total prescription cost and the cost to the patient as these have been demonstrated to affect patient adherence.32 Although in Quebec the standard supply of medications is for 30 days, patients at risk for non-adherence can receive weekly medication supply and patients stabilised on chronic therapies can receive 90-day supplies. Adherence was, therefore, expected to be worse for patients receiving less than 30 days’ supply and better for patients receiving more than 30 days’ supply.

Patient-level characteristics were those known to affect adherence such as sex, age and income, with older men and patients with higher income anticipated to be more compliant.20 33 As our previous work indicated that patients within their first 6 months of antihypertensive therapy are less compliant as are those on single drug therapy or receiving their antihypertensive medications from more than one physician or pharmacy, these variables were also included.12 34–36

Pharmacy-level characteristics included workload as higher numbers of prescriptions dispensed has been identified as a factor limiting community pharmacists’ ability to provide professional services 37 and predisposing to dispensing errors.15 17 38 Workload has been reported variously as prescriptions dispensed per year, which can readily be determined from administrative claims data, to prescriptions per pharmacist per hour, which has only been reported using self-reported estimates.17 We received from RAMQ the total number of billings and open days for each pharmacy over the 8-month study period and used the administrative claims data to calculate for each pharmacy the average number of: open hours per day, pharmacists billing per hour, prescriptions dispensed per hour and prescriptions dispensed per pharmacist per hour. Related to workload, as medication dispensing errors occur more frequently when only one pharmacist is working, there have been calls for mandatory overlapping of pharmacists’ schedules to allow one pharmacist to focus uninterruptedly on prescription verification while a second pharmacist provides professional services.16 39 To measure pharmacist overlap for each pharmacy, we created a matrix of the number of pharmacists billing each open hour over each open day during the 8-month study period. From this, we calculated the average per cent of each pharmacy’s open hours where more than one pharmacist was billing (Pharmacist Overlap Index). Finally, although continuity of care measuring whether patients received all antihypertensive medications from a single pharmacy was included as a patient-level variable, based on evidence from other health professions that care from the same healthcare professional is important in creating trust, professional relationships, we determined the likelihood that a patient would be cared for by the same pharmacist on multiple visits (Within-Pharmacy Continuity of Care Index).34 We calculated, for each pharmacy, the total number of pharmacists working over the 8-month study period (weighted to emphasise differences in high and low numbers of pharmacists) and divided this by the average number of pharmacists working per day at that pharmacy. The lowest value of the index is 1, representing the best within-pharmacy continuity of care when there is only one single pharmacist working in the pharmacy over the 8 months. Increasing indices indicate a lower chance that the patient would be cared for by the same pharmacist at multiple visits. To determine the culture within the pharmacy, we calculated the total number of pharmacists’ professional services billed per 100 prescriptions dispensed over the 8-month period, including refusals to dispense, pharmaceutical opinions, transmission of medication profiles and emergency contraception. We also counted the number of professional services billed specifically for management of underuse of antihypertensive medications.

Data sources/measurement

Baseline community pharmacy claims data for all dispensings of antihypertensive medications and pharmacist services were received from RAMQ for all Quebec community pharmacies for the period of 1 October 2008 to 30 June 2010.14 Patient, pharmacy, pharmacy chain/banner group, pharmacist and prescriber identifiers were anonymised by RAMQ prior to data transfer. Data for the 8-month period of 1 November 2009 to 30 June 2010 were used to calculate dispensing, patient and pharmacy-level characteristics and estimate determinants of non-adherence.

Statistical methods

Descriptive statistics summarised the characteristics of the dispensings, patients and pharmacies including the incidence of dispensing to non-adherent patients by type of antihypertensive, patient sex and age. Multivariate alternating logistic regression (ALR) estimated the association among the dispensing, patient and pharmacy-level characteristics and non-adherence. ALR allows analysis of dichotomous outcomes when observations have more than one level of clustering.40 For our results, ALR first measured the extent of clustering of non-adherence among multiple dispensings within the same patient and then for multiple patients receiving their medications from the same pharmacy. All analyses were completed using SAS V.9.4, with ALR using PROC GENMOD.

Where multiple measures could be calculated to reflect a single construct, results for each measure were first compared with previously reported estimates (if available) to test the accuracy of the calculations. Next, each measure was tested individually for association with non-adherence. A single measure of each construct was selected for inclusion based on the accuracy of the calculation, the strength of evidence supporting its use and the strength of association. Collinearity was evaluated for all variables considered for the final analysis using the variance inflation factor. When collinearity was present, variables that were calculated as interim steps were considered for exclusion and the variables retained were those most directly measuring the constructs of interest. To account for interactions between patient income and the cost of the medication to the patient, we divided both variables into low, medium and high categories and created dummy variables for each of the nine possible interactions, setting low income and low cost to the patient as the reference.41


Study participants

A total of 1872 pharmacies were enrolled in the study, after 19 (1%) opted out of the previous trial (figure 1). Ninety-one pharmacies open for <61 days and 39 additional pharmacies dispensing >165 317 prescriptions over the 8-month period were removed from the analysis. Evaluation was carried out for 8 655 348 dispensings of antihypertensive medications to 760 700 patients in 1742 pharmacies

Figure 1

Consort diagram. *Dispensings that were provided to patients who had been either adherent or non-adherent with their antihypertensive medication over the previous 90 days. **Patients with at least one adherent dispensing over the 8-month study period. ***Patients with at least one non-adherent dispensing over the 8-month study period. As patients received multiple dispensings, they could be counted as both adherent and non-adherent; therefore, the total of adherent and non-adherent patients is more than 760 700.

Population characteristics

Angiotensin receptor blockers (ARBs) were the most commonly dispensed antihypertensive medications (23.2% of dispensings) with <1% of dispensings for each of alpha agonists, alpha blockers, potassium sparing diuretics and vasodilators (table 1). Most prescriptions were dispensed in the morning and were for an approximate 1-month duration. A single physician prescribed antihypertensive medications to 74.1% of patients and 86.0% went to a single pharmacy for all of their antihypertensive medications over the previous 6 months. Most patients had been taking antihypertensive medications for more than 6 months (98.5%) and were on multiple antihypertensive medications (79.4%). The majority of pharmacies were either chains or banners (89.9%). Pharmacists dispensed an average 18.4 prescriptions per pharmacist per hour, billing for 0.18 professional services for every 100 prescriptions dispensed. Most pharmacies did not have any billings for pharmacists’ services for antihypertensive non-adherence, leading to an average of less than one billing over the 8 months (0.35±1.8). Pharmacies had more than one pharmacist billing for 15.5% of their open hours and an average of nine different pharmacists worked in each pharmacy over the 8-month study period.

Table 1

Characteristics of prescriptions dispensed, patients and their pharmacies


Over 8 months, 9.2% of all dispensings of antihypertensive medications were provided to non-adherent patients (795 031 of 8 655 348 dispensings) (table 2). Antihypertensive dispensings were provided to 760 700 distinct patients, 31% of whom were non-adherent to their antihypertensive medication at least once over the study period (235 885 of 760 700). The highest incidence of non-adherence occurred with alpha agonists (21.49%) and for dispensings provided in the evening (12.03%). The incidence of non-adherence was also higher if the patient was <65 years old (12.41%), new to therapy (18.29%) or on a single antihypertensive medication (12.47%).

Table 2

Dispensed prescription, patient and pharmacy characteristics associated with risk of non-adherence with antihypertensive medications

When adjusted for the three levels of variables and clustering, the odds of non-adherence were significantly greater for medications supplied for less than 10 days and for medications dispensed at times other than morning (p<0.05)(table 2). Relative to beta-blockers, the odds of dispensing an ARB or ACE inhibitor to a non-adherent patient were decreased by 17% (OR: 0.83; 95% CI: 0.82 to 0.84).

Older, female patients were less likely to be non-adherent at the time of receiving an antihypertensive medication, with a 41% decrease in the odds for patients ≥80 years relative to patients <65 years old (OR: 0.59; 95% CI: 0.58 to 0.60). Patients newly started on their antihypertensive medication within the past 6 months experienced a 27% increase in odds of non-adherence at the time of dispensing. Patients with decreased continuity of care were also more likely to be non-adherent at the time of dispensing, with the odds of non-adherence increased by 10% if the patient had used multiple pharmacies and 16% if she/he had used multiple physicians for their antihypertensive medications over the past 6 months. The impact of cost of the medication to the patient was modified by the patient’s income and, in contrast to the unadjusted incidence of non-adherence where increasing out-of-pocket costs lead to higher non-adherence, when adjusted for all three levels of characteristics, higher out-of-pocket costs resulted in a decreased odds of non-adherence within all of low-income, middle-income and high-income patients. High-income patients with low out-of-pocket medication costs were 15% more likely to be non-adherent at the time of dispensing as compared with low-income patients with low medication costs (OR: 1.15, 95% CI: 1.12 to 1.18).

At the pharmacy level, the odds of non-adherence decreased by 40% per 1 increase in the number of professional services billed per 100 prescriptions dispensed (OR: 0.60; 95% CI: 0.57 to 0.62). Neither the number of billings for pharmacists’ services targeted at managing non-adherence with antihypertensive medications nor the percentage of open hours with overlapping pharmacists influenced non-adherence. However, pharmacist overlap was highly correlated with dispensing volume (Pearson correlation coefficient 0.51, p<0.0001). Higher workload decreased the odds of non-adherence by 4% per 10 prescription increase in number of prescriptions dispensed per pharmacist per hour (OR: 0.96; 95% CI: 0.96 to 0.97). Higher scores on the Within-Pharmacy Continuity Care Index, indicating a decreased chance of patients being cared for by the same pharmacist, slightly but significantly increased the odds of non-adherence (OR: 1.003; 95% CI: 1.001 to 1.005). There was significant variability in the odds of non-adherence among pharmacies belonging to various banners or chains and the odds of non-adherence were significantly higher for chains/banners relative to independent pharmacies (OR: 1.02; 95% CI: 1.00 to 1.05).


Statement of principal findings

This study is the first to document that linked community pharmacy claims and health administrative data can be used to directly measure a range of pharmacy-level characteristics and quality measures. It is also the first study that investigated the association between the provision of pharmacists’ professional services and better within-pharmacy continuity of care with adherence, showing that each of these pharmacists’ practices are associated with a decreased odds of dispensing antihypertensive medications to non-adherent patients.

Strengths and limitations

The main strengths of this study are the direct measurement of pharmacy characteristics from administrative claims data and the use of an objective, validated quality of care measure of adherence.10–12 As significant variability in results has been reported from studies using differing measures of adherence, use of standardised methods for measuring adherence is particularly important in determining predictors of non-adherence.10 As only 1% of community pharmacies in Quebec did not consent to participate (18 of 1891), a second strength is that the sample approximated a population-based cohort and selection bias was minimised. Limitations include that we evaluated performance on only one quality of medication-use measure and results cannot be generalised to other measures of pharmacists’ quality of care. Although underuse measures of other therapeutic categories such as lipid-lowering or diabetes may show similar results, determinants of performance on quality indicators measuring medication overuse (eg, rescue inhalers for asthma) or unsafe dispensing may differ as the professional services pharmacists provide to detect and manage these medication-use problems differ from those provided for medication underuse. Evaluation of performance on additional quality indicators measuring both adherence and unsafe dispensing is required to determine if results are generalizable. In addition, our methodology for calculating adherence did not allow for detection of primary non-adherence or non-adherence/non-persistence within the first 90 days of therapy. As these types of non-adherence are problematic with antihypertensive medications, our results may have underestimated non-adherence and measures of these additional types of non-adherence should be evaluated. Finally, administrative claims data are limited in the extent to which they can measure whether pharmacists provided a service but did not bill for it.42–44


Our overall rate of non-adherence is consistent with previous reports that use community pharmacy administrative claims data and similar measures of non-adherence.10 45 Calculation of pharmacy-level characteristics required multiple steps and complex analysis and for characteristics that had previously been estimated via self-report, such as prescriptions per pharmacist per hour, our results were higher (18.4±10.5 our study vs 14.1±4.9).17 This is consistent with national reports documenting higher total prescriptions dispensed in Quebec relative to other provinces.37 Results of the drug and patient characteristics affecting non-adherence agree with previous research documenting that there is higher adherence to antihypertensive medications with fewer side effects, such as ARB and ACE, and that increasing age is associated with increased adherence to antihypertensive medications.32 46 However, given the variability in results of non-adherence rates and predictors from studies that used differing measurement methodologies, our results should be compared with studies using pharmacy administrative claims data and standardised methods for measuring non-adherence.10 To our knowledge, this literature is limited to the study that used an ecological approach to measuring pharmacy and patient characteristics.20 Our results differ from this ecological study for the impact of patient sex and income, and independent pharmacy ownership on the odds of dispensing to a non-adherent patient. Our results demonstrate the impact of measuring these characteristics directly for each dispensing and adjusting for clustering. When only considering whether the pharmacy is independent versus a chain/banner, the incidence of non-adherence is higher in independent pharmacies. However, when adjusted for clustering and the remaining dispensing, patient and pharmacy characteristics, this association reverses with chain/banner pharmacies demonstrating a greater odds of non-adherence. The same is true for the impact of patient costs relative to income. Without adjustment, the incidence of non-adherence increases as cost to the patient increases. However, when adjusted for all characteristics, this relationship reverses. As higher patient cost typically occurs with second-line treatments for hypertension, this may represent patients who required switches or additions to their therapies due to side or insufficient effects from their initial treatments, which has been shown to increase adherence.47

The most striking results of our analysis are the reductions in the odds of non-adherence with both an increasing rate of provision of pharmacists’ professional services and improved within-pharmacy continuity of care. It is hypothesised that the relationship between the rate of provision of these services and lower non-adherence indicates that improved quality of care is provided at pharmacies where pharmacists prioritise provision professional services versus involvement in technical distributive functions.48 49 The relationship between improved within-pharmacy continuity of care and decreased odds of non-adherence supports such a hypothesis as patients can more easily develop trusting relationships with their pharmacist when continuity of care is improved. Our findings that increased workload is associated with lower odds of non-adherence would not appear to support that increased workload challenges pharmacists’ provision of quality care. However, we had removed very high-volume pharmacies so we did not see the previously reported results of lower quality of care in pharmacies with both very low and very high dispensing volumes.15 The strong positive correlation between workload and pharmacist overlap suggests that pharmacists are not being scheduled to provide professional services but to enable increased number of prescriptions to be processed. As both culture and workflow are determined predominantly by the pharmacist owner, greater freedom to emphasise professional pharmacists’ practice by owners of independent pharmacies could account for their lower odds of non-adherence relative to chains/banners.50 Similarly, differences in practical philosophy among the chains/banners could account for the variability in performance among the different banners and chains.

Implications and future research

Our results indicate that emphasis on the caring role of pharmacists both during dispensing and via provision of professional services appears key to improving patients’ use of medications. Results also support policies that encourage continuity of care and that focus adherence strategies on younger men, new to treatment and taking single antihypertensive therapy. Pharmacy administrative claims data can be used to directly measure dispensing, patient and pharmacy characteristics, thereby increasing the range and accuracy of pharmacy-level characteristics evaluated. Evaluation of additional measures both of non-adherence and dispensing of contraindicated medications is needed to determine if there is consistency across the measures of pharmacy-level characteristics identified in our study as being related to pharmacists’ quality of care.


We thank Danielle Fagnan, from the Quebec Order of Pharmacists, for the support provided throughout the project. We thank Sherry Shi for her analysis of the data. The critical review of the manuscript and recommendations for analysis by Drs Cees van der Vleuten and Lambert Schuwirth are gratefully acknowledged. We acknowledge the financial support of Green Shield Canada Foundation and the financial and administrative support of the Quebec Order of Pharmacists.


  1. 1.
  2. 2.
  3. 3.
  4. 4.
  5. 5.
  6. 6.
  7. 7.
  8. 8.
  9. 9.
  10. 10.
  11. 11.
  12. 12.
  13. 13.
  14. 14.
  15. 15.
  16. 16.
  17. 17.
  18. 18.
  19. 19.
  20. 20.
  21. 21.
  22. 22.
  23. 23.
  24. 24.
  25. 25.
  26. 26.
  27. 27.
  28. 28.
  29. 29.
  30. 30.
  31. 31.
  32. 32.
  33. 33.
  34. 34.
  35. 35.
  36. 36.
  37. 37.
  38. 38.
  39. 39.
  40. 40.
  41. 41.
  42. 42.
  43. 43.
  44. 44.
  45. 45.
  46. 46.
  47. 47.
  48. 48.
  49. 49.
  50. 50.
View Abstract


  • Contributors NW and RT conceived and designed the project. NW reviewed literature to determine relevant potential predictors. RT provided expertise on the statistical analysis. NW completed the basic analysis and drafted the manuscript. RT contributed substantially to the manuscript revision. NW redrafted and finalised the manuscript following reviewers’ recommendations. RT contributed substantially to the final revisions.

  • Funding Green Shield Foundation Canada provided a grant for the completion of this project and was both informed of the decision to submit this paper and received a copy of the manuscript. The Quebec Order of Pharmacists provided a matching grant in the form of salary support, facilitated and paid the required data access fees to RAMQ.

  • Competing interests NW reports grants from Green Shield Foundation Canada, personal fees and non-financial support from Quebec Order of Pharmacists during the conduct of the study. RT reports grants from Green Shield Foundation Canada and non-financial support from Quebec Order of Pharmacists during the conduct of the study.

  • Ethics approval McGill University’s Faculty of Medicine Institutional Review Board.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data sharing statement No additional unpublished data are available from this study.

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.