1,720
Views
28
CrossRef citations to date
0
Altmetric
Research Paper

Validation of infant immunization billing codes in administrative data

, , , , , , , , & show all
Pages 1840-1847 | Received 13 Jan 2015, Accepted 16 Apr 2015, Published online: 06 Jul 2015

Abstract

Ontario has a single payer provincial health insurance program. Administrative data may provide a potentially robust source of information for post-marketing vaccine studies. Vaccine-specific immunization billing codes were introduced in 2011. Our objective was to validate Ontario's universal health care administrative datasets to assess infant immunization status. Electronic medical record data from the Electronic Medical Record Administrative data Linked Database (EMRALD) was used as the reference standard to calculate performance characteristics of the Ontario Health Insurance Plan (OHIP) database vaccine-specific and general immunization codes for 4 primary infant immunizations: diphtheria, tetanus, acellular pertussis, inactivated polio, Haemophilus influenzae type B (DTaP-IPV-Hib) combination vaccine, pneumococcal conjugate vaccine, measles, mumps, rubella (MMR) vaccine, and meningococcal conjugate serogroup C vaccine. OHIP billing claims had specificity ranging from 81% to 92%, sensitivity 70% to 83%, positive predictive value (PPV) 97% to 99%, and negative predictive value (NPV) 13% to 46% for identifying the various specific vaccines in administrative data. For cohorts vaccinated in the new code introduction phase, using both the vaccine-specific and general codes had higher sensitivity than the vaccine-specific codes alone. In conclusion, immunization billing claims from administrative data in Ontario had high specificity and PPV, moderate sensitivity, and low NPV. This study identifies some of the applications of utilizing administrative data for post-marketing vaccine studies. However, limitations of these data decrease their utility for measuring vaccine coverage and effectiveness. Therefore, the establishment of a comprehensive and linkable immunization registry should be a provincial priority.

Introduction

Vaccines are one of the greatest public health achievements of the last century.Citation1 Public confidence in vaccines is important for ensuring continued success of this achievement. Therefore, enhanced post-marketing vaccine effectiveness, coverage, and safety surveillance is a critical component of any immunization program.Citation2-4

In Ontario, Canada's largest province (population 13.5 million with 142,448 births in 2013),Citation5 infant and toddler immunizations are almost exclusively administered through physician offices and funded under the Ontario Health Insurance Plan (OHIP).Citation6 Measles, mumps, rubella (MMR) and meningococcal conjugate serogroup C (MenC) vaccines are recommended at 12 months of age, while diphtheria, tetanus, acellular pertussis, inactivated polio, Haemophilus influenzae type B (DTaP-IPV-Hib) combination vaccine and pneumococcal conjugate vaccine (PC) are given at 2, 4, 6, and 18 months, and at 2, 4, and 12 months, respectively. In addition, rotavirus vaccine is recommended at 2 and 4 months of age, varicella at 15 months of age, and booster doses of MMR-varicella vaccine and DTaP-IPV are currently recommended between 4 to 6 y of age.

Physicians submit billing claims to OHIP for providing immunizations, which are captured in the OHIP database. Prior to 2011, physicians submitted claims using general immunization codes for any vaccine administered. In September 2011, 9 vaccine-specific codes were introduced. Due to the lack of a comprehensive centralized immunization registry that includes pre-school aged children, administrative data sets are a potential data source to assess vaccine coverage as they contain individual-level data spanning the whole province. However, the validity of using both the newer vaccine-specific codes and the earlier general immunization codes recorded in the administrative data is unknown.

Across Canada, the use of electronic medical record (EMR) systems by family physicians has been increasing.Citation7 The objective of this study was to utilize primary care EMR data recorded in family physician offices to validate OHIP physician billing claims for general and vaccine-specific immunization codes.

Results

We initially identified 2204 children for the MMR/MenC 2009 birth cohort, 1823 for the DTaP-IPV-Hib/PC 2010 birth cohort, 3329 children for the MMR/MenC 2011 birth cohort, and 2661 for the DTaP-IPV-Hib/PC 2012 birth cohort. Application of the exclusion criteria left final cohorts of 1604, 1557, 2404, and 2194, respectively (). These cohorts were drawn from the practices of 147, 167, 246, and 257 family physicians, respectively.

