Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Influence of peer networks on physician adoption of new drugs

  • Julie M. Donohue ,

    Roles Conceptualization, Funding acquisition, Investigation, Methodology, Supervision, Writing – original draft, Writing – review & editing

    jdonohue@pitt.edu

    Affiliation Department of Health Policy and Management, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America

  • Hasan Guclu,

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliations Department of Statistics, School of Engineering and Natural Sciences, Istanbul Medeniyet University, Istanbul, Turkey, Department of Biostatistics and Medical Informatics, School of Medicine, Istanbul Medeniyet University, Istanbul, Turkey

  • Walid F. Gellad,

    Roles Conceptualization, Funding acquisition, Investigation, Writing – review & editing

    Affiliations Center for Pharmaceutical, Policy and Prescribing, Health Policy Institute, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America, Department of Medicine, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America, Center for Health Equity Research and Promotion, Veterans Affairs Pittsburgh Healthcare System, Pittsburgh, Pennsylvania, United States of America

  • Chung-Chou H. Chang,

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Validation, Writing – review & editing

    Affiliations Department of Medicine, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America, Department of Biostatistics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America

  • Haiden A. Huskamp,

    Roles Conceptualization, Funding acquisition, Investigation, Writing – review & editing

    Affiliation Department of Health Care Policy, Harvard Medical School, Harvard University, Boston, Massachusetts, United States of America

  • Niteesh K. Choudhry,

    Roles Conceptualization, Funding acquisition, Investigation, Writing – review & editing

    Affiliation Department of Medicine, Division of Pharmacoepidemiology and Pharmacoeconomics, Brigham and Women’s Hospital and Harvard Medical School, Boston, Massachusetts, United States of America

  • Ruoxin Zhang,

    Roles Data curation, Formal analysis, Methodology, Validation, Writing – review & editing

    Affiliation Quality and Operations Support, The Permanente Medical Group, Inc., Oakland, California, United States of America

  • Wei-Hsuan Lo-Ciganic,

    Roles Data curation, Formal analysis, Investigation, Methodology, Writing – review & editing

    Affiliation Department of Pharmacy, Practice and Science, College of Pharmacy, University of Arizona, Tucson, Arizona, United States of America

  • Stefanie P. Junker,

    Roles Project administration, Resources, Writing – review & editing

    Affiliation Department of Health Policy and Management, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America

  • Timothy Anderson,

    Roles Formal analysis, Investigation, Methodology, Writing – review & editing

    Affiliation Department of Medicine, University of California San Francisco, San Francisco, California, United States of America

  • Seth Richards-Shubik

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Validation, Writing – original draft, Writing – review & editing

    Affiliation Department of Economics, College of Business and Economics, Lehigh University, Bethlehem, Pennsylvania, United States of America

Abstract

Although physicians learn about new medical technologies from their peers, the magnitude and source of peer influence is unknown. We estimate the effect of peer adoption of three first-in-class medications (dabigatran, sitigliptin, and aliskiren) on physicians’ own adoption of those medications. We included 11,958 physicians in Pennsylvania prescribing anticoagulant, antidiabetic, and antihypertensive medications. We constructed 4 types of peer networks based on shared Medicare and Medicaid patients, medical group affiliation, hospital affiliation, and medical school/residency training. Instrumental variables analysis was used to estimate the causal effect of peer adoption (fraction of peers in each network adopting the new drug) on physician adoption (prescribing at least the median number prescriptions within 15 months of the new drug’s introduction). We illustrate how physician network position can inform targeting of interventions to physicians by computing a social multiplier. Dabigatran was adopted by 25.2%, sitagliptin by 24.5% and aliskiren by 8.3% of physicians. A 10-percentage point increase in peer adoption in the patient-sharing network led to a 5.90% (SE = 1.50%, p<0.001) increase in physician adoption of dabigatran, 8.32% (SE = 1.51%, p<0.001) increase in sitagliptin, and 7.84% increase in aliskiren adoption (SE = 2.93%, p<0.001). Peer effects through shared hospital affiliation were positive but not significant, and medical group and training network effects were not reliably estimated. Physicians in the top decile of patient-sharing network peers were estimated to have nearly 2-fold stronger influence on their peers’ adoption compared to physicians in the top decile of prescribing volume. Limitations include lack of detailed clinical information and pharmaceutical promotion, variables which may influence physician adoption but which are unlikely to bias our peer effect estimates. Peer adoption, especially by those with whom physicians share patients, strongly influenced physician adoption of new drugs. Our study shows the potential for using information on physician peer networks to improve technology diffusion.

Introduction

Diffusion of technology in US healthcare, while influenced partly by payer policies regarding coverage and reimbursement, is to a large extent driven by the individual decisions of practicing physicians. The relative absence of centralized technology assessment in the US reflects the importance of physician autonomy, allows for flexible responses to rapidly changing evidence, and creates opportunities to tailor decisions to individual patients. However, it comes at the cost of wide geographic variation,[1, 2] excess spending, and sluggish translation of evidence into practice.[3, 4] Recognizing that no US regions provide uniformly better care, an Institute of Medicine report recommended that efforts to achieve high-value healthcare target the loci of decision-making–hospitals, physician groups, and individual providers[5]. Yet the scale of changing provider behavior at the individual-level is daunting.

To make the most efficient use of resources for educating the workforce,[68] providers may be viewed as embedded in social systems.[911] While it is known that physicians learn from each other, the magnitude of peer influence is poorly understood and is a largely untapped resource. The tools of social network analysis can be used to map the connections among providers, and identify those playing a central role among their peers.[12] The extant studies applying social network methods to technology diffusion among physicians,[11, 1322] have been limited to uptake of a single technology, and small physician samples.[16, 21, 23, 24] Prior studies have relied primarily on physicians’ self-reported information on peer connections—information that, while informative, would be cost-prohibitive to collect on a large scale.

We take advantage of increasingly large and detailed healthcare databases to examine the value of harnessing social network information to drive physician adoption of evidence-based technologies. First, we constructed peer networks among nearly 12,000 physicians drawing on multiple sources of information to form peer networks based on shared patients, practice settings, and training. Second, we estimated the magnitude of peer influence on technology adoption using as natural experiments the introductions of three new prescription drugs varying in novelty, clinical indication, number of competitors, and the specialties of physicians prescribing them. Third, we illustrate how simple information on physicians’ positions in their peer networks can be used to target interventions to change physician behavior.

