This article has Open Peer Review reports available.

# Use of re-randomized data in meta-analysis

- Iztok Hozo
^{1}Email author, - Benjamin Djulbegovic
^{2}, - Otavio Clark
^{3}and - Gary H Lyman
^{4}

**5**:17

https://doi.org/10.1186/1471-2288-5-17

© Hozo et al; licensee BioMed Central Ltd. 2005

**Received: **26 September 2004

**Accepted: **10 May 2005

**Published: **10 May 2005

## Abstract

### Background

Outcomes collected in randomized clinical trials are observations of random variables that should be independent and identically distributed. However, in some trials, the patients are randomized more than once thus violating both of these assumptions. The probability of an event is not always the same when a patient is re-randomized; there is probably a non-zero covariance coming from observations on the same patient. This is of particular importance to the meta-analysts.

### Methods

We developed a method to estimate the relative error in the risk differences with and without re-randomization of the patients. The relative error can be estimated by an expression depending on the percentage of the patients who were re-randomized, multipliers (how many times more likely it is to repeat an event) for the probability of reoccurrences, and the ratio of the total events reported and the initial number of patients entering the trial.

### Results

We illustrate our methods using two randomized trials testing growth factors in febrile neutropenia. We showed that under some circumstances the relative error of taking into account re-randomized patients was sufficiently small to allow using the results in the meta-analysis. Our findings indicate that if the study in question is of similar size to other studies included in the meta-analysis, the error introduced by re-randomization will only minimally affect meta-analytic summary point estimate.

We also show that in our model the risk ratio remains constant during the re-randomization, and therefore, if a meta-analyst is concerned about the effect of re-randomization on the meta-analysis, one way to sidestep the issue and still obtain reliable results is to use risk ratio as the measure of interest.

### Conclusion

Our method should be helpful in the understanding of the results of clinical trials and particularly helpful to the meta-analysts to assess if re-randomized patient data can be used in their analyses.

## Background

Statistical tests performed to evaluate differences in the outcomes collected in randomized clinical trials (RCTs) are based on the assumption that all the data come from random variables that are independent of each other. However, in some trials, the patients are randomized more than once, usually twice. Investigators often pool these events and report data for all randomized patient episodes instead of episode per patient. Therefore, the independence assumption can be violated.

Concerns about the re-randomization of patients on febrile neutropenia trials have been raised by some authors [1] due to a possible violation of the principle of independency. However, this practice is accepted by the Immunocompromised host society as a valid one [2] and has been used in some of the trials performed in neutropenic patients.

We have recently encountered this problem during our conduct of meta-analysis (MA) on the role of colony-stimulating factors in the treatment of febrile neutropenia [3] where two out of the 12 trials reported analysis according to febrile neutropenia episodes instead of episodes per patients. In this paper, we address the issue whether the combined (pooled) risk difference (RD) after re-randomization is considerably different from the RD if the patients were not re-randomized.

## Methods

### Assumptions

Assuming that once an event (FN episode) occurred in a patient the probability of another occurrence is likely to be different from its initial value, we introduce a parameter *d* measuring this change:

Previously, Chouaid et al. [4] provided the data on the probability of FN after each randomization. Data were available from 39 patients, of which seven had an event; six of these patients were re-randomized and four of them had another occurrence of event.

Thus, accordingly an estimate for *d* can be obtained by
. To avoid messy notation, we will omit the 'hat' from the estimates of these parameters. This indicates that the patient who develops an FN episode after first cycle will have about 3.71 times higher likelihood of development of FN in the second chemotherapy cycle. Note that our calculations also work for values *d* < 1.

On the other hand, the patients who did not experience an episode of FN are usually less likely to experience it after re-randomization. The probability of an FN episode in a patient who didn't experience an episode of FN in the first round is *c* times smaller (or larger, if *c >* 1) in the second round of randomization. In this example,
.

## Results

### Calculations

