 Research article
 Open Access
 Published:
Multiple imputation for estimation of an occurrence rate in cohorts with attrition and discrete followup time points: a simulation study
BMC Medical Research Methodology volume 10, Article number: 79 (2010)
Abstract
Background
In longitudinal cohort studies, subjects may be lost to followup at any time during the study. This leads to attrition and thus to a risk of inaccurate and biased estimations. The purpose of this paper is to show how multiple imputation can take advantage of all the information collected during followup in order to estimate the cumulative probability P(E) of an event E, when the first occurrence of this event is observed at t successive time points of a longitudinal study with attrition.
Methods
We compared the performance of multiple imputation with that of KaplanMeier estimation in several simulated attrition scenarios.
Results
In missingcompletelyatrandom scenarios, the multiple imputation and KaplanMeier methods performed well in terms of bias (less than 1%) and coverage rate (range = [94.4%; 95.8%]). In missingatrandom scenarios, the KaplanMeier method was associated with a bias ranging from 5.1% to 7.0% and with a very poor coverage rate (as low as 0.2%). Multiple imputation performed much better in this situation (bias <2%, coverage rate >83.4%).
Conclusions
Multiple imputation shows promise for estimation of an occurrence rate in cohorts with attrition. This study is a first step towards defining appropriate use of multiple imputation in longitudinal studies.
Background
Longitudinal studies are the most appropriate way of analysing the occurrence of an event E and of estimating its cumulative probability during followup P(E). One of the most important drawbacks of such studies is that subjects may be lost to followup at any time point. These dropouts accumulate over time and lead to attrition of the initial cohort, which has fewer and fewer participants as time goes by. Attrition can lead to biased and inaccurate estimations of P(E) [1], especially in classical completecase analyses [2, 3]. For example, in a UK cohort study of ocular outcome after premature birth, which involved 558 children born before 32 weeks of gestation in 19901991, abnormalities were more frequent among dropouts, and completecase analysis thus underestimated the prevalence of ocular abnormalities in the study population [4].
Survival analysis and KaplanMeier estimation are commonly used to examine timetoevent measurements [5, 6]. This approach takes into account the fact that subjects are followed for different lengths of time. However, it assumes that censored patients (including patients lost to followup) would have the same probability of experiencing a subsequent event as noncensored patients. This assumption has been challenged in various fields [4, 7, 8]. One specific example, which motivated the present study, is the issue of in vitro fertilization (IVF) success rates [9]. In IVF cohorts, dropout corresponds to treatment interruption, which is partly linked to a poor likelihood of success [10]. Couples who drop out therefore have a lower chance of success than couples who persevere [11]. The KaplanMeier approach tends to overestimate the IVF success rate, and alternative methods have recently been proposed: dropouts are divided into two groups depending on the chances of success (poor/not poor) [12], or, equivalently, according to whether IVF treatment is interrupted for medical causes or not [13]. Couples with a poor chance of success will be considered as having a zero probability of subsequent success, whereas those with a good prognosis will have the same chance of subsequent success as those who persevere. This method makes it possible to take into account auxiliary information on the probability of success among dropouts, through individual prognostic factors, considered in two groups. This idea may be developed by looking for ways of taking into account more of this auxiliary information. Multiple imputation is a good candidate approach in this setting.
Multiple imputation was developed thirty years ago [14] and is now used in observational as well as in randomized studies [15, 16]. However, epidemiological applications have been limited [17]. Previous works have addressed the use of multiple imputation in longitudinal studies [18–20], but few focused on estimation of the occurrence rate at the end of the followup. This is a special issue, as it implies imputation of the missing covariates and outcomes. Hsu et al. developed a weighted KaplanMeier method to take into account dependant interval censoring in the estimation of a recurrence rate [21, 22]. Our objective is similar. However, we preferred to use multiple imputation, which consists in replacing missing values with a set of plausible values, based on auxiliary information and which has the advantage of being routinely implemented and easily applicable. Moreover, IVF studies have a particular way of defining time, as each time point of the study is an IVF attempt. The length of time elapsing between the attempts is not considered. The method of Hsu et al. implied imputing midpoints of intervals and is therefore not appropriate here. KaplanMeier estimation was chosen for comparison with multiple imputation, because this method is often used to estimate cumulative success rates in IVF.
The purpose of this paper is to show how multiple imputation can take advantage of all information collected during followup. It focuses specifically on the estimation of a cumulative probability P(E) of an event E, when the first occurrence of this event is observed at t successive time points of a longitudinal study. The performance of multiple imputation was compared with that of the KaplanMeier approach by simulating several attrition scenarios.
Methods
Attrition types
Three types of missingness mechanisms can be distinguished: missing completely at random (MCAR), missing at random (MAR) and not missing at random (NMAR) [23]. If the probability of dropout depends neither on missing data nor on observed data, then the data are MCAR. In this case attrition implies a loss of power, but no bias occurs in completecase estimations. An example of MCAR situation in the context of IVF is the accidental loss of patient records. If, given the observed data, the probability of dropout does not depend on missing data, then the data are MAR. Bias may arise but it can be minimized with appropriate methods, such as multiple imputation, that take into account available information explaining the missingness mechanism. As an example, suppose that records of patients older than 35 years are recorded separately (file A for patients younger than 35, and file B for patients older than 35). The situation is MAR if all the file A is available but a part of file B is accidentally lost. Finally, if the probability of dropout depends on missing data or on other unobserved variables, then the data are NMAR. Bias due to attrition cannot be reduced in this case, and sensitivity analyses have to be conducted under various NMAR assumptions. An example is the situation in which patients are recorded differentially according to the outcome of the attempt: all patients with a delivery would be recorded, but only a part of those with pregnancy loss would. The present study focuses on MAR attrition.
Simulation study
Data were simulated using SAS Software (SAS Institute Inc. 2004. SAS/STAT^{®} 9.1 User's Guide. Cary, NC: SAS Institute Inc.).
The simulation process took place in three steps detailed below: first, a complete cohort was simulated, then attrition was added and samples were randomly drawn, and finally the different methods of estimation were applied to the samples with attrition.
Complete cohort generation
A cohort of 100 000 subjects with four followup time points was simulated. Simulated data consisted in a binary outcome E representing occurrence of the event, and a continuous covariate X predictive of both occurrence and dropout. Measures for subject i at time point t were named E _{ it } and X _{ it } . In our IVF motivating example, the outcome is delivery and the covariate is the age of the woman.
The covariate X was simulated Gaussian and the link between the successive time points was the following:
For simplicity, coefficients were taken as being equal to one: α _{ it } = β _{ it } = 1 and ${s}_{it}^{2}$ = 10.
The outcome E _{ it } was drawn from a Bernoulli distribution with a probability of occurrence of the event equal to p _{ it } , so that p _{ it } = P(E _{ it } = 1). In order to make X a predictor of E, the value of p _{ it } was linked to the value of X _{ it } . The supposition was that two subpopulations with specific probabilities of occurrence of the event coexist in the population. To achieve this, subjects were divided into two categories C _{ 1 } and C _{ 2 } according to their X _{ it } values. The probabilities of occurrence of the event p _{ it } could differ according to the category. The threshold delimiting the two categories was the median of X _{ 1 } (0.014). The same threshold was kept for all time points, i.e. C _{ 1 } was defined as X _{ it } < 0.014 whatever the value of t. For a given category of subjects, the probabilities of occurrence p _{ it } were taken as identical at each time point. The different values considered for p _{ it } were 10% and 60%.
Only the first occurrence of the event was of interest. Thus, the final endpoint E _{ i } (occurrence/nonoccurrence) was either nonoccurrence of the event in subject i if all E _{ it } values (1≤ t ≤4) were zeros, or occurrence of the event in subject i if one of the E _{ it } values was one. In other words, the final endpoint E _{ i } was a binary variable indicating whether a subject i experienced the event during followup.
The final true occurrence rate P(E) represented the cumulative occurrence rate of E _{ i } at the end of the followup, calculated in this complete cohort of 100 000 subjects.
Attrition generation
The binary variable D _{ it } indicates whether subject i dropped out immediately after time t. If D _{ ik } = 1, then X _{ ij } and E _{ ij } are missing for all j >k. Note that dropout can only occur after nonoccurrence of the event at the previous time point(s), and that all variables were measured in all subjects at the first time point (i.e. no missing data on X _{ i1 } and E _{ i1 } ).
The variable D _{ it } (interrupt/continue) was drawn from a Bernoulli distribution with a probability of dropout of r _{ it } , so that r _{ it } = P(D _{ it } = 1). In order to make X a predictor of D, the value of r _{ it } was linked to the value of X _{ it } in the same way as for p _{ it } . Indeed, the probabilities of dropout r _{ it } could differ according to the category of subjects, thus being linked to X _{ it } . This corresponds to the assumption that the probability of dropout immediately after time point t depended on the observed covariate at time point t. For a given category of subjects, the probabilities of dropout r _{ it } were taken as identical at each time point. The different values considered for r _{ it } were 10%, 30% and 60% (see Table 1).
The total dropout rate represented the cumulative dropout rate at the end of the followup, calculated in the cohort of 100 000 subjects with dropouts.
Methods
The KaplanMeier method [5, 6] was applied considering dropout as censor. Subjects encountering four nonevents were also censored at the end of the study. The event occurrence rate was estimated as [1S(t)] at time t = 4, S(t) being the survival function.
Multiple imputation consists of replacing missing values with a set of plausible values, based on auxiliary information. In this way, all subjects can be included in the analysis. Multiple imputation was implemented with the SAS MI procedure (SAS Institute Inc. 2004. SAS/STAT^{®} 9.1 User's Guide. Cary, NC: SAS Institute Inc.). Multiple imputation was performed at each time point (except the first, where all subjects were observed) for the subjects who failed and dropped out at a previous time point. For instance at time point 2, values were imputed for subjects i with E _{ i1 } = 0 (ie subjects with failure at time point 1); at period 3, values were imputed for subjects i with E _{ i2 } = 0 (ie subjects with failure at time point 2, regardless this failure is observed or imputed at the previous step). A specific SAS program executed the successive imputations one after another. The covariate X _{ it } was first imputed linearly according to the covariate and the outcome value at the previous time (X _{ it1 } and E _{ it1 } respectively). The response E _{ it } was then imputed by using logistic regression according to observed and imputed values of this covariate X _{ it } . The covariate X _{ it } was included linearly in the logistic regression to enable an evaluation of the results of the multiple imputation procedure as close as possible to its practical use, in which the threshold 0.014 is unknown. As only the first occurrence of the event was of interest, the response E _{ it } was imputed only when E _{ ij } = 0 for all j <t. Thus, the covariate X _{ it } and the intermediate outcome E _{ it } were imputed at each time point following dropout of subject i, until the event occurred or until the end of followup. Then the final endpoint E _{ i } was calculated. Following Rubin's recommendations, five imputations were made at each time point [24].
Design of the evaluation
Table 1 presents the eight scenarios that were tested and gives the final dropout and true occurrence rates resulting from each of the eight scenarios. A scenario was determined by a probability of occurrence p _{ it } and a probability of dropout r _{ it } for each category of subjects C _{ 1 } and C _{ 2 } . Scenarios explored different contrasts between the categories C _{ 1 } and C _{ 2 } , starting from scenario 1 which was the reference scenario where all probabilities equalled 10%. The robustness of the KaplanMeier method and the performance of multiple imputation were examined as the difference between the categories C _{ 1 } and C _{ 2 } was amplified.
In scenarios 2 and 3, the difference between the categories C _{ 1 } and C _{ 2 } lay only in the probabilities of occurrence p _{ it } . The probabilities of dropout r _{ it } were equal in both categories (r _{ it } = 10% for all subjects). This means that dropout did not depend on the observed value of the covariate X (nor on any other observed or unobserved value), and so that the data were missing completely at random (MCAR). In these cases, KaplanMeier and multiple imputation estimates were expected to be unbiased.
In scenarios 4 to 8, the probabilities of occurrence p _{ it } and the probabilities of dropout r _{ it } differed according to the category (p _{ it } = 10%, 30%, 60%; r _{ it } = 10%, 60%). The difference increased to peak with the last scenario, in which probabilities of occurrence and dropout corresponded to inverse situations between the categories. The difference between probabilities of dropout means that dropout depended on the observed value of the covariate X, and so that the data were missing at random (MAR). Scenario 4 was a particular case, where both categories had the same probability of occurrence of the event (p _{ it } = 10%). Consequently, even if probabilities of dropout differed, both methods were expected to be unbiased. Otherwise, multiple imputation was expected to handle MAR attrition, whereas the KaplanMeier method was not. Scenarios 4 to 6 corresponded to the situation where the second category of subjects C _{ 2 } had a higher chance of occurrence of the event p _{ it } and also a higher chance of dropout r _{ it } . This could happen when the event is negative, such as a disease: patients more likely to be ill are more likely to drop out of the cohort. Scenarios 7 and 8 corresponded to the inverse situation where the second category of subjects C _{ 2 } had a higher chance of occurrence of the event p _{ it } but a lower chance of dropout r _{ it } . This could reflect the IVF situation, where women more likely to have a child are less likely to discontinue.
For each scenario, a cohort of 100 000 subjects was created with corresponding p _{ it } and r _{ it } . For the simulation study, 500 replications of samples of 2000 subjects were drawn from each simulated cohort of 100 000 subjects. The cumulative probability P(E) of the event E during followup was estimated in each sample using both methods (KaplanMeier and multiple imputation). For each scenario, mean value of the estimate and of the standard error of the 500 replications were computed. Bias and the coverage rate for the 95% confidence interval of P(E) were also provided. Bias was defined as the mean difference over the 500 replications between the estimation obtained by the method and the true occurrence rate of the event E in the complete simulated cohort (100 000 subjects with no dropouts). The coverage rate was the percentage of times when the 95% confidence interval contained the true value of P(E).
Results
Results of the simulations for the 8 scenarios are displayed in Table 2.
As expected, in MCAR scenarios 1, 2 and 3, both methods (KaplanMeier and multiple imputation) were unbiased. Their coverage rates were close to the nominal level of 95% (range = [94.6%; 95.8%]).
In scenario 4, data were MAR but probabilities of occurrence of the event were equal in both categories of subjects. Both methods were as efficient as in MCAR scenarios: estimates were unbiased and coverage rates close to 95%. This was expected because both categories had the same probability of occurrence of the event (= 10%), so the difference in the probabilities of dropout did not impact the estimation.
In the MAR scenarios 5 to 8, the amplitude of bias obtained with the KaplanMeier method was approximately 10% of P(E) (bias≈5%) and the coverage rate was very poor (less than 15%). The bias tended towards an underestimation of the true occurrence rate of the event E when subjects in the category C _{ 2 } with a higher probability of the event had a higher probability of dropout (scenarios 5 and 6), and towards an overestimation when subjects in the category C _{ 2 } with a lower probability of the event had a higher probability of dropout (scenarios 7 and 8). The bias was worse in the last scenario where the difference between the categories was maximised.
In these MAR scenarios 5 to 8, multiple imputation yielded a very small bias (less than 2%) and a very good coverage rate (above 83%). The method seemed to be less efficient in scenarios 7 and 8 when the category C _{ 2 } with the higher probability of occurrence had the lower probability of dropout.
To point up the importance of the difference between the estimations distributions in MAR scenarios, histograms of all estimates according to the method used are provided in Figure 1 for scenarios 3 (MCAR) and 5 (MAR).
Simulations were done with a smaller sample size of 500 (results not shown). This led to larger standard errors, that mechanically increased the coverage rate. In scenarios 5 to 8, the coverage rates ranged from 41.6% to 69.4% for the KaplanMeier estimate and from 88.4% to 94.8% for multiple imputation estimations. Bias was unchanged. Conclusions remained the same concerning the performance of the methods.
Discussion
Multiple imputation has been developed to deal with missing data, which is a particular problem in longitudinal studies which suffer from attrition. However, little attention has been paid to how multiple imputation might be applied to studies with successive data collection points. The estimation of an occurrence rate in a cohort with attrition is a recurrent issue in reproductive health in general, and in IVF studies in particular. These studies have the characteristics of having discrete time points corresponding to attempts. Intervals between time points are not taken into account, thus leading to no distribution of the event times. Methods developed in previous publications thus cannot be applied here [25, 26] We therefore used a simulation study to compare the performance of multiple imputation with that of the widely used KaplanMeier method, for the estimation of an occurrence rate in a cohort with attrition.
We found that in MAR scenarios, KaplanMeier estimation performed fairly well with respect to bias, but not with respect to the coverage rate. The direction of the biases was as expected, owing to the differences in the probabilities of occurrence of the event and of dropout between the two subject categories. Indeed, when subjects dropping out tended to have a higher chance of experiencing the event, the KaplanMeier assumption that subjects dropping out have the same chance of occurrence of the event as other subjects led to an underestimation, and viceversa. The amplitude of bias was quite small. However, the KaplanMeier estimate was very precise around this biased estimation. This implies that the confidence interval rarely included the true value, and therefore that the coverage rate was poor. Thus, even if the bias is considered acceptable, one must be aware that the confidence interval will almost never contain the true value.
In contrast, multiple imputation was virtually biasfree and gave robust estimates in all the scenarios. This approach enables auxiliary information to be incorporated in the estimation. This makes it possible to take into account the fact that subjects' chances of experiencing the event may differ. This was the case in the simulated cohort in which two categories of subjects coexisted and were defined by their covariate values. Thus, multiple imputation succeeded in correcting the bias inherent in the KaplanMeier method. By replacing each missing value with several values, multiple imputation takes into account the uncertainty of the imputed value [23, 27]. Therefore, the standard errors of multiple imputation were larger than those of the KaplanMeier estimation, resulting in larger confidence intervals and better coverage rates.
In comparison with the method of Hsu et al. [22], multiple imputation makes it possible to include the information contained in auxiliary variables without summarizing them in risk scores. Moreover, it gives a good estimation of variance, as it takes into account the uncertainty of the imputed values.
Multiple imputation was evaluated here under the assumption that the probability of dropout immediately after time point t depended on the observed covariate at time point t. This type of attrition is likely to occur in most epidemiological cohorts. For example, a cohort of patients attending regular medical appointments would fit this attrition mechanism if the probability of dropping out after a given visit was likely to depend on the results of the same visit. This is the case in IVF, where failure of a procedure may indicate a poor chance of subsequent success and lead a couple to discontinue treatment. Multiple imputation at each time point is compatible with this attrition mechanism, as imputation at time point t is performed according to the observed or imputed covariate values at the same time point. This highlights the fact that data have to be investigated in order to apply a multiple imputation strategy that fits their pattern. The part of the auxiliary information that has to be included in the imputation model would be of particular interest: in some situations, baseline data could be sufficient, whereas in other situations, inclusion of intermediate measures could be advantageous. Imputation appears a preferable method than complete case analysis that assumes stronger hypotheses (MCAR data) and induces an important loss of statistical power.
The multiple imputation method could be generalized without theoretical problems to situations with a larger number of time intervals. However, it would rise practical difficulties such as a (much) longer computational time and a decreasing accuracy of the estimates, especially with limited sample size, due to categories of subjects with dropout rates almost equal to 100% as the time goes by.
Data not missing at random were not tested in this work. In practice, implementations of multiple imputation in statistical software are made under the MAR assumption. Attrition was fully explained in the simulation study, as all the variables explaining attrition were taken into account. This may not be the case with observational data. Sensitivity analyses are necessary to determine whether the results are robust to deviations from this assumption.
The second difference between simulated and observational data is that, with simulated data, the model underlying the response and the covariates is known. The relations are then correctly used for imputation. With observational data,matters can be much more difficult. The relations used for imputation rarely reproduce exactly the relations underlying the data. With observational data, multiple imputation could lose precision and eventually generate bias. The consequences of misspecification of the imputation model were not explored here and require more research.
Conclusions
Although our findings have the limitations of any simulation study, they point to a new way of estimating an occurrence rate in cohorts with attrition, using auxiliary information collected during followup. The robust estimations of multiple imputation in this simulation study are very satisfactory and promising. Our results suggest that in an observational context, comparison of the results of KaplanMeier estimation and of multiple imputation could provide a measure of the impact of attrition on the estimation. If this impact appears small, further analyses could be made using the easiertoimplement KaplanMeier approach.
This study is a first step towards defining the appropriate use of multiple imputation in longitudinal studies.
Abbreviations
 IVF:

in vitro fertilization
 MAR:

missing at random
 NMAR:

not missing at random
References
 1.
Deeg DJ: Attrition in longitudinal population studies: does it affect the generalizability of the findings?. J Clin Epidemiol. 2002, 55 (3): 213215. 10.1016/S08954356(01)004723.
 2.
Crawford SL, Tennstedt SL, McKinlay JB: A comparison of analytic methods for nonrandom missingness of outcome data. J Clin Epidemiol. 1995, 48 (2): 209219. 10.1016/08954356(94)001249.
 3.
Greenland S: Response and followup bias in cohort studies. Am J Epidemiol. 1977, 106 (3): 184187.
 4.
Pennefather PM, Tin W, Clarke MP, Dutton J, Fritz S, Hey EN: Bias due to incomplete follow up in a cohort study. Br J Ophthalmol. 1999, 83 (6): 643645. 10.1136/bjo.83.6.643.
 5.
Kaplan EL, Meier P: Nonparametric estimation from incomplete observations. J Am Stat Assoc. 1958, 53 (282): 457481. 10.2307/2281868.
 6.
Altman DG: Practical statistics for medical research. 1991, London: Chapman & Hall
 7.
Van Beijsterveldt CE, van Boxtel MP, Bosma H, Houx PJ, Buntinx F, Jolles J: Predictors of attrition in a longitudinal cognitive aging study: the Maastricht Aging Study (MAAS). J Clin Epidemiol. 2002, 55 (3): 216223. 10.1016/S08954356(01)004735.
 8.
McLean RR, Hannan MT, Epstein BE, Bouxsein ML, Cupples LA, Murabito J, Kiel DP: Elderly cohort study subjects unable to return for followup have lower bone mass than those who can return. Am J Epidemiol. 2000, 151 (7): 689692.
 9.
Daya S: Life table (survival) analysis to generate cumulative pregnancy rates in assisted reproduction: are we overestimating our success rates?. Hum Reprod. 2005, 20 (5): 11351143. 10.1093/humrep/deh889.
 10.
Olivius C, Friden B, Borg G, Bergh C: Why do couples discontinue in vitro fertilization treatment? A cohort study. Fertil Steril. 2004, 81 (2): 258261. 10.1016/j.fertnstert.2003.06.029.
 11.
Sharma V, Allgar V, Rajkhowa M: Factors influencing the cumulative conception rate and discontinuation of in vitro fertilization treatment for infertility. Fertil Steril. 2002, 78 (1): 4046. 10.1016/S00150282(02)031606.
 12.
Olivius K, Friden B, Lundin K, Bergh C: Cumulative probability of live birth after three in vitro fertilization/intracytoplasmic sperm injection cycles. Fertil Steril. 2002, 77 (3): 505510. 10.1016/S00150282(01)032174.
 13.
Stolwijk AM, Wetzels AM, Braat DD: Cumulative probability of achieving an ongoing pregnancy after invitro fertilization and intracytoplasmic sperm injection according to a woman's age, subfertility diagnosis and primary or secondary subfertility. Hum Reprod. 2000, 15 (1): 203209. 10.1093/humrep/15.1.203.
 14.
Rubin DB: Multiple imputations in sample surveys  a phenomenological Bayesian approach to nonresponse. Proc Survey Res Meth Sec (American Statistical Association). 1978, 2034.
 15.
Barzi F, Woodward M: Imputations of missing values in practice: results from imputations of serum cholesterol in 28 cohort studies. Am J Epidemiol. 2004, 160 (1): 3445. 10.1093/aje/kwh175.
 16.
Wood AM, White IR, Thompson SG: Are missing outcome data adequately handled? A review of published randomized controlled trials in major medical journals. Clin Trials. 2004, 1 (4): 368376. 10.1191/1740774504cn032oa.
 17.
Klebanoff MA, Cole SR: Use of multiple imputation in the epidemiologic literature. Am J Epidemiol. 2008, 168 (4): 355357. 10.1093/aje/kwn071.
 18.
Kristman VL, Manno M, Côté P: Methods to account for attrition in longitudinal data: do they work? A simulation study. Eur J Epidemiol. 2005, 20 (8): 657662. 10.1007/s1065400579197.
 19.
Engels JM, Diehr P: Imputation of missing longitudinal data: a comparison of methods. J Clin Epidemiol. 2003, 56 (10): 968976. 10.1016/S08954356(03)001707.
 20.
Mazumdar S, Tang G, Houck PR, Dew MA, Begley AE, Scott J, Mulsant BH, Reynolds CF: Statistical analysis of longitudinal psychiatric data with dropouts. J Psychiatr Res. 2007, 41 (12): 10321041. 10.1016/j.jpsychires.2006.09.007.
 21.
Hsu CH, Taylor JM, Murray S, Commenges D: Survival analysis using auxiliary variables via nonparametric multiple imputation. Stat Med. 2006, 25 (20): 35033517. 10.1002/sim.2452.
 22.
Hsu CH, Long Q, Alberts D: Estimation of colorectal adenoma recurrence with dependent censoring. BMC Med Res Methodol. 2009, 9 (1): 6610.1186/14712288966.
 23.
Little RJ, Rubin DB: Statistical Analysis with Missing Data. 2002, Hoboken, NJ: John Wiley & Sons Inc, 2
 24.
Rubin DB: Multiple imputation after 18+ years. J Am Stat Assoc. 1996, 91 (434): 473489. 10.2307/2291635.
 25.
Taylor JMG, Murray S, Hsu CH: Survival estimation and testing via multiple imputation. Statistics & Probability Letters. 2002, 58 (3): 221232.
 26.
Jan B, Shah S, Shah S, Fazli Qadir M: Weighted Kaplan Meier estimation of survival function in heavy censoring. Pak J Statist. 2005, 21 (1): 5563.
 27.
Schafer JL: Analysis of Incomplete Multivariate Data. 1997, Boca Raton, FL: Chapman & Hall/CRC, 72:
Prepublication history
The prepublication history for this paper can be accessed here:http://www.biomedcentral.com/14712288/10/79/prepub
Acknowledgements
This work was supported by the French Institute for Public Health Research (IReSP); and the IledeFrance Region [doctoral grant to NS].
Author information
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
NS designed and implemented the simulation study, conducted the analysis and drafted the article. ELR and JB helped to design the study and interpret the results and they critically revised the article. All the 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
About this article
Cite this article
Soullier, N., de La Rochebrochard, E. & Bouyer, J. Multiple imputation for estimation of an occurrence rate in cohorts with attrition and discrete followup time points: a simulation study. BMC Med Res Methodol 10, 79 (2010). https://doi.org/10.1186/147122881079
Received:
Accepted:
Published:
Keywords
 Multiple Imputation
 Coverage Rate
 Occurrence Rate
 Auxiliary Information
 Simulated Cohort