Open Access Research article

Estimating time-to-onset of adverse drug reactions from spontaneous reporting databases

Fanny Leroy12*, Jean-Yves Dauxois3, Hélène Théophile45, Françoise Haramburu45 and Pascale Tubert-Bitter12

Author Affiliations

1 Inserm, CESP Centre for research in Epidemiology and Population Health, U1018, Biostatistics Team, F-94807 Villejuif, France

2 Univ Paris-Sud, UMRS1018, F-94807 Villejuif, France

3 Université de Toulouse-INSA, IMT UMR CNRS 5219, Toulouse, France

4 Département de pharmacologie, Centre de pharmacovigilance, CHU de Bordeaux, Bordeaux, France

5 INSERM U657, Bordeaux, France

For all author emails, please log on.

BMC Medical Research Methodology 2014, 14:17  doi:10.1186/1471-2288-14-17

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2288/14/17


Received:10 October 2013
Accepted:22 January 2014
Published:3 February 2014

© 2014 Leroy et al.; licensee 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.

Abstract

Background

Analyzing time-to-onset of adverse drug reactions from treatment exposure contributes to meeting pharmacovigilance objectives, i.e. identification and prevention. Post-marketing data are available from reporting systems. Times-to-onset from such databases are right-truncated because some patients who were exposed to the drug and who will eventually develop the adverse drug reaction may do it after the time of analysis and thus are not included in the data. Acknowledgment of the developments adapted to right-truncated data is not widespread and these methods have never been used in pharmacovigilance. We assess the use of appropriate methods as well as the consequences of not taking right truncation into account (naive approach) on parametric maximum likelihood estimation of time-to-onset distribution.

Methods

Both approaches, naive or taking right truncation into account, were compared with a simulation study. We used twelve scenarios for the exponential distribution and twenty-four for the Weibull and log-logistic distributions. These scenarios are defined by a set of parameters: the parameters of the time-to-onset distribution, the probability of this distribution falling within an observable values interval and the sample size. An application to reported lymphoma after anti TNF- α treatment from the French pharmacovigilance is presented.

Results

The simulation study shows that the bias and the mean squared error might in some instances be unacceptably large when right truncation is not considered while the truncation-based estimator shows always better and often satisfactory performances and the gap may be large. For the real dataset, the estimated expected time-to-onset leads to a minimum difference of 58 weeks between both approaches, which is not negligible. This difference is obtained for the Weibull model, under which the estimated probability of this distribution falling within an observable values interval is not far from 1.

Conclusions

It is necessary to take right truncation into account for estimating time-to-onset of adverse drug reactions from spontaneous reporting databases.

Keywords:
Pharmacovigilance; Reporting databases; Right truncation; Parametric estimation; Maximum likelihood estimation; Bias; Simulation study

Background

Identifying and preventing adverse drug reactions are major objectives of pharmacovigilance. Owing to design constraints, pre-marketing clinical trials fail to identify rare events, which lead in the last decades to an increased focus placed on the development of post-marketing surveillance methods [1-11]. Post-marketing spontaneous reporting of suspected adverse drug reactions has proved a valuable resource for signal detection [12-17]. It has recently been suggested that the modeling of the time-to-onset of adverse drug reactions could be a useful adjunct to signal detection methods, either from spontaneous reports [18,19] or longitudinal observational data [20]. Timely acquiring knowledge with respect to the time-to-onset distribution of adverse drug reactions contributes to meeting pharmacovigilance objectives. Early estimation procedures tailored to available pharmacovigilance data, i.e. spontaneous reporting data, should be sought.

The data consisting of the time-to-onset among patients who were reported to have potentially developed an adverse drug reaction are right-truncated. Truncation arises because some patients who were exposed to the drug and who will eventually develop the adverse drug reaction may do it after the time of analysis (Figure 1). Among patients exposed to the drug, only those who experienced adverse reactions before time of analysis are included in the database. No information is available for the other patients. If all the patients begin their treatment at the same time, the data are right-truncated with a single truncation time. If they do not all begin their treatment at the same time, the data are right-truncated with different truncation times. In spontaneous reporting, data are right-truncated with different truncation times and they require appropriate statistical methods.

