This article has Open Peer Review reports available.

# Comparison of the risk factors effects between two populations: two alternative approaches illustrated by the analysis of first and second kidney transplant recipients

- Katy Trébern-Launay
^{1, 2}, - Magali Giral
^{2}, - Jacques Dantal
^{2}and - Yohann Foucher
^{1}Email author

**13**:102

https://doi.org/10.1186/1471-2288-13-102

© Trébern-Launay et al.; licensee BioMed Central Ltd. 2013

**Received: **8 April 2013

**Accepted: **1 August 2013

**Published: **6 August 2013

## Abstract

### Background

Whereas the prognosis of second kidney transplant recipients (STR) compared to the first ones has been frequently analyzed, no study has addressed the issue of comparing the risk factor effects on graft failure between both groups.

### Methods

Here, we propose two alternative strategies to study the heterogeneity of risk factors between two groups of patients: *(i)* a multiplicative-regression model for relative survival (MRS) and *(ii)* a stratified Cox model (SCM) specifying the graft rank as strata and assuming subvectors of the explicatives variables. These developments were motivated by the analysis of factors associated with time to graft failure (return-to-dialysis or patient death) in second kidney transplant recipients (STR) compared to the first ones. Estimation of the parameters was based on partial likelihood maximization. Monte-Carlo simulations associated with bootstrap re-sampling was performed to calculate the standard deviations for the MRS.

### Results

We demonstrate, for the first time in renal transplantation, that: *(i)* male donor gender is a specific risk factor for STR, *(ii)* the adverse effect of recipient age is enhanced for STR and *(iii)* the graft failure risk related to donor age is attenuated for STR.

### Conclusion

While the traditional Cox model did not provide original results based on the renal transplantation literature, the proposed relative and stratified models revealed new findings that are useful for clinicians. These methodologies may be of interest in other medical fields when the principal objective is the comparison of risk factors between two populations.

## Keywords

## Background

In patients facing a first allograft loss, repeat kidney transplantation provides a better chance for both long-term survival and quality of life than a return to dialysis [1, 2]. The prognosis of second kidney transplant recipients (STR) compared to first kidney transplant recipients (FTR) has been frequently studied. The older literature tends to conclude that STR have a worse prognosis than FTR [3, 4]. However, recent analyses with adjustments for confounding factors have challenged this generally accepted idea [5, 6], with the exception of one study [7]. By modelling the time-dependent hazard between FTR and STR, we recently demonstrated that STR have a higher risk of graft failure than FTR, but this excess risk appears several years after transplantation [8]. According to this literature, one can accept that the excess risk for STR compared to FTR is negligible considering the improvements in life expectancy and quality of life compared to dialysis therapy. Nevertheless, as the demand for kidney transplants largely exceeds the supply [9], it is necessary to evaluate the differences in risk factors between STR and FTR so as to improve graft allocation.

For this purpose, traditional survival models can be used by merging STR and FTR. Nevertheless, one can notice two important limitations of these traditional approaches. Firstly, the comparisons of risk factors between both groups would imply testing all the interactions with the graft rank. Secondly, STR-specific explicative variables (survival time of the first transplant, transplantectomy or time in dialysis before re-transplantation) cannot be included despite the knowledge that their use would improve risk evaluation [5, 6, 10, 11].

To overcome these limits, in this paper we propose two alternative strategies. The first one is an adaptation of a multiplicative-regression model for relative survival (MRS). This type of relative approach is often used to study the net survival of patients with cancer, i.e. the survival of patients if the only cause of death is related to the disease [12–17]. The principle of such additive-regression models is to introduce expected mortality rates by using life tables adjusted for gender, age and calendar year. Still using life tables, Andersen et al. [18] proposed a multiplicative-regression model. To our knowledge, the development of such methodology to endpoints other than mortality and with a reference group without a life time, has never been explored.

Moreover, we propose a second method by adapting a stratified Cox model (SCM) specifying the graft rank as strata and assuming a vector of explicative variables decomposed into subvectors of variables that enter either in the reference hazard only, or in the relative hazard only, or in both groups but with common or separate effects.

## Methods

### Study population