Methods

Study setting and data sources

Our study setting was Pennsylvania, the 5th largest US state, the population of which mirrors national averages in socio-demographic characteristics and on measures of health care utilization.[25] The study period during which the three drugs of interest were introduced and over which we measured their adoption was 2007–2011. We obtained 5 data sources all of which contained and were linked by National Provider Identifier (NPI). First, physician-level prescribing data were obtained from QuintilesIMS’s Xponent database which directly captures >70% of all US prescriptions and uses a projection method to represent all retail prescriptions filled. Second, from the American Medical Association (AMA) Masterfile we obtained information on physician characteristics including demographics (age, sex), and training (e.g., medical school, graduation year). Third, we obtained information on physician’s organizational affiliations from QuintilesIMS’s Healthcare Organizational Services (HCOS) database which captures provider affiliations with >29,000 US practices, clinics, and hospitals. Fourth, we obtained administrative claims data on all (fee-for-service and managed care) enrollees in PA Medicaid from the PA Department of Human Services. Fifth, Medicare data for all fee-for-service enrollees with Part D pharmacy benefits were obtained from the Centers for Medicare and Medicaid Services (CMS). Medicaid and Medicare cover approximately 1/3 of Pennsylvanians.[26]

New medications of interest

We measured physician adoption of dabigatran (an oral anticoagulant initially approved to treat atrial fibrillation on 10/19/2010), sitagliptin (an oral dipeptidyl peptidase-4 inhibitor approved to treat diabetes on 10/16/2006), and aliskiren (an oral direct renin inhibitor approved to treat hypertension on 3/05/2007) (S1 Supporting Information). All three medications were first-in-class, with a novel mechanism of action, although they varied in the extent to which they were viewed as superior within the broader therapeutic class,[2730] the relevant patient populations, and the availability of substitutes.

Physician cohorts

We constructed three cohorts of physicians applying 4 broad inclusion criteria. We required physicians to: a) prescribe medications in one or more drug classes of interest (oral anticoagulant, antidiabetic or antihypertensive medications) during the study period (Parts a-c of S1 Table); b) have an AMA Masterfile record and a Pennsylvania practice address; c) have a record in HCOS database; and d) demonstrate some minimal prescribing in the therapeutic class in the first 15 months after the new drug was introduced (minimal defined as ≥1 prescription/quarter) (S1S3 Figs for cohort construction).

Peer network construction

Physicians form relationships with peers whom they meet during training, and in office- and hospital-based practice settings, and form referral and information networks with peers both within and outside of their own health systems.[11, 14, 15, 3135] To capture these rich, and potentially overlapping peer networks we formed 4 types of physician social networks illustrated in Fig 1. We used the network analysis library Igraph in python[36] to construct the networks, and measure network characteristics. Network construction is briefly summarized here with additional information provided in S1 Supporting Information.

thumbnail
Fig 1. Illustration of physician peer networks.

The figure shows the affiliation networks of a typical physician in our study cohorts (bottom left). Within the cohort of physicians prescribing the drug class of interest, she is connected to peers (shown in pink) with whom she attended the same medical school (one year +/-) or with whom she completed the same residency program (one year +/-). She is connected to peers (shown in yellow) through the medical group where she has an outpatient practice and to peers (shown in blue) through the hospital where she admits patients. In addition, she shares Medicare and Medicaid patients with several physicians. The patient-sharing network is represented by the lines in the figure; line thickness corresponds to the number of patients shared between physicians. Connections shown in orange are affiliated with the physician in this illustration through shared training institution and medical group. Connections shown in green share a medical group and hospital affiliation in common with the physician. Connections shown in purple only have shared patients with the physician.

https://doi.org/10.1371/journal.pone.0204826.g001

First, we constructed patient-sharing networks (P) (depicted by black lines in Fig 1) using a previously published approach developed using patient-sharing in Medicare data.[3739] Barnett and colleagues used self-reported connections to validate claims-based connections, reporting that two physicians billing Medicare for at least 9–10 patients in common during the same year were highly likely to self-report having a relationship based on referrals and/or advice.[37] We extend prior studies using data from a single payer to construct physician networks[4043] by combining Medicare and Medicaid claims for unique patients to form the patient-sharing network. Second, we constructed a medical group network (G) using data from HCOS. We identified as medical group peers all physicians in the prescribing cohort with whom a physician shared a medical group or clinic affiliation. Third, a hospital network (H) was constructed similarly using data from HCOS on shared hospital affiliation (e.g., attending or admitting); we include in the hospital network all physicians in the cohort with a shared hospital affiliation. Last, a training network (T) was constructed using AMA Masterfile data on institutions attended and dates of graduation. Two physicians were connected if they attended the same medical school or the same residency program within +/-1 year of each other.

Measuring adoption

We defined adoption as writing at least the median number of prescriptions among physicians prescribing the new drug at least once in its first 15 months on the market. The medians were 7, 13, and 7 prescriptions for dabigatran, sitagliptin, and aliskiren, respectively. We tested the robustness of our findings to alternative specifications of the adoption measure (e.g., > = 1, > = 15) (S12 Table).

Peer adoption in each network was the key independent variable in our analysis and was measured as the fraction of a physician’s peers who adopted the new drug (i.e., who had at least the median number of prescriptions in the first 15 months the new drug was on the market). This peer adoption measure (the “peer adoption rate”) was computed in each of the four networks separately as we hypothesized each to have a distinct influence. Hence, there are four separate adoption variables among peers in the: (1) patient-sharing network, (2) medical group network, (3) hospital network, and (4) training network. Peers in the training, medical group and hospital networks were assumed to have equal influence, while peers in the patient-sharing network were weighted based on the number of patients shared (i.e., a weighted average was used to compute the adoption rate in the patient-sharing network) (S1 Supporting Information).

Statistical analyses

