Variation in plasma calcium analysis in primary care in Sweden - a multilevel analysis
- Sofia Dalemo†1, 2Email author,
- Per Hjerpe†2, 3,
- Henrik Ohlsson3,
- Robert Eggertsen1,
- Juan Merlo3 and
- Kristina Bengtsson Boström2
© Dalemo et al; licensee BioMed Central Ltd. 2010
Received: 9 November 2009
Accepted: 30 May 2010
Published: 30 May 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
Examples of diagnoses with strong influence of the chance of having a plasma calcium analyses the risk score equation. Total number of patients 154 629.
All the positive diagnoses in the stepwise regression
Number of P-Ca analysis
Other disorders of thyroid
Nonorganic sleeping disorders
Diseases of the circulatory system
Atrial fibrillation and flutter
Diseases of the digestive system
Diseases of the musculoskeletal system
Other artritis and rheumatism unspecified
Abnormal blood-pressure reading, without diagnosis
Malaise and fatigue
Contact with health services
General medical examination
Calcium and vitamin D supplements
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 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 .
Staffing of physicians at health care centres, characteristics of physicians and number of patients visits and frequency of plasma calcium analyses per title in the county of Skaraborg during 2005.
Mean age (year)
Median age (year)
Range age (year)
With P-Ca test %
Preregistration house officer
GP < 46 year
GP ≥ 46 year
Locum < 46 year
Locum ≥ 46 year
The multilevel model
Multi-level logistic regression analysis of plasma calcium analyses in primary care in the county of Skaraborg, Sweden
P-Ca test 2004
Preregistration house officer
GP < 46 year
GP ≥ 46 year
Locum < 46 year
Locum ≥ 46 year
Number groups include P-Ca
HCC and Physician (Intercept)
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.
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.
- 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: SPRIGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Goldstein H: Multilevel Statistical Models. 2003, London: Hodder ArnoldGoogle Scholar
- Snijders T, Bosker R: Multilevel Analysis An Intro- duction to Basic and Advanced Multilevel Modeling. 1999, London: SageGoogle Scholar
- Fielding A, Goldstein H: Cross-classified and Multiple Membership Structures in Multilevel Models: An Introduction and Review. 2006, Birmingham: University of BirminghamGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Gilks WR, Richardson S, Spiegelhalter DJ: Markov Chain Monte Carlo in Pratice. 1996, London, Chapman and HallGoogle Scholar
- 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 LondonGoogle Scholar
- 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-Google Scholar
- 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.PubMedPubMed CentralGoogle Scholar
- Salloum S, Franssen E: Laboratory investigations in general practice. Can Fam Physician. 1993, 39: 1055-1061.PubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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,Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2296/11/43/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.