thumbnailFigure 1. Right truncation and data on time-to-onset of adverse drug reactions from spontaneous reporting databases. Some patients who were exposed to the drug and who will eventually develop the adverse drug reaction may do it after the time of analysis. Here, in these hypothetical examples, the patient on the top line is included in the database because he experienced the adverse drug reaction before the time of analysis, i.e.x1 ≤ t1. The patient on the bottom line is not included in the database because he has not yet experienced the adverse drug reaction, i.e. <a onClick="popup('http://www.biomedcentral.com/1471-2288/14/17/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/14/17/mathml/M1">View MathML</a>, when data are analyzed.

This paper investigates parametric maximum likelihood estimation of the time-to-onset distribution of adverse drug reactions from spontaneous reporting data for different types of hazard functions likely to be encountered in pharmacovigilance. Acknowledgment of the developments adapted to right-truncated data is not widespread and these methods have never been used in pharmacovigilance. No simulation studies are available on the accuracy of their estimates. Furthermore, a naive approach that does not take into account right truncation features of spontaneous reports and uses classical parametric methods instead of appropriate methods may lead to misleading estimates. We consider the two approaches, i.e. taking or not taking right truncation into account, and the corresponding parametric maximum likelihood estimators. Both approaches are compared with a simulation study conducted to evaluate the consequences, notably in terms of bias, of not considering right truncation on the maximum likelihood estimates, as well as assessing the performances of the right truncation-based estimation. We also apply these methods to a set of 64 cases of lymphoma occurring after anti TNF- α treatment from the French pharmacovigilance.

Methods

Proper estimation of the time-to-onset distribution

We consider a given time of analysis and the population of exposed patients who will eventually experience the adverse drug reaction before they die. Let X be the time-to-onset of the adverse drug reaction of interest in that population and F its cumulative distribution function one is willing to estimate. Observations arising from n reported cases are (x1,t1),(x2,t2),…,(xn,tn), where xi is the time-to-onset calculated as the lag between the time of the occurrence of the reaction and the time of initiation of treatment, and ti is the truncation time calculated as the lag between the time of analysis and the time of initiation of treatment. Let t be the maximum of the observed truncation times. All observed data meet the condition xi ≤ ti.

We consider a parametric model for the time-to-onset X, with cumulative distribution function F (x; θ) and density f(x; θ), and derive the following maximum likelihood estimations of θ.

When right truncation, i.e. the condition xi ≤ ti, is ignored, the likelihood of the sample is written as:

<a onClick="popup('http://www.biomedcentral.com/1471-2288/14/17/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/14/17/mathml/M2">View MathML</a>

maximizing this likelihood yields the naive estimator of θ.

When right truncation is considered, the likelihood is modified. Observed times-to-onset consist of n independent realizations of random variables with respective distribution the conditional distribution of Xi given {Xi ≤ ti}, that is with cumulative distribution function <a onClick="popup('http://www.biomedcentral.com/1471-2288/14/17/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/14/17/mathml/M3">View MathML</a> and density <a onClick="popup('http://www.biomedcentral.com/1471-2288/14/17/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/14/17/mathml/M4">View MathML</a>. The likelihood is now written as:

<a onClick="popup('http://www.biomedcentral.com/1471-2288/14/17/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/14/17/mathml/M5">View MathML</a>

the maximum likelihood estimator from this likelihood, <a onClick="popup('http://www.biomedcentral.com/1471-2288/14/17/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/14/17/mathml/M6">View MathML</a>, is the proper estimation of θ and is called the truncation-based estimator (TBE).

The non-parametric maximum likelihood estimation for right-truncated data was developed and used to estimate the incubation period distribution for AIDS [21,22]. However, in a non-parametric setting, one can only estimate the distribution function conditional on the time to event as being less than t :

<a onClick="popup('http://www.biomedcentral.com/1471-2288/14/17/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/14/17/mathml/M7">View MathML</a>

where the vj’s are the m distinct values of the xi’s, i =1,…,n, taken by <a onClick="popup('http://www.biomedcentral.com/1471-2288/14/17/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/14/17/mathml/M8">View MathML</a> patients and <a onClick="popup('http://www.biomedcentral.com/1471-2288/14/17/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/14/17/mathml/M9">View MathML</a> for 1 ≤ j ≤ m, I denoting the indicator function. The unconditional distribution function is not identifiable, as F (t) is not known and cannot be estimated from the data.