Second transplant recipients (STR) constituted the relative group of interest. Recipients older than 18 years at the date of transplantation between 1996 and 2010 were selected from the French DIVAT (Données Informatisées et VAlidées en Transplantation - http://www.divat.fr/en) multicentric prospective cohort [19]. Codes were used to assure donor and recipient anonymity and blind assay. The ’Comité National Informatique et Liberté’ approved the study (N° CNIL 891735) and written informed consent was obtained from the participants. Only recipients with a maintenance therapy with calcineurin inhibitors, mammalian target of rapamycin inhibitors or belatacept, in addition to mycophenolic acid and steroids, were included. Simultaneous transplantations were excluded. Recipients with at least one missing data for all the variables taken into account in the expected hazard were excluded. The same criteria were applied to the reference group composed of first transplant recipients (FTR). The principal outcome was the time between transplantation and graft failure, which was the first event between return to dialysis and patient death with a functioning graft.

### The multiplicative-regression model for relative survival (MRS)

*j*(

*j*= 1,...,

*n*

_{ e }) in the reference group and by

*i*(

*i*= 1,...,

*n*

_{ r }) in the relative group.

*n*

_{ e }and

*n*

_{ r }represent the sample sizes of the respective groups. We note

*h*

^{ o }(

*t*

_{ i }|

*z*

_{ i }) the observed instantaneous hazard function at time

*t*

_{ i }for the

*i*th individual, where

*z*

_{ i }is the vector of explicative variables. The observed hazard of the

*i*th subject of the relative group can be decomposed in the multiplication of two hazards [18, 20]

where ${h}^{e}\left({t}_{i}\right|{z}_{i}^{e})$ is the expected hazard for an individual in the reference group with similar characteristics to the *i*th individual. ${z}_{i}^{e}$ is a subset of *z*
_{
i
} and represents these common characteristics. ${h}^{r}\left({t}_{i}\right|{z}_{i}^{r})$ is the relative hazard with ${z}_{i}^{r}$ being a subset of *z*
_{
i
}.

*j*th individual (

*j*= 1,...,

*n*

_{ e })

*β*

^{ e }is the vector of regression coefficients associated with ${z}_{j}^{e}$. The estimations ${\widehat{\beta}}^{e}$ are obtained by maximizing the partial log-likelihood among the

*n*

_{ e }patients of the reference group

where *δ*
_{
j
} equals 1 if the failure was observed for the *j*th subject and 0 otherwise. The estimation of the variance-covariance matrix, $\widehat{V}\left({\widehat{\beta}}^{e}\right)$, is obtained via the corresponding information matrix.

*i*th individual (

*i*= 1,...,

*n*

_{ r }), the instantaneous hazard is defined as

*β*

^{ r }is the vector of regression coefficients associated with ${z}_{i}^{r}$. By adapting the partial likelihood function (3) and assuming the previous estimations of expected parameters as constants, the regression coefficients

*β*

^{ r }are estimated by maximizing

Of note, for explicative variables not taken into account in the expected hazard in the model (1), exp(*β*
^{
r
}) represents the observed hazard ratios (HR) in the reference group, just as exp(*β*
^{
r
}) estimated from the PH model (4) among the relative group. In contrast, for explicative variables taken into account in the expected hazard model in the model (1), exp(*β*
^{
r
}) represents the weighting factors between the expected HR, i.e. $\text{exp}\left({\widehat{\beta}}^{e}\right)$, and the observed HR in the relative group, i.e. $\text{exp}\left({\widehat{\beta}}^{e}\right)\times \text{exp}\left({\beta}^{r}\right)$. In other words, for explicative variables involved in both models (2) and (4), *β*
^{
r
} = 0 means that the variable has the same effect in both groups. If *β*
^{
r
} > 0, the hazard ratio increases in the relative group compared to the reference group. If *β*
^{
r
} < 0, the hazard ratio decreases.

Nevertheless, in contrast to the traditional relative survival models based on life tables, the expected hazards cannot be reasonably assumed as constants since the corresponding parameters were estimated from the reference sample. To take into account the variability associated with the expected model (2) in the estimation of the relative model (4), we used Monte-Carlo simulations associated with bootstrap resampling [22]. At each of the *B* iterations ($b=1,\dots ,B$), this procedure can be divided into the following steps

- (a)
Generation of a vector of parameters ${\widehat{\beta}}_{b}^{e\ast}$ using the multivariate normal distribution $\mathcal{N}\left({\widehat{\beta}}^{e},\widehat{V}\left({\widehat{\beta}}^{e}\right)\right)$ obtained from the maximisation of the partial likelihood (3). This first step takes into account the variance of the expected hazard.

- (b)
Generation of a bootstrap sample from the relative sample comprising

*n*_{ r }subjects. This second step takes into account the variance due to sample-to-sample fluctuation. - (c)
From this bootstrap sample, the model (4) is estimated by maximizing (5) in which the simulated parameters ${\widehat{\beta}}_{b}^{e\ast}$ are used instead of ${\widehat{\beta}}^{e}$. ${\widehat{\beta}}_{b}^{r}$ is the resulting estimation of the relative regression coefficients.

Means, standard deviations and 95% confidence intervals can be calculated from the *B* estimations of ${\widehat{\beta}}_{b}^{e\ast}$.

### The stratified Cox model (SCM)

*i*th individual of the overall sample (

*i*= 1,...,

*n*), i.e. the reference and the relative group together (

*n*=

*n*

_{ e }+

*n*

_{ r }), let

*z*

_{ i }be the vector of covariates that are applied to the model for either first or second transplant recipients. Let

*k*be the indicator of the strata with

*k*=

*e*for the reference group and

*k*=

*r*for the relative group. The stratified Cox model is given by

*h*

_{ k,0}(

*t*

_{ i }) the baseline hazard function in the strata

*k*. The vector

*z*

_{ i }can be decomposed into four different subvectors:

*(i)*the vector ${z}_{i}^{e}$ of explicative variables that enter in the reference hazard only (their values equal 0 if

*k*=

*r*);

*(ii)*the vector ${z}_{i}^{r}$ of explicative variables that enter in the relative hazard only (their values equal 0 if

*k*=

*e*);

*(iii)*the vector ${z}_{i}^{c}$ of explicative variables associated to both groups; and

*(iv)*the vector ${z}_{i}^{s}$ a subvector of explicative variables included in ${z}_{i}^{c}$ but with separate effects. The model (6) can be developed as follows

*β*= (

*β*

^{ e },

*β*

^{ c },

*β*

^{ s },

*β*

^{ r }) the vector of regression coefficients and

*δ*

_{ i r }equals 1 if the subject

*i*belongs to the strata

*k*=

*r*and 0 otherwise. Then, for the reference group, we obtain

Therefore, *e* *x* *p*(*β*
^{
e
}) represents the HR associated with specific variables for the reference group. And *e* *x* *p*(*β*
^{
r
}) represents the HR associated with specific variables for the relative group. For variables only included in *z*
^{
c
} (and not in *z*
^{
s
}), *e* *x* *p*(*β*
^{
c
}) represents the common HR in both groups associated with *z*
^{
c
}. For variables included in both *z*
^{
c
} and *z*
^{
s
}, the HR associated with *z*
^{
c
} equals *e* *x* *p*(*β*
^{
c
}) in the reference group and *e* *x* *p*(*β*
^{
c
}+*β*
^{
s
}) in the relative group. Thus, *e* *x* *p*(*β*
^{
s
}) represents the weighting factors between the expected HR, i.e. *e* *x* *p*(*β*
^{
c
}), and the observed HR in the relative group, i.e. *e* *x* *p*(*β*
^{
c
}) × *e* *x* *p*(*β*
^{
s
}).

### Evaluation of the proportional hazards assumption

In models (2), (4) and (7), hazard proportionality was checked for each explicative variable by plotting log-minus-log survival curves obtained by the Kaplan and Meier estimator [23] and by testing the scaled Schoenfeld residuals [24] separately in the reference and relative samples. If the observed hazard ratios are constant regardless of time in both groups, the ratio between both observed hazard ratios, i.e. the weighting factor *e* *x* *p*(*β*
^{
r
}), will also be constant.

### Software

All statistical analyses were performed using R version 2.15.1 [25]. The proposed multiplicative-regression model for relative survival was implemented in an R package MRsurv available at http://www.divat.fr/en/softwares/mrsurv. The adapted stratified Cox model was implemented by using the R package survival (function coxph, option strata).

## Results

### Description of the cohort

641 STR potentially made up the relative group of interest, but 75 STR (11.7%) with missing data for explicative variables of the expected hazard were excluded. Finally, 566 STR were included in the group of interest. The mean follow-up was 3.1 years with a maximum of 13.1 years. During the observation period, 72 returns to dialysis and 34 deaths were observed. We identified 2462 FTR who met the inclusion criteria. We excluded 256 FTR (10.3%) with one missing data for at least one of the variables taken into account in the expected hazard model. Finally, 2206 FTR made up the reference group. The mean follow-up was 3.4 years with a maximum of 13.7 years. During the observation period, 191 returns to dialysis and 109 deaths were observed.

**Demographic characteristics at the date of transplantation for**
(i)
**whole cohort and**
(ii)
**FTR and STR separately; the last three rows of the table concern covariates specific for STR**

All (N = 2772) | FTR (N = 2206) | STR (N = 566) | |||||
---|---|---|---|---|---|---|---|

Demographic characteristics | Number (%) | Number (%) | Number (%) | p-value | |||

Transplantation period < 2005 | 594 | (21.4) | 457 | (20.7) | 137 | (24.2) | 0.0806 |

Recipient ≥ 55 years of age | 1175 | (42.4) | 994 | (45.1) | 181 | (32.0) | <0.0001 |

Male recipient | 1705 | (61.5) | 1362 | (61.7) | 343 | (60.6) | 0.6536 |

Recurrent causal nephropathy | 906 | (32.7) | 666 | (30.2) | 240 | (42.4) | <0.0001 |

History of diabetes | 306 | (11.0) | 269 | (12.2) | 37 | (6.5) | 0.0002 |

History of hypertension | 2263 | (81.6) | 1804 | (81.8) | 459 | (81.1) | 0.7545 |

History of vascular disease | 352 | (12.7) | 272 | (12.3) | 80 | (14.1) | 0.2804 |

History of cardiac disease | 903 | (32.6) | 686 | (31.1) | 217 | (38.3) | 0.0012 |

History of dyslipemia | 799 | (28.8) | 661 | (30.0) | 138 | (24.4) | 0.0104 |

History of malignancy | 228 | (8.2) | 147 | (6.7) | 81 | (14.3) | <0.0001 |

History of hepatitis B or C | 168 | (6.1) | 96 | (4.4) | 72 | (12.7) | <0.0001 |

Recipient BMI ≥ 30 kg.m | 263 | (9.5) | 235 | (10.7) | 28 | (4.9) | <0.0001 |

Positive anti-class I PRA | 706 | (25.5) | 355 | (16.1) | 351 | (62.0) | <0.0001 |

Positive anti-class II PRA | 733 | (26.4) | 319 | (14.5) | 414 | (73.1) | <0.0001 |

Donor ≥ 55 years of age | 1172 | (42.3) | 973 | (44.1) | 199 | (35.2) | 0.0002 |

Deceased donor | 2470 | (89.1) | 1940 | (87.9) | 530 | (93.6) | 0.0002 |

Donor serum creatinine ≥ 133 | 342 | (12.5) | 279 | (12.8) | 63 | (11.4) | 0.3807 |

Positive donor EBV serology | 2613 | (94.3) | 2087 | (94.6) | 526 | (92.9) | 0.1540 |

HLA-A-B-DR incompatibilities > 4 | 365 | (13.2) | 326 | (14.8) | 39 | (6.9) | <0.0001 |

Cold ischemia time ≥ 24h | 754 | (27.2) | 552 | (25.0) | 202 | (35.7) | <0.0001 |

Lymphocyte-depleting induction | 1223 | (44.1) | 793 | (35.9) | 430 | (76.0) | <0.0001 |

First graft survival < 1 year | - | - | - | - | 131 | (24.1) | - |

Waiting time before regraft ≥ 3 years | - | - | - | - | 272 | (49.8) | - |

Transplantectomy of the first graft | - | - | - | - | 220 | (38.9) | - |

Among FTR meeting the inclusion criteria, some patients were also part of the STR group as they had received two transplants during the observation period. These 37 patients, who were included in both cohorts, represented 2% and 7% of the FTR and STR groups respectively. Given the large number of explicative variables, it seemed reasonable to assume conditional independence of the two transplantations of a given patient. In order to validate this assumption, we performed a frailty Cox model [26] based on the 37 individuals who were included in both groups. The frailty term was assumed to be Gamma distributed. The variance of the random variable was estimated at 5.10 ^{-9} (p = 0.9948). Therefore, no intra-individual dependency was demonstrated. In order to validate the robustness of the results, we also performed both models after exclusion of the 37 STR also included in FTR. These results are presented in Additional files 1 and 2.

### Analysis of risk factors in the FTR sample

*(i)*for all of the possible pre- or per-transplant immunological and non-immunological confounding factors according to experts and

*(ii)*for all the baseline parameters differentially distributed between FTR and STR. All together, the expected hazard of graft failure was estimated according to recipient age and gender, causal nephropathy, comorbidities (including history of diabetes, hypertension, cardiac or vascular disease, dyslipemia, hepatitis B or C, and malignancy), obesity, pre-transplant immunization (panel reactive antibody, PRA) against class I and class II antigens), donor age, deceased or living donor status, Epstein-Barr Virus (EBV) serology, period of transplantation, level of HLA-A-B-DR mismatches, induction therapy and cold ischemia time. This modelling is explained in detail in the paper by Trébern-Launay et al. [8]. The final multivariate model in the reference group of FTR is presented in the first three columns of Table 2.

