Research article | Open | Open Peer Review | Published:
Variation in plasma calcium analysis in primary care in Sweden - a multilevel analysis
BMC Family Practicevolume 11, Article number: 43 (2010)
Primary hyperparathyroidism (pHPT) is a common disease that often remains undetected and causes severe disturbance especially in postmenopausal women. Therefore, national recommendations promoting early pHPT detection by plasma calcium (P-Ca) have been issued in Sweden. In this study we aimed to investigate variation of P-Ca analysis between physicians and health care centres (HCCs) in primary care in county of Skaraborg, Sweden.
In this cross sectional study of patients' records during 2005 we analysed records from 154 629 patients attending 457 physicians at 24 HCCs. We used multilevel logistic regression analysis (MLRA) and adjusted for patient, physician and HCC characteristics. Differences were expressed as median odds ratio (MOR).
There was a substantial variation in number of P-Ca analyses between both HCCs (MORHCC 1.65 [1.44-2.07]) and physicians (MORphysician 1.95 [1.85-2.08]). The odds for a P-Ca analysis were lower for male patients (OR 0.80 [0.77-0.83]) and increased with the number of diagnoses (OR 25.8 [23.5-28.5]). Sex of the physician had no influence on P-Ca test ordering (OR 0.93 [0.78-1.09]). Physicians under education ordered most P-Ca analyses (OR 1.69 [1.35-2.24]) and locum least (OR 0.73 [0.57-0.94]). More of the variance was attributed to the physician level than the HCC level. Different mix of patients did not explain this variance between physicians. Theoretically, if a patient were able to change both GP and HCC, the odds of a P-Ca analysis would in median increase by 2.45. Including characteristics of the patients, physicians and HCCs in the MLRA model did not explain the variance.
The physician level was more important than the HCC level for the variation in P-Ca analysis, but further exploration of unidentified contextual factors is crucial for future monitoring of practice variation.
Primary hyperparathyroidism (pHPT) is a common disease that often remains undetected and causes severe disturbance especially in postmenopausal women. Therefore, national recommendations promoting early pHPT detection by plasma calcium (P-Ca) have been issued in Sweden [1, 2]. In this study we aimed to investigate variation of P-Ca analysis between physicians and health care centres (HCCs).
pHPT is a potentially serious condition leading to increased morbidity and mortality from cardiovascular disease  and cancer . Although mild disease might not increase the risk . PHPT gives raised plasma calcium (P-Ca) and because of the vague symptoms, pHPT is difficult to detect without an analysis of P-Ca.
Even though previous studies indicate that the frequency of P-Ca analyses differs between health care centres (HCC)  the understanding of the relative importance of the different levels (patients, physicians, HCCs) for these differences is limited. However, in a study from New Zealand where a defined clinical situation was presented to GPs, it was shown that inherent characteristics of the physicians more than the patients clinical situation determined which laboratory tests were ordered . A study from the Netherlands, without patient characteristics, found a regional variation in laboratory testing and that factors at both the physician and HCC level influenced the inclination to order tests .
The aim of this study was to investigate the relative importance of the different levels in the health care organization for P-Ca analyses using the Skaraborg Primary Care Database (SPCD). Identification of factors contributing to the variation can be of relevance for planning interventions for an optimal frequency of P-Ca analyses and for evaluating the national recommendations.
Skaraborg is a rural area in Sweden and comprised 255 758 inhabitants in 2005. The public primary care is a part of the Västra Götaland region and serves 97% of the population (n = 247 985). All the HCCs (n = 24) use the same computerised medical record, ProfDoc Journal III (PDIII ProfDoc AB: Uppsala) facilitating data extraction. SPCD has been created containing encrypted data from patients and caregivers from all HCCs. The database contains patients' age, sex, diagnoses, laboratory analyses, and drug prescriptions. The HCCs' laboratory facilities are accredited by SWEDAC (the Swedish Board for Accreditation and Conformity Assessment). The validity of the information in the database has recently been audited and judged to be mostly appropriate but varying with type of diagnosis . All 154 629 individuals that attended any of the 24 HCCs during 2005 are included in the analysis. The local ethics committee at Gothenburg University approved the study (255-09).
Study procedure and assessment of variables
The outcome variable was P-Ca analyses during 2005 (yes/no). Sex of the patient and P-Ca analyses during 2004, were included as explanatory variables. We also selected ICD-10 coded diagnoses and symptoms associated with pHPT  . A risk score for a P-Ca analysis was created with stepwise logistic regression  based on age, concomitant diagnosis and drug treatment, in order to control for confounding factors. The risk score was divided in quintiles, patients with the lowest risk of P-Ca analyses (group 1) were used as reference. The main characteristics included in the risk score are listed in Table 1 A more detailed description can be found in additional file 1.
The physicians were categorised according to sex and title. GP and locum were also dichotomised at 46 year. GPs, 46 years or older, were used as reference in the analysis. As only six doctors among house officers and preregistration house officers were above 45 years, they were not dichotomised.
The HCCs had different standardised group analyses, for instance analyses of electrolytes, hypertension check ups and diagnosing dementia, in which P-Ca was included. We categorised HCCs as having none, 1-2, and ≥ 3 standardised groups including P-Ca. The HCCs having no group analyses were used as reference.
We used multilevel logistic regression analysis (MLRA) to estimate the odds of patients being ordered a P-Ca analysis, as the data had a hierarchical structure (i.e., patients nested within physicians nested within HCCs) [12, 13]. As one patient could attend several physicians and several HCCs, we used a multiple membership model (Figure 1)  . The weights were constructed according to number of visits to a certain physician/HCC during our study period.
We developed four consecutive models. Model A included the random parameters (physicians and HCCs), in order to partition the variance at different levels. Model B included the patient characteristics, model C the patient and physician characteristics and model D the patient, physician and HCC characteristics. In this way we could investigate whether the contextual characteristics explained the residual variation at the physician and HCC levels.
In the fixed-effects part of the MLRA, we calculated odds ratios (OR) and their 95% credible intervals (95% CI). In the random-effects part of the MLRA, we obtained the variance at the physician and HCC levels. In order to quantify the importance of the different levels in the analysis we calculated the median odds ratio (MOR) [15, 16]. The MOR translates the variance into the widely used OR scale, and can thereby be directly compared with the ORs of individual or area variables. In very simple terms, the MOR could be interpreted as how much a patient's odds of being ordered a P-Ca analysis would (in median) increase if this patient moved to a physician/HCC with higher odds of ordering a P-Ca analysis. A MOR of one indicates that there are no differences between physicians/HCCs in their odds of ordering P-Ca analysis. The larger the differences between physicians/HCCs are, the larger the MOR will be. The MORphysician+HCC is applicable to patients that visited only one physician during the study period. Parameters were estimated by MCMC methods  in the MLwiN 2.02 software .
Overall 5.8% of the inhabitants in Skaraborg and 9% of the patients (11% of the women and 8% of the men) attending the HCCs had a P-Ca analysis. The mean age of the patients with P-Ca analysis was 62 years compared to 45 years for patients with no P-Ca analyses. At the different HCCs the number of standardised group analyses including P-Ca analyses varied from zero to seven. The locums were most numerous shown in Table 2, but had short periods of attendance.
The multilevel model
There was a substantial variation in number of P-Ca analyses between HCCs and physicians. The four models used in the analyses are shown in Table 3. In model A the MORphysician+HCC indicated that for a patient changing both GP and HCC, to a GP and HCC with higher odds for a P-Ca analysis, the odds would in median increase by 2.31. The physician level, MORphysician = 1. 95 (95% CI: 1. 85-2.08) contributed more than the HCC level, MORHCC = 1.65 (95% CI: 1.44-2.07). Figure 2 shows the residuals for physicians (Panel A) and for HCCs (Panel B) from the multilevel analysis.
Model B, C and D
Model B illustrates that female sex and increased number of diagnoses in the risk score was associated with higher propensity of a P-Ca analysis However, the inclusion of other patient characteristics did not explain the variance at the physician or HCC level. Model C illustrates that house officers, preregistration house officer and younger GPs ordered more P-Ca analysis compared to older GPs. Locums, regardless of age, ordered fewer P-Ca analyses. There were no differences between male and female physicians. Inclusion of the physician characteristics did not explain the variance in model A. Model D illustrates that a high number of standardised group analyses were associated with a high number of P-Ca analyses. The inclusion of HCC characteristics did not explain the variance at the higher levels.
Including all explanatory variables and controlling for confounders, a patient changing both GP and HCC, from low to high odds for P-Ca analysis, the odds for a P-Ca analysis would in median increase by 2.5 times, MORphysician+HCC 2.45.
The main finding of this study was that the ordering of P-Ca analyses was influenced by factors both at the physician and at the HCC level, with the physician level being more important than the HCC level. Theoretically, if a patient were able to change both GP and HCC, the odds of undergoing a P-Ca analysis would in median increase by 2.45 times. Including compositional and contextual characteristics in the model did not explain the variance at the higher levels.
Overall 5,8% of the inhabitants underwent a P-Ca analysis, which is comparable with an earlier study from Skaraborg (6,1%)  and two-fold compared with a study from primary care in southern Stockholm  Female patients and patients with previous P-Ca analysis were more likely to have a P-Ca analysis, which could be explained by women's greater risk of pHPT and recurrent check-ups of patients with chronic diseases.
In order to control for compositional confounding at the patient level we included an individual risk score for P-Ca-analysis. The inclusion of this variable did not explain the variation between physicians and between HCCs. Further, our empirical analysis found that the sex of the physician had no influence on P-Ca test ordering, in contrast to a study from Israel where female physicians ordered more test . Older and more experienced physicians were less likely to order a P-Ca-test, which is in line with previous studies indicating that test ordering behaviour of GPs was influenced by years of experience . P-Ca analyses done as part of group analyses used in surveillance of different chronic conditions may inflate the number of P-Ca analyses . However, even though the number of group analyses was associated with higher frequency of P-Ca analysis, it could not explain the variation at the HCC level.
As explained in previous studies ; the measures of variation (e.g. median odds ratio) should be interpreted only for the specific time and place of the study, as there may be pattern of variance produced by different conditions. The associations, however, between characteristics of, on the one hand, the patients, physicians, and HCCs and on the other the frequency of P-Ca analysis, intend to provide information that can be generalised and applied to contexts beyond the one where the study was performed.
The risk for selection bias is low since this study is based on a large sample from a primary care area serving 97% of the population. Moreover, as this study is a retrospective database study, the ordering of analyses is not influenced by the study. A limitation of the study is that the frequency of ICD coded patient visits varies both between HCCs and according to diagnosis . This might affect the risk score calculation.
Different views of the reason for screening could also affect the result. However national recommendations are well known in Swedish primary care [1, 2] thus the risk for bias is minor. Due to regional variation in laboratory testing  the results from this study might not be applicable in all regions in Sweden.
In this study only the variables available in the SPCD database were included. In previous studies, other characteristics of the physician, such as attitude to risk taking and involvement in development of guidelines, explained parts of the higher level variance .
We found that there was variation between physicians and between HCC in ordering of P-Ca analysis, which is in line with previous studies . However, in this study we also tried to quantify the contribution of each level by using the median odds ratio. Even though our multilevel approach identified factors, both at the physician and HCC level, which are important to consider for understanding the inclination to order a P-Ca test, none of the included variables could explain the variation at the higher level. The identification of yet unidentified factors that contribute to the variation is needed for monitoring of practice variation and quality assessment and for applying appropriate interventions to achieve optimal frequency of P-Ca analyses.
National recommendations in Sweden have been issued to increase the frequency of P-Ca analyses to detect more patients with pHPT. There is a substantial variation in number of P-Ca analyses primarily between physicians but also between Health Care Centres. Female sex of the patient and increasing number of diagnoses is associated with higher propensity of P-Ca analysis. Physicians under education order most P-Ca analyses and locum least, but sex of the physician has no influence.
Tryding N: Laborera rätt och lagom i primärvården- klinisk kemi i primärvården (Correct and Optimal use of Clinical Chemistry in Primary Care, English summary). 1999, Stockholm: SPRI
Laborera rätt och lagom - klinisk kemi i primärvården, analysdatabasen, (Correct and Optimal use of Clinical Chemistry in Primary Care, English version. [http://www3.svls.se/sektioner/sfkk/nyckel/]
Ogard CG, Engholm G, Almdal TP, Vestergaard H: Increased mortality in patients hospitalized with primary hyperparathyroidism during the period 1977-1993 in Denmark. World J Surg. 2004, 28: 108-111. 10.1007/s00268-003-7046-0.
Nilsson IL, Zedenius J, Yin L, Ekbom A: The association between primary hyperparathyroidism and malignancy: nationwide cohort analysis on cancer incidence after parathyroidectomy. Endocr Relat Cancer. 2007, 14: 135-140. 10.1677/erc.1.01261.
Wermers RA, Khosla S, Atkinson EJ, Grant CS, Hodgson SF, O'Fallon WM, Melton LJ: Survival after the diagnosis of hyperparathyroidism: a population-based study. Am J Med. 1998, 104: 115-122. 10.1016/S0002-9343(97)00270-2.
Dalemo S, Hjerpe P, Bostrom Bengtsson K: Diagnosis of patients with raised serum calcium level in primary care, Sweden. Scand J Prim Health Care. 2006, 24: 160-165. 10.1080/02813430600784443.
Malcolm L, Wright L, Seers M, Davies L, Guthrie J: Laboratory expenditure in Pegasus Medical Group: a comparison of high and low users of laboratory tests with academics. N Z Med J. 2000, 113: 79-81.
Verstappen WH, ter Riet G, Dubois WI, Winkens R, Grol RP, Weijden van der T: Variation in test ordering behaviour of GPs: professional or context-related factors?. Fam Pract. 2004, 21: 387-395. 10.1093/fampra/cmh408.
Hjerpe P, Merlo J, Ohlsson H, Bengtsson Boström K, Lindblad U: Validity of registration of ICD codes and prescriptions in a research database in Swedish primary care - A cross-sectional study in Skaraborg primary care database. BMC Med Inform Decis Mak. 10: 23-10.1186/1472-6947-10-23.
Bilezikian JP, Brandi ML, Rubin M, Silverberg SJ: Primary hyperparathyroidism: new concepts in clinical, densitometric and biochemical features. J Intern Med. 2005, 257: 6-17. 10.1111/j.1365-2796.2004.01422.x.
Merlo J, Gerdtham UG, Eckerlund I, Hakansson S, Otterblad-Olausson P, Pakkanen M, Lindqvist PG: Hospital level of care and neonatal mortality in low- and high-risk deliveries: reassessing the question in Sweden by multilevel analysis. Med Care. 2005, 43: 1092-1100. 10.1097/01.mlr.0000182484.14608.b9.
Goldstein H: Multilevel Statistical Models. 2003, London: Hodder Arnold
Snijders T, Bosker R: Multilevel Analysis An Intro- duction to Basic and Advanced Multilevel Modeling. 1999, London: Sage
Fielding A, Goldstein H: Cross-classified and Multiple Membership Structures in Multilevel Models: An Introduction and Review. 2006, Birmingham: University of Birmingham
Larsen K, Merlo J: Appropriate assessment of neighborhood effects on individual health: integrating random and fixed effects in multilevel logistic regression. Am J Epidemiol. 2005, 161: 81-88. 10.1093/aje/kwi017.
Larsen K, Petersen JH, Budtz-Jorgensen E, Endahl L: Interpreting parameters in the logistic regression model with random effects. Biometrics. 2000, 56: 909-914. 10.1111/j.0006-341X.2000.00909.x.
Gilks WR, Richardson S, Spiegelhalter DJ: Markov Chain Monte Carlo in Pratice. 1996, London, Chapman and Hall
Rasbash JSF, Browne W: A User's Guide to MLwiN, Version 2.0. Documentation Version 2.1e. 2003, London: Centre for Multilevel Modelling, Institute of Education, University of London
Borgquist L, Gustafsson S, Hultén G, Jansson U, Paulsson E, Tryding N: Klinisk kemi i primärvården (Clinical chemistry in primary care, English summary). 1996, Stockholm: SPRI, 422-
Vinker S, Kvint I, Erez R, Elhayany A, Kahan E: Effect of the characteristics of family physicians on their utilisation of laboratory tests. Br J Gen Pract. 2007, 57: 377-382.
Salloum S, Franssen E: Laboratory investigations in general practice. Can Fam Physician. 1993, 39: 1055-1061.
Larsson A, Biom S, Wernroth ML, Hulten G, Tryding N: Effects of an education programme to change clinical laboratory testing habits in primary care. Scand J Prim Health Care. 1999, 17: 238-243. 10.1080/028134399750002476.
Merlo J, Ohlsson H, Lynch KF, Chaix B, Subramanian S: Individual and collective bodies: using measures of variance and association in contextual epidemiology. Journal of epidemiology and community health. 2009,
Larsson A, Palmer M, Hulten G, Tryding N: Large differences in laboratory utilisation between hospitals in Sweden. Clin Chem Lab Med. 2000, 38: 383-389. 10.1515/CCLM.2000.056.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2296/11/43/prepub
This research was funded by Health and Medical Care Executive Board of the Region Västra Götaland, the Swedish Research Council and ALF Government Research Grant. We are most grateful to Birgitta Lindberg for excellent technical assistance in preparation of the manuscript. This paper was presented in part at the Nordic Congresses of General Practice in June 2007.
The authors declare that they have no competing interests.
SD conceived the study, drafted the manuscript, responded to the reviewer comments and critically revised the manuscript. PH conceived the study, participated on the design of the study, performed multilevel analyses and interpretation of data responded to the reviewer comments and critically revised the manuscript. HO participated on the design of the study, supported PH in the performance of the multilevel analyses and interpretation of data. RE conceived the study and critically revised the manuscript for important intellectual content. JM participated in the design of the study and interpretation of data and critically revised the manuscript for important intellectual content. KBB Conceived the study, drafted the manuscript, responded to the reviewer comments and critically revised the manuscript. All authors read and approved the final manuscript.
Sofia Dalemo, Per Hjerpe contributed equally to this work.