We fit both linear probability models and logistic regression models to estimate the influence of these peer adoption rates and other factors, on the individual adoption outcome for each physician. Comparing the estimated outcomes obtained from linear and logistic models, we found that the mean and median differences of the estimated values were ≤0.1% and the first and third quartiles of the difference were within 3%. Accordingly, and to allow us to easily implement instrumental variables methods, we used linear probability models specified as follows: where yi is the binary indicator of adoption. Covariates included: characteristics of the individual physician (xi) such as demographics, training, specialty, prescribing volume, and location; the age distribution and payer mix of the physician’s patients filling prescriptions (also included in xi); and the four variables for peer adoption rates in each network ( to )–the coefficients on these variables being key estimates of interest (S1 Supporting Information for additional details). We checked the overlap in connections among these 4 networks and determined that it was low (<20%) allowing us to estimate the effect of all 4 simultaneously in the same model.

The estimation of peer effects with observational data has many challenges.[4446] The most basic problem is that each physician influences her peers just as they influence her (simultaneity). In addition, physicians may choose peers who are similar to them (homophily), and there may be other, unobserved factors that are common among groups of peers. To address these sources of bias we use an instrumental variables approach to estimation that is common in the econometric literature on peer effects.[47] A set of exogenous characteristics of peers, such as their sex, age, and location, serve as “instruments” which predict peer adoption rates in the absence of confounding factors. Estimation proceeds via two-stage least squares. In the first stage, the means of exogenous peer characteristics were used to predict the adoption rates in each peer network (S1 Supporting Information). In the second stage, the predicted adoption rates were then used in place of the observed adoption rates to estimate our main model.

This approach yields consistent estimates as long as peer characteristics meet two criteria. First, they must have no direct influence on an individual’s adoption outcome and must not be correlated with unobserved confounders. For example, the proportion of peers who are female should not directly affect whether a physician adopts the drug, net of the adoption rate among those peers. By contrast, we believe the proportion of peers who are relevant specialists or high-volume prescribers may have a direct influence on physician adoption. For example, a primary care physician may be less likely to adopt a new drug if he can easily refer complex patients to specialists in the same hospital. Hence, neither the specialty mix nor the prescribing volume of peers are used as instruments. Instead, we include eight variables for the proportions of peers in each network who are relevant specialists ( to ) or are high-volume prescribers ( to ) in the main model. The assumption that the other peer characteristics are valid instruments can be partially assessed,[48] and this assessment is presented in detail in the S1 Supporting Information. Overall the results indicate that our instruments are not correlated with unobserved confounders: the null hypothesis of no confounding is not rejected for two out of the three drugs (p = 0.213 for sitagliptin, p = 0.500 for sitagliptin). The diagnostic testing for dabigatran suggests that unmeasured confounding could be an issue (p = 0.031), but given that a rejection of the null can also occur as a consequence of functional misspecification (e.g., variable definitions, presence or absence of interactions, etc.), and because there is no known reason why this drug would be different from the others, we maintain dabigatran in the present analysis.

Second, the means of peer characteristics must be sufficiently predictive of the peer adoption rates. We assessed this with two sets of tests for weak instruments commonly used in the literature.[49] For all three drugs the peer characteristics were sufficiently predictive in the patient-sharing and hospital networks, with first-stage F-statistics above 10, but not in the medical group and training networks. We then tested whether the instruments were sufficiently predictive for the patient-sharing and hospital networks jointly, using the minimum eigenvalue statistic, which confirmed that the instruments have sufficient power in these two networks (S8 Table). As a consequence of this assessment, we consider the estimated peer effects in the patient-sharing and hospital networks to be statistically reliable, but not those in the medical group or training networks.

After quantifying the magnitude of peer influence in each network, we used our estimates to compute an aggregate “social multiplier” on adoption[50] (S1 Supporting Information). Conceptually, the social multiplier represents the number of other physicians expected to adopt a new drug following adoption by a given individual physician. This reflects both the direct influence a physician has on her own peers as well as her indirect influences on others throughout the network. The formula for the social multiplier in the patient-sharing network, our focal network, is as follows: where wji is the weight on the link between physicians j and i based on the number of shared patients, and γP is the coefficient on the peer adoption rate from the linear probability model of adoption. The social multiplier is an eigenvector centrality equivalent to a type of Bonacich power centrality[51] in a weighted, directed graph, where γP serves as the attenuation parameter. We calculated the social multiplier for each physician in the network, and then used this to conduct a simulation using the patient-sharing network and adoption of dabigatran to illustrate. Imagine a health system wishing to achieve rapid diffusion of a new technology seen as cost-effective for treating a condition. In order to target scarce resources efficiently that system could target interventions to physicians seeing a high volume of the relevant patient population. Alternatively, the payer could target interventions to physicians with many peer connections observed in administrative claims data. To illustrate the value of targeting physician interventions based on social network vs. other physician characteristics, we present the average multiplier (i.e., number of other physicians who would adopt a drug) by physician prescribing volume, and a simple measure of network centrality: the number of direct peer connections a physician has (through patient-sharing), also known as degree.[52]

We used SAS 9.4 (Cary, NC) for data management and variable construction, the network analysis library Igraph in python[36] to construct the peer networks, STATA V14 to estimate the instrumental variable models, and R x64 3.3.2 to compute the multiplier described above. This study was approved by the University of Pittsburgh Institutional Review Board.

Results

Descriptive characteristics of physicians

We examined adoption of new drugs among 7,785 physicians prescribing anticoagulants, 8,257 prescribing antidiabetic medications, and 9,974 prescribing antihypertensives. Physicians were, on average, 49–50 years old and roughly one-quarter were female (Table 1). Most (61–72%) physicians prescribing these medications were in primary care specialties (e.g., internal medicine). The medical sub-specialties of interest–cardiology, nephrology, and endocrinology–made up 3–14% of physicians depending on the cohort. Most physicians were affiliated with at least one medical group (69–75%) and hospital (90%). The proportion of physicians adopting the new drugs varied; one-quarter (25.2%) of anticoagulant prescribers adopted dabigatran, a similar fraction of antidiabetic prescribers adopted sitagliptin (24.5%), but fewer antihypertensive prescribers adopted aliskiren (8.3%) (Table 1). Several characteristics differed significantly in bivariate analyses between adopters and non-adopters including the percent who were specialists and had high prescribing volumes (S3S5 Tables).