Assume that *N* patients in a clinical control trial are being randomized (for the first time) in such way that *N*
_{
C
}patients are randomized into the control group and *N*
_{
A
}patients are selected into the experimental (alternate) group. In order to simplify our calculations, we have set *N*
_{
A
}≈ *N*
_{
C
}≈ *N* / 2. The percentage of people who were re-randomized (randomized for the second time) is *x*. After re-randomization, the total number of events is reported as *E*
_{
CT
}for the control group and *E*
_{
AT
}for the experimental group. The total number of events reported is given by *E*
_{
T
}= *E*
_{
AT
}+ *E*
_{
CT
}.

Consider the following two variables:

*E*
_{
A
}– the number of events in experimental group in the initial randomization;

*E*
_{
C
}– the number of events in control group in the initial randomization;

*E* = *E*
_{
A
}+ *E*
_{
C
}– the total number of events in the initial randomization;

*NE* = *N* - *E* – the number of patients without an event in the initial randomization;

Assuming that re-randomization is approximately evenly divided (if not, we can change the fraction 1/2 with appropriate value), the re-randomized control group and re-randomized experimental group will have the following sizes:

Having in mind that the reoccurrence of an event is *d* times more likely (or less likely if *d* < 1) if the patient comes from the *E* – group (patients with events), and *c* times less likely (or more likely if *c* > 1) if it comes from NE-group (patients without events), and that the events occur at the initial risks *p*
_{
A
}and *p*
_{
C
}for the experimental and control groups, respectively, we have:

Therefore, the risk difference of the re-randomization alone, can be calculated as

The overall risk difference for both randomizations, the initial one and the re-randomization, can then be calculated as

The relative error between the overall risk difference and the risk difference resulting from considering the patients only from the initial randomization is then estimated by the expression

and

Adding these two equations, we get the quadratic equation in terms of *p*
, which can be solved as:

*E*

_{ T }=

*E*

_{ AT }+

*E*

_{ CT }. Replacing this resulting estimate into the relative difference equation above, leads to

which after some simplifying becomes our main formula:

Therefore, the relative error of the reported risk difference and the actual risk difference before re-randomization can be estimated by the expression depending on *x* – the percentage of the patients who were re-randomized, *d* and *c* – multipliers for the probability of reoccurrence of the event, and the ratio
of the total events reported and the initial number of patients entering the trial.

An Excel file performing all of these calculations and allowing the readers to enter their own data is enclosed as an Additional File 1.

*x*– the percent of the patients who were re-randomized. The horizontal axis is

*d*– the multiplier by which the probability of an event changed after an initial occurrence of an event. Figure helps to determine under which circumstances the error in risk difference is less than 5%. The curves represent 5% relative error level curves for different values of the ratio of the total events reported and the initial number of patients entering the trial. The values of the ratio are indicated at the top of the graph.

### Risk ratio

*N*

_{ A }≈ (

*N*/ 2) ≈

*N*

_{ C }and . Then, using the definitions of

*p*

_{ A }and

*p*

_{ C }, we obtain the following equation:

which is exactly the risk ratio of the initial randomization. Therefore, the overall risk ratio after the re-randomization is equal to the initial risk ratio after the first randomization. Interestingly enough, the risk ratio is not affected by the re-randomization.

### Odds ratio

### Data from CSF trials

To illustrate our method, we will use examples from two papers that studied use of CSFs in the treatment of chemotherapy-related febrile neutropenia (FN) and reported the number of patients who were hospitalized and remained in the hospital longer than 10 days [5, 6]. In both of these trials some of the patients were re-randomized after first FN episode (hospitalization).

Mitchell et al[5] studied 112 patients in the first randomization for the effect of G-CSF in FN. Seventy-four of these patients were re-randomized (*x* = 66.07%). From a graph in their paper we deduced that 18% and 20% of the patients on GCSF and placebo, respectively, remained in the hospital longer than 10 days. This translates into 17 and 18 events in these two treatment arms
. From the graph in the Figure 3, and assuming [4] that *d* = 3.71 and *c* = 0.40, we estimate the relative difference between the reported risk difference evaluated for all the patients, and the risk difference resulting from considering the patients only in the first initial randomization to be only 0.70%. The point corresponding to this example is marked by the purple circle in the Figure 3. In our case, the point is safely between the two identical
level curves (the first and the third curve from the right) indicating that the relative error is smaller than 5%.