In a parametric framework, the unconditional distribution is completely specified by a parameter θ of finite dimension. Maximum likelihood estimation of the parameter of interest can be conducted with the conditional distributions that describe the observations and the unconditional distribution can be estimated secondarily by <a onClick="popup('http://www.biomedcentral.com/1471-2288/14/17/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/14/17/mathml/M10">View MathML</a>. Hence parametric maximum likelihood estimation is potentially more useful than non-parametric estimation since the unconditional distribution is of interest for pharmacovigilance purposes [18,20].

Simulation study

Some adverse reactions have a very short time-to-onset, from several minutes to several hours after the beginning of treatment. Others occur only after several days, weeks, months or even years of exposure. This variation depends on numerous factors such as the pharmacokinetics of the drug and its metabolites, or the pathophysiological mechanism of the effect. The multiplicity of the underlying mechanisms results in a range of possible hazard functions that can be observed in pharmacovigilance [23]. The simplest model is given by a constant hazard function of time; the corresponding distribution is the exponential distribution with a rate parameter λ. Effects may also have an early or a late onset, the latter being the case for instance, when the rate of occurrence of the adverse reaction depends on the duration of exposure. Two distribution families among others make it possible to handle a wide range of hazard functions: the Weibull distributions and the log-logistic distributions (Table 1). Both are defined with two scalar parameters (λ,β); λ is the scale parameter and β is the shape parameter. The hazard function for the Weibull model is increasing if β > 1, decreasing if β < 1 and constant if β = 1 where it reduces to the exponential distribution. The hazard function for the log-logistic model is decreasing if β < 1 and has a single maximum if β > 1. We therefore consider the families of the exponential, Weibull and log-logistic distributions.

Table 1. Exponential, Weibull and log-logistic distributions

The times-to-onset were generated from these three distributions. Two values of λ were considered for the exponential distribution: 0.05 and 1. The same values were used for the scale parameter λ of the Weibull and log-logistic distributions. For the shape parameter β, the values 0.5 and 2 were chosen. The truncation times were uniformly distributed in [0,τ]. Survival and truncation times were independently generated. For a chosen value of p, with p representing the probability of X falling within the observable values interval [0, τ], the parameter τ was determined as P (X< τ) = p. The probability 1 - p is also a lower bound of the actual proportion of truncated data P (X> T), the truncation time T being randomly generated. The probability p was chosen in {0.25, 0.50, 0.80}. The sample size n was chosen in {100, 500}. For each drawn pair (X,T), if the time-to-onset was shorter than the truncation time, then the pair was included in the data. If not, another pair (X,T) was generated. Pairs were generated until the sample size of observations included was equal to n.

Parametric likelihood maximization with and without considering right truncation were performed for each generated sample. An iterative algorithm is necessary to solve this optimization problem except for the naive exponential estimation. Calculations were made with the R [24] function maxLik from the package maxLik. For each set of simulation parameters, 1000 replications were run.

Application study

We analyzed 64 French cases of lymphoma that occurred after anti TNF- α treatment using the national pharmacovigilance database at the date of February 1, 2010 [25]. The population included patients suffering from rheumatoid arthritis, Crohn’s disease, ankylosing spondylitis, psoriatic arthritis, psoriasis, Sjögren’s syndrome, dermatomyositis, polymyositis or polyarthropathy and exposed to one or (successively) more of the three anti TNF- α available at the study date: etanercept, adalimumab and infliximab. The occurrence of a malignant lymphoma was confirmed by histopathological analysis. Marketing authorization was obtained in August 1999 for infliximab, in September 2002 for etanercept and in September 2003 for adalimumab. These 64 adverse effects occurred between July 2001 and October 2009. None of the survival or truncation times was missing in the database. The observed maximum truncation time was 529 weeks.

All anti TNF-agents taken together, we derived the parametric maximum likelihood estimates and secondarily corresponding estimated mean times, with and without considering right truncation, for the exponential, Weibull and log-logistic distributions. For completeness, we also derived the non-parametric maximum likelihood estimation.

The French pharmacovigilance database is developed by the French drug agency (Agence Nationale de Sécurité du Médicament et des produits de santé, ANSM) and is not publicly available. It is build up and used on an ongoing basis by the network of regional pharmacovigilance centres, which have a direct access to the data. This set of data has already been extracted for another study [25] with the authorization of the ANSM and the network of regional centres, according to the internal rule.