thumbnail
Table 1. Characteristics of prescriber cohorts for each class of new drug.

https://doi.org/10.1371/journal.pone.0204826.t001

In each cohort, the patient-sharing network had the greatest number of peers, ranging from 200–344 depending on the cohort, with a small share of providers (2–6%) having no peers in the patient-sharing network (S7 Table). Physicians in each cohort shared at least one hospital affiliation with an average of 166–204 peers. Physicians had fewer peers in the medical group (9–12) and training networks (39–46).

Unadjusted association between peer adoption and own adoption

In unadjusted analyses, a physician’s likelihood of adopting a drug was strongly associated with the likelihood that her peers had also adopted the drug. For example, a physician for whom 2/3 of her peers had adopted dabigatran was twice as likely to adopt the drug compared to a physician for whom one-third of her peers had adopted the same drug (Fig 2a).

thumbnail
Fig 2. a-c: Unadjusted relationship between physicians’ own adoption rate and that of their peers in the patient-sharing network.

The graphs show the average physicians’ own adoption rate for a given fraction of peers in the patient-sharing network adopting the drug. The size of the bubble corresponds to the number of physicians with a particular peer adoption rate. The graphs show a positive association between the fraction of peers adopting the new drug and the physicians’ own likelihood of adopting. For example, with dabigatran, physicians in networks in which the fraction of peers adopting is 0.65 have a 0.38 probability of adopting compared to physicians whose peers’ adoption rate is 0.35 who have only a 0.19 probability of adopting.

https://doi.org/10.1371/journal.pone.0204826.g002

Adjusted estimates of peer effects on adoption

In the fully-adjusted instrumental variables analysis, peer effects on adoption were largest in the patient-sharing network for all three drugs (Fig 3). For example, among anticoagulant prescribers, the coefficient for peer effects in the patient-sharing network was 0.590 (standard error (SE) = 0.150, p<0.001). This implies that for every 10 percentage point increase in the fraction of peers in the patient-sharing network adopting dabigatran, a physician’s own adoption probability increased by 5.90%. The coefficients on peer adoption in the patient-sharing network were of similar magnitude for the other two new medications. A 10 percentage point increase in peer adoption corresponded to a 8.32% (SE = 1.51%, p<0.001) increase in sitagliptin adoption, and a 7.84% increase in aliskiren adoption (SE = 2.93%, p<0.001) (S10 Table displays full model estimates).

thumbnail
Fig 3. Estimated peer effects on drug adoption with 95% confidence intervals.