**Multivariate Cox model for FTR and results of the MRS in the STR group :**
(i)
**the first three columns provides the results of the multivariate Cox model analysis of graft failure risk factors for FTR (N = 2206);**
(ii)
**the next three columns provide the results of the relative survival model based on 540 STR (26 recipients presenting missing data for the waiting time before re-transplantation were excluded)**

Cox model in | MRS in | |||||
---|---|---|---|---|---|---|

the FTR group | the STR group | |||||

HR | 95% CI | p-value | HR | 95% CI | p-value | |

| ||||||

Causal nephropathy (recurrent / non recurrent) | 1.24 | 0.96-1.59 | 0.0987 | - | - | - |

History of diabetes (positive / negative) | 1.34 | 0.96-1.85 | 0.0819 | - | - | - |

History of hypertension (positive / negative) | 0.77 | 0.57-1.05 | 0.0986 | - | - | - |

History of cardiac disease (positive / negative) | 1.41 | 1.11-1.79 | 0.0051 | - | - | - |

History of vascular disease (positive / negative) | 1.10 | 0.81-1.51 | 0.5351 | - | - | - |

History of dyslipemia (positive / negative) | 1.12 | 0.87-1.45 | 0.3828 | - | - | - |

History of hepatitis B/C (positive / negative) | 0.82 | 0.45-1.47 | 0.4969 | - | - | - |