Figure 1. Flow diagram describing the application of the exclusion criteria to the 4 cohorts used in the analysis. MMR = measles, mumps, rubella vaccine; MenC = Meningococcal C vaccine; DTaP-IPV-Hib = Diptheria, tetanus, acellular pertussis, inactivated polio, Haemophilus influenzae type b vaccine; PC = pneumococcal conjugate vaccine; BC = birth cohort; ID = Identification number *No billed primary care assessments or fewer than 4 visits to their family physician during the first year of life.

Figure 1. Flow diagram describing the application of the exclusion criteria to the 4 cohorts used in the analysis. MMR = measles, mumps, rubella vaccine; MenC = Meningococcal C vaccine; DTaP-IPV-Hib = Diptheria, tetanus, acellular pertussis, inactivated polio, Haemophilus influenzae type b vaccine; PC = pneumococcal conjugate vaccine; BC = birth cohort; ID = Identification number *No billed primary care assessments or fewer than 4 visits to their family physician during the first year of life.

Characteristics of the study children, mothers, and physicians were summarized and compared to the entire Ontario 2011 birth cohort in . Compared to the Ontario population, rural patients are overrepresented in the Electronic Medical Record Administrative data Linked Database (EMRALD); 10% Ontario infants live in a rural setting compared to 21% in EMRALD. The immigrant population is underrepresented in EMRALD and physicians using this EMR system are more likely to have graduated in the last 1–2 decades and less likely to be foreign trained, compared to all Ontario physicians ().

Table 1. Child, maternal, and physician characteristics of the combined four cohorts used for validation analysis compared to the 2011 Ontario population

The performance measures of OHIP immunization codes were relatively similar across all 4 vaccines (). For the general immunization codes, sensitivity was 79.4%–83.2%, specificity 81.0%–92.2%, PPV 96.8%–99.5%, and NPV 18.9%–46.4% for the various vaccines. For the vaccine-specific codes, the sensitivity was 70.1%–72.4%, specificity 88.5%-91.5%, PPV 96.9–99.4%, and NPV 12.9%–38.8%. Comparing the performance of the vaccine-specific codes (post-2011) with the general vaccine codes (pre-2011), overall sensitivity declined. Specificity decreased for DTaP-IPV-Hib and PC but increased for MMR and MenC. PPV was unchanged, and NPV decreased.

Table 2. Evaluating Ontario Health Insurance Plan physician billing claims performance characteristics with 95% confidence intervals compared with electronic medical records as the reference standard

In the sensitivity analysis that included both the general and vaccine-specific immunization codes post-2011, we observed a marginal improvement in sensitivity for all vaccines, with modest reductions in specificity compared to using the vaccine-specific codes alone ().

We identified predictors of discordance between immunization billing claims and documentation in the EMR (). Children with lower continuity of care scores tended to have higher rates of discordant results (p < 0.05). Physicians who saw proportionally more young children in their practice had higher concordance (p < 0.05). Maternal influenza immunization status also correlated with concordance (p < 0.05); however, the magnitude was small. Predictors of discordance for each vaccine and cohort can be seen in Appendix .

Table 3. Predictors of discordant results for the 2-month dose of diphtheria, tetanus, acellular pertussis, inactivated polio, and Haemophilus influenzae type b vaccine using Ontario Health Insurance Plan vaccine-specific billing claims for the 2012 birth cohort

We evaluated performance measures stratifying by COC and pediatric roster size (). As an example, the sensitivity of MMR in the 2009 birth cohort with high COC was 85.4% compared to 69.1% with low COC. Specificity was 92.4% and 43.8% in the high and low COC groups, respectively.

Table 4. Performance measures for select birth cohorts stratified by continuity of care and physician pediatric roster size

Discussion

By utilizing EMR patient charts as the reference standard for immunization status, we identified that OHIP billing claims have high specificity and PPV, moderate sensitivity, and low NPV in establishing immunization status of children in Ontario. Newer, vaccine-specific codes have worse sensitivity than previous general immunization codes.

