2,619
Views
18
CrossRef citations to date
0
Altmetric
Articles

Reduced survival of total knee arthroplasty after previous unicompartmental knee arthroplasty compared with previous high tibial osteotomy: a propensity-score weighted mid-term cohort study based on 2,133 observations from the Danish Knee Arthroplasty Registry

, , &

Abstract

Background and purpose — Both medial unicompartmental knee arthroplasties (UKA) and high tibial osteotomies (HTO) are reliable treatments for isolated medial knee osteoarthritis. However, both may with time need conversion to a total knee arthroplasty (TKA). We conducted the largest nationwide registry comparison of the survival of TKA following UKA with TKA following HTO.

Patients and methods — From the Danish Knee Arthroplasty Registry, aseptic conversions to TKA from UKA and TKA converted from HTO within the period of 1997–2018 were retrieved. The Kaplan–Meier method and the Cox proportional hazards regression were used to estimate the survival and hazard ratio (HR) for revision, considering confounding by indication utilizing propensity-score based inverse probability of treatment weighting (PS-IPTW).

Results — PS-IPTW yielded a well-balanced pseudo-cohort (standard mean difference (SMD) < 0.1 for all covariates, except implant supplementation) of 963.8 TKAs following UKA and 1139.1 TKAs following HTO. The survival of TKA following UKA was significantly less than that of TKA following HTO with a 5-year estimated survival of 0.88 (95% confidence interval (CI) 0.85–0.90) and 0.94 (CI 0.93–0.96), respectively. The differences in survival corresponded to an implant-supplementation adjusted HR of 2.7 (CI 2.4–3.1) for TKA following UKA compared with TKA following HTO.

Interpretation — Previous UKA more than doubled the revision risk of a subsequent TKA compared with previous HTO. This potential risk should be considered in the shared treatment decision of patients who are candidates for both UKA and HTO.

In isolated osteoarthritis of the medial knee compartment, both medial unicompartmental knee arthroplasties (UKA) and high tibial osteotomies (HTO) are solutions with reliable clinical outcomes (Cao et al. Citation2018). The survival of UKA is secondary to that of total knee arthroplasties (TKA) with a recent meta-analysis reporting 15-year survival of 76% and 93%, respectively (Evans et al. Citation2019). The long-term survival of HTO seems inferior to both UKA and TKA with a declining survival from 75% at 10 years to 55% at 15 years (van Wulfften Palthe et al. Citation2018). When UKA or HTO fail, conversion to TKA is a common solution (Lee et al. Citation2019). Nationwide registry studies have investigated the survival of either TKA following UKA or TKA following HTO compared with primary or revision TKA. They have reported an increased risk of revision in TKA following UKA (Robertsson and W-Dahl Citation2015, Leta et al. Citation2016, Lewis et al. Citation2018, El-Galaly et al. Citation2019) while no consensus regarding the influence of HTO on the survival of a subsequent TKA has been reached (Niinimäki et al. Citation2014, Badawy et al. Citation2015, Robertsson and W-Dahl Citation2015, El-Galaly et al. Citation2018). However, a direct comparison of the survival estimates from these studies is prone to confounding by indication due to a range of unadjusted baseline characteristics associated with the survival of TKA, such as implant constraints and hospital volume of arthroplasty surgeries (Jasper et al. Citation2016). This concern is further encouraged by a recent single-center study reporting similar short-term survival of TKA following UKA and TKA following HTO (Lim et al. Citation2017). Based on the Danish Knee Arthroplasty Registry (DKR), our study compares the survival of TKA converted from UKA with TKA converted from HTO with consideration for confounding by indication utilizing propensity-score based inverse probability of treatment weighting (PS-IPTW) (Inacio et al. Citation2015).  

Patients and methods

Data source

Since 1997, the DKR has prospectively collected information on Danish knee arthroplasties through standardized forms completed by surgeons. Since 2007, the registration of arthroplasties has been mandatory for all hospitals in Denmark leading to a registry completeness above 90% for primary arthroplasties and 80% for revision arthroplasties (Danish Knee Arthroplasty Registry Citation2019). The DKR is reported suitable for epidemiological studies and is crosslinked with the Danish Civil Registration System (DCRS) which contains vital and emigration status for all Danish citizens (Pedersen et al. Citation2012, Schmidt et al. Citation2014). Mandatory registration and linkage to the DCRS enable complete follow-up in a population-based cohort (Schmidt et al. Citation2019).