History of malignancy (positive / negative) | 1.25 | 0.84-1.86 | 0.2698 | - | - | - |

Body mass index (≥ 30 kg.m-2 / < 30 kg.m-2) | 1.58 | 1.12-2.14 | 0.0084 | - | - | - |

Anti-class I PRA (positive / negative) | 1.45 | 1.07-1.97 | 0.0182 | - | - | - |

Anti-class II PRA (positive / negative) | 1.09 | 0.78-1.52 | 0.6299 | - | - | - |

Donor status (deceased/living) | 2.50 | 1.41-4.43 | 0.0016 | - | - | - |

Donor EBV serology (positive / negative) | 1.65 | 0.98-2.78 | 0.0606 | - | - | - |

Number of HLA-A-B-DR mismatches (> 4 / ≤ 4) | 1.30 | 0.97-1.76 | 0.0824 | - | - | - |

Induction therapy (depleting / non depleting) | 0.79 | 0.60-1.05 | 0.1091 | - | - | - |

Cold ischemia time (≥ 24 h / < 24 h) | 1.29 | 1.01-1.66 | 0.0441 | - | - | - |

| ||||||

Transplantation period (< 2005 / ≥ 2005) | 1.33 | 0.97-1.82 | 0.0693 | 0.97 | 0.55-1.74 | 0.9360 |