The most likely explanation for the moderate sensitivity of OHIP claims is that physicians do not bill because remuneration is small. Ontario's schedule of benefits reimburses CDN $4.50 for each immunization (the fee is independent of the vaccine given), however depending on the practice remuneration model, a physician may receive as little as $0.68 per vaccine billing claim. Furthermore, primary care reform has been a major initiative in Ontario that has led to the creation of multidisciplinary team practices.Citation8 In this model, more vaccines are administered by nursing staff, so physicians may not always submit billing claims for these services. This factor may partially explain the drop in sensitivity from 2009 and 2010 to 2011 and 2012 that persisted in the sensitivity analysis, using both general and vaccine-specific codes. This sensitivity analysis also revealed minimal change from the primary vaccine-specific codes analysis, confirming that most EMRALD physicians have adopted these new codes. However, multiple vaccines are administered at the 2-month visit. Documentation of the vaccine of interest (e.g., DTaP-IPV-Hib) corresponding to a general OHIP billing code resulted in classification as a true positive. But since the general code can be used for any vaccine, there is a small chance of misclassification and this may have inflated the sensitivity and PPV.

Immunizations of young children are mostly administered through physician offices in Ontario, and it is unlikely that children who are seeing an EMRALD-participating physician on a regular basis (at least 4 times during the first year of life) would receive their vaccines in a walk-in-clinic or other office besides their primary care physician. This assumption is exemplified by the very high specificity; however a small percentage of children may have been vaccinated elsewhere or by salaried physicians.

Important limitations of using an EMR as the reference standard for immunizations should be noted. EMR records are dependent on the clinician documenting administration of the vaccine. There are multiple areas of free text within the EMR, and while we attempted to capture all relevant text entries, missed recordings are possible. The data available through EMRALD are a voluntary sample of Ontarian physicians who all use one type of EMR system and practice under some type of primary care reform model of care, and therefore may not be entirely representative of all physicians in the province. Patients captured by EMRALD are more likely to live in rural areas and to be cared for by younger physicians who have adopted EMR systems.

This study provides important information on the validity of using province-wide administrative data for assessing immunization status, particularly given the absence of a comprehensive all-of-life vaccine registry in Ontario. These findings will help to support a number of future applications. First, immunization status in the OHIP database can be linked to health services utilization to study potential adverse effects of specific vaccines better.Citation2 The excellent PPV makes this an optimal data source for such studies, because they often use self-controlled designs that require inclusion only of individuals who have had both the exposure (i.e., immunization) and the outcome of interest. The high specificity and PPV values for these immunization codes indicate that they lend themselves well to these analyses.

The second application is to assess immunization coverage across different jurisdictions and populations, in particular coverage in early infancy and at the crucial milestone of 2 y of age. Ontario's immunization coverage estimates, derived from its centralized repository, are limited to immunization records received following school entry (generally at age 4 to 6 years). However, the moderate sensitivity seen in this study suggests that administrative data will underestimate vaccine coverage in the population.

Finally, by linking immunization data to microbiologic results, we can evaluate vaccine effectiveness on a population level. This analysis is becoming increasingly important with outbreaks of vaccine-preventable diseases such as measlesCitation9 and pertussis.Citation10 This application may however underestimate vaccine effectiveness because of the non-differential misclassification bias in immunization status that we have demonstrated using administrative data.Citation11

EMR data contained more immunizations than administrative data, although some records observed in the administrative data were not recorded in the EMR. Therefore a provincial immunization registry drawn from combining EMR data with administrative data may provide a more complete picture of immunization coverage. Several approaches could be undertaken to improve the quality of EMRALD data. ICES currently provides semi-annual feedback to physicians using EMRALD in the form of quality indicators related to several chronic diseases, including diabetes and ischemic heart disease.Citation12 Expanding this feedback to include up-to-date immunization status for each vaccine preventable disease may provide additional incentive to physicians using EMRs to improve data quality and/or coverage. Providing continuing education to clinicians on methods to improve the completeness and accuracy of this data source may enhance the utility of these data for research purposes. We strongly encourage further expansion of efforts to link EMR data to administrative data, in order to increase our knowledge about coverage of the population until a comprehensive provincial immunization registry is established.

In summary, OHIP's new vaccine-specific (and older general immunization) billing codes have high specificity and PPV, but only moderate sensitivity and low NPV for predicting immunization status of children in Ontario. OHIP billing claims are a potentially rich source of information for post-marketing vaccine safety surveillance. We plan to use the results of this analysis to evaluate waning immunity from pertussis immunization as well as measles vaccine effectiveness. These results will also be used in ongoing studies of vaccine safety and to evaluate immunization coverage in infants. However, establishment of a comprehensive vaccine registry that can be linked to administrative data should be a priority, given the limitations of this dataset.

