PoS {LongCART} | R Documentation |
Probability of trial ans clinical success for a prospective trial using normal-normal approximation
Description
This function can be used to determine probability of trial success and clinical success based on the prior distribution for each of continuous, binary and time-to-event endpoints. The calculation is carried out assuming normal distribution for estimated parameter and normal prior distribution.
Usage
PoS(type, nsamples, null.value = NULL, alternative = "greater",
N = NULL, D = NULL, a = 1,
succ.crit = "trial", Z.crit.final = 1.96, alpha.final = 0.025,
clin.succ.threshold = NULL, se.exp = NULL,
meandiff.prior = NULL, mean.prior = NULL, sd.prior = NULL,
propdiff.prior = NULL, prop.prior = NULL, hr.prior = NULL, D.prior = NULL)
Arguments
type |
Type of the endpoint. It could be |
nsamples |
Number of samples. For continuous and binary case, it can be 1 or 2. For survival endpoint, it can be only 2. |
null.value |
The specified value under null hypothesis. Default is 0 for continuous and binomial case and 1 for survival case. |
alternative |
Direction of alternate hypothesis. Can be "greater" or "less". Default is "less" for test of HR and "greater" otherwise. |
N |
Total sample size at final analysis. Cannot be missing for continuous and binary endpoint. |
D |
Total number of events at final analysis. Cannot be missing for survival endpoint. |
a |
Allocation ratio in two sample case. |
succ.crit |
Specify "trial" for trial success (i.e., null hypothesis is rejected at final analysis) or "clinical" for clinical success (i.e., estimated value at the final analysis is greater than clinically meaningful value as specified under |
Z.crit.final |
The rejection boundary at final analysis in Z-value scale. Either |
alpha.final |
The rejection boundary at final analysis in alpha (1-sided) scale (e.g., 0.025). Either |
clin.succ.threshold |
Clinically meaningful value. Required when |
se.exp |
Expected standard error to be observed in the study. Must be specified in continuous case and two-sample binary case. |
meandiff.prior |
Mean value of prior distribution for mean difference. Relevant for two-sample continuous case. |
mean.prior |
Mean value of prior distribution for mean. Relevant for one-sample continuous case. |
sd.prior |
Standard deviation of prior distribution for mean difference (2-sample continuous case) or mean (1-sample continuous case) or prop (2-sample binary case) or difference of proportion (1-sample binary case) or log(HR) (2 sample survival case). |
propdiff.prior |
Mean value of prior distribution for difference in proportion. Relevant for two-sample binomial case. |
prop.prior |
Mean value of prior distribution for proportion. Relevant for one-sample binomial case. |
hr.prior |
Mean value of prior distribution for hazards ratio (HR). Relevant for two-sample survival case. |
D.prior |
Ignored if |
Details
This function can be used to determine probability of success (PoS) for a prospective trial for each of continuous (one-sample or two-samples), binary (one-sample or two-samples) and time-to-event endpoints (two-samples). The PoS is calculated based on the prior distribution and expected standard error in the estimate in trial. The calculation PoS is carried out assuming normal distribution for estimated parameter and normal prior distribution. This function can be used to determine clinical success (succ.crit="clinical"
) and trial success (succ.crit="trial"
). For clinical success, clin.succ.threshold
must be specified. For trial success, Z.crit.final or alpha.final
must be specified.
Author(s)
Madan Gopal Kundu <madan_g.kundu@yahoo.com>
References
Kundu, M. G., Samanta, S., and Mondal, S. (2021). An introduction to the determination of the probability of a successful trial: Frequentist and Bayesian approaches. arXiv preprint arXiv:2102.13550.
See Also
succ_ia_betabinom_one
, succ_ia_betabinom_two
, succ_ia
Examples
#--- Example 1
PoS(type="cont", nsamples=2, null.value=-0.05, alternative="greater",
N=1552, a=1,
succ.crit="trial", Z.crit.final=1.97,
se.exp=0.12*sqrt(1/776 + 1/776),
meandiff.prior=0, sd.prior=0.02)
#--- Example 2
PoS(type="bin", nsamples=2, null.value=0, alternative="greater",
N=210, a=2,
succ.crit="trial", Z.crit.final=2.012,
se.exp=0.5*sqrt(1/140 + 1/70),
propdiff.prior=0.20, sd.prior=sqrt(0.06))
PoS(type="bin", nsamples=2, null.value=0, alternative="greater",
N=210, a=2,
succ.crit="clinical", clin.succ.threshold =0.15,
se.exp=0.5*sqrt(1/140 + 1/70),
propdiff.prior=0.20, sd.prior=sqrt(0.06))
#--- Example 4
PoS(type="surv", nsamples=2, null.value=1, alternative="less",
D=441,
succ.crit="trial", Z.crit.final=1.96,
hr.prior=0.71, D.prior=133)
PoS(type="surv", nsamples=2, null.value=1, alternative="less",
D=441,
succ.crit="clinical", clin.succ.threshold =0.8,
hr.prior=0.71, D.prior=133)