 Research article
 Open Access
 Published:
Network metaanalysis of (individual patient) time to event data alongside (aggregate) count data
BMC Medical Research Methodology volume 14, Article number: 105 (2014)
Abstract
Background
Network metaanalysis methods extend the standard pairwise framework to allow simultaneous comparison of multiple interventions in a single statistical model. Despite published work on network metaanalysis mainly focussing on the synthesis of aggregate data, methods have been developed that allow the use of individual patientlevel data specifically when outcomes are dichotomous or continuous. This paper focuses on the synthesis of individual patientlevel and summary time to event data, motivated by a real data example looking at the effectiveness of high compression treatments on the healing of venous leg ulcers.
Methods
This paper introduces a novel network metaanalysis modelling approach that allows individual patientlevel (time to event with censoring) and summarylevel data (event count for a given followup time) to be synthesised jointly by assuming an underlying, common, distribution of time to healing. Alternative model assumptions were tested within the motivating example. Model fit and adequacy measures were used to compare and select models.
Results
Due to the availability of individual patientlevel data in our example we were able to use a Weibull distribution to describe time to healing; otherwise, we would have been limited to specifying a uniparametric distribution. Absolute effectiveness estimates were more sensitive than relative effectiveness estimates to a range of alternative specifications for the model.
Conclusions
The synthesis of time to event data considering individual patientlevel data provides modelling flexibility, and can be particularly important when absolute effectiveness estimates, and not just relative effect estimates, are of interest.
Background
Treatment decisions in medicine, whether at the patient or policy level, should consider all relevant alternative health care technologies potentially capable of delivering the change being sought. Such informed decision making on the use of competing treatments requires evidence of relative effects from randomised controlled trials (RCTs) (which may be further used to inform estimates of costeffectiveness). However, trying to synthesise evidence from several different pairwise comparisons relevant to a decision problem can be challenging. A piecemeal approach can be avoided if all RCTs evaluating interventions relevant to the treatment decision are considered collectively, for example, with the use of network metaanalysis (NMA). NMA is a wellestablished statistical technique that extends standard pairwise metaanalysis framework to allow simultaneous comparison of multiple interventions in a single statistical model [1, 2]. This approach then produces relative effect estimates (and associated descriptions of uncertainty) for all treatments connected by the network of evidence [3] – even where headtohead trials for comparisons do not exist.
NMA using individual patient data
Published work on NMA mainly focuses on the synthesis of aggregate data (AD) (sometimes called summary data, e.g. group means and standard errors available from study reports) [4, 5]; however, methods have been developed that allow use of individual patientlevel data (IPD) in NMA, specifically for binary and continuous data [6–8]. The appeal of including IPD in a NMA is that it is likely to reduce statistical heterogeneity across the network (and in this way help resolve possible inconsistencies); it may also allow subgroup effects to be estimated which in turn could guide more personalised treatment decisions [6]. The use of IPD, alone or in combination with AD, has been shown to improve inference in NMAs where the outcome of interest is binary (dichotomous) by aiding convergence, and by providing unbiased treatment–covariate interactions (that would otherwise be affected by ecological bias [9]). For continuous outcomes, IPD is likely to produce more precise estimates of treatment effects, even in the absence of treatment–covariate interactions [10].
NMA using time to event related outcome data
When undertaking an NMA, aggregate time to event data presented as hazard ratios (HR), and some measure of uncertainty around these estimates, can be pooled directly using standard methods (analogous to pooling odds ratios or relative risks) [11]. However, other AD outputs such as median/mean time to event [12] and cumulative counts of patients having the outcome event in a period of time [13] have also been metaanalysed in a network. In modelling these data, an underlying time to event distribution can be specified which will then allow HRs to be generated from the original AD [5, 14, 15]. To date, evidence synthesis that includes IPD on time to event has been limited to pairwise analysis; additionally there has been little methodological exploration how outcomes in the form of IPD and AD might be considered together in the same NMA [13].
Developing a NMA combining AD and IPD data to synthesise time to event related outcomes
This paper describes a modelling framework that combines AD and IPD in the synthesis of time to event related outcomes within an NMA. It extends the work of Saramago et al. [6] and Sutton et al. [16] on the synthesis of AD and IPD for binary data to consider time to event outcomes. And extends the work on NMAs for time to event AD outcomes by Soares et al. [15] and Woods et al. [14] (the general modelling framework is also described in Dias et al. [5]). Our work was motivated by an NMA on treatments for healing venous leg ulcers, for which we had data from multiple RCTs (more information provide in the next section). A proportion of the RCT data on ulcer healing was available in IPD format (time to healing and time to censoring), with the remaining data available as AD, i.e. count data (proportion of patients healed at end of followup which differs across studies). To maximally draw from available data we aimed to jointly synthesise the available AD and IPD, and, in this way, generating better estimates and providing fuller characterisations of uncertainty to best inform decisions on the use of the treatments of interest.
Motivating example: high compression treatments for venous leg ulcers
The case study relates to compression systems aiming to deliver high compression (classed as ≥40 mmHg compression at the ankle) to promote venous leg ulcer healing. Available standardised systems are: two layer hosiery (HH), the four layer bandage (4LB), the short stretch bandage (SSB), the zinc paste bandage (ZINC), and the two layer bandage system (2LB). A detailed description is provided in Additional file 1 with further details of these systems presented elsewhere [17].
Effectiveness evidence from RCTs was obtained from the most recent update of the relevant Cochrane review available to us, and from a recent multicentre RCT which compared 4LB with HH. All available RCT evidence was assessed for inclusion in the current NMA: a detailed process that has been reported elsewhere [17]. The final NMA contained data from 16 RCTs on the relative effectiveness of high compression systems for the treatment of venous leg ulcers. Data for two of the 16 included RCTs (VenUS I and VenUS IV, hereby denominated studies 1 and 2) had full IPD data available (841 participants) which included time to healing or censoring for each participant, together with other individuallevel characteristics such as treatment centre, ulcer duration and size and also patient mobility. For the remaining RCTs (1105 participants), aggregate data on the number of healed ulcers were extracted from the source review alongside information regarding treatment type, number of participants allocated to each treatment group, mean duration of followup (if this was not stated, trial duration was used), mean ulcer duration and size.
The 16 included RCTs described nine unique high compression treatments: the five standard treatments (4LB, SSB, ZINC, HH and the 2LB) and four ad hoc systems [17]. The ad hoc group consisted of treatments deemed irrelevant to current clinical practice, and are not reported further (results can be provided upon request). These studies were, however, included in the NMA as their data was potentially relevant, for example, in describing determinants of healing.
Table 1 describes the data available and Figure 1 presents the treatment network formed by the evidence [17–32]. The most populated comparison was the 4LB vs. SSB, informed by seven RCTs: six with healing data available as AD [18–22, 31] and one as IPD [32]. The link between the 2LB and 4LB was informed by two RCTs and each of the remaining six comparisons in the NMA were informed by AD extracted from one RCT for each comparison (Table 1).
Methods
We first describe the modelling framework for our main analysis, model A. This model brings together IPD and AD into the same synthesis model, with both types of evidence contributing to the estimation of all key models parameters (i.e. absolute and relative effectiveness estimates). We then detail the evaluation of alternative assumptions applied to model A, thus highlighting and challenging specific assumptions of the modelling framework proposed. All synthesis was conducted in a Bayesian framework.
Suppose there are two sets of studies, one set for which IPD are available and one set for which AD are available. Consider the set of studies to be T = {1,…, number of IPD studies (NS_{IPD}),…, NS_{IPD} + number of AD studies (NS_{AD})} where the total number of studies is J; consider the set of treatments to be {1,2,3,…}, where the total number of treatments is K, and treatment 1 the reference.
Statistical model for the data
We describe model A in two interrelated parts: part I describes the modelling of the IPD and part II the modelling of the AD.
Time to ulcer healing (t _{ ijk }) of the i ^{th} participant in the j ^{th} study (where j =1,…, NS _{ IPD }) and in the k ^{th} treatment arm was assumed to be Weibull distributed [33] with shape^{a.} parameter, s, and scale, λ _{ ijk }. For some participants, time to event was not observed, and these observations were censored at the time the participant last had trial data recorded, {t}_{\mathit{\text{ijk}}}^{c}. The linear predictor log(λ _{ ijk }), was modelled as a function of the loghazard of an event for the baseline treatment b, {\mu}_{j}^{\mathit{\text{IPD}}}, of the log of the treatment effects, d _{ bk } and of a studyspecific individuallevel regression term, β _{0j } x _{ ijk }, where β _{0j } exemplifies a covariate effect, i.e. the difference in the log hazard ratio per unit increase in the covariate x _{ ijk } a patientlevel covariate for the i ^{th} patient in the j ^{th} trial on treatment k available in the IPD data sets [17]). The effect of each covariate on the hazard of healing was assumed to be equal in both IPD studies. Due to the possibility of missing covariate information existing for some individuals, x _{ ijk } was represented as a Normally distributed random variable with mean m and precision p, common across all IPD studies. This is essentially a multiple imputation technique through MCMC^{b.} and assumes that the missingness mechanism was at random^{c.}. Additionally, to account for centre variability within each IPD study, {\gamma}_{j}^{c} was defined for each centre, c, in the j ^{th} study, these were combined using a common frailty effect described by a normal distribution with mean zero and precision τ.
The treatment effects, d _{ bk }, were loghazard ratios for treatment k relative to the studyspecific baseline treatment b, partitioned further as d _{1k }  d _{1b }. This equation defines the set of functional parameters, that is the log effects of each treatment in relation to the reference treatment 1. Note that d _{11} = 0, where treatment 1 was 4LB, arbitrarily chosen as reference treatment.
Vague prior distributions were specified for \left[{\mu}_{b}^{\mathit{\text{IPD}}}\left(\sim N\left(0,{10}^{6}\right)\right)\right] for the regression coefficient β _{0j } (∼N(0, 10^{6})) for m(∼N(0, 10^{6})) and p(∼Gamma(0.01, 0.01)) for the shared centre effect τ(∼Gamma(0.01, 0.01)), for s(∼Gamma(0.01, 0.01)) and for d _{1k } ∼ N(0, 10^{6}).
Within the AD studies, the observed number of participants with a healed ulcer within each study, r _{ jk }, out of the total number of individuals in the j ^{th} trial [(where j = NS _{ IPD } +1,…, (NS _{ IPD } + NS _{ AD })] and k ^{th} treatment (intention to treat), n _{ jk }, was assumed to be Binomially distributed. The underlying probabilities of an event for each arm and in each trial were represented by p _{ jk }. In turn, p _{ jk } was expressed as a function of the scale parameter, {\lambda}_{\mathit{\text{jk}}}^{\mathit{\text{AD}}}, of followup time, {t}_{\mathit{\text{jk}}}^{\mathit{\text{AD}}}, and the shape parameter, s, of the Weibull distribution. The linear predictor, log\left({\lambda}_{\mathit{\text{jk}}}^{\mathit{\text{AD}}}\right) was a function of the baseline loghazard of an event for treatment b in study j, {\mu}_{\mathit{\text{jb}}}^{\mathit{\text{AD}}}, and by the loghazard ratio for treatment k and baseline treatment b, d _{ bk }(=d _{1k }  d _{1b }) Note that there are parameters common to both model parts (equations A1 and A2), namely the loghazard ratios and the shape parameter of the time to healing distribution. Vague prior distributions were specified for {\mu}_{\mathit{\text{jb}}}^{\mathit{\text{AD}}}\left(\sim N\left(0,{10}^{6}\right)\right).
Alternative modelling assumptions
Underlying model A are a set of assumptions that are detailed below. Such assumptions have been relaxed in models B to D.
Exploring betweenstudy variation
Model A assumed that each included RCT aimed to measure a common treatment effect (fixedeffect); however, it is likely that there was betweenstudy variation. Model B included a random effect to characterise betweenstudy heterogeneity, where d _{ bk } was replaced by a study specific δ _{ jbk }, that is then described by a Normal distribution with mean d _{ bk } and precision, prec – this is common to both parts I (eq. A1) and II (eq. A2).
Time to healing distributions
Model A used the Weibull distribution to describe time to healing. Our choice of survival distribution was limited as distributions such as the LogLogistic or the LogNormal do not allow the probability of healing over time to be expressed in a closed form, and hence impede the approach proposed here for the joint synthesis of IPD and AD. Other distributions, such as the Gompertz, were not readily defined within the software used in this work (WinBugs/OpenBugs), specifically under censoring. Nonetheless, the goodness of fit could still be assessed in each IPD data source individually. To do this, we applied parametric regression survivaltime models [33] to both IPD data sources [17, 32] independently (covariates and frailty effect considered, as in model A).
Distributional shape parameter
Model A assumes that the Weibull shape parameter of the hazard of healing was common to both IPD data sources. It is possible that this parameter differed between studies, in which case HRs could be affected. Thus, we implemented two alternative NMA models to ascertain the impact of this assumption on the relative effectiveness estimates: model C1 used the shape parameter from the first IPD study to describe the AD studies and model C2 used the shape parameter from the second IPD study to describe these same studies. Because models C1 and C2 represent simple modifications of model A we do not present these algebraically.
Treatmentcovariate associations
Model A uses baseline covariates to adjust for clinical heterogeneity in the IPD. To further explore the impact of covariates on the relative treatment effects (i.e. whether they were effect modifiers), and potentially help explain betweenstudy heterogeneity, we also included interaction terms between alternative treatments and baseline ulcer area and duration— as described by Cooper et al. 2009 [34] and Saramago et al. 2012 [6]. Model D assumed a regression (slope) coefficient for the interaction terms, this effect was common across treatments and thus common to parts I (eq. D1) and II (eq. D2). This assumption was data driven, as this was the only option we were able to implement with the data available (i.e. compared to assuming ‘exchangeability’ or ‘independence’). Note that interaction estimates obtained were influenced by the full evidence base for which study mean covariate(s) values were available, including trials considering ad hoc treatments.
Model D, part I modelling the IPD studies, considering common treatmentbycovariate interactions
Time to ulcer healing was modelled in the same way as in model A (equation A1). A treatmentbycovariate interaction regression term, βx _{ ijk }, was defined, where β is the association effect, assumed common across studies and the same regardless of treatment (excluding control). We included no interaction term for each comparison of k versus b when b ≠ 1 because the common regression coefficient cancels out; but, as β _{11} = 0, we included an interaction term for each comparison of k versus treatment 1. For the remainder of the parameters of interest, vague prior distributions were assigned as in model A.
In part two of model D, {\overline{x}}_{j} represented the mean covariate value for the j ^{th} trial ((where j = NS _{ IPD } +1,…, (NS _{ IPD } + NS _{ AD })). Both the IPD and AD contributed to the estimation of this interaction term (i.e. same β term in equations D1 and D2). All other components of the model were as described for model A.
Model selection and implementation
The NMA analyses were undertaken in the WinBUGs software [35]. In all models the MCMC sampler was run for 10 000 iterations and these were discarded as ‘burnin’. Models were run for a further 5000 iterations, on which inferences were based. Chain convergence was assessed by running the model for two disparate sets of initials values, and by then checking that the BrooksGelmanRubin diagnostic converges to 1.0. The WinBUGS code is included for reference in Additional file 2. Within the NMA, goodness of fit was assessed using the deviance information criterion (DIC) [36]. Results were presented using hazard ratio estimates (and associated credibility intervals, CrIs) and also using the probability of each compression system being the ‘best’ treatment in terms of being the most clinically effective [37].
The statistical software STATA [38] was used to fit alternative time to event distributions to the IPD datasets individually. Goodness of fit was assessed with the Akaike Information Criteria (AIC) statistic [39].
Results
Table 2 shows parameter estimates obtained for model A (first column) and alternative models relaxing its assumptions (models B to D, second to fifth columns). Model A considers the baseline information for the following covariates: the log of the ulcer area (in cm^{2}) and duration (in months)  both centred to its sample mean; and dummies generated from a categorical variable on participants’ mobility with three categories (‘immobile’, ‘walks with difficulty’ and ‘walks freely’, the latter being the reference category). The results for model A highlight that the modelling framework proposed is feasible. A graphical representation of model A results is shown in Figure 2. This plots the probability of healing over time for the five main high compression ulcer treatments, but considers only the uncertainty over relative treatment effects. The results of testing the assumptions are described next.
Exploring betweenstudy variation
Despite estimates of HRs from the random effect model (model B) being associated with wider CrIs than those from model A (as expected), point estimates were found to be fairly similar except for the comparison between HH vs. 4LB: HH is estimated to be more effective in model B (HR 1.63, 95% CrI 0.763.53) compared to model A (HR: 1.05, 95% CrI 0.85 to 1.29), although the CrI of the former includes the later. The treatment with the greatest estimated probability of healing was HH in model B (59%), rather than 2LB (72%) as in model A. Differences may be explained by any existing variation between studies of SSB vs. 4LB indirectly impacting on the evidence loop 4LB vs. SSB vs. HH. Baseline covariate effect estimates remained similar. However, note that the gain in quality fitting of the randomeffects model compared to the fixedeffects is null (DIC: 5396.21 and 5396.22, respectively). Previous published work assessing evidence on the SSB vs. 4LB comparison [40] similarly found no evidence of betweenstudy heterogeneity.
Time to healing distributions
The Weibull was used to describe time to healing in model A. Whilst we were limited in the use of other distributions, goodness of fit was explored by applying alternative time to event distributions to the IPD studies individually. Table 3 shows results of such analysis (AIC statistic) and Additional file 3 graphically presents the fitted curves and the Kaplan Meier counterpart for control arms. The best fitting distributions for both studies were the LogLogistic and LogNormal. Of the remaining, the Weibull and Gompertz distributions provided better fit than the Exponential; this was expected given the flexibility of these distributions in assuming increasing, decreasing or constant hazards over time. The Weibull was best in IPD study 2 and the Gompertz best in IPD study 1.
Distributional Weibull shape parameter of healing hazard
The Weibull shape parameters estimated within models C1 and C2 indicate that in IPD study 1 [32] the hazard of healing was expected to decrease over time (s _{ 1 } = 0.93, 95% CrI 0.861.01), while in IPD study 2 it is expected to increase (s _{ 2 } = 1.27, 95% CrI 1.171.38). Note that there is no overlap in the CrIs. However, results show that relative effectiveness estimates are robust to the range of assumptions tested: the estimated HRs did not differ between models C1 and C2, and did not substantially differ from model A.
Treatmentcovariate associations
Model D tested the inclusion of interaction terms. Two covariates on ulceration area (in cm^{2}) and duration (in months) were considered and their association with treatment evaluated. Results (column 5 of Table 2) show that the covariates included did not appear to be treatment effect modifiers in this case study. However, estimating these two additional regression terms increased uncertainty in relative treatment effects estimates, specifically for ZINC and 2LB.
Discussion
This paper introduces a novel NMA modelling approach that allows IPD on time to event (with censoring) and AD on event count (for a given followup time) to be synthesised jointly, by assuming an underlying, common, distribution of time to healing. Available IPD is used directly to inform this distribution (likelihood). Studies reporting the number of participants healed (AD) are used to inform a probability parameter, and a Binomial likelihood was defined for this subset of the evidenceset. The probability of healing was then related (algebraically) to the common distribution of time to healing, by taking the duration of followup in each AD study into account. This modelling framework extends approaches in the literature [5, 14, 15] and is also a natural extension of previously published methodologies of synthesising IPD and AD jointly [6–8, 16]. This work was motivated by a real data example looking at the effectiveness of high compression treatments on the healing of venous leg ulcers.
We found that the key strength of the use of IPD in this context (additional to the known advantages described in the introduction) was the flexibility in modelling these data allowed. For example, had all evidence been available as AD, the modelling process would have been limited to the specification of uniparametric distributions for time to healing i.e. the Exponential, with constant healing hazard over time [5, 14, 15]. In our motivating example, the Exponential distribution was shown to be less adequate than other distributions in describing the time to event data in the studies for which IPD were available. The availability of IPD allowed a more complex distribution for time to event outcomes to be implemented, in this case the Weibull. This may be of particular importance when absolute effectiveness estimates, and not just relative effect estimates, are of interest – and especially where results may need to be extrapolated beyond the followup time horizon.
We note that even with the flexibility offered by the use of IPD we were, in practice, limited to using the Weibull distribution. Semiparametric methods are also unmanageable as the evidence in aggregate form consists of a single data point and thus does not allow defining a nonparametric distribution of hazard of healing. Given this limitation, the synthesis of time to event data will still often require the use of potentially suboptimal distributional assumptions, in which case estimates obtained may be biased. We suggest further research could focuses on using numerical analysis techniques within NMA, to try and resolve this issue.
This work was also relevant in again highlighting the importance of including IPD when wanting to consider patient characteristics in the model. In the presence of heterogeneity, incorporating information on either a) patients’ baseline characteristics or b) control for treatmenteffect modifiers within the synthesis model has been shown to improve estimates and, by doing so, possibly resolve evidence inconsistencies [6, 34]. The first relates to potential heterogeneity in the baseline hazard, which cannot be explicitly explored with AD alone (this is important when analyses aim to explore determinants of baseline hazard, for example). In this study we assumed a common effect of baseline covariates on the hazard of healing across IPD studies (as commonly assumed in IPD metaanalysis). The second relates to treatmentcovariate interactions, which are generally acknowledged to be best estimated using IPD, as ecological bias can be avoided [9]. For the proportion of evidence only available as AD, the model here implemented considered study level mean covariate values. Nonetheless, not all studies provided information for these, and imputation was undertaken (imputation is naturally done through the MCMC, and assumes values are ‘missing at random’).
In our work we assumed the shape parameter of the Weibull time to event distribution to be common across studies. However, testing proved this assumption was not valid, thus highlighting the importance of evaluating any assumptions of similarity imposed across studies. Despite relative effectiveness estimates being mainly unaffected, such potential heterogeneity between studies should be explored and accounted for in analyses. Such assumptions of commonality also mean that information may be shared throughout the network, in which case evidence on treatments other than those on our decision set (the five treatments of interest for which results were reported). This is the case of model D that makes use of all evidence (including ad hoc treatments) to estimate treatmentcovariate interaction, which may indirectly affect the relative effectiveness estimates of interest.
In the network of evidence there was one closed loop where both direct and indirect data informed relative treatment effect estimates. The existence of inconsistency was explored elsewhere [17], showing no evidence of statistically significant discrepancies between the direct and the indirect data. We note, however, that given the fairly high uncertainty in the evidence base, only large differences in direct and indirect data within the loop would have returned a statistically significant result.
The evidence included as AD in our case study was limited to the primary outcome data on proportion of patients healed at the end of followup. However, we could have extended the modelling framework to use any other summary data reported, such as number of patients healed at different time points or KaplanMeier (KM) curves. Information conveyed in Kaplan Meier plots can be reconstructed [41], and further research could focus on including this information in our NMA modelling framework to strengthen inferences.
Conclusions
The use of IPD for a time to event outcome is particularly useful in guiding HTA decision making by allowing flexibility in the specification of more appropriate survival distributions and in dealing with potential existing study heterogeneity [42]. Increasingly, data sharing is promoted in research [43, 44] and this example highlights how use of IPD allows the development of more informative and flexible models that are better able to summarise existing evidence. However, it is important to acknowledge that accessing and analysing IPD can be time consuming and may cause delay [4]. The process needs to be well planned and implemented. In our case both sources of IPD were easily accessed and that directly facilitated the conduct of these analyses and the associated methodological work presented here.
Endnotes
^{a.}The shape parameter of the Weibull distribution, s, can be interpreted directly as follows: i) if 0 < s < 1, hazard rate decreases over time; ii) if s = 1, hazard rate is constant over time (hazard exponentially distributed); and if s > 1, it indicates that the hazard increases with time.
^{b.}When imputing missing information, MCMC generates independent draws of the missing data from its predictive distribution. Multiple imputation through MCMC techniques is attractive for exploratory or multipurpose analyses involving a large number of estimands.
^{c.}The missingatrandom assumption (sometimes called the ignorability assumption) considers that the probability that an observation is missing may depend on the observed values but not the missing values, as sufficient data has already been collected.
Abbreviations
 2LB:

Two layer bandage
 4LB:

Four layer bandage
 AIC:

Akaike information criteria
 AD:

Aggregate data
 DIC:

Deviance information criteria
 HH:

Two layer hosiery
 HTA:

Health technology assessment
 IPD:

Individual patientlevel data
 MCMC:

Markov chain Monte Carlo
 NMA:

Network metaanalysis
 RCT:

Randomised controlled trial
 SSB:

Short stretch bandage
 ZINC:

Zinc paste bandage.
References
Bucher H, Guyatt G, Griffith L, Walter S: The results of direct and indirect treatment comparisons in metaanalysis of randomized controlled trials. J Clin Epidemiol. 1997, 50 (6): 683691.
Sutton A, Ades A, Cooper N, Abrams K: Use of indirect and mixed treatment comparisons for technology assessment. Pharmacoeconomics. 2008, 26 (9): 753767.
Dias S, Welton NJ, Sutton AJ, Ades AE: Evidence synthesis for decision making 1: introduction. Med Decis Making. 2013, 33 (5): 597606.
Saramago P, Manca A, Sutton AJ: Deriving input parameters for costeffectiveness modelling: taxonomy of data types and approaches to their statistical synthesis. Value Health. 2012, 15 (5): 639649.
Dias S, Sutton AJ, Ades AE, Welton NJ: Evidence synthesis for decision making 2: a generalized linear modeling framework for pairwise and network metaanalysis of randomized controlled trials. Med Decis Making. 2013, 33 (5): 607617.
Saramago P, Sutton AJ, Cooper NJ, Manca A: Mixed treatment comparisons using aggregate and individualparticipant level data. Stat Med. 2012, 10; 31 (28): 35163536.
Jansen J, Cope S: Network metaanalysis of individual and aggregate level data. Value Health. 2012, 15 (4): A159
Donegan S, Williamson P, D’Alessandro U, Garner P, Smith C: Combining individual patient data and aggregate data in mixed treatment comparison metaanalysis: individual patient data may be beneficial if only for a subset of trials. Stat Med. 2013, 32 (6): 914930.
Berlin J, Santanna J, Schmid C, Szczech L, Feldman H: Individual patient versus grouplevel data metaregressions for the investigation of treatment effect modifiers: ecological bias rears its ugly head. Stat Med. 2002, 21 (3): 371387.
Riley R, Lambert P, Staessen J, Wang J, Gueyffier F, Thijs L, Boutitie F: Metaanalysis of continuous outcomes combining individual patient data and aggregate data. Stat Med. 2008, 27 (11): 18701893.
Sutton AJ: Methods for MetaAnalysis in Medical Research. 2000, Chichester ; New York: J. Wiley
Welton NJ, Willis SR, Ades AE: Synthesis of survival and disease progression outcomes for health technology assessment of cancer therapies. Res Synth Method. 2010, 1: 239257. doi: 10.1002/jrsm.21
Smith C, Williamson P, Marson A: An overview of methods and empirical comparison of aggregate data and individual patient data results for investigating heterogeneity in metaanalysis of timetoevent outcomes. J Eval Clin Pract. 2005, 11 (5): 468478.
Woods B, Hawkins N, Scott D: Network metaanalysis on the loghazard scale, combining count and hazard ratio statistics accounting for multiarm trials: a tutorial. BMC Med Res Methodol. 2010, 10: 54
Soares MO, Dumville JC, Ades AE, Welton NJ: Treatment comparisons for decision making: facing the problems of sparse and few data. J R Stat Soc A Stat Soc. 2014, 177: 259279. doi: 10.1111/rssa.12010
Sutton A, Kendrick D, Coupland C: Metaanalysis of individual and aggregatelevel data. Stat Med. 2008, 27 (5): 651669.
Ashby RL, Gabe R, Ali S, Saramago P, Chuang LH, Adderley U, Bland JM, Cullum NA, Dumville JC, Iglesias CP, Soares MO, Stubbs NC, Torgerson DJ: VenUS IV (venous leg ulcer study IV): a randomised controlled trial of compression hosiery versus compression bandaging in the treatment of venous leg ulcers. Health Technology Assessment (forthcoming). 2014
Duby T, Hoffman D, Cameron J, Doblhoffbrown D, Cherry G, Ryan T: A randomized trial in the treatment of venous Leg ulcers comparing short stretch bandages, 4 layer bandage system, and a long stretch paste bandage system. Woundsa Compendium Clin Res Pract. 1993, 5 (6): 276279.
Scriven J, Taylor L, Wood Rgn A, Bell P, Naylor A, London N: A prospective randomised trial of fourlayer versus short stretch compression bandages for the treatment of venous leg ulcers. Ann R Coll Surg Engl. 1998, 80 (3): 215220.
Partsch H, Damstra R, Tazelaar D, SchullerPetrovic S, Velders A, De Rooij M, Tjon Lim Sang R, Quinlan D: Multicentre, randomised controlled trial of fourlayer bandaging versus shortstretch bandaging in the treatment of venous leg ulcers. Vasa. 2001, 30 (2): 108113.
Ukat A, Konig M, Vanscheidt W, Munter K: Shortstretch versus multilayer compression for venous leg ulcers: a comparison of healing rates. J Wound Care. 2003, 12 (4): 139143.
Franks P, Moody M, Moffatt C, Martin R, Blewett R, Seymour E, Hildreth A, Hourican C, Collins J, Heron A: Randomized trial of cohesive shortstretch versus fourlayer bandaging in the management of venous ulceration. Wound Repair Regen. 2004, 12 (2): 157162.
Junger M, Wollina U, Kohnen R, Rabe E: Efficacy and tolerability of an ulcer compression stocking for therapy of chronic venous ulcer compared with a belowknee compression bandage: results from a prospective, randomized, multicentre trial. Curr Med Res Opin. 2004, 20 (10): 16131623.
Kralj B, Kosicek M: Randomised Comparative Trial of SingleLayer and MultiLayer Bandages in the Treatment of Venous leg Ulcer. 1995, Harrogate, UK: 6th European Conference on Advances in Wound Management: 1995, 158160.
Polignano R, Bonadeo P, Gasbarro S, Allegra C: A randomised controlled study of fourlayer compression versus Unna’s boot for venous ulcers. J Wound Care. 2004, 13 (1): 2124.
Wilkinson E, Buttfield S, Cooper S, Young E: Trial of two bandaging systems for chronic venous leg ulcers. J Wound Care. 1997, 6 (7): 339340.
Colgan MP, Teevan M, McBride C, O’Sullivan L, Moore D, Shanik G: Cost Comparisons in the Management of Venous Ulceration. 1995, Harrogate, UK: 5th European Conference on Advances in Wound Management: 21–24 November 1995
Blecken S, Villavicencio J, Kao T: Comparison of elastic versus nonelastic compression in bilateral venous ulcers: a randomized trial. J Vasc Surg. 2005, 42 (6): 11501155.
Moffatt C, Edwards L, Collier M, Treadwell T, Miller M, Shafer L, Sibbald G, Brassard A, McIntosh A, Reyzelman A, Price P, Kraus SM, Walters SA, Harding K: A randomised controlled 8week crossover clinical evaluation of the 3M coban 2 layer compression system versus profore to evaluate the product performance in patients with venous leg ulcers. Int Wound J. 2008, 5 (2): 267279.
Szewczyk M, Jawie A, Cierzniakowska K, CwajdaBialasik J, Moscicka P: Comparison of the effectiveness of compression stockings and layer compression systems in venous ulceration treatment. Archives of Medical Science. 2010, 6 (5): 793799.
Wong I, Andriessen A, Charles H, Thompson D, Lee D, So W, Abel M: Randomized controlled trial comparing treatment outcome of two compression bandaging systems and standard care without compression in patients with venous leg ulcers. J Eur Acad Dermatol Venereol. 2012, 26 (1): 102110.
Iglesias C, Nelson E, Cullum N, Torgerson D: VenUS I: a randomised controlled trial of two types of bandage for treating venous leg ulcers. Health Technol Assess. 2004, 8 (29): iii, 1iii, 105.
Collett D: Modelling Survival Data in Medical Research. 2003, Boca Raton, Fla: Chapman & Hall/CRC, 2
Cooper N, Sutton A, Morris D, Ades A, Welton N: Addressing betweenstudy heterogeneity and inconsistency in mixed treatment comparisons: application to stroke prevention treatments in individuals with nonrheumatic atrial fibrillation. Stat Med. 2009, 28 (14): 18611881.
Lunn DJ, Thomas A, Best N, Spiegelhalter D: WinBUGS  a Bayesian modelling framework: concepts, structure, and extensibility. Statistics and Computing. 2000, 10 (4): 325337.
Spiegelhalter DJ, Best NG, Carlin BR, van der Linde A: Bayesian measures of model complexity and fit. J R Stat Soc B Stat Soc. 2002, 64: 583616.
Ades A, Sculpher M, Sutton A, Abrams K, Cooper N, Welton N, Lu G: Bayesian methods for evidence synthesis in costeffectiveness analysis. Pharmacoeconomics. 2006, 24 (1): 119.
StataCorp A: Stata Statistical Software: Release 12. 2011, College Station, TX: StataCorp LP
Akaike H: Information Theory and the Maximum Likelihood Principle. 1973, 2nd International Symposium in Information Theory: 1973, 267281.
O’Meara S, Tierney J, Cullum N, Bland J, Franks P, Mole T, Scriven M: Four layer bandage compared with short stretch bandage for venous leg ulcers: systematic review and metaanalysis of randomised controlled trials with data from individual patients. BMJ (Clinical research ed). 2009, 338 (7702): 10541057.
Guyot P, Ades AE, Ouwens MJ, Welton NJ: Enhanced secondary analysis of survival data: reconstructing the data from published KaplanMeier survival curves. BMC Med Res Methodol. 2012, 12: 9
Espinoza M: Heterogeneity in costeffectiveness analysis: methods to explore the value of subgroups and individualised care in a collectively funded health system. 2012, University of York
Medical Research Council Data Sharing. http://www.mrc.ac.uk/research/researchpolicyethics/datasharing/ (accessed in September 2014)
National Institute for Health Research  Health Technology Assessment Data Sharing. http://www.journalslibrary.nihr.ac.uk/informationforauthors/datasharing (accessed in September 2014)
Prepublication history
The prepublication history for this paper can be accessed here:http://www.biomedcentral.com/14712288/14/105/prepub
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
PS developed and implemented the novel synthesis models, conducted the analysis and prepared the final manuscript. LHC structured the evidence base and conducted preliminary data analysis using existing published synthesis methodology. MS oversaw and advised on all elements of the work undertaken and prepared the final manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Saramago, P., Chuang, LH. & Soares, M.O. Network metaanalysis of (individual patient) time to event data alongside (aggregate) count data. BMC Med Res Methodol 14, 105 (2014). https://doi.org/10.1186/1471228814105
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/1471228814105
Keywords
 Evidence synthesis
 Time to event data
 Survival data
 Individual patient data
 Aggregate data
 Metaanalysis
 Mixed treatment comparisons
 Networkmetaanalysis
 Treatmenteffect modifiers