Materials and Methods

We conducted a validation study of vaccine billing codes submitted by physicians compared to the reference standard of documentation in a primary care EMR.

Study population and setting

We included pre-specified birth cohorts of Ontario children under the care of family physicians who share their practices' EMR data with EMRALD, a centralized repository of EMR data used for research and evaluation.Citation13 EMRALD has been previously used for administrative data validation studies to identify patients with a variety of disease conditions.Citation14-16

We created 4 cohorts of children. For the MMR and MenC vaccines, we enrolled children born between January 1 2009 to December 31 2009 and January 1 2011 to December 31 2011 to evaluate the general and vaccine-specific immunization codes, respectively. Similarly, for the first DTaP-IPV-Hib and PC vaccine doses, we enrolled children born between January 1 2010 to December 31 2010 and January 1 2012 to December 31 2012 to evaluate the general and vaccine-specific immunization codes, respectively. We selected these study periods to validate the general codes before the introduction of the vaccine-specific codes in September 2011, and the vaccine-specific codes after that date, allowing for a period of several months for physicians to adjust to using the newer codes.

We excluded children who had duplicate identification numbers (preventing 1:1 linkage between EMRALD and the OHIP database) and those who died during the observation period. We also excluded children with no billed primary care assessments, and those with fewer than 4 visits to their family physician during the first year of life. These children were likely to have resided in Ontario intermittently or received their primary care from a salaried physician who did not submit billing claims to OHIP.

Data source for immunization reference standard

EMRALD currently comprises EMR data from over 350 Ontario family physicians who use Practice Solutions Suite® software, the most widely used EMR in Ontario.Citation17 Individual-level data from EMRALD are collected semi-annually and linked to other administrative databases at the Institute for Clinical Evaluative Sciences (ICES). All clinically relevant information is collected in EMRALD, including clinical encounters, the cumulative patient profile, family history, allergies, immunizations, diagnostic tests, prescriptions, discharge summaries, and consultation notes. Participating physicians contribute to EMRALD on a voluntary basis. This data set was used as the reference standard for immunization status.

We conducted text searches for each vaccine of interest in the prescriptions and immunizations fields in EMRALD using a series of keywords (e.g., “MMR,” “measles, mumps, and rubella vaccine;” full list in Supplemental ). We also searched EMRALD for records of immunization billing claims submitted through the EMR software.

We used different observation windows depending on the vaccine of interest. Since both MMR and MenC vaccines are scheduled at 12 months of age, the observation window spanned 335 d to 455 d after birth (approximately 11–15 months of age) to capture children vaccinated before or after the recommended age. We limited the validation of the DTaP-IPV-Hib and PC vaccines to the 2-month dose (i.e., the first dose in the series). The observation window was 53 d (2 months less 1 week) to 112 d (4 months less 1 week) after birth. We used the same observation windows to evaluate the general immunization codes.

Data source for immunization codes

The OHIP database contains physician billing claims paid for by OHIP, which covers virtually all of Ontario's approximately 13.5 million residents, except for recent immigrants and migrants (i.e., residing in Ontario for <3 months) as well as a very small percentage of patients who see non-billing physicians.Citation18 We used these data to identify infant immunization codes. We searched for physician billing claims for MMR vaccine (G845), MenC vaccine (G844), DTaP-IPV-Hib vaccine (G841), PC vaccine (G846), and general immunizations with (G538) and without (G539) physician consultation.

Data sources for covariates

The data were linked to administrative datasets housed at ICES to identify factors related to concordance between EMR and OHIP records.

The Registered Persons Database (RPDB) was used to identify the child's sex, vital status, and socio-demographic data. The child's postal code was linked to Canadian census data to attribute mean household income quintile of their neighborhood and rural residence (community size under 10 000).

The MOMBABY database is an administrative data set maintained at ICES comprised of admission records of delivering mothers and their newborn babies, which are linked through a unique matching number on each hospitalization record. We used this dataset and the Canadian Institute of Health Information's Discharge Abstract Database (CIHI-DAD) to identify the mothers of the children in our cohort and to determine whether each mother had a previous delivery (primiparous vs. not), to measure the maternal age at first parity, and to identify infants of low birth weight. The CIHI-DAD was also used to identify chronic medical conditions among children, as described by Feudtner et al.,Citation19 within the first year of life.

