The function SampleSize.CondPoisson obtains the required sample size (length of surveillance) needed
to guarantee a desired statistical power for a pre-specified relative risk, when doing continuous sequential CMaxSPRT, using a Wald-type upper boundary,
which is flat with respect to the likelihood ratio function.
SampleSize.CondPoisson(cc,D=0,M=1,alpha=0.05,power=0.9,RR=2,
Tailed="upper",NegBinApprox= TRUE)A table containing the main performance measures associated to the required samples sizes, expressed in the scale of the number of events in the surveillance period, for each combination of RR and power.
The total number of observed adverse events in the historical data. There is no default.
The minium number for the ratio \(P_k/V\) before the null hypothesis can be rejected. The default value is \(D=0\). A delayed start with \(D>0\) is recommended to avoid signaling very early on such that very little information would be available to judge whether the signal is more likely to be a true signal or chance finding.
The minimum number of events needed before the null hypothesis can be rejected. The default value is M=1. A delayed start with \(M>1\) is recommended to avoid signaling very early on such that very little information would be available to judge whether the signal is more likely to be a true signal or chance finding.
The significance level. It must be in the range (0,0.5]. The default value is alpha=0.05.
The target vector of overall statistical powers to detect an increased relative risk (RR). The default value is power=0.90.
The target vector of relative risks to be detected with the requested vector of statistical powers. The default value is RR=2.
Tailed="upper" (default) for H0:RR<=1, and Tailed="lower" for H0:RR>=1 or Tailed="two" for H0:RR=1.
NegBinApprox=TRUE (default) for a fast calculation based on the negative binomial approximation. See details.
Ivair Ramos Silva, Lingling Li
Development of the SampleSize.CondPoisson function was funded by:
- National Institute of General Medical Sciences, NIH, USA, through grant number R01GM108999 (v2.0.1, v2.0.2).
- Foundation for Research Support of Minas Gerais State (FAPEMIG), MG, Brazil, through the grant Demanda Universal.
When using the CMaxSPRT (Li and Kulldorff, 2010) and the CV.CondPoisson function to conduct continuous sequential analysis
for Poisson data and limited historical data, the null hypothesis is rejected when the log likelihood ratio
exceeds the predetermined critical value calculated by CV.CondPoisson. The sequential analysis ends
without rejecting the null hypothesis when a predetermined upper limit on the sample size is reached, expressed
either in terms of the ratio of the cumulative person-time in the surveillance population divided by the total
cumulative person-time in historical data (StopType="Tal"), or in terms of the observed number of events in the
surveillance population (StopType="Cases"). For example, the sequential analysis may end as soon as the sample size
is such that the cumulative person-time in the surveillance population is twice the cumulative person-time in historical data,
or there are 50 observed events in the surveillance population. The function SampleSize.CondPoisson calculates what
the upper limit on the sample size (length of surveillance) that is required for the continuous CMaxSPRT to achieve the desired
statistical power for a pre-specified relative risk RR. It outputs the upper limit on sample size for both definitions of the
surveillance length, one expressed in terms of the ratio of the cumulative person-time in the surveillance population divided by
the total cumulative person-time in historical data (\(T\)), and the other one expressed in terms of the observed number of
events in the surveillance population (\(K\)). To save computing time, the negative binomial approximation proposed by
Silva and Montalban (2023) is used.
Li L, Kulldorff M. (2010). A conditional maximized sequential probability ratio test for Pharmacovigilance. Statistics in Medicine, 29 (2), 284--295.
Silva IR, Li L, Kulldorff M. (2019). Exact Conditional Sequential Testing for Poisson Data. Sequential Analysis, in press.
Silva IR., Lopes LM., Dias P., Yih WK. (2019). Alpha Spending for Historical Versus Surveillance Poisson Data With CMaxSPRT. Statistics in Medicine, DOI: 10.1002/sim.8097, 1--13.
Silva IR, Montalban, J. (2023), The person-time ratio distribution for the exact monitoring of adverse events: Historical vs surveillance Poisson data, Statistics in Medicine, 42(18), 3283--3301.
CV.CondPoisson: calculating the critical value for continuous CMaxSPRT.
Performance.CondPoisson: calculating the statistical power, expected time to signal and expected time of analysis for continuous CMaxSPRT.
# Sample size required to obtain a power of 90%, for a relative risk of 2,
# no delay for starting the surveilla# nce (D=0), under an alpha level of 5%,
# with 50 events in the historical data.
# res<- SampleSize.CondPoisson(cc=50,D=0,M=1,alpha=0.05,power=0.9,RR=2)
Run the code above in your browser using DataLab