Estimates show the effect of a 1% absolute change in the adoption rate of a physician’s peers on the likelihood of a physician’s decision to adopt the new drug. For example, a 1% increase in patient-sharing network peer adoption of dabigatran corresponded to a 0.59% increase in the probability of own adoption. Estimates are from a two-stage least squares regression model including physician-level characteristics: sex, medical school graduation year, US vs. non-US medical school, Top 20 U.S. medical school vs. not, geographic indicators (hospital referral region and metropolitan vs. non-metropolitan, prescription share paid for by Medicare, Medicaid fee-for-service, and cash, and age of patients filling prescriptions (see S10 Table for full model estimates). The means of peer characteristics serve as instruments for peer adoption rates. Included in the second stage were: physician specialty (primary care physician, relevant sub-specialty; e.g., cardiologist, endocrinologist, nephrologist vs. other), proportion of peers in each network who are in each specialty group, an indicator for whether the physician was a high-volume prescriber (> = median), the proportion of peers in each network who were high-volume prescribers, and an indicator for physicians who do not have peers in a particular type of network in the cohort. Peer effects estimates in the medical group and training networks are not reliable due to poor predictive power of the instruments.

https://doi.org/10.1371/journal.pone.0204826.g003

After adjusting for other peer effects and covariates, own adoption was positively influenced by peer adoption in the same hospital although the effect was not statistically significant for any of the drugs at p<0.05 (Fig 3). The coefficient for hospital peer adoption was 0.319, SE = 0.188, p = 0.09 for dabigatran; 0.376, SE = 0.211, p = 0.08 for sitagliptin; and 0.367, SE = 0.261, p = 0.16 for aliskiren. Estimates of medical group and training network effects were not reliable for any of the three drugs due to weak instruments.

Illustration of use of social networks to target interventions

Fig 4 shows the social multipliers using the adoption of dabigatran in the patient-sharing network to illustrate the potential for using social network information to target interventions. The figure presents the average multiplier by decile of peer connections (i.e., number of physicians with whom a physician shares patients) and by decile of prescribing volume and shows the value of targeting interventions using peer network information. Adoption by physicians in the top decile of connections is projected to induce 28 times as many other physicians to adopt compared with physicians in the bottom decile (4.81 vs. 0.16 physicians adopting). By contrast, targeting physicians in the top volume decile is projected to induce only twice as many adoptions as targeting physicians in the bottom decile (2.64 vs. 1.33 physicians adopting), and little more than half as many as targeting the top decile of connections (2.64 vs. 4.81).

thumbnail
Fig 4. Physicians’ average social multiplier by prescribing volume and number of connections in the patient-sharing network.

The social multiplier captures the number of other physicians who might adopt dabigatran following adoption by an individual physician in their peer network. Thus, the figure shows the simulated number of physicians who would adopt dabigatran per physician in each decile of prescribing volume (dark blue bars) and number of connections in the patient-sharing network (striped bars) adopting dabigatran. Physicians in the top decile of patient-sharing who had the most connections with other physicians who adopted dabigatran appear to have a nearly 2-fold stronger influence on adoption by other physicians compared to physicians in the top decile of prescribing volume (4.81 vs. 2.64).

https://doi.org/10.1371/journal.pone.0204826.g004

Discussion

We find that physician adoption of new drugs is heavily influenced by the extent to which their peers have adopted those new drugs. These effects were particularly large for peers with whom physicians share patients. Our study points to a potential mechanism underlying the tremendous variation in US medical care and to the importance of viewing physicians as part of a larger social system.

We sought to measure and test the effects of a rich set of peer influences derived from the institutional affiliations that physicians have with medical group practices, hospitals, and health systems. We also examined the influence of peers in the informal networks physicians form as they develop referral relationships, and interact with other physicians in the management of shared patients. Notably, peers in these informal patient-sharing networks had the largest effects on adoption of all three drugs after adjusting for the influence of peers in the same medical group and hospital setting. The patient-sharing network may exert the most influence on physician adoption because it captures more active connections over which physicians exercise the most control. A physician likely has more discretion, for example, over whom he refers his patients to than he does over the physicians with admitting privileges to the same hospital. This element of choice in the patient-sharing network raises the possibility that our estimates may be influenced in part by homophily[53, 54] although our use of instrumental variables reduces this as a concern. A key advantage of using patient-sharing information to measure peer networks is the routine availability of claims data to payers, the stakeholders best-positioned to invest in large-scale interventions to improve the quality of care.

The magnitude of the patient-sharing peer effects on physician adoption was broadly consistent across three drugs introduced over a 4-year time period that were prescribed to different patient populations, and by physicians with primary care and medical sub-specialty training. We studied the introduction of three newly-approved drugs that were first-in-class, each with a novel mechanism of action. All three medications have a place in the treatment armamentarium alongside older, existing therapies but varied with respect to whether they are considered first- or second-line treatments. That the magnitude of peer influence on adoption was comparable for all three of these drugs adds to the generalizability of our findings.

Meltzer and colleagues describe an approach to using social network analysis to form quality improvement teams to maximize the reach of interventions.[11] Similarly, our social multiplier exercise illustrates that targeting interventions to physicians who are well-connected in patient-sharing networks may be a more efficient way of improving the diffusion of evidence-based therapies. For example, health systems routinely adopt intensive educational interventions such as academic detailing—face-to-face educational programs borrowing principles from pharmaceutical promotional efforts. These educational interventions may reduce over-prescribing of ineffective medications, or increase use of highly effective agents.[7] Because academic detailing visits are meant to be in-depth and frequent[55] they are difficult to deliver on a large scale. However, systematic reviews suggest that academic detailing is seldom targeted.[56] Our study relates to an extensive literature on speeding diffusion of innovations in healthcare and other sectors by targeting ‘key opinion leaders’ who by virtue of their technical competence, credibility, and/or social acceptability are seen as influential with their peers.[5760] What our study adds to this body of work is information on the potential for using patient-sharing based measures of network centrality to target influential physicians at a larger scale than a single health care setting or community.[61] We focus on the most basic centrality measure, e.g., degree centrality, because it is simple to construct facilitating its widespread use and does not require complete observation of the network. Other centrality measures such as eigenvector centrality, Bonacich centrality, PageRank centrality, density centrality have been used to characterize other technological diffusion processes[51, 6264] and may be useful in this context.

Our study improves on prior work by including a larger physician sample, examining the diffusion of multiple new drugs, drawing on several information sources to form 4 types of physician networks, and using an instrumental variables approach to overcome some challenges in estimating network effects. Nevertheless, we note some limitations. First, although Pennsylvania resembles national averages on most measures of healthcare utilization,[25] our findings are limited to a single state and estimates of peer effects on physician adoption may not generalize to other geographic areas with different physician network structures[38] or adoption patterns. Second, while our study of adoption of three new prescription drugs improves on prior studies of uptake of a single technology, our estimates of peer effects may not generalize to all medical technologies. Third, while we adjusted for several physician characteristics associated with adoption we lack information on important sources of influence on physician behavior, namely patient clinical characteristics, pharmaceutical company promotion, and payer formularies. Only having access to information on the age and source of payment for a physician’s patients who filled their prescriptions we are unable to adjust for differences in adoption due to patient health state. We note that all three of the new drugs were first-in-class and likely to be heavily promoted by manufacturers; however, our study predates availability of physician-level data on exposure to industry promotion. We were not able to adjust for differences in formulary coverage of the new drugs for a physician’s patient panel. Fourth, as with any study that uses instrumental variables, our estimates capture the effects of peer adoption rates that were driven by the peer mean characteristics we use as instruments, not by other sources of variation, which potentially limits their generalizability.

Using multiple data sources to measure the rich set of peer relationships formed by physicians we find that peers can exert significant influence over physician technology adoption decisions. Our study shows the potential for using information on physician social networks routinely available to health systems to improve the targeting of interventions to speed the diffusion of evidence-based health care technologies.

Supporting information

S1 Fig. Sample size flow chart for anticoagulant prescribers.

Sources: QuintilesIMS’s XPonent, QuintilesIMS’s HCOS, and AMA Masterfile Abbreviations: AC: Anticoagulant; PA: Pennsylvania. Notes: Dabigratran was introduced in October 2010 and our XPonent data end in December 2011. We included physicians with at least some minimal anticoagulant prescribing, although not necessarily the drug of interest, and an AMA Masterfile record and HCOS record for the maximum duration given data availability (15 months post-market introduction).

https://doi.org/10.1371/journal.pone.0204826.s001

(DOCX)

S2 Fig. Sample size flow chart for antidiabetic prescribers.

Sources: QuintilesIMS ‘s XPonent, QuintilesIMS’s HCOS, and AMA Masterfile Abbreviations: AD: Antidiabetic; PA: Pennsylvania. Notes: Sitagliptin was introduced in October 2006 a few months before our data was available. We included physicians with at least some minimal antidiabetic prescribing, although not necessarily the drug of interest, and an AMA Masterfile and HCOS record from months 3 to 15 following sitagliptin’s introduction, the period during which adoption is measured.

https://doi.org/10.1371/journal.pone.0204826.s002

(DOCX)

S3 Fig. Sample size flow chart for antihypertensive prescribers.

Sources: QuintilesIMS’s XPonent, QuintilesIMS’s HCOS, and AMA Masterfile Abbreviations: AH: Antihypertensive; PA: Pennsylvania. Notes: Aliskiren was introduced in March 2007. To provide a comparable period over which adoption is measured for all classes we included physicians with at least some minimal antihypertensive prescribing, although not necessarily the drug of interest, and an AMA Masterfile record and HCOS record for 15 months following aliskiren’s introduction.

https://doi.org/10.1371/journal.pone.0204826.s003

(DOCX)

S1 Table. A) Oral Anti-coagulant products; B) Antihypertensive products in ACEI, ARB, and RAAS classes; C) Oral anti-diabetic medications.