Results

Simulation study

For each set of simulations parameters, for both approaches and for both parameters, the bias and the mean squared error of the parametric maximum likelihood estimator, based on the 1000 replications, were calculated as well as the proportion of replications where the estimate is larger than the true value. As the iterative algorithm may fail to find a maximum, those three quantities were actually calculated on the replications where there was no problem of maximization. The mean squared error is a measure of the dispersion of the estimator around the true value of the parameter - the smaller the better - and is used for global comparative purposes between two estimation procedures, as it incorporates both the variance of the estimator and its bias. The proportion of replications where the estimate is larger than the true value makes it possible to know if the estimators tend to overestimate or underestimate systematically the true value of the parameter.

Bias and mean squared error

For both approaches, for all distributions and for both parameters, the smaller is p, the larger are the bias and the mean squared error (Tables 2, 3 and 4). This increase with p is smaller for the parameter β than for the parameter λ. These estimators tend to be positively biased. However, the bias might be almost naught for the TBE. The bias and the mean squared error of the naive estimator are always larger than the bias and the mean squared error of the TBE, but to a lesser extent for the parameter β. When the sample size n increases, the bias and the mean squared error are almost constant for the naive estimator, while for the TBE, they decrease clearly (Tables 2, 3 and 4). The naive estimator might be unacceptably large whatever the value of p, whereas the TBE shows good performances when p is equal to 0.8, and often even less according to the distribution.

Table 2. Simulation results: estimations of bias and mean squared error for the exponential model

Table 3. Simulation results: estimations of bias and mean squared error for the Weibull model

Table 4. Simulation results: estimations of bias and mean squared error for the log-logistic model

Proportion of replications where the estimator is larger than the true value

For both approaches, for all distributions and for both parameters, Tables 5, 6 and 7 show that the naive estimator of λ appears to be almost always larger than the theoretical value λ, and that this is not far from being true for the naive estimator of β. This suggests that the naive estimator of λ might be almost surely larger than the true value of the parameter, which would be a - non desirable - statistical feature of the naive estimator.

Table 5. Simulation results: proportion of replications where the maximum likelihood estimator is larger than the true value of the parameter for the exponential model

Table 6. Simulation results: proportion of replications where the maximum likelihood estimator is larger than the true value of the parameter for the Weibull model

Table 7. Simulation results: proportion of replications where the maximum likelihood estimator is larger than the true value of the parameter for the log-logistic model

Application study

Table 8 presents the estimates of the parameters for the three models and both approaches. There was no problem of maximization. The naive estimates are always larger than the truncation-based estimates. From the simulation results, it might be thought that the naive estimator overestimates the true values of parameters λ and β, and that the size of the bias is related to the unknown probability p. Estimations of the parameters for the truncation-based approach make it possible to estimate p by calculating <a onClick="popup('http://www.biomedcentral.com/1471-2288/14/17/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/14/17/mathml/M36">View MathML</a>. However, estimates of p are different according to the model (Table 8). In particular, for the Weibull model, the estimate is large (<a onClick="popup('http://www.biomedcentral.com/1471-2288/14/17/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/14/17/mathml/M37">View MathML</a> = 0.98). The larger is <a onClick="popup('http://www.biomedcentral.com/1471-2288/14/17/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/14/17/mathml/M38">View MathML</a>, the closer are the naive and the truncation-based estimates.

Table 8. Parameter estimation and estimated mean time-to-onset for 64 cases of lymphoma that occurred after anti TNF- α treatment

Figure 2 shows the non-parametric maximum likelihood estimation of the conditional survival function, <a onClick="popup('http://www.biomedcentral.com/1471-2288/14/17/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/14/17/mathml/M45">View MathML</a>, and the parametric maximum likelihood estimation of the conditional, <a onClick="popup('http://www.biomedcentral.com/1471-2288/14/17/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/14/17/mathml/M46">View MathML</a>, and unconditional, <a onClick="popup('http://www.biomedcentral.com/1471-2288/14/17/mathml/M47','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/14/17/mathml/M47">View MathML</a>, survival functions for the truncation-based approach for these data. The estimations of the conditional survival functions are always closer to the non-parametric estimation than the estimations of the unconditional survival functions. The conditional and unconditional estimations of the Weibull survival functions are almost similar because the estimate of p is about 1. This figure shows that the estimation of the conditional Weibull survival function is closer to the non-parametric maximum likelihood estimation of the conditional survival function than the estimations of the conditional exponential and conditional log-logistic survival functions. Thus, Weibull could be a reasonable candidate model to describe the data.