Study cohort

In the registry, each patient is identified by a unique code, and the side of surgery is denoted. Therefore, each knee can be considered a unique observation, which has been reported to provide unbiased results in large arthroplasty studies (Robertsson and Ranstam Citation2003). All UKAs indicated by osteoarthritis in knees without prior surgery from January 1, 1997 until December 31, 2018 were retrieved. We identified later revisions conducted on the same knee and excluded revisions to UKA and conversions from UKA to TKA due to infections (). In the same timeframe, we retrieved all TKAs indicated by osteoarthritis in knees previously treated with HTO. The validity of the registration of a previous HTO was evaluated in a recent study and confirmed in 96% of the cases (El-Galaly et al. Citation2018). However, HTO is not divided in open- and closed-wedge osteotomies, and thus the methods are considered as one in the DKR.

Figure 1. Flowchart depicting the formation of the study cohort.

Figure 1. Flowchart depicting the formation of the study cohort.

Covariates

At the time of TKA surgery, several patients and surgical characteristics are registered in the DKR (). Comorbidity is depicted by Charnley class sorted into class A (unilateral arthritis), class B1 (bilateral arthritis), class B2 (arthroplasty in the opposite knee), and class C (other condition affecting walking capacity) (Bjorgul et al. Citation2010). Knee function at the time of surgery is registered using American Knee Society Score (KKS) divided into clinical and functional sub-scores, both ranging from 0 to 100 (Insall et al. Citation1989). Level of implant constraint is divided into cruciate retaining (CR), posterior stabilized (PS), constraint condylar (CCR), and hinged. Perioperative supplementation (stems, augments, or cones) and perioperative complications (e.g. fractures, rupture of the patellar tendon or ligament injuries) are registered. We defined the hospital volume of arthroplasties as the mean annual volume during the study period and divided it into 4 groups (< 100, 100–249, 250–449, > 449). Year of surgery was classified into 2 periods (1997–2007 and 2008–2018). For TKAs following UKAs, we retrieved type of UKA bearings and indications of conversion to TKA.

Table 1. Baseline characteristics for original cohort and the PS-IPTW cohort at the time of conversion to TKA

Outcome

The outcome was TKA revision of any indication with revision defined in accordance with the DKR as removal, exchange, or addition of an implant. The indications for TKA revision have recently been thoroughly evaluated in both groups, and thus are not presented in this study (El-Galaly et al. Citation2018, Citation2019).

Missing values

Missing values existed in height (n = 1,214), weight (n = 67), KSS clinical sub-score (n = 55), KSS functional sub-score (n = 39), Charnley class (n = 12), fixation (n = 11), duration of surgery (n = 8), and patella replacement (n = 3). Missing values in height were deemed too high for meaningful imputation and discarded. The remaining missing values were estimated by multiple imputation with chain equation (MICE), generating 5 datasets under the assumption of missing at random (Azur et al. Citation2011).

Statistics

PS-IPTW

This study is subjected to confounding due to the non-random assignment of prior UKA or HTO. Therefore, PS-IPTW was utilized to account for confounding by indication. PS were estimated with logistic regression and applied by IPTW with stabilized weights aiming to estimate the average effect of treatment (Austin Citation2014). Based on the considerations depicted in the directed acyclic graph (Williams et al. Citation2018) (Figure 2, see Supplementary data), the following covariates were included in the model: sex, age (quantiles), weight (quantiles), KSSs, Charnley class, level of constraint, patella resurfacing, fixation, hospitals annual arthroplasty volume, and period of surgery. Implant supplementation was rare in TKA following HTO and therefore omitted from the PS estimation to avoid overweighting rare observations. The balance of the baseline characteristics was evaluated graphically and by standardized mean differences (SMD) with an SMD of 0 indicating perfect balance and SMD < 0.1 deemed an acceptable balance between the groups (Austin Citation2009).

Survival analyses

The Kaplan–Meier method was used to estimate the survival with revision as the primary endpoint. Unrevised knees were censored by death, emigration, or end of study period at December 31, 2018. The risk of revision was estimated by Cox regression with robust variance estimator to account for dependencies in the PS-IPTW cohort. The assumption of proportional hazards was evaluated by Schoenfeld’s plots and Schoenfeld’s residual test. Implant supplementation was included as covariate in the Cox regression to account for remaining imbalance following the PS-IPTW.