Anaissie et al. [6] studied effect of GM-CSF on FN in 92 patients, of which 15.9% were re-randomized (*x* = 0.159). Overall, 18 and 26 patients in GM-CSF and placebo arm remained in the hospital after 10 days
. The difference is now 8.23%. The point corresponding to this example is indicated by the blue square in the Figure 3. Note that this time the point is not between two
-level curves, indicating that the relative error is larger than 5%.

The odds ratio for this data also changes from *OR* = 0.554 before the re-randomization to the overall value after re-randomization *OR* = 0.535, a change of 3.49%.

### Estimating the risks in the initial randomization

Meta-analysts would also like to know at what risks the events occur in the experimental and control groups before any re-randomization is performed. For example, while in the first example the relative difference between the reported risk difference and the risk difference before the re-randomization was not considerable; in the second example we found a considerable difference.

Some meta-analysts might prefer to use the risks of the events from the initial randomization, i.e., before any re-randomization was performed. In this section we will try to estimate these risks.

*p*, we have , where

*E*

_{ T }=

*E*

_{ AT }+

*E*

_{ CT }. Using this parameter, the formulas for

*p*

_{ A }and

*p*

_{ C }are easily obtained as:

and

Using these estimates, we can get the estimates of *E*
_{
A
}and *E*
_{
C
}– the number of events in the initial randomization in experimental and control arms of the trial, respectively.

If the number of patients in each arm of the trial in the initial randomization is given, the events are estimated as *E*
_{
A
}= *p*
_{
A
}·*N*
_{
A
}and *E*
_{
C
}= *p*
_{
C
}·*N*
_{
C
}. If the only number we have is the total number of patients entering the initial randomization, *N*, the events must be estimated more crudely as
and
.

In the case of Anaissie et al. example, we obtain the following estimates for the number of events: *E*
_{
A
}= 14 and *E*
_{
C
}= 21, and since we are not given the number of patients in each arm, we have to estimate each to 46. Therefore, the numbers we could use for a meta-analysis from this particular trial would be

Experimental 14 / 46

Control 21 / 46

## Discussion

Our method has been developed using the example of febrile neutropenia (FN) in which the patients who once developed FN are considered to be at increased risk for subsequent event [7]. We assume that it is not possible to reliably predict which patients will develop the event of interest from the outset of the trial, but that once the event occurred the risk for subsequent event will be higher in those patients who had developed the first event. Because of this we adopted the values of *d* > 1 and *c* < 1. However, it is conceivable that biology of the process may differ in different circumstances. A meta-analyst is well advised to work with the content-specific experts to address this issue. The same formulas can be used to reflect any combination of the values for *d* and *c* parameters, except when *c* = *d*.

A reader should note that in our model the parameters *c* and *d* work in "conditional probability" framework, i.e., it represents a multiplicative rather than additive approach to asses the effect of previous occurrence of the event on the next one. This makes sense from purely clinical point of view.

### Error analysis

A natural question to ask at this point is "What is the extent of error in the estimation of the parameters *c* and *d* in these formulas, and in particular, in the formula (1.1)."

We first attempted to determine confidence intervals for the parameters *c* and *d*. Unfortunately, since the proportions on top and bottom of the formulas for evaluation of *c* and *d* are obviously dependent random variables, none of the classical examples for evaluation of a confidence interval of a ratio of two proportions can be applied in this case. For example, Sutton et al [8] recommend that a confidence interval for a ratio of two risks
can be calculated using the formula
, where
. However, this formula, as well as most other formulas, are derived from Taylor's expansion of a transformation of random variables, and depend on the assumption that the ratios *a* / *b* and *c* / *d* are independent (have zero covariance). In our case this assumption can not be applied.

Error analysis using Anaissie et al. data

≈ c | ||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|

≈ d | 0.0 | 0.1 | 0.2 | 0.3 | 0.4 | 0.7 | 0.9 | 1.2 | 1.5 | 1.7 | 2.0 | |