thumbnailFigure 2. Right truncation-based estimations of time-to-onset of lymphoma that occurred after anti TNF-α treatment. Data include 64 cases. Three models are fitted: exponential, Weibull and log-logistic. Estimations of the conditional survival function (C), estimations of the unconditional survival function (U) and the non-parametric maximum likelihood estimation of the survival function (NPMLE) are displayed.

Figure 3 shows the parametric maximum likelihood estimation of the unconditional survival function for both approaches. The distance between both survivals, naive and truncation-based, decreases with the estimated probability <a onClick="popup('http://www.biomedcentral.com/1471-2288/14/17/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/14/17/mathml/M48">View MathML</a> (in the order: exponential, log-logistic and Weibull). Furthermore, the survival functions from the truncation-based estimates are always above the survival functions from the naive estimates, which is consistent with the naive estimator overestimating the true values of the parameters λ and β. Even for the Weibull model, i.e. the model with the largest <a onClick="popup('http://www.biomedcentral.com/1471-2288/14/17/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2288/14/17/mathml/M49">View MathML</a>, the estimated expected time-to-onset would be 135 weeks with the naive approach and 193 weeks with the truncation-based estimates, which corresponds to a markedly large gap (Table 8). For completeness, we also calculated the 95% simple bootstrap confidence intervals of the expected time (BCa method) [26,27] based on 5000 bootstrap samples, for the truncation-based approach. They do not include the naive estimated mean time, whatever the fitted model, and even though these confidence intervals are extremely wide.

thumbnailFigure 3. Naive and right truncation-based estimations of time-to-onset of lymphoma that occurred after anti TNF-α treatment. Data include 64 cases. Three models are fitted: exponential, Weibull and log-logistic. Estimations of the unconditional survival function for the naive approach (Naive) and for the truncation-based approach (TBE) are displayed.

Discussion and conclusions

In drug safety assessment, the temporal relationship between drug administration and time-to-onset is of utmost relevance. A better understanding of the underlying mechanism of the occurrence of an adverse effect is crucial, as it could allow the identification of particular groups of patients at risk and of particular risk time-windows in the course of a treatment and lead to preventing or diagnosing earlier the occurrence of adverse reactions. In this framework, the time-to-onset of an adverse drug reaction constitutes an essential feature to be analyzed. Its accurate estimation and modeling could help in understanding the mechanism of a drug’s action.

As rare adverse effects are not generally identified by cohort studies of exposed patients but from spontaneous reporting systems, we investigated with a simulation study the accuracy of estimates that can be obtained from these data in a parametric framework. As one can only estimate a conditional distribution function in a non-parametric setting, the non-parametric maximum likelihood estimator is of rather little interest for pharmacovigilance people. For a finite sample size, the simulations show that, whatever the approach, naive or truncation-based, the parametric maximum likelihood estimator may be positively biased and that this bias and the corresponding mean squared error increase when the theoretical probability p for the time-to-onset to fall within the observable values interval decreases. However, for a fixed value of p, the bias and the mean squared error are always larger when the right truncation is not considered than when it is, and the gap may be large. In addition, bias and mean squared error might in some instances (Weibull, log-logistic) be unacceptably large for the naive approach, even for a large value of p, while with a probability p of 0.8, or sometime even less, the TBE shows good performances. Asymptotically, the naive estimator may not be unbiased because the bias and the mean squared error seem to be constant with the sample size and the maximization is based on a misleading likelihood, while the bias and the mean squared error for the TBE decrease as the sample size increases. Therefore, even if the sample size is large, the gap between both estimators does not disappear and the truncation-based approach should be used.

The probability p plays an important role in the estimation of the distribution of the time-to-onset of adverse reaction for right-truncated data. Knowledge exists on a range of possible pharmacological mechanisms. It is thus possible to get a rough idea of the fraction of potentially missed cases (the adverse reactions of treated patients that have yet to occur) and then to decide on the relevance of the time of analysis. Spontaneous reports result from three processes: the occurrence case process, its diagnosis and the reporting process. It is well known that under-reporting is widespread, even for serious events. In addition, factors of under-reporting include the seriousness of the effect, the age of the patient and the novelty of the effect, but also time-related variables such as the length of marketing or the time since exposure [28-33]. In the approach proposed here, it is assumed that the under-reporting is uniform. Such a hypothesis might not always be acceptable. However, with long-term effects such as lymphoma and a homogeneous observation period within the marketing life of the product, non-stationarity of reporting is unlikely.