E-value

The robustness of the estimated hazard ratios (HR) was evaluated by calculating their E-values, which estimates the magnitude of association unmeasured confounders must have with both the exposure and outcome to negate the observed HRs (Van Der Weele and Ding Citation2017).

Significance

Means are presented with absolute range, medians with interquartile range (IQR), and SMD are calculated to assess balance between the groups. Estimates from the imputed datasets were combined by Rubin’s rule (White et al. Citation2011), and all estimates are presented with 95% confidence interval (CI) to address their significance (Ranstam Citation2019).

Statistical programs

Data were sorted in STATA 15 (StataCorp, College Station, TX, USA) and all analyses were conducted in R© Version 3.5.1 (R Foundation for Statistical Computing, Vienna, Austria).

Ethics, funding, and potential conflict of interests

The study was approved by the Danish Data Protection Agency (entry number: 2008-58-0028) and financed by the Orthopaedic Research Unit at Aalborg University Hospital. No conflict of interest is present among the authors.

Results

Original cohort

Characteristics

The baseline covariates differed in sex, age, Charnley class, KKSs, level of constraint, patella resurfacing, fixation, implant supplementation, annual arthroplasty volume, and period of surgery as depicted in . The median duration of surgery in TKA following UKA was 90 minutes (IQR 75–115) compared with 80 minutes (IQR 65–100) in TKA following HTO (SMD = 0.35). In 12 of the TKAs following UKAs perioperative complications (8 fractures, 4 other) were registered compared with 18 registered complications (5 fractures, 5 ligament/tendon rupture, 8 other) in TKAs following HTOs (SMD = 0.11).

Survival

Of the 978 TKAs following UKAs, 121 (12%) were revised within the study period and 93 (10%) were censored due to either death or emigration. In comparison, 101 (9%) and 234 (20%) of the 1,155 TKAs following HTOs were revised or censored, respectively. The median follow-up in TKA following UKA was 4.7 years (IQR 1.9–7.7) compared with 9.3 years (IQR 5.0–13) for TKA following HTO (SMD = 0.87). The 1st, 5th, and 10th year survival estimates were 0.97 (CI 0.96–0.98), 0.88 (CI 0.86–0.91), and 0. 82 (CI 0.78–0.85) for TKA following UKA compared with 0.98 (CI 0.97–0.99), 0.95 (CI 0.93–0.96), and 0.92 (CI 0.90–0.94) for TKA following HTO, which corresponds to an HR of 2.3 (CI 2.1–2.6) associated with TKA following UKA ().

Table 3. Survival estimates, hazard ratios (HR), and E-value for the original cohort and PS-IPTW cohort

PS-IPTW cohort

Characteristics

Following PS-IPTW, all covariates included in the estimation of the PS were well balanced between TKAs following UKAs and TKAs following HTOs ( and Figure 3, see Supplementary data). However, the difference in implant supplementation was still unbalanced following PS-IPTW (SMD = 0.58). depicts the distribution of indication of UKA conversion and type of UKA-bearing in TKA following UKA, which was clinically comparable before and after PS-IPTW. The imbalance in duration of surgery was unchanged by PS-IPTW with a median duration of 90 minutes (IQR 75–120) in TKA following UKA and 80 minutes in TKA following HTO (IQR 66–100) (SMD = 0.35). PS-IPTW did not balance the difference in perioperative complications with registered complications in 7.7 of the TKAs following UKAs and 24.1 in TKAs following HTOs (SMD = 0.13).

Table 2. Baseline characteristics for original cohort and PS-IPTW cohort at the time of conversion to TKA. Values are counts/weighted counts (%)

Survival

Of the 963.8 TKAs following UKAs, 169.1 (17%) were revised within the study period and 117.5 (12%) were censored due to death or emigration. Similarly, 89.4 (8%) and 187.8 (16%) of the 1139.1 TKAs following HTOs were revised or censored, respectively. This corresponded to a significantly inferior survival of TKA following UKA compared with TKA following HTO ( and ), with an implant-supplementation adjusted HR of 2.7 (CI 2.4–3.1) associated with TKA following UKA (). 

Figure 4. Kaplan–Meier survival estimates for the PS-IPTW cohort with confidence interval and weighted knees at risk.

Figure 4. Kaplan–Meier survival estimates for the PS-IPTW cohort with confidence interval and weighted knees at risk.

Discussion