https://doi.org/10.1371/journal.pone.0204826.s004

(DOCX)

S2 Table. Distribution of number of prescriptions for the new drugs among physicians adopting them.

https://doi.org/10.1371/journal.pone.0204826.s005

(DOCX)

S3 Table. Characteristics of anticoagulant prescriber cohort and comparisons by adoption vs. non-adoption of new drug (dabigatran).

Data sources: QuintilesIMS, HCOS; XPonent; AMA Masterfile 1 age of physician in 2010.

https://doi.org/10.1371/journal.pone.0204826.s006

(DOCX)

S4 Table. Characteristics of anti-diabetes medication prescriber cohort and comparisons by adoption vs. non-adoption of new drug (sitagliptin).

Data sources: QuintilesIMS, HCOS; XPonent; AMA Masterfile * Age in 2007.

https://doi.org/10.1371/journal.pone.0204826.s007

(DOCX)

S5 Table. Characteristics of antihypertensive prescriber cohort and comparisons by adoption vs. non-adoption of new drug (aliskiren).

Data sources: QuintilesIMS, HCOS; XPonent; AMA Masterfile. *Age in 2007.

https://doi.org/10.1371/journal.pone.0204826.s008

(DOCX)

S6 Table. Sample size of physicians in each prescribing cohort with Medicare and/or Medicaid claims during adoption measurement period.

Sources: Medicare data were obtained from CMS. Medicaid data were obtained from the Pennsylvania Department of Human Services. Notes: In the anticoagulant cohort, there are 7,785 physicians meeting inclusion criteria of whom 7,522 (96.6%) had Medicare claims and 6,680 (85.8%) had Medicaid claims. We included claims submitted by those physicians to Medicare and Medicaid with dates of service between 10/1/2010 and 12/31/2011 to match the period over which we measured adoption of dabigatran. In the antidiabetic cohort, 8,257 physicians met inclusion criteria. We included claims from Medicare and Medicaid submitted by those physicians with dates of service between 1/1/2007 and 1/31/2008 (the measurement period for adoption of sitagliptin). There are 9,974 physicians meeting inclusion criteria for the antihypertensive prescriber cohort. We included claims from Medicare or Medicaid with dates of service between 3/1/2007 and 5/31/2008, the measurement period for adoption of aliskiren.

https://doi.org/10.1371/journal.pone.0204826.s009

(DOCX)

S7 Table. Number of physician peers to which physicians are connected in patient-sharing, medical group, hospital and training networks.

Data sources: QuintilesIMS, HCOS; XPonent; AMA Masterfile 1 Column displays the number of physicians who are not connected to any peers in a particular network. For example, physicians may lack peers in the patient-sharing network because they do not see patients with Medicaid and Medicare patients or because they did not share any patients with that source of coverage with other physicians in the prescribing cohort.

https://doi.org/10.1371/journal.pone.0204826.s010

(DOCX)

S8 Table. Assessments of instrument exogeneity and relevance.

https://doi.org/10.1371/journal.pone.0204826.s011

(DOCX)

S9 Table. Peer characteristics used as instruments in each network.

Data sources: QuintilesIMS, HCOS; XPonent; AMA Masterfile *In the analytical dataset, the unit is physician. For categorical variables, we obtain the proportion of each type (e.g., proportion of peers who are female) for each physician. For the continuous variables, we use the average percentage of pay type (e.g., Medicare) of all peers for each physician. Table shows the proportion of physicians in each network with that characteristic.

https://doi.org/10.1371/journal.pone.0204826.s012

(DOCX)

S10 Table. Estimates of peer effects and covariates on adoption of new drugs from instrumental variables regression.

Primary analysis showing all covariates. Data sources: QuintilesIMS, HCOS; XPonent; AMA Masterfile Notes: Robust standard errors in parentheses *** p<0.001, ** p<0.01, * p<0.05 ahigh-volume = total prescribing volume in the category above the median.

https://doi.org/10.1371/journal.pone.0204826.s013

(DOCX)

S11 Table. Estimates of peer effects on adoption of new drugs from linear model without instrumental variables.

Data sources: QuintilesIMS, HCOS; XPonent; AMA Masterfile Notes: Robust standard errors in parentheses *** p<0.001, ** p<0.01, * p<0.05.

https://doi.org/10.1371/journal.pone.0204826.s014

(DOCX)

S12 Table. Estimates of peer effects and covariates on adoption of new drugs.

Alternate specification of own and peer adoption (using> = 1 and 15). Data sources: QuintilesIMS, HCOS; XPonent; AMA Masterfile Notes: Robust standard errors in parentheses *** p<0.001, ** p<0.01, * p<0.05.

https://doi.org/10.1371/journal.pone.0204826.s015

(DOCX)

S13 Table. Distribution of aggregate social multipliers on individual adoption by drug and network.

https://doi.org/10.1371/journal.pone.0204826.s016

(DOCX)

Acknowledgments

We gratefully acknowledge technical support from the Health Services Research Data Center at the University of Pittsburgh, and Aiju Men, Jie Li and Joo Yeon Kim at the University of Pittsburgh for expert programming.