Recipient gender (male / female) | 1.17 | 0.91-1.51 | 0.2186 | 0.61 | 0.38-1.05 | 0.0720 |

Recipient age (≥ 55 years / < 55 years) | 1.39 | 1.05-1.83 | 0.0204 | 1.65 | 1.01-2.72 | 0.0480 |

Donor age (≥55 years / <55 years) | 1.34 | 1.03-1.74 | 0.0313 | 0.59 | 0.33-0.99 | 0.0440 |

| ||||||

Donor gender (male / female) | - | - | - | 1.53 | 1.03-2.48 | 0.0320 |

Waiting time before regraft ≥ 3 years | - | - | - | 1.92 | 1.22-3.00 | <0.0001 |

### Relative hazard modelling in the STR group using the MRS

A first selection of variables was performed (p < 0.20) followed by a step-by-step descending procedure (Wald test with p < 0.05). In line with the requirements of additive-regression models, adjustments were forced for recipient gender and age and transplantation period. All the variables were categorized in order to avoid any log-linearity assumption and to obtain interpretable results.

The final relative model is presented in the last three columns of Table 2. Expected HR previously estimated in FTR are presented in the first columns to enable a direct comparison between FTR (Cox model) and STR (relative model). Donor gender and waiting time before retransplantation were not taken into account in the expected hazards for FTR. Donor gender was not a significant risk factor for FTR and waiting time is by definition a specific factor for STR. More precisely, we estimated a 1.5-fold increase in risk of graft failure for STR with grafts from males compared to STR with grafts from females (p = 0.0320). Moreover, STR who waited more than 3 years in dialysis before retransplantation had a 1.9-fold increased risk compared to STR with a shorter waiting time (p < 0.0001).

In contrast, the effect of recipient age and donor age seemed significantly different between FTR and STR (p < 0.05). More precisely, if we assumed a similar effect of recipient age between both groups, the expected HR associated with recipient age ≥ 55 years would be 1.39 in the STR group, regarding the HR observed in the FTR group. In fact, the relative model showed that this HR was 1.6-fold higher for STR compared to FTR (CI95% = [1.01-2.72], p = 0.0480). Similarly, the effect of donor age ≥ 55 years was nearly two fold lower for STR than for FTR (CI95% = [0.33-0.99], p = 0.0440), while it was identified as a significant risk factor for FTR (HR = 1.34, p = 0.0313). Of note, the relationship between the recipient gender and the risk of graft failure was not found to be significantly different between FTR and STR (p = 0.0720).

### Relative hazard modelling in the STR group using the SCM

*z*

^{ r }. The other four variables (transplantation period, recipient gender and age, and donor age) were included in

*z*

^{ c }and

*z*

^{ s }to evaluate the difference in their effect between both groups. The results are presented in Table 3.

**Results of the stratified Cox model based on 2746 patients with 2206 FTR and 540 STR (26 STR presenting missing data for the waiting time before re-transplantation were excluded)**

FTR strata | STR strata | |||||
---|---|---|---|---|---|---|

HR | 95% CI | p-value | HR | 95% CI | p-value | |

| ||||||

Causal nephropathy (recurrent / non recurrent) | 1.12 | 0.90-1.39 | 0.3031 | - | - | - |

History of diabetes (positive / negative) | 1.28 | 0.95-1.72 | 0.1001 | - | - | - |

History of hypertension (positive / negative) | 0.88 | 0.68-1.14 | 0.3270 | - | - | - |

History of cardiac disease (positive / negative) | 1.35 | 1.10-1.66 | 0.0042 | - | - | - |

History of vascular disease (positive / negative) | 1.11 | 0.85-1.47 | 0.4433 | - | - | - |

History of dyslipemia (positive / negative) | 1.19 | 0.95-1.49 | 0.1263 | - | - | - |

History of hepatitis B/C (positive / negative) | 0.97 | 0.65-1.46 | 0.9091 | - | - | - |

History of malignancy (positive / negative) | 1.21 | 0.87-1.67 | 0.2646 | - | - | - |

Body mass index (≥ 30 kg.m-2 / < 30 kg.m-2) | 1.54 | 1.13-2.08 | 0.0057 | - | - | - |

Anti-class I PRA (positive / negative) | 1.39 | 1.07-1.82 | 0.0153 | - | - | - |