The OHIP database was used to identify the number of primary care visits in the first year of life and to calculate a continuity of care (COC) score. COC is defined as the number of visits to an individual's primary care physician divided by the total number of physician visits during the first year of life. We defined a COC of less than 50% as low.Citation20,21 We also used the OHIP database to identify maternal influenza immunization status during the year following delivery using the influenza-specific billing codes G590 (influenza immunization with a physician consultation) and G591 (influenza immunization without a physician consultation).

The Citizen and Immigration Canada (CIC) database contains information on individuals who have landed in Ontario since 1985, and we used this data set to ascertain maternal immigration status. Mothers were considered recent immigrants if they arrived to Canada within the previous 5 y

The ICES Physician Database contains information on physician demographics and specialization. We obtained the following physician covariates: sex, rural practice, decade of graduation, and place of medical training (foreign vs. domestic).

The Client Agency Program Enrolment (CAPE) and OHIP databases were used to determine the volume of patients aged 6 y or younger on a physician's roster.

Analysis

EMRALD immunization status was set as the reference standard for immunization, and we calculated performance measures of sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) for OHIP immunization codes for each vaccine, during each time period (i.e., before and after introduction of the vaccine-specific codes). In a sensitivity analysis, we analyzed the 2011 birth cohort for the MMR and MenC vaccines and the 2012 birth cohort for DTaP-IPV-Hib and PC vaccines to include both general and vaccine-specific immunization codes. We also examined discordant results between EMRALD and OHIP to ascertain the nature of the discordance. We stratified selected variables to demonstrate performance measures within each stratum.

Disclosure of Potential Conflicts of Interest

No potential conflicts of interest were disclosed.

Disclaimer

The opinions, results and conclusions reported in this paper are those of the authors and are independent from the funding sources. No endorsement by ICES, PHO, or MOHLTC is intended or should be inferred.

Supplemental Material

Supplemental data for this article can be accessed on the publisher's website

Supplemental material

Supplementary_Tables.zip

Download Zip (30.5 KB)

Funding

This work was supported by the Ontario Ministry of Health and Long-Term Care (MOHLTC). This study was also supported by the Institute for Clinical Evaluative Sciences (ICES) and Public Health Ontario (PHO), which are funded by annual grants from the MOHLTC. KS is supported by a Fellowship Award from the Canadian Institutes of Health Research (CIHR). KT is supported by a Clinician Scientist Award from the University of Toronto Department of Family and Community Medicine (DFCM). JCK is supported by a New Investigator Award from CIHR and a Clinician Scientist Award from the University of Toronto DFCM. The study sponsors did not participate in: the design and conduct of the study; collection, management, analysis and interpretation of the data; preparation, review or approval of the manuscript; or the decision to submit the manuscript for publication.