Our study showed that in a cohort from the Danish Knee Arthroplasty Registry, with well-balanced baseline covariates, the survival of TKA following UKA was lower than the survival of TKA following HTO. More specifically, the risk of revision more than doubled when TKA was preceded by UKA compared with HTO.

During recent years, UKA has gained popularity while the use of HTO has decreased, indicating a trend towards treating patients with UKA instead of HTO (Niinimäki et al. Citation2012, Henkel et al. Citation2019). While both UKAs and HTOs relieve pain, the procedures differ. UKAs replace the diseased compartment with an implant and thus preserve the mechanical axis of the knee. In contrast, HTOs shift the mechanical axis laterally, unloading the diseased compartment while increasing the load on the lateral compartment. Therefore, a subsequent conversion from HTO to TKA due to progression of osteoarthritis might be expected, whereas progression of arthritis or implant failure are considered an adverse event in UKA surgery. To our knowledge, this study constitutes the largest direct comparison of TKA following UKA with TKA following HTO based on a nationwide registry (Pearse et al. Citation2012) and expands upon a range of recent nationwide registry studies comparing the survival of either TKA following UKA or TKA following HTO with primary TKA (Niinimäki et al. Citation2014, Badawy et al. Citation2015, Robertsson and W-Dahl Citation2015, Leta et al. Citation2016, El-Galaly et al. Citation2018, Citation2019, Lewis et al. Citation2018). Due to inconsistent adjustment for confounding, a direct comparison of the results in the current literature might be affected by residual confounding, which could result in acceptance of a false causal relationship (Kyriacou and Lewis Citation2016). Our study compliments current literature by directly comparing TKA following UKA with TKA following HTO while expanding the statistical adjustment for confounding using PS-IPTW. The PS is the probability of an observation receiving a treatment given a set of baseline covariates and, thus, dependence on the PS creates balance in the included covariates between the groups. Dependence on the PS can be achieved by matching, weighting, adjusting, or stratifying (Austin Citation2011). We used IPTW and, thus, weighted the observation based on their inverse probability of treatment (i.e., PS) to create a pseudo-cohort with comparable baseline characteristics between the groups. As depicted in Figure 3 (see Supplementary data), this approach eliminated imbalances in a range of baseline covariates, and thus diminished the influence of the confounders presented in Figure 2 (see Supplementary data) except implant supplementation, which was including in the Cox regression. In this pseudo-cohort, TKA following UKA was associated with a 2.7-fold increase in the risk of revision compared with TKA following HTO.

Limitations

The study has some limitations. First, nationwide registries are prone to misclassifications. However, as the data are collected prospectively by the surgeon on a standardized form, the misclassifications are assumed to be non-differential and thus bias the results towards no difference between the groups. Second, even though the PS-IPTW successfully balances a wide range of covariates, residual confounding is unavoidable in non-randomized studies. We calculated the E-value for the presented HRs to elucidate which magnitude unmeasured confounders must have to negate the presented HRs (). The E-value indicated that unmeasured confounders must be associated with both TKA following UKA (exposure) and subsequent revisions (outcome) by a ratio of at least 4.3 (lower CI) to move the HR’s CI to include 1. In comparison, diabetes has recently been associated with a risk ratio of revision at 1.3 (CI 1.02–1.6) in a large retrospective study of both TKAs and total hip arthroplasties (Maradit Kremers et al. Citation2017). Therefore, the presented HRs seemed robust for residual confounding. Third, the completeness in the DKR has increased from 69% of primary arthroplasties in 1997 to above 91% since 2008, with a similar evolution in revision arthroplasties with a completeness from 54% in 1997 to above 87% since 2008 (Danish Knee Arthroplasty Registry Citation2019). The overall completeness of TKA following HTO might be less than the overall completeness of TKA following UKA, as more HTOs were converted before 2008. This imbalance might have overestimated the risk of revision associated with TKA following UKA compared with TKA following HTO. We included the period of surgery in the PS estimation to contain the bias induced by the difference in completeness.

Conclusion

In this propensity-score weighted cohort study, TKA following UKA was associated with a more than 2-fold increased risk of revision compared with TKA following HTO. This potential risk emphasized that UKA should be considered a definitive treatment in line with TKA rather than a temporary treatment to postpone TKA.

Supplementary data

Figures 2 and 3 are available as supplementary data in the online version of this article, http://dx.doi.org/10.1080/17453674.2019.1709711