| -.286 | -.272 | -.259 | -.245 | -.232 | -.197 | -.162 | -.128 | -.096 | -.064 | .033 | |

| -.198 | -.186 | -.175 | -.163 | -.152 | -.122 | -.093 | -.064 | -.036 | -.009 | .018 | |

| -.131 | -.121 | -.111 | -.101 | -.091 | -.065 | -.039 | -.014 | .011 | .035 | .059 | |

| -.077 | -.068 | -.059 | -.050 | -.041 | -.018 | .005 | .028 | .050 | .072 | .093 | |

| -.032 | -.024 | -.016 | -.008 | .021 | .042 | .063 | .083 | .103 | .123 | ||

| .062 | .069 | .075 | .082 | .089 | .106 | .123 | .140 | .156 | .173 | .189 | |

| .131 | .136 | .142 | .147 | .153 | .168 | .182 | .196 | .211 | .225 | .238 | |

| .183 | .188 | .193 | .198 | .203 | .229 | .241 | .254 | .266 | .278 | ||

| .226 | .231 | .235 | .239 | .244 | .255 | .267 | .278 | .289 | .300 | .311 | |

| .262 | .266 | .270 | .274 | .278 | .288 | .299 | .309 | .319 | .329 | .339 | |

| .293 | .296 | .300 | .303 | .307 | .317 | .326 | .335 | .345 | .354 | .363 |

The top row of Table 1 represents the estimate of the parameter *c* as it varies from 0 to 2, while the leftmost column indicates value of the estimate of the parameter *d* as it varies from 1 to 15. The numbers in the table show the difference between the real relative difference from formula (1.1) and the relative difference calculated with these alternate values for *c* and *d*. The value .000 (boxed) in the table corresponds to the difference from the RD calculated with *c* = 0.4 and *d* = 3.7 in our example (Anaissie et al. [6]). However, if we assume that *c* = 0.7, and that *d* = 9.4, the value of RD (circled) will be larger by 21.6%, i.e., instead of 8.23%, the value of relative Risk Difference given by the formula (1.1) will be 29.83%. In this case we estimated that the original number of patients who underwent first randomization was 12/46 in experimental group, and 17/46 in the control group, respectively (see above). Note that the while some variation exists, the relative risk difference is not extraordinarily sensitive to the variation in the parameters *c* and *d*.

## Conclusion

In general, it is considered that meta-analyses represent the best method to provide reliable assessment of the effectiveness of health care interventions. However, reliability of meta-analyses depends on the validity of underlying methods used to combine data from the individual studies. One of the key assumptions for the valid MA is the independence of the data being analyzed. As illustrated here, in practice this assumption is often violated. The question for meta-analysts then becomes under which circumstances data can still be properly analyzed even if the independence axiom was violated. In the case of FN, for example, Immunocompromised Host Society issued a guideline in which they considered it acceptable to have multiple randomization of the same patients and reintroduce fully recovered patients into clinical trials [2, 9]. However, we believe that the answer regarding the acceptability of re-randomized data is not a simple categorical "yes or no" answer, but it will depend on the amount of error introduced into calculation. If the relative error of risk differences in the results with and without re-randomization is considered small, then such calculations can be safely accepted by meta-analysts. On the other hand, if such an error is considered large, then the analysts should not use the re-randomized data. Instead, the estimates of initial randomization from the method illustrated in the previous section should be used.

It is very intriguing that if the summary effect measure selected is the risk ratio, under the assumptions of our model there is no difference in the final outcome before and after the re-randomization. Therefore, if a meta-analyst is concerned about the effect of re-randomization on the meta-analysis, one way to sidestep the issue and still obtain reliable results is to use risk ratio as the measure of interest.

In this particular example – change in the input for the meta-analysis did not change the outcome and only marginally changed the pooled risk difference from -9% to -8%. In general, the effect of a given error on the results of meta-analysis can be approximated by multiplying the weight of a given study by the estimate of its relative error. In our example, weight was ~8% and relative error was also ~8% resulting in the change of the summary point estimate in meta-analysis less than 1%. Even when the relative error is high as 29.83% (see above), the summary point estimate would change only by 2.3%. However, should this study have contributed to the meta-analysis with the weight of 40%, the overall change in the summary point estimate would be about 12% (0.4 × 0.3 = 0.12). In general, however, it appears that if the study in question is of small size relative to the rest of the studies included in the meta-analysis, the error introduced by re-randomization will only minimally affect meta-analytic summary point estimate.