Problems of maximization may arise when right truncation is taken into account. The smaller is p, the more the iterative algorithm is likely to fail. Some papers mentioned the existence of a problem in the parametric likelihood maximization and explained that, because of right truncation, the likelihood may be flat and the maximum may be difficult to find [21,34-36].

For the 64 cases of lymphoma after anti TNF- α treatment, there was no problem of convergence of the iterative algorithm. Both estimates, naive and truncation-based, were available for each fitted model. From the truncation-based estimates, it is possible to estimate p. Here it ranges from 0.98 (Weibull) to 0.60 (exponential). Since this probability is unknown, the non-parametric maximum likelihood estimation estimates only the distribution function conditional on the time-to-event being less than the maximum observed truncation time. However, although conditional, the non-parametric estimate is a reference that provides an idea of how the data fit a given model. We followed the graphical procedure for checking goodness-of-fit for right-truncated data suggested by Lawless (2003) that is based on the non-parametric maximum likelihood estimator and consists in plotting the conditional fitted parametric survivals together with the non-parametric estimation [36]. Here, the conditional Weibull survival function seems the closest to the non-parametric estimation. This finding underlines the interest for developing goodness-of-fit tests adapted to right-truncated data. While only three families of distributions were considered for the present simulation study, other families could be explored such as the gamma or the log-normal families or mixture models. For instance, in more complex situations, the treatment might be a combination of drugs, each of them inducing the effect but in a different time window. In that case, the hazard function may vary several times and a family of more complex distributions could be of greater interest. Additionally, we chose to consider the truncation times as deterministic, which is equivalent to working on conditional distributions for the likelihood. However, another possible approach is to consider the truncation time as a random variable and to study the random pair (X,T) where X is the survival time and T is the truncation time [37-39].

Finally, improvement of time-to-onset distribution assessment could make it possible to compare two drug profiles or more generally to assess risk factors with regression models.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

FL, JYD and PTB conceived and designed the work. FL implemented the simulations, performed data analysis and wrote the initial draft of the manuscript. HT and FH made the extraction of the data from the national pharmacovigilance database. All authors contributed to the interpretation of the results of the data analysis. All authors reviewed and revised the draft version of the manuscript. All authors read and approved the final version of the manuscript.

Acknowledgements

This work was supported by the Fondation ARC (fellowship DOC20121206119 to Fanny Leroy).