Anti-class II PRA (positive / negative) | 0.98 | 0.73-1.31 | 0.8857 | - | - | - |

Donor status (deceased/living) | 2.06 | 1.27-3.36 | 0.0036 | - | - | - |

Donor EBV serology (positive / negative) | 1.65 | 1.07-2.54 | 0.0235 | - | - | - |

Number of HLA-A-B-DR mismatches (> 4 / ≤ 4) | 1.33 | 1.01-1.75 | 0.0397 | - | - | - |

Induction therapy (depleting / non depleting) | 0.88 | 0.70-1.11 | 0.2742 | - | - | - |

Cold ischemia time (≥ 24 h / < 24 h) | 1.20 | 0.97-1.49 | 0.0894 | - | - | - |

| ||||||

Transplantation period (< 2005 / ≥ 2005) | 1.42 | 1.09-1.86 | 0.0099 | 0.94 | 0.54-1.64 | 0.8295 |

Recipient gender (male / female) | 1.17 | 0.91-1.50 | 0.2200 | 0.63 | 0.40-1.02 | 0.0581 |

Recipient age (≥ 55 years / < 55 years) | 1.36 | 1.03-1.78 | 0.0274 | 1.60 | 0.95-2.72 | 0.0785 |

Donor age (≥55 years / <55 years) | 1.36 | 1.04-1.77 | 0.0238 | 0.60 | 0.35-1.05 | 0.0725 |

| ||||||

Donor gender (male / female) | - | - | - | 1.51 | 0.97-2.36 | 0.0674 |

Waiting time before regraft ≥ 3 years | - | - | - | 1.99 | 1.29-3.07 | 0.0019 |

Estimations and corresponding 95% confidence intervals were very similar to those obtained in the MRS. Indeed, as in the previous model, we estimated a 2-fold increase in risk of graft failure for STR who waited more than 3 years in dialysis before retransplantation compared to STR who waited less than 3 years (p = 0.0019). The relationship between the donor gender and the risk of graft failure among STR was similar to that obtained in the MRS but was not found to be significant (HR = 1.51, p = 0.0674).

Transplantation period, recipient gender, recipient age and donor age were included in variables applied in both models. Results were also concordant with the MRS. For the four explicative variables, estimations and corresponding 95% confidence intervals were similar to those obtained in the MRS. However, conversely to the MRS, recipient age and donor age were not found to be significantly differently associated with the risk of graft failure between the two groups.

## Discussion

Although the comparison of survival between first and second kidney transplants has been frequently performed, no study has addressed the issue of comparing the risk factors associated with the time to graft failure between both groups. Understanding the factors influencing the long-term evolution of STR compared to FTR would benefit the medical management of graft attribution by identifying patients with the best chances.

The absence of literature focusing on this question may be partially explained by the methodological issues associated with such studies. Indeed, the Cox model is classically used to explore risk factors influencing graft survival and interactions can be included to evaluate risk factor differences between FTR and STR. However, this approach has several limitations. Firstly, it implies testing interactions between the graft rank and each explicative variable, increasing the number of parameters and making interpretations difficult. Secondly and certainly more importantly, only covariates common to both groups can be taken into account. This excludes explicative variables specific for one group. Concerning our application, this constitutes a limitation as several STR-specific explicative variables are known to be associated with second graft prognosis: the first graft transplantectomy [10], the first graft survival duration [6, 11] or the time in dialysis before re-transplantation [5].

This paper describes two alternative models to overcome these difficulties. Firstly, the adaption of a multiplicative-regression model for relative survival allows a direct comparison of risk factors between two groups of patients without presupposing the role of each variable, i.e. common, different or specific relationships. The corresponding semi-parametric models are the Cox model for the expected hazard and the multiplicative-regression for the relative hazard. The main difficulty and limit of these models is the estimation of standard deviations which were obtained by Monte-Carlo simulations associated with bootstrap re-sampling. In this multiplicative modelling, the regression coefficients are straightforward to interpret in terms of their interactions. We propose an R package for a simple way of using the model.

Secondly, we demonstrated that a stratified Cox model specifying the graft rank as strata may be fitted to take into account STR-specific variables as a subvector of variables that enter in the model for STR only. In addition, some variables would enter either in the model for FTR only or in both models (with common or separate effects). The main limit to this approach is that the corresponding structure presupposes knowledge of variables potentially applicable to both models (in contrast with the relative model) unless testing a very large number of models. Indeed, whereas explicative variables entering in a single model (for FTR or STR) would easily be clinically assumed, those applicable to both models and with common or separate effects are not known in advance. Nevertheless, the SCM can be simply estimated by maximising a single partial likelihood function.

As expected, the results were concordant between both approaches. Regression coefficients were similar while standard deviations appeared a little smaller with the MRS approach. The results showed that male donor gender and a long waiting time before retransplantation were two specific-STR risk factors: donor gender was not significantly associated with the risk of graft failure in the FTR population and the waiting time before retransplantation was only related to STR by definition. The interpretations were similar to hazard ratios from a Cox model preformed on the STR group.