Although we developed our method to help with meta-analyses, it should be stressed that our method is also relevant to the analysis of the original trials, and may help assess impact of re-randomization on the results of the trial. We strongly suggest that meta-analysts always compare overall results of their meta-analysis to results based on first randomizations only, as a general control of the reliability of the meta-analysis. In order to make that possible researchers conducting randomized clinical trials would have to report a much more detailed information on their methodology, specifying the numbers of patients being re-randomized, experiencing only one occurrence of an event, experiencing two occurrences of an event, or dropped out of the study.

We believe that the method provided in this paper represents a valuable contribution to the improvement of interpretation of results of clinical trials and in particular their use in meta-analysis. Although the development of this method was stimulated by a concrete problem in the field of febrile neutropenia, practice of re-randomization is used in other areas of medicine such as testing of treatments for epilepsy, recurrent headaches and anti-emetics in cancer chemotherapy.

## Declarations

## Authors’ Affiliations

## References

- Elting LS PM: Controversies and quandaries: The design of clinical trials in fever and neutropenia. Textbook of febrile neutropenia. Edited by: Rolston KVI, Rubenstein EB. 2001, London , Martin Dunitz, p. 27-44.Google Scholar
- Immunocompromised Host Society. The design, analysis, and reporting of clinical trials on the empirical antibiotic management of the neutropenic patient. Report of a consensus panel. J Infect Dis. 1990, 161 (3): 397-401.Google Scholar
- Clark OA, Lyman G, Castro AA, Clark LG, Djulbegovic B: Colony stimulating factors for chemotherapy induced febrile neutropenia. Cochrane Database Syst Rev. 2003, CD003039-Google Scholar
- Chouaid C, Bassinet L, Fuhrman C, Monnet I, Housset B: Routine use of granulocyte colony-stimulating factor is not cost-effective and does not increase patient comfort in the treatment of small-cell lung cancer: an analysis using a Markov model. J Clin Oncol. 1998, 16 (8): 2700-2707.PubMedGoogle Scholar
- Mitchell PL, Morland B, Stevens MC, Dick G, Easlea D, Meyer LC, Pinkerton CR: Granulocyte colony-stimulating factor in established febrile neutropenia: a randomized study of pediatric patients. J Clin Oncol. 1997, 15 (3): 1163-1170.PubMedGoogle Scholar
- Anaissie EJ, Vartivarian S, Bodey GP, Legrand C, Kantarjian H, Abi-Said D, Karl C, Vadhan-Raj S: Randomized comparison between antibiotics alone and antibiotics plus granulocyte-macrophage colony-stimulating factor (Escherichia coli-derived in cancer patients with fever and neutropenia. Am J Med. 1996, 100 (1): 17-23. 10.1016/S0002-9343(96)90006-6.View ArticlePubMedGoogle Scholar
- Ozer H, Armitage JO, Bennett CL, Crawford J, Demetri GD, Pizzo PA, Schiffer CA, Smith TJ, Somlo G, Wade JC, Wade JL, Winn RJ, Wozniak AJ, Somerfield MR: 2000 update of recommendations for the use of hematopoietic colony-stimulating factors: evidence-based, clinical practice guidelines. American Society of Clinical Oncology Growth Factors Expert Panel. J Clin Oncol. 2000, 18 (20): 3558-3585.PubMedGoogle Scholar
- Sutton AJ: Methods for meta-analysis in medical research. Wiley series in probability and statistics. 2000, Chichester, West Sussex, England ; New York , John Wiley, xvii, 317-Google Scholar
- Paesmans M: Statistical considerations in clinical trials testing empiric antibiotic regimens in patients with febrile neutropenia. Support Care Cancer. 1998, 6 (5): 438-443. 10.1007/s005200050191.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2288/5/17/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.