Abstract
PURPOSE Evidence that fewer children are being seen at family physician (FP) practices has not been confirmed using population-level data. This study examines the proportion of children seen at FP and pediatrician practices over time and the influence of patient demographics and rurality on this trend.
METHODS We conducted a retrospective longitudinal analysis of Vermont all-payer claims (2009-2016) for children aged 0 to 21 years. The sample included 184,794 children with 2 or more claims over 8 years. Generalized estimating equations modeled the outcome of child attribution to a FP practice annually, with covariates for calendar year, child age, sex, insurance, and child Rural Urban Commuting Area (RUCA) category.
RESULTS Over time, controlling for other covariates, children were 5% less likely to be attributed to a FP practice (P <.001). Children had greater odds of attribution to a FP practice as they aged (odds ratio (OR) = 1.11, 95% CI, 1.10-1.11), if they were female (OR = 1.05, 95% CI, 1.03-1.07) or had Medicaid (OR = 1.09, 95% CI, 1.07-1.10). Compared with urban children, those from large rural cities (OR = 1.54, 95% CI, 1.51-1.57), small rural towns (OR = 1.45, 95% CI, 1.42-1.48), or isolated/small rural towns (OR = 1.96, 95% CI, 1.93-2.00) had greater odds of FP attribution. When stratified by RUCA, however, children had 3% lower odds of attending a FP practice in urban areas and 8% lower odds in isolated/small rural towns.
CONCLUSIONS The declining proportion of children attending FP practices, confirmed in this population-based analysis and more pronounced in rural areas, represents a continuing challenge.
INTRODUCTION
The overwhelming majority of US children are seen for primary care by family physicians (FPs) and pediatricians (PEDs).1 Studies suggest, however, that the proportion of children seeing FPs has been decreasing.2–5 This trend is significant for several reasons. First, family medicine and primary care pediatrics represent different approaches to care, offering families different areas of expertise. Second, if FPs are seeing fewer children over time, the concept of family medicine as a discipline serving entire families becomes problematic, calling for careful reconsideration of how the discipline is described.1,3 Third, insofar as FPs and PEDs are not evenly spread geographically,6 with marked disparities in access within states, such a trend could be worsening health care disparities. Fourth, fewer children seen in FP practices would mean more children would need to transition from pediatric-only to adult-oriented care. Such transitions are challenging, especially for youth with chronic health problems.7 Finally, to the extent that workforce policy affects the distribution of clinicians, policy makers’ understanding of this shift has implications for future physician training and practice sustainability.
Our understanding of trends in the decreasing proportion of children seeing FPs has not been confirmed on a population basis and is grounded solely in surveys: serial cross-sectional surveys of patient visits,8 and physician-reported data from the American Board of Family Medicine recertification surveys. It is important, therefore, to confirm the shift in proportion of children seeking primary care from FPs to PEDs at the population level. Vermont represents an excellent testing ground for exploring FP/PED trends—especially with respect to the geographic distribution of clincians. The Vermont Health Care Uniform Reporting & Evaluation System9 is an all-payer claims database. It is not subject to the sampling and incomplete response bias present in the National Ambulatory Medical Care Survey methodology or the self-report bias of the American Board of Family Medicine recertification surveys, and it offers patient-level longitudinal data of utilization. Also, the highly rural state of Vermont has an ample per capita supply of both FPs and PEDs,6,10,11 but also has 25% of its child population living in rural areas with a low supply of physicians trained to care for children.6,12 This demographic situation provides the opportunity to address rural residence as a factor in the site of primary care for children.
The objectives of this study were to test, based on Vermont claims data, (1) whether the proportion of children at the population level seen at FP and PED practices is changing over time, (2) how this change is influenced by patient demographics (age, sex, rurality, insurance), and (3) whether these associations vary by the rurality of patient residence.
METHODS
Data and Population
The study team accessed the Vermont Health Care Uniform Reporting & Evaluation System for services provided from January 1, 2009 through December 31, 2016. Our sample contained claims for visits for all children (aged 0 to 21 years) with a reported service, physician, and home zip code in Vermont enrolled in any health plan including commercial health plans, Medicaid, and Medicare. Our study population began with 309,722 children with at least 1 claim during 2009-2016. Of those, 198,562 children were attributed to a FP or PEDS practice. The final sample was limited to 184,794 children having at least 2 years with claims during the 8 year period.
Attribution Methodology
We used a hierarchical approach to attribute children to primary care clinicians (physicians, nurse practitioners, or physician assistants) each year based on services children received, similar to previous work by Christensen and colleagues.13 Visits were classified using Current Procedural Terminology codes and International Classification of Diseases Clinical Modification codes (version 9 before October 2015, and version 10 thereafter; see Supplemental Table at http://www.annfammed.org/content/17/5/390/suppl/ DC1). Children were attributed to primary care clinicians who provided all or the majority of the well-child visits. If a child had no well-child visits, they were attributed to the clinician for all or the majority of specific ambulatory visits including evaluation and management visits and other preventive visits, such as vision and hearing screening. If a child was still unattributed at this point, they were attributed to the PED or FP primary care clinician who provided the majority of the child’s care for any reason (ie, visits could have any of the Current Procedural Terminology or International Classification of Diseases Clinical Modification codes not used in the previous steps). Any remaining children were attributed to the clinician they had seen most often, most recently.
Because the Vermont Health Care Uniform Reporting & Evaluation System does not contain fields linking all clinicians with their practice or organization, we developed procedures for matching the clinician on a child’s medical claim with a practice or organization historically, for each year. We used various state and public resources including Vermont’s health care reform program, the Vermont Blueprint for Health, and the Vermont Area Health Education Centers. The Vermont Department of Health also conducts a census of selected health care professionals that includes practice locations and time spent by clinicians at each location. The Department of Vermont Health Access (Medicaid) maintains a list of active Medicaid providers. In addition to these resources, we contacted the largest private insurance provider in Vermont for in-network clinicians and their practice information, and consulted practice websites, practice and organization press releases and newsletters, and commercial websites with clinician information to confirm when clinicians may have joined or left a practice. Finally, we contacted practices by e-mail and telephone to verify dates of service for some clinicians. When possible, clinicians were associated with the practice where they spent the majority of their time in a calendar year. In cases where a clinician appeared in medical claims in 1 year but not in following years, only in later years, or intermittently over many years, the clinician’s history of practice in Vermont was verified.
Statistical Analyses
Analyses included the years 2009-2016, child age modeled as a continuous variable, and sex as binary female vs male, and children with any Medicaid product were classified as having Medicaid as a binary variable. Rural geographic location was categorized by linking children’s home zip code to the modified Rural Urban Commuting Area (RUCA) classification system. RUCA categories were assigned to zip codes and incorporated commuting data from census tracts into their criteria.14 The analysis used the 4-tier consolidation (urban, large town, small town, and isolated rural) of the original RUCA classification system.15,16
To address study objectives 1 and 2, longitudinal multiple regression analyses including 8 years of data were conducted using a generalized estimating equation statistical model. The generalized estimating equation model was expanded to examine the moderating effect of rurality (RUCA category) on the associations between all covariates and the site of primary care practice type to address objective 3. Interaction terms between RUCA category and covariates were added to the above generalized estimating equation model, first to test for significant interactions, and then the model was stratified into the 4 RUCA categories to explore effects within geographic regions. All statistical analyses were conducted in Stata SE version 15 (StataCorp LLC) statistical software and associations were considered statistically significant at P <.05.
The scope of work for this study was approved by the Institutional Review Board at the University of Vermont and by the Green Mountain Care Board of Vermont.
RESULTS
Table 1 shows the number and percentage of children (aged 0-21 years) with a medical claim attributed to a FP practice by year, child age (categories derived from Bright Futures Guidelines),17 sex, RUCA category for home location, and type of insurance. Percentages are in relation to 100% of children with claims to either a FP or PED primary care practice. The data suggest that the percentage of children attributed to a FP practice decreased from 2009 to 2016. The data in Table 1 show similar trends over time within categories of the child demographic variables as well. The changes over the years were tested in multiple regression models controlling for the other covariates. Over the study period, we found the average distribution of children to be: urban (39%), large town (20%), small town (16%), and isolated rural (25%). This is consistent with the US Census identification of most of Vermont residents living in rural areas.18
All covariates included in the generalized estimating equation multiple regression model were significantly associated with attribution to a FP practice (P <.001). In any given year, children had 5% less odds of attending a FP than PEDS practice compared with the year before, controlling for child demographics (age, sex, RUCA category, and type of insurance), and accounting for the correlation among repeated child observations over 8 years (OR = 0.95, 95% CI, 0.94-0.95). In this adjusted model, the child demographic variables independently predicted attending a FP practice. Females, compared with males, had 5% greater odds of attending a FP practice (OR = 1.05, 95% CI, 1.03-1.07). Children with Medicaid, compared with children with other insurance, had 9% greater odds of attending a FP practice (OR = 1.09, 95% CI, 1.07-1.10). For each year increase in age, children had 11% greater odds of attending a FP practice (OR = 1.11, 95% CI, 1.10-1.11). We wondered if these findings were influenced by including young adults (aged 19 to 21 years) in the sample. Thus, we conducted a sensitivity analysis of children aged 0 to 18 years and found results of similar magnitude, direction, and significance. Children living in nonurban areas, compared with urban areas, had greater odds of attending a FP practice when living in large towns (OR = 1.54, 95% CI, 1.51-1.57), small towns (OR = 1.45, 95% CI, 1.42-1.48), and isolated rural areas (OR = 1.96, 95% CI, 1.93-2.00).
In investigating the moderating effects of rurality on our primary associations of interest, we found significant interactions between RUCA and all covariates in the model (P <.02). Detailed results from the same generalized estimating equation model stratified into 4 RUCA categories are shown in Table 2. The effect of calendar year on attending a FP practice was largest in the isolated rural RUCA areas, where for every 1 year increase, children living in isolated rural areas had an 8% lower odds of attending a FP practice. Regarding RUCA interactions among the demographic variables, the odds for age over all 4 RUCA categories were significant, with a somewhat more pronounced association for residence in urban areas or large towns. For sex, the association between sex and attending a FP practice was only significant for residence in the urban areas, with females from urban areas having a 13% greater odds of attending a FP than a PED practice. The odds of children with Medicaid, as opposed to private insurance coverage, attending a FP over a PED practice was also highly associated with living in urban areas, although significant smaller associations in either direction also were present in large towns and small towns.
DISCUSSION
This longitudinal, population-based analysis of the child population in Vermont for the years 2009-2016 adds substantially to our understanding of the state of FP vs PED practice for child health care. This study confirms the findings of previous studies and yields new information about the influence of the rurality of the child’s residence and of child’s age, sex, and insurance status on the type of practice they are receiving primary preventative medical care from.
This study’s finding of a 5% overall lower odds of attending a FP practice over time, controlling for confounders, verifies previous findings4,19 of a temporal trend of child visits moving from FP to PED practices. Others have delved into the reasons for this overall trend1,3,4,20,21 and its implications, and we will not add to those speculations. The evidence shown in Table 2 that this trend varied according to the rurality of children’s address, with decreasing odds of attending a FP practice associated with increasing rurality of residence, runs counter to the American Board of Family Medicine survey finding20 that rural FPs were more likely than counterparts to see pediatric patients. However, that latter finding was based on the location of the FP, whereas the present findings are based on the location of children’s residence.
There are several possibilities for the shift in practice type of children’s primary care. The number of children initially establishing care with family physicians may be a reflection of the field’s national decline in offering obstetric services.22,23 Analysis of the National Ambulatory Medical Care Survey data, however, suggested the decreased role of FPs providing prenatal care was not accompanied by a similar drop in provision of well-child care to children younger than 2 years.5 It is possible children and families new to Vermont and looking to establish primary care preferentially choose PED practices because they are accepting children as new patients. In 2010, Vermont had the second highest state median age,24 and FP practices may have a higher caseload of older adults and thus not room in their practice for younger children. Projected shortfalls in the primary care workforce to address more complex, chronic ambulatory care needs of older individuals may divert resources from children’s primary care25 and is worthy of further study.
A recent report on Vermont’s health care workforce offers insight into this shift from FPs to PEDs. There has been a steady decrease in FP numbers, and fewer FPs (81%) than PEDs (97%) are accepting new patients.26 A shortage in adult primary care clinicians in rural areas may decrease access for pediatric patients in those FP practices. In our study, we found that for every year of increase in age children were 11% to 13% more likely to attend an FP practice in more urban areas, whereas in more rural areas they were only 9% more likely. This particular finding may be most generalizable to states where the numbers of FPs and PEDs are adequate.
The association of increasing child age with likelihood of attending a FP practice has been described previously in the National Ambulatory Medical Care Survey analyses and is readily understood, as research on adolescents suggests that they are more comfortable in settings appropriate for their age.27 The finding that females are more likely to receive care in a FP practice, and that this finding is largely driven by an association with patients’ urban residence, is novel. Unfortunately, our data set did not allow us to query for the sex of the practitioners within FP or PED practices to investigate whether this was driven in some way by the sex of the clinicians within practices. This seems unlikely, however, as in the years of our study, the number of female primary care clincians in Vermont grew by 69%, compared with 9% for males,26 and other data suggest that this has occurred equally among FPs and PEDs. The National Ambulatory Medical Care Survey data could be analyzed to assess the association between patient sex and type of primary care physician to assess whether this finding is applicable more widely.
The association of FP practice over PED practice among patients covered by Medicaid, largely attributable to families living in more urban areas, is challenging to explain, and runs counter to the higher concentration of PED practices in urban areas. This finding contradicts a separate data source of self-reported capacity in Vermont to accept new patients when considering insurance status, with FPs accepting smaller proportion (73%) with Medicaid compared with PEDs (97%).26 One possible explanation, although beyond the scope of this analysis, is that FP practices are also more prevalent in more urban areas, and that access to care for parents who themselves are covered by Medicaid is easier, translating to more pediatric visits to those practices. Beyond the potential convenience of parent and child being seen the same day, it is possible that FP practices that maintain acceptable access to parents with Medicaid are also those that maintain their pediatric clientele over time.
Limitations
In spite of the unique opportunities and advantages offered by Vermont as a model for investigating our research question, it is also possible that what happens in a rural New England state with a small population is not generalizable to the rest of the nation. Similar studies using the all-payer claims database, however, could run similar analyses based on the method of attribution described herein. Research using the all-payer claims database is known for potential problems, including mislabeling of activity.28 Our algorithm for practice attribution relied on codes selected by clinicians, who may have coded erroneously when multiple clinical issues arose in one visit. Primary care sought outside Vermont, or for which there is no medical claim submitted, is not part of this analysis. Finally, a drop in the number of claims submitted to the Vermont Health Care Uniform Reporting and Evaluation System occurred in the final year of analysis, which coincided with a decision by some payers not to submit data, although no evidence suggests that payers preferentially reported on patients of distinct age, sex, or home location. Across the claim years examined, the average number of children we attributed to clincians was within 2% of the 2010 US Census.29
CONCLUSIONS
In this study, we used medical claims data to determine whether the proportion of children at the population level seen at FP and PED practices was changing over time, how patient demographics influenced any change, and whether identified associations varied by the rurality of patient residence. Our analysis confirmed the temporal trend from FP to PED practices noted in survey studies, and suggested a new finding related to an enhanced trend toward PED practices among children living in rural areas. This will need to be confirmed by research in other settings.
Footnotes
Conflicts of interest: authors report none.
To read or post commentaries in response to this article, see it online at http://www.AnnFamMed.org/content/17/5/390.
Previous presentations: Academic Pediatric Association Region I Meeting; March 26, 2018; Waltham, Massachusetts; and the Pediatric Academic Societies Meeting; May 5-8, 2018; Toronto, Ontario, Canada.
Disclaimer: The analyses, conclusions, and recommendations from these data are solely those of the study authors, and are not necessarily those of the Green Mountain Care Board (GMCB). The GMCB had no input into the study design, implementation, or interpretation of the findings. The conclusions are those of the authors alone, and no official endorsement by the GMCB was intended or should be inferred.
Supplemental materials: Available at http://www.AnnFamMed.org/content/17/5/390/suppl/DC1/.
- Received for publication October 22, 2018.
- Revision received March 1, 2019.
- Accepted for publication March 26, 2019.
- © 2019 Annals of Family Medicine, Inc.