References

  1. Fourrier A, Bégaud B, Alpérovitch A, Verdier-Taillefer M-H, Decker N, Imbs J-L, Touzé E: Hepatitis B vaccine and first episodes of central nervous system demyelinating disorders: a comparison between reported and expected number of cases.

    Br J Clin Pharmacol 2001, 51(5):489-490. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Tubert P, Bégaud B, Haramburu F, Péré JC: Spontaneous reporting: how many cases are required to trigger a warning?

    Br J Clin Pharmacol 1991, 32(4):407-408. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  3. Moore N, Kreft-Jais C, Haramburu F, Noblet C, Andrejak M, Ollagnier M, Bégaud B: Reports of hypoglycaemia associated with the use of ACE inhibitors and other drugs: a case/non-case study in the French pharmacovigilance system database.

    Br J Clin Pharmacol 1997, 44(5):513-518. PubMed Abstract | PubMed Central Full Text OpenURL

  4. Tubert-Bitter P, Bégaud B, Moride Y, Chaslerie A, Haramburu F: Comparing the toxicity of two drugs in the framework of spontaneous reporting: a confidence interval approach.

    J Clin Epidemiol 1996, 49(1):121-123. PubMed Abstract | Publisher Full Text OpenURL

  5. van der Heijden PG, van Puijenbroek EP, van Buuren S, van der Hofstede JW: On the assessment of adverse drug reactions from spontaneous reporting systems: the influence of under-reporting on odds ratios.

    Stat Med 2002, 21(14):2027-2044. PubMed Abstract | Publisher Full Text OpenURL

  6. Bate A, Lindquist M, Edwards IR, Olsson S, Orre R, Lansner A, De Freitas RM: A bayesian neural network method for adverse drug reaction signal generation.

    Eur J Clin Pharmacol 1998, 54(4):315-321. PubMed Abstract | Publisher Full Text OpenURL

  7. DuMouchel W: Bayesian data mining in large frequency tables, with an application to the FDA spontaneous reporting system.

    Am Stat 1999, 53(3):177-190. OpenURL

  8. Szarfman A, Machado SG, O’Neill RT: Use of screening algorithms and computer systems to efficiently signal higher-than-expected combinations of drugs and events in the US FDA’s spontaneous reports database.

    Drug Saf 2002, 25(6):381-392. PubMed Abstract | Publisher Full Text OpenURL

  9. Evans SJW, Waller PC, Davis S: Use of proportional reporting ratios (PRRs) for signal generation from spontaneous adverse drug reaction reports.

    Pharmacoepidemiol Drug Saf 2001, 10(6):483-486. PubMed Abstract | Publisher Full Text OpenURL

  10. Ahmed I, Haramburu F, Fourrier-Réglat A, Thiessard F, Kreft-Jais C, Bégaud B, Tubert-Bitter P, Miremont-Salamé G: Bayesian pharmacovigilance signal detection methods revisited in a multiple comparison setting.

    Stat Med 2009, 28(13):1774-1792. PubMed Abstract | Publisher Full Text OpenURL

  11. Ahmed I, Dalmasso C, Haramburu F, Thiessard F, Broët P, Tubert-Bitter P: False discovery rate estimation for frequentist pharmacovigilance signal detection methods.

    Biometrics 2010, 66(1):301-309. PubMed Abstract | Publisher Full Text OpenURL

  12. Roux E, Thiessard F, Fourrier A, Bégaud B, Tubert-Bitter P: Evaluation of statistical association measures for the automatic signal generation in pharmacovigilance.

    IEEE Trans Inf Technol Biomed 2005, 9(4):518-527. PubMed Abstract OpenURL

  13. Ahmed I, Thiessard F, Bégaud B, Tubert-Bitter P, Miremont-Salamé G: Pharmacovigilance data mining with methods based on false discovery rates: a comparative simulation study.

    Clin Pharmacol Ther 2010, 88(4):492-498. PubMed Abstract | Publisher Full Text OpenURL

  14. Bate A, Evans SJW: Quantitative signal detection using spontaneous ADR reporting.

    Pharmacoepidemiol Drug Saf 2009, 18(6):427-436. PubMed Abstract | Publisher Full Text OpenURL

  15. Alvarez Y, Hidalgo A, Maignen F, Slattery J: Validation of statistical signal detection procedures in eudravigilance post-authorization data: a retrospective evaluation of the potential for earlier signalling.

    Drug Saf 2010, 33(6):475-487. PubMed Abstract | Publisher Full Text OpenURL

  16. Hochberg AM, Hauben M: Time-to-signal comparison for drug safety data-mining algorithms vs. traditional signaling criteria.

    Clin Pharmacol Ther 2009, 85(6):600-606. PubMed Abstract | Publisher Full Text OpenURL

  17. Ahmed I, Thiessard F, Haramburu F, Kreft-Jais C, Bégaud B, Tubert-Bitter P, Miremont-Salamé G: Early detection of pharmacovigilance signals with automated methods based on false discovery rates: a comparative study.

    Drug Saf 2012, 35(6):495-506. PubMed Abstract | Publisher Full Text OpenURL

  18. Maignen F, Hauben M, Tsintis P: Modelling the time to onset of adverse reactions with parametric survival distributions.

    Drug Saf 2010, 33(5):417-434. PubMed Abstract | Publisher Full Text OpenURL

  19. Van Holle L, Zeinoun Z, Bauchau V, Verstraeten T: Using time-to-onset for detecting safety signals in spontaneous reports of adverse events following immunization: a proof of concept study.

    Pharmacoepidemiol Drug Saf 2012, 21(6):603-610. PubMed Abstract | Publisher Full Text OpenURL

  20. Cornelius VR, Sauzet O, Evans SJW: A signal detection method to detect adverse drug reactions using a parametric time-to-event model in simulated cohort data.

    Drug Saf 2012, 35(7):599-610. PubMed Abstract | Publisher Full Text OpenURL

  21. Lagakos SW, Barraj LM, De Gruttola V: Nonparametric analysis of truncated survival data, with application to aids.

    Biometrika 1988, 75(3):515-523. Publisher Full Text OpenURL

  22. Kalbfleisch JD, Lawless JF: Regression models for right truncated data with applications to AIDS incubation times and reporting lags.

    Stat Sin 1991, 1:19-32. OpenURL

  23. Bégaud B, Miremont G, Péré JC: Estimation of the denominator in spontaneous reporting. In Methodological Approaches in Pharmacoepidemiology: Application to Spontaneous Reporting. Amsterdam: Elsevier; 1993:51-70. OpenURL

  24. R Development Core Team R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. Vienna, Austria;

    [http://cran.r-project.org/ webcite]

    OpenURL

  25. Théophile H, Schaeverbeke T, Abouelfath A, Kahn V, Haramburu F, Bégaud B, Miremont-Salamé G: Sources of information on lymphoma associated with anti-tumour necrosis factor agents.

    Drug Saf 2011, 34(7):577-585. PubMed Abstract | Publisher Full Text OpenURL

  26. Efron B, Tibshirani RJ: An Introduction to the Bootstrap. New York: Chapman & Hall; 1993. OpenURL

  27. Gross ST, Lai TL: Bootstrap methods for truncated and censored data.

    Stat Sin 1996, 6:509-530. OpenURL

  28. Weber JCP: Mathematical models in adverse drug reaction assessment. In Iatrogenic Diseases. 3rd Ed. Edited by Arcy PF, Griffin JP. Oxford: Oxford University Press; 1986. OpenURL

  29. Tubert-Bitter P, Haramburu F, Bégaud B, Chaslerie A, Abraham E, Hagry C: Spontaneous reporting of adverse drug reactions: who reports and what?

    Pharmacoepidemiol Drug Saf 1998, 7(5):323-329. PubMed Abstract | Publisher Full Text OpenURL

  30. Haramburu F Bégaud B, Moride Y: Temporal trends in spontaneous reporting of unlabelled adverse drug reactions.

    Br J Clin Pharmacol 1997, 44(3):299-301. PubMed Abstract | PubMed Central Full Text OpenURL

  31. Moride Y, Haramburu F, Requejo AA, Bégaud B: Under-reporting of adverse drug reactions in general practice.

    Br J of Clin Pharmacol 1997, 43(2):177-181. OpenURL

  32. Bégaud B, Martin K, Haramburu F, Moore N: Rates of spontaneous reporting of adverse drug reactions in France (letter).

    JAMA 2002, 288(13):1588-1588. PubMed Abstract | Publisher Full Text OpenURL

  33. Tubert P, Bégaud B, Haramburu F, Lellouch J, Péré J-C: Power and weakness of spontaneous reporting: a probabilistic approach.

    J Clin Epidemiol 1992, 45(3):283-286. PubMed Abstract | Publisher Full Text OpenURL

  34. Kalbfleisch JD, Lawless JF: Inference based on retrospective ascertainment: an analysis of the data on transfusion-related AIDS.

    J Am Stat Assoc 1989, 84(406):360-372. Publisher Full Text OpenURL

  35. Colton T: Biased Sampling of Cohorts in Epidemiology. In Encyclopedia of Biostatistics, Vol 1. Edited by Armitage P, Colton T. Chichester: Wiley; 1998:338-350. OpenURL

  36. Lawless JF: Statistical Models and Methods for Lifetime Data, 2nd Ed. Hokoben, New Jersey: Wiley; 2003. OpenURL

  37. Keiding N: Nonparametric estimation under truncation. In Encyclopedia of Statistical Sciences, Vol 14, 2nd Ed. Hokoben, New Jersey: Wiley; 2006:8775-8777. OpenURL

  38. Gürler Ü: Bivariate estimation with right-truncated data.

    J Am Stat Assoc 1996, 91(435):1152-1165. OpenURL

  39. Gross ST, Huber-Carol C: Regression models for truncated survival data.

    Scandinavian J Stat 1992, 193-213. OpenURL

Pre-publication history

The pre-publication history for this paper can be accessed here:

http://www.biomedcentral.com/1471-2288/14/17/prepub