References

  • The Canadian Public Health Association: 12 great achievements. Accessed January 12, 2015. Available from: http://www.cpha.ca/en/programs/history/achievements.aspx.
  • Wilson K, Hawken S, Kwong JC, Deeks SL, Crowcroft NS, Manuel D. Vaccine and immunization surveillance in ontario (VISION) - using linked health administrative databases to monitor vaccine safety. Vaccine 2012; 30:6115-20; PMID:22709951; http://dx.doi.org/10.1016/j.vaccine.2012.06.004
  • Chen RT, Glasser JW, Rhodes PH, Davis RL, Barlow WE, Thompson RS, Mullooly JP, Black SB, Shinefield HR, Vadheim CM, et al. Vaccine Safety Datalink project: a new tool for improving vaccine safety monitoring in the United States. The Vaccine Safety Datalink team. Pediatrics 1997; 99:765-73; PMID:9164767; http://dx.doi.org/10.1542/peds.99.6.765
  • Lankinen KS, Pastila S, Kilpi T, Nohynek H, Makela PH, Olin P. Vaccinovigilance in Europe–need for timeliness, standardization and resources. Bull W H O 2004; 82:828-35; PMID:15640918
  • Statistics Canada: Accessed November 18, 2014. Available from: http://www.statcan.gc.ca/tables-tableaux/sum-som/l01/cst01/demo02a-eng.htm.
  • Publically funded immunization schedule in Ontario (2011): Accessed May 30 2014. Available from: http://www.health.gov.on.ca/en/public/programs/immunization/docs/schedule.pdf.
  • College of Family Physicians of Canada CMA, Royal College of Physicians and Surgeons of Canada. National physician survey 2014. Results 2014. [accessed 2014 December 11] Available from: http://nationalphysiciansurvey.ca/surveys/2014-survey/survey-results-2/
  • Hutchison B, Glazier R. Ontario's primary care reforms have transformed the local care landscape, but a plan is needed for ongoing improvement. Health Aff (Millwood) 2013; 32:695-703; PMID:23569049; http://dx.doi.org/10.1377/hlthaff.2012.1087
  • Bird C. Measles outbreaks threaten those averse to vaccine. Canad Med Associat J 2013; 185:E393-E4; PMID:23649422; http://dx.doi.org/10.1503/cmaj.109-4480
  • Waters V, Jamieson F, Richardson SE, Finkelstein M, Wormsbecker A, Halperin SA. Outbreak of atypical pertussis detected by polymerase chain reaction in immunized preschool-aged children. Pediatr Infect Dis J 2009; 28:582-7; PMID:19561423; http://dx.doi.org/10.1097/INF.0b013e318197fac1
  • Ozasa K. The effect of misclassification on evaluating the effectiveness of influenza vaccines. Vaccine 2008; 26:6462-5; PMID:18573297; http://dx.doi.org/10.1016/j.vaccine.2008.06.039
  • Ivers NM, Tu K, Young J, Francis JJ, Barnsley J, Shah BR, Upshur RE, Moineddin R, Grimshaw JM, Zwarenstein M. Feedback GAP: pragmatic, cluster-randomized trial of goal setting and action plans to increase the effectiveness of audit and feedback interventions in primary care. Implement Sci 2013; 8:142; PMID:24341511; http://dx.doi.org/10.1186/1748-5908-8-142
  • Tu K, Mitiku TF, Ivers NM, Guo H, Lu H, Jaakkimainen L, Kavanagh DG, Lee DS, Tu JV. Evaluation of electronic medical record administrative data linked database (EMRALD). Am J Manag Care 2014; 20:e15–21; PMID:24669409.
  • Tu K, Mitiku T, Lee DS, Guo H, Tu JV. Validation of physician billing and hospitalization data to identify patients with ischemic heart disease using data from the Electronic Medical Record Administrative data Linked Database (EMRALD). Can J Cardiol 2010; 26:e225-8; PMID:20847968; http://dx.doi.org/10.1016/S0828-282X(10)70412-8
  • Tu K, Wang M, Young J, Green D, Ivers NM, Butt D, Jaakkimainen L, Kapral MK. Validity of administrative data for identifying patients who have had a stroke or transient ischemic attack using EMRALD as a reference standard. Can J Cardiol 2013; 29:1388-94; PMID:24075778; http://dx.doi.org/10.1016/j.cjca.2013.07.676
  • Tu K, Wang M, Jaakkimainen RL, Butt D, Ivers NM, Young J, Green D, Jette N. Assessing the validity of using administrative data to identify patients with epilepsy. Epilepsia 2014; 55:335-43; PMID:24417710; http://dx.doi.org/10.1111/epi.12506
  • Ontario MD. EMR Enhanced Use Program. [accessed 2014 December 11]; Available from: https://http://www.ontariomd.ca/portal/server.pt/community/emr_funding/emr_enhanced_use_program/.
  • Chan B. Supply of physicians' services in ontario. Hosp Q 1999; 3:17; PMID:11570359
  • Feudtner C, Christakis DA, Connell FA. Pediatric deaths attributable to complex chronic conditions: a population-based study of Washington State, 1980–1997. Pediatrics 2000; 106:205-9; PMID:10888693
  • Campitelli MA, Inoue M, Calzavara AJ, Kwong JC, Guttmann A. Low rates of influenza immunization in young children under Ontario's universal influenza immunization program. Pediatrics 2012; 129:e1421–30; PMID:22585770; http://dx.doi.org/10.1542/peds.2011-2441
  • Guttmann A, Manuel D, Dick PT, To T, Lam K, Stukel TA. Volume matters: physician practice characteristics and immunization coverage among young children insured through a universal health plan. Pediatrics 2006; 117:595-602; PMID:16510636; http://dx.doi.org/10.1542/peds.2004-2784