Conversely, two explicative variables appeared to be differently associated with the risk of graft failure between STR and FTR. More precisely, we showed for the first time that the adverse effect of recipient age was enhanced for STR as compared to FTR. The main clinical explanation is a cumulative effect of the risk factors for STR, in particular because of the cumulative exposure to immunosuppressive drugs during the first transplantation period. From a clinical point of view, this result may imply that clinicians should pay particular attention to recipient age in second kidney transplantations. Also, for the first time to our knowledge, this study identified an attenuation of the risk factor related to older transplants for STR as compared to FTR. Two explanations are: *(i)* an indication bias with only high-quality donors (without diabetes, hypertension or cardiovascular disease) proposed to STR; *(ii)* a higher non-HLA immunization in STR, explaining why graft failure is generally due to immunological phenomena rather than transplant quality.

Although we illustrated the advantages of both alternative approaches in renal transplantation, this methodology may be useful in number of other clinical and epidemiological applications. For practical use, we propose an R package to compute the MRS. The adaptation of the SCM can be computed by using many statistical software. Of course, the aim of such models is not to replace traditional survival models, but rather to provide a more suitable alternative when the main objective is to compare risk factors between two populations, in particular when population-specific covariates need to be included.

As always, there are several avenues worth exploring from this work. First, both models can be generalized for time-dependent explicative variables by adapting the likelihood functions as proposed by Therneau and Grambsch [27] (chapter 5, pages 111-115). Second, both models assumed the independence of FTR and STR. While this assumption was evaluated by using a frailty model among the 37 individuals common in both groups, a low statistical power may explain the non-rejection of this independence hypothesis. To ensure the validity of our results, we reperformed both models after exclusion of the corresponding 37 STR also included in FTR. The results presented in additional files showed the robustness of the results. Third, other strategies for variable selection can be adapted, such as partial likelihood generalization. Finally, further work is needed to develop a Goodness-of-fit statistic for the MRS approach, in particular concerning the proportional hazards assumption.

## Conclusions

MRS and SCM consitute two original approaches to compare risk factors between two populations. The advantage of MRS is to allow a direct modelling strategy but it is not straightforward to estimate the standard deviations. In contrast, SCM allows an overall estimation of parameters and standard deviations but its structure presupposes knowledge of the role of each explicative variable. This study also highlighted novel risk factor differences between first and second kidney transplant recipients. These results could help improve the management of patients waiting for a second graft. They may also encourage the widespread use of this original methodology in other medical fields.

## Declarations

### Acknowledgements

We thank Dr. J.Ashton-Chess for editing the manuscript. We also wish to thank the DIVAT scientific council and the members of the clinical research assistant team (M. Kessler, M. Ladrière, JP. Soulillou, C. Legendre, H. Kreis, G. Mourad, V. Garrigue, L. Rostaing, N. Kamar, E. Morelon, F. Buron, S. Le Floch, C. Scellier, V. Eschbach, P. Przednowed, K. Zurbonsen, V. Godel, K. Zurbonsen, X. Longy, C. Dagot, F. M’Raiagh) for the collection of the data in the DIVAT cohort and finally Dr. M. Pohar Perme for her helpful advice. This work was partly supported by the RTRS, the Fondation de Coopération Scientifique CENTAURE and Roche Laboratory.

## Authors’ Affiliations

## References