The authors thank the Danish knee surgeons for a thorough registration of their procedures, the steering committee of the Danish Knee Arthroplasty Registry for their goodwill in relation to data acquisition, and Angélica Meleñdez-Muñoz for her linguistic contribution to the manuscript.

All authors contributed to the study design. AEG received and analyzed the data with supervision from all authors. AEG wrote the initial draft, which was revised and accepted by all authors.

Acta thanks Leif Ryd and Annette W-Dahl for help with peer review of this study.

Supplemental material

Supplemental Material

Download PDF (270.5 KB)

  • Austin P C. Balance diagnostics for comparing the distribution of baseline covariates between treatment groups in propensity-score matched samples. Stat Med 2009; 28(25): 3083–107.
  • Austin P C. An introduction to propensity score methods for reducing the effects of confounding in observational studies. Multivariate Behav Res 2011; 46(3): 399–424.
  • Austin P C. The use of propensity score methods with survival or time-to-event outcomes: reporting measures of effect similar to those used in randomized experiments. Stat Med 2014; 33(7): 1242–58.
  • Azur M J, Stuart E A, Frangakis C, Leaf P J. Multiple imputation by chained equations: what is it and how does it work? Int J Methods Psychiatr Res 2011; 20(1): 40–9.
  • Badawy M, Fenstad A M, Indrekvam K, Havelin L I, Furnes O. The risk of revision in total knee arthroplasty is not affected by previous high tibial osteotomy: a 15-year follow-up of 32,476 total knee arthroplasties in the Norwegian Arthroplasty Register 2015; 86(6): 734–9.
  • Bjorgul K, Novicoff W M, Saleh K J. Evaluating comorbidities in total hip and knee arthroplasty: available instruments. J Orthop Traumatol 2010; 11(4): 203–9.
  • Cao Z, Mai X, Wang J, Feng E, Huang Y. Unicompartmental knee arthroplasty vs high tibial osteotomy for knee osteoarthritis: a systematic review and meta-analysis. J Arthroplasty 2018; 33: 952–9.
  • Danish Knee Arthroplasty Registry. Danish Knee Arthroplasty Registry—Annual Report 2019; 2019.
  • El-Galaly A, Nielsen P T, Jensen S L, Kappel A. Prior high tibial osteotomy does not affect the survival of total knee arthroplasties: results from the Danish Knee Arthroplasty Registry. J Arthroplasty 2018; 33(7): 2131–5.e1.
  • El-Galaly A, Kappel A, Nielsen P T, Jensen S L. Revision risk for total knee arthroplasty converted from medial unicompartmental knee arthroplasty. Comparison with primary and revision arthroplasties, Based on mid-term results from the Danish Knee Arthroplasty Registry. J Bone Joint Surg 2019; 101(22): 1999–2006.
  • Evans J T, Walker R W, Evans J P, Blom A W, Sayers A, Whitehouse M R. How long does a knee replacement last? A systematic review and meta-analysis of case series and national registry reports with more than 15 years of follow-up. Lancet 2019; 393(10172): 655–63.
  • Henkel C, Mikkelsen M, Pedersen A B, Rasmussen L E, Gromov K, Price A, Troelsen A. Medial unicompartmental knee arthroplasty: increasingly uniform patient demographics despite differences in surgical volume and usage—a descriptive study of 8,501 cases from the Danish Knee Arthroplasty Registry. Acta Orthop 2019; 90(4): 354–9 .
  • Inacio M C S, Chen Y, Paxton E W, Namba R S, Kurtz S M, Cafri G. Statistics in brief: an introduction to the use of propensity scores. Clin Orthop Relat Res 2015; 473(8): 2722–6.
  • Insall J N, Dorr L D, Scott R D, Scott W N. Rationale of the Knee Society clinical rating system. Clin Orthop Relat Res 1989; (248): 13–14.
  • Jasper L L, Jones C A, Mollins J, Pohar S L, Beaupre L A. Risk factors for revision of total knee arthroplasty: a scoping review. BMC Musculoskelet Disord 2016; 17(1): 1–9.
  • Kyriacou D N, Lewis R J. Confounding by indication in clinical research. JAMA 2016; 316(17): 1818–19.
  • Lee Y, Kim H, Mok S, Lee O-S. Similar outcome, but different surgical requirement in conversion total knee arthroplasty following high tibial osteotomy and unicompartmental knee arthroplasty: a meta-analysis. J Knee Surg 2019; 32(7): 686–700.
  • Leta T H, Lygre S H L, Skredderstuen A, Hallan G, Gjertsen J E, Rokne B, Furnes O. Outcomes of unicompartmental knee arthroplasty after aseptic revision to total knee arthroplasty: a comparative study of 768 TKAs and 578 UKAs revised to TKAs from the Norwegian Arthroplasty Register (1994 to 2011). J Bone Joint Surg Am 2016; 98(6): 431–40.
  • Lewis P L, Davidson D C, Graves S E, De Steiger R N, Donnelly W, Cuthbert A. Unicompartmental knee arthroplasty revision to TKA: are tibial stems and augments associated with improved survivorship? Clin Orthop Relat Res 2018; 476(4): 854–62.
  • Lim J B T, Chong H C, Pang H N, Tay K J D, Chia S L, Lo N N, Yeo SJ, Lim J B T, Joint B. Revision total knee arthroplasty for failed high tibial osteotomy and unicompartmental knee arthroplasty have similar patient-reported outcome measures in a two-year follow-up study. Bone Joint J 2017; 99-B(10): 1329–34 .
  • Maradit Kremers H, Schleck C D, Lewallen E A, Larson D R, Van Wijnen A J, Lewallen D G. Diabetes mellitus and hyperglycemia and the risk of aseptic loosening in total joint arthroplasty. J Arthroplasty 2017; 32(9): S251–S3.
  • Niinimäki T T, Eskelinen A, Ohtonen P, Junnila M, Leppilahti J. Incidence of osteotomies around the knee for the treatment of knee osteoarthritis: a 22-year population-based study. Int Orthop 2012; 36(7): 1399–402.
  • Niinimäki T, Eskelinen A, Ohtonen P, Puhto A-P, Mann B S, Leppilahti J. Total knee arthroplasty after high tibial osteotomy: a registry-based case-control study of 1,036 knees. Arch Orthop Trauma Surg 2014; 134(1): 73–7.
  • Pearse A J, Hooper G J, Rothwell A G, Frampton C. Osteotomy and unicompartmental knee arthroplasty converted to total knee arthroplasty: data from the New Zealand Joint Registry. J Arthroplasty 2012; 27(10): 1827–31.
  • Pedersen A B, Mehnert F, Odgaard A, Schroder H M. Existing data sources for clinical epidemiology: the Danish Knee Arthroplasty Register. Clin Epidemiol 2012; 4:125–35.
  • Ranstam J. Time to restrict the use of p-values in Acta Orthopaedica. Acta Orthop 2019; 90(1): 1–2.
  • Robertsson O, Ranstam J. No bias of ignored bilaterality when analysing the revision risk of knee prostheses: analysis of a population based sample of 44,590 patients with 55,298 knee prostheses from the national Swedish Knee Arthroplasty Register. BMC Musculoskelet Disord 2003; 4: 1.
  • Robertsson O, W-Dahl A. The risk of revision After TKA is affected by previous HTO or UKA. Clin Orthop Relat Res 2015; 473(1): 90–3.
  • Schmidt M, Pedersen L, Sorensen H T. The Danish Civil Registration System as a tool in epidemiology. Eur J Epidemiol 2014; 29(8): 541–9.
  • Schmidt M, Schmidt S A J, Adelborg K, Sundbøll J, Laugesen K, Ehrenstein V, Sørensen HT. The Danish health care system and epidemiological research: from health care contacts to database records. Clin Epidemiol 2019; 11: 563–91.
  • Van Der Weele T J, Ding P. Sensitivity analysis in observational research: introducing the E-value. Ann Intern Med 2017; 167(4): 268–74.
  • White I R, Royston P, Wood A M. Multiple imputation using chained equations: issues and guidance for practice. Stat Med 2011; 30(4): 377–99.
  • Williams T C, Bach C C, Matthiesen N B, Henriksen T B, Gagliardi L. Directed acyclic graphs: a tool for causal studies in paediatrics. Pediatr Res 2018; 84(4): 487–93.
  • van Wulfften Palthe A F Y, Clement N D, Temmerman O P P, Burger B J. Survival and functional outcome of high tibial osteotomy for medial knee osteoarthritis: a 10–20-year cohort study. Eur J Orthop Surg Traumatol 2018; 28(7): 1381–9.