References

  1. 1. Neumann PJ, Rosen AB, Weinstein MC. Medicare and cost-effectiveness analysis. N Engl J Med. 2005;353:1516–22. pmid:16207857
  2. 2. Fisher ES, Wennberg DE, Stukel TA, Gottlieb DJ, Lucas FL, Pinder EL. The Implications of Regional Variations in Medicare Spending. Part 1: The Content, Quality, and Accessibility of Care. Annals of Internal Medicine. 2003;138(4):273–87. pmid:12585825
  3. 3. Morris ZS, Wooding S, Grant J. The answer is 17 years, what is the question: understanding time lags in translational research. J R Soc Med. 2011;104(510–520). pmid:22179294
  4. 4. Institute of Medicine. Best care at lower cost: The path to continuously learning health care in America. Washington, DC.: 2013.
  5. 5. Institute of Medicine. Variation in health care spending: target decision making, not geography. Washington DC: Institute of Medicine, 2013. https://doi.org/10.17226/18393 pmid:24851301
  6. 6. Davis D, O’Brien MA, Freemantle N, Wolf FM, Mazmanian P, Taylor-Vaisey A. Impact of formal continuing medical education: do conferences, workshops, rounds and other traditional continuing education activities change physician behavior or health outcomes? JAMA. 1999;282(9):867–74. pmid:10478694
  7. 7. Soumerai SB, Avorn J. Principles of educational outreach ('academic detailing') to improve clinical decision making. JAMA. 1990;263(4):549–56. pmid:2104640.
  8. 8. Grimshaw JM, Shirran L, Thomas R, Mowatt G, Fraser C, Bero LA, et al. Changing provider behavior: an overview of systematic reviews of interventions. Medical Care. 2001;39(8, supplement 2):II2–II45.
  9. 9. McGettigan P, Golden J, Chan R, Feely J. Prescribers prefer people: the sources of information used by doctors for prescribing suggest that the medium is more than the message. Br J Clin Pharmacol. 2001;51(2):184–9. pmid:11259994
  10. 10. Keating N. Peer Influence and Opportunities for Physician Behavior Change. J Natl Cancer Instit. 2017;109(8). Epub March 14, 2017.
  11. 11. Meltzer D, Chung J, Khalili P, Marlow E, Arora V, Schumock G, et al. Exploring the use of social network methods in designing healthcare quality improvement teams. Soc Sci Med 2010;71:1119–30. pmid:20674116
  12. 12. Chambers D, Wilson P, Thompson C, Harden M. Social network analysis in healthcare settings: a systematic scoping review. PLoS One. 2012;7(8):e41911. Epub Aug 3. pmid:22870261
  13. 13. Wasserman S, Faust K. Social Network Analysis: Methods and Applications. Cambridge, MA.: Cambridge University Press; 1994.
  14. 14. Keating N, Ayanian J, Cleary P, Marsden P. Factors Affecting Influential Discussions Among Physicians: A Social Network Analysis of a Primary Care Practice. Journal of General Internal Medicine. 2007;22(6):794–98. pmid:17404798
  15. 15. Fattore G, Frostini F, Salvatore D, Tozzi V. Social network analysis in primary care: the impact of interventions on prescribing behavior. Health Policy 2009;92:141–8. pmid:19356822
  16. 16. Iyengar R, Van de Bulte C, Valente TW. Opinion leadership and social contagion in new product diffusion Marketing Science. 2011;30(2):195–212.
  17. 17. Coleman J, Katz E, Menzel H. The diffusion of an innovation among physicians. Sociometry. 1957;20(4):253–70.
  18. 18. Bhatia T, Wang L. Identifying physician peer-to-peer effects using patient movement data. International Journal of Research in Marketing. 2011;28:51–61.
  19. 19. Van den Bulte C, Lilien G. Medical innovation revisited: social contagion versus marketing effort. American Journal of Sociology. 2001;106(5):1409–35.
  20. 20. Manchanda P, Xie Y, Youn N. The role of targeted communication and contagion in product adoption. Marketing Science. 2008;27(6):961–76.
  21. 21. Nair HS, Manchanda P, Bhatia T. Asymmetric social interactions in physician prescription behavior: the role of opinion leaders. Journal of Marketing Research 2010;XLVII:883–95.
  22. 22. Pollack C, Soulos P, Herrin J, Xu X, Christakis N, Forman H, et al. The Impact of Social Contagion on Physician Adoption of Advanced Imaging Tests in Breast Cancer. J Natl Cancer Instit. 2017;109(8). Epub March 14, 2017.
  23. 23. Zheng K, Padman R, Krackhardt D, Johnson MP, Diamond HS. Social networks and physician adoption of electronic health records: insights from an empirical study. J Am Med Inform Assoc. 2010;17(3):328–36. pmid:20442152
  24. 24. Mascia D, Cicchetti A. Physician social capital and the reported adoption of evidence-based medicine: exploring the role of structural holes. Soc Sci Med. 2011;72(5):798–805. pmid:21306807
  25. 25. Dartmouth Atlas of Health Care. [cited 2010]. http://www.dartmouthatlas.org/Downloads.aspx
  26. 26. Kaiser Family Foundation. Health Insurance Coverage of the Total Population 2015 [cited 2017 August]. http://www.kff.org/other/state-indicator/total-population/?currentTimeframe=0&sortModel=%7B%22colId%22:%22Location%22,%22sort%22:%22asc%22%7D.
  27. 27. Connolly SJ, Ezekowitz MD, Yusuf S, Eikelboom J, Oldgren J, Parekh A, et al. Dabigatran versus Warfarin in Patients with Atrial Fibrillation. NEJM. 2009;361:1139–51. pmid:19717844
  28. 28. Yusuf S, Teo KK, Pogue J, Dyal L, Copland I, Schumacher H, et al. Telmisartan, ramipril, or both in patients at high risk for vascular events. N Engl J Med. 2008;358:1547–59. pmid:18378520
  29. 29. Desai NR, Shrank WH, Fischer MA, Avorn J, Liberman JN, Schneeweiss S, et al. Patterns of medication initiation in newly diagnosed diabetes mellitus: quality and cost implications. Am J Med. 2012;125(3):302.e1–7.
  30. 30. Norris S, Weinstein J, Peterson K, Thakurta S. Drug class review: direct renin inhibitors, angiotensin converting enzyme inhibitors, and angiotensin II receptor blockers: Drug Effectiveness Review Project, Oregon Evidence-based Practice Center, Oregon Health & Science University; 2010.
  31. 31. Menchik D, Meltzer D. The Cultivation of Esteem and Retrieval of Scientific Knowledge in Physician Networks. Journal of Health and Social Behavior. 2010;51:137–52. pmid:20617755
  32. 32. Scott J, Talia A, Crosson JC, Orzana AJ, Stroebel C, DiCicco-Bloom B, et al. Social Network Analysis as an Analytic Tool for Interaction Patterns in Primary Care Practices. Ann Fam Med. 2005;3:443–8. pmid:16189061
  33. 33. Lewis JM, Baeza JI, Alexander D. Partnerships in primary care in Australia: Network structure, dynamics and sustainability. Social Science & Medicine. 2008;67(2):280–91.
  34. 34. Prosser H, Walley T. New drug prescribing by hospital doctors: the nature and meaning of knowledge. Soc Sci Med. 2006;62(1565–1578). pmid:16199121
  35. 35. Prosser H, Almond S, Walley T. Influences on GP’s decision to prescribe new drugs—the importance of who says what. Family Practice. 2003;20(1):61–8. pmid:12509373
  36. 36. Csardi G, Nepusz T. The igraph software package for complex network research. Inter-Journal, Complex Systems. 2006.
  37. 37. Barnett ML, Landon BE, O’Malley AJ, Keating NL, Christakis NA. Mapping physician networks with self-reported and administrative data. Health Services Research 2011;46 1592–609. pmid:21521213
  38. 38. Landon BE, Keating NL, Barnett ML, Jukka-Pekka O, Paul S, O’Malley AJ, et al. Variation in patient-sharing networks of physicians across the United States. JAMA. 2012;308(3):265–73. pmid:22797644
  39. 39. Barnett ML, Christakis NA, O’Malley J, Onnela JP, Keating NL, Landon BE. Physician patient-sharing networks and the cost and intensity of care in US hospitals. Med Care 2012;50(2):152–60. pmid:22249922
  40. 40. Casalino LP, Pesko MF, Ryan AM, Nyweide DJ, Iwahyna TJ, Sun X, et al. Physician Networks and Ambulatory Care-Sensitive Admissions. Med Care. 2015;53(6):534–41. pmid:25906013
  41. 41. Pollack CE, Weissman G, Bekelman J, Liao K, Armstrong K. Physician social networks and variation in prostate cancer treatment in three cities. Health Services Research. 2012;47(1 Pt 2):380–403. pmid:22092259
  42. 42. Moen EL, Austin AM, Bynum JP, Skinner JS, O’Malley AJ. An analysis of patient-sharing physician networks and implantable cardioverter defibrillator therapy. Health Serv Outcomes Res Methodol. 2016;16(3):132–53. pmid:27597812
  43. 43. Altalib HH, Fenton BT, Cheung KH, Pugh MJV, Bates J, Valente TW, et al. Care coordination in epilepsy: Measuring neurologists’ connectivity using social network analysis. Epilepsy Behav. 2017;73:31–5. pmid:28605631
  44. 44. Manski CF. Identification of Endogenous Social Effects: The Reflection Problem. Review of Economic Studies. 1993;60(3):531–42.
  45. 45. Moffitt RA. Policy Interventions, Low-Level Equilibria, and Social Interactions. In: Durlauf S, Young H, editors. Social Dynamics 2001. p. 45–82.
  46. 46. Blume LE, Brock WA, Durlauf SN, Loannides YM. Identification of Social Interactions. In: Benhabib J, Bisin A, Jackson MO, editors. Handbook of Social Economics 2011
  47. 47. Gaviria A, Raphael S. School-based peer effects and juvenile behavior. Review of Economics and Statistics. 2001;83(2):257–68.
  48. 48. Wooldridge JM. Econometric analysis of cross section and panel data: MIT press; 2010.
  49. 49. Stock JH, Yogo M. Testing for weak intruments in linear IV regression. In: Donald WA, Stock JH, editors. Identification and inference for econometric models: Essays in honor of Thomas Rothenberg: Cambridge University Press; 2005.
  50. 50. Becker GS, Murphy K. Social Economics: Harvard University Press; 2000.
  51. 51. Bonacich P. Power and centrality: a family of measures. American Journal of Sociology. 1987;92(5):1170–82.
  52. 52. Christakis NA, Fowler JH. Contagion in prescribing behavior among networks of doctors. Marketing Science. 2010;epub ahead of print.
  53. 53. McPherson M, Smith-Lovin L, Cook JM. Birds of a feather: homophily in social networks Annu Rev Sociol. 2001;27:415–44.
  54. 54. Centola D. An experimental study of homophily in the adoption of health behavior. Science. 2011;334:1269–72. pmid:22144624
  55. 55. Yeh JS, Van Hoof TJ, Fischer MA. Key features of academic detailing: Development of an expert consensus using the Delphi method. Am Drug Health Benefits. 2016;9(1):42–50.
  56. 56. Van Hoof TJ, Harrison LG, Miller NE, Pappas MS, Fischer MA. Characteristics of academic detailing: Results of a literature review. Am Health Drug Benefits. 2015;8(8):414–22. pmid:26702333
  57. 57. Valente TW. Network interventions. Science. 2012;337(6090):49–53. pmid:22767921
  58. 58. Valente TW, Davis RL. Accelerating the diffusion of innovations using opinion leaders. The Annals of the American Academy of Political and Social Science. 1999;566(1):55–67.
  59. 59. Flodgren G, Parmelli E, Doumit G, Gattellari M, O’Brian MA, Grimshaw J, et al. Local opinion leaders: effects on professional practice and health care outcomes. Cochrane Database Syst Rev. 2011;8.
  60. 60. Liu QH, Lu FM, Tang M, Zhou T. Impacts of opinion leaders on social contagion. Chaos: An Interdisciplinary Journal of Nonlinear Science. 2018;28(5).
  61. 61. Lomas J, Enkin M, Anderson GM, Hannah WJ, Vayda E, Singer J. Opinion leaders vs audit and feedback to implement practice guidelines: delivery after previous cesarean section. JAMA. 1991;265(17):2202–7. pmid:2013952
  62. 62. Brin S, Page L. The anatomy of a large-scale hypertextual Wed search engine. Computer Networks and ISDN Systems. 1998;30(107–117).
  63. 63. Kang C, Kraus S, Molinaro C, Spezzzano F, Subrahamanian VS. Diffusion centrality: a paradigm to maximize spread in social networks Artificial Intelligence. 2016;239(70–96).
  64. 64. Banerjee A, Chandrasekhar AG, Duflo E, Jackson MO. The diffusion of microfinance. Science. 2013;341(6144).