- Ojo A, Wolfe R, Agodoa L, Held P, Port F, Leavey S, Callard S, Dickinson D, Schmouder R, Leichtman A: Prognosis after primary renal transplant failure and the beneficial effects of repeat transplantation: multivariate analyses from the United States Renal data system. Transplantation. 1998, 66 (12): 1651-1659. 10.1097/00007890-199812270-00014.View ArticlePubMedGoogle Scholar
- Rao P, Schaubel D, Wei G, Fenton S: Evaluating the survival benefit of kidney retransplantation. Transplantation. 2006, 82 (5): 669-674. 10.1097/01.tp.0000235434.13327.11.View ArticlePubMedGoogle Scholar
- Stratta R, Oh C, Sollinger H, Pirsch J, Kalayoglu M, Belzer F: Kidney retransplantation in the cyclosporine era. Transplantation. 1988, 45: 40-45. 10.1097/00007890-198801000-00010.View ArticlePubMedGoogle Scholar
- Almond P, Matas A, Gillingham K, Troppmann C, Payne W, Dunn D, Sutherland D, Najarian J: Risk factors for second renal allografts immunosuppressed with cyclosporine. Transplantation. 1991, 52 (2): 253-258. 10.1097/00007890-199108000-00013.View ArticlePubMedGoogle Scholar
- Arnol M, Prather J, Mittalhenkle A, Barry J, Norman D: Long-term kidney regraft survival from deceased donors: risk factors and outcomes in a single center. Transplantation. 2008, 86 (8): 1084-1089. 10.1097/TP.0b013e318187ba5c.View ArticlePubMedGoogle Scholar
- Wang D, Xu T, Chen J, Wu W, Yang S, Lin W, Cai J, Tan J: Factors influencing second renal allograft survival: a single center experience in China. Transpl Immunol. 2009, 20 (3): 150-154. 10.1016/j.trim.2008.09.010.View ArticlePubMedGoogle Scholar
- Magee J, Barr M, Basadonna G, Johnson M, Mahadevan S, McBride M, Schaubel D, Leichtman A: Repeat organ transplantation in the United States, 1996-2005. Am J Transplant. 2007, 7 (5 Pt 2): 1424-1433.View ArticlePubMedGoogle Scholar
- Trébern-Launay K, Foucher Y, Giral M, Legendre C, Kreis H, Kessler M, Ladrière M, Kamar N, Rostaing L, Garrigue V, Mourad G, Morelon E, Soulillou JP, Dantal J: Poor long-term outcome in second kidney transplantation: a delayed event. PLoS ONE. 2012, 7 (10): e47915-10.1371/journal.pone.0047915.View ArticlePubMedPubMed CentralGoogle Scholar
- Wolfe R, McCullough K, Leichtman A: Predictability of survival models for waiting list and transplant patients: calculating LYFT. Am J Transplant. 2009, 9 (7): 1523-1527. 10.1111/j.1600-6143.2009.02708.x.View ArticlePubMedGoogle Scholar
- Abouljoud M, Deierhoi M, Hudson S, Diethelm A: Risk factors affecting second renal transplant outcome, with special reference to primary allograft nephrectomy. Transplantation. 1995, 60 (2): 138-144.View ArticlePubMedGoogle Scholar
- Gjertson D: A multi-factor analysis of kidney regraft outcomes. Clin Transpl. 2002, 335-349.Google Scholar
- Ederer F, Axtell L, Cutler S: The relative survival rate: a statistical methodology. Natl Cancer Inst Monogr. 1961, 6: 101-121.PubMedGoogle Scholar
- Hakulinen T, Tenkanen L: Regression analysis of relative survival rates. J R Stat Soc Series C (Appl Stat). 1987, 36 (3): 309-317.Google Scholar
- Dickman P, Sloggett A, Hills M, Hakulinen T: Regression models for relative survival. Stat Med. 2004, 23: 51-64. 10.1002/sim.1597.View ArticlePubMedGoogle Scholar
- Perme M, Henderson R, Stare J: An approach to estimation in relative survival regression. Biostatistics. 2009, 10: 136-146.View ArticlePubMedGoogle Scholar
- Pohar Perme M, Stare J, Estève J: On estimation in relative survival. Biometrics. 2012, 68: 113-120. 10.1111/j.1541-0420.2011.01640.x.View ArticleGoogle Scholar
- Eloranta S, Lambert P, Andersson T, Czene K, Hall P, Björkholm M, Dickman P: Partitioning of excess mortality in population-based cancer patient survival studies using flexible parametric survival models. BMC Med Res Methodol. 2012, 12: 86-10.1186/1471-2288-12-86.View ArticlePubMedPubMed CentralGoogle Scholar
- Andersen P, Borch-Johnsen K, Deckert T, Green A, Hougaard P, Keiding N, Kreiner S: A Cox regression model for the relative mortality and its application to diabetes mellitus survival data. Biometrics. 1985, 41 (4): 921-932. 10.2307/2530964.View ArticlePubMedGoogle Scholar
- Ladrière M, Foucher Y, Legendre C, Kamar N, Garrigue V, Morélon E, Kessler M, Soulillou J: M G: The western europe cohort of kidney transplanted recipients - the DIVAT network. Clin Transpl. 2010, 460-461.Google Scholar
- Buckley J: Additive and multiplicative models for relative survival rates. Biometrics. 1984, 40: 51-62. 10.2307/2530743.View ArticlePubMedGoogle Scholar
- Cox D: Regression models and life-tables. J R Stat Soc Series B (Methodol). 1972, 34 (2): 187-220.Google Scholar
- Efron B: Bootstrap methods: another look at the jackknife. Ann Stat. 1979, 7: 1-26. 10.1214/aos/1176344552.View ArticleGoogle Scholar
- Kaplan E, Meier P: Nonparametric estimation from incomplete observations. J Am Stat Assoc. 1958, 53 (282): 457-481. 10.1080/01621459.1958.10501452.View ArticleGoogle Scholar
- Grambsch P, Therneau T: Proportional hazards tests and diagnostics based on weighted residuals. Biometrika. 1994, 81: 515-526. 10.1093/biomet/81.3.515.View ArticleGoogle Scholar
- Team RDC: R: A Language and Environment for Statistical Computing. 2010, Vienna: R Foundation for Statistical Computing, http://www.R-project.org/,Google Scholar
- Ripatti S, Palmgren J: Estimation of multivariate frailty models using penalized partial likelihood. Biometrics. 2004, 56 (4): 1016-1022.View ArticleGoogle Scholar
- Therneau TM, Grambsch PM: Modeling Survival Data: Extending the Cox Model (Statistics for Biology and Health). 2000, New York: Springer VerlagView ArticleGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2288/13/102/prepub

### Pre-publication history

## Copyright

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.