pgreater_NIX {RARtrials} | R Documentation |
Calculate the Futility Stopping Probability for Continuous Endpoint with Unknown Variances Using a Normal-Inverse-Chi-Squared Distribution
Description
Calculate the futility stopping probability in Bayesian response-adaptive randomization with
a control group using the Thall \&
Wathen method for continuous outcomes with unknown variances. The prior distributions
follow Normal-Inverse-Chi-Squared (NIX) distributions and can be specified individually for each treatment group.
Usage
pgreater_NIX(par1, par2, delta = 0, side, ...)
Arguments
par1 |
current parameters including mu, kappa, nu, sigsq of a Normal-Inverse-Chi-Squared distribution from the control group. |
par2 |
current parameters including mu, kappa, nu, sigsq of a Normal-Inverse-Chi-Squared distribution from the compared treatment group. |
delta |
pre-specified minimal effect size expected to be observed between the control group and the compared treatment group. |
side |
direction of a one-sided test, with values 'upper' or 'lower'. |
... |
additional arguments to be passed to stats::integrate() (such as rel.tol) from this function. |
Details
This function calculates the results of Pr(\mu_k>\mu_{control}+\delta|data)
for side
equals to
'upper' and the results of Pr(\mu_{control}>\mu_k+\delta|data)
for side
equals to 'lower'.
The result indicates the posterior probability of stopping a treatment group due to futility around 1\%
in Bayesian
response-adaptive randomization with a control arm using Thall \&
Wathen method, with accumulated results
during the conduct of trials. Parameters used in a Normal-Inverse-Gamma ((\mu,\sigma^2) \sim NIG(mean=m,variance=V \times \sigma^2,shape=a,rate=b)
)
distribution should be converted to parameters equivalent in a Normal-Inverse-Chi-Squared
((\mu,\sigma^2) \sim NIX(mean=\mu,effective sample size=\kappa,degrees of freedom=\nu,variance=\sigma^2/\kappa)
)
distribution using convert_gamma_to_chisq
before applying this function.
Value
a posterior probability of Pr(\mu_k>\mu_{control}+\delta|data)
with side
equals to 'upper';
a posterior probability of Pr(\mu_{control}>\mu_k+\delta|data)
with side
equals to 'lower'.
References
Wathen J, Thall P (2017). “A simulation study of outcome adaptive randomization in multi-arm clinical trials.” Clinical Trials, 14, 174077451769230. doi:10.1177/1740774517692302. Murphy K (2007). “Conjugate Bayesian analysis of the Gaussian distribution.” University of British Columbia. https://www.cs.ubc.ca/~murphyk/Papers/bayesGauss.pdf.
Examples
para<-list(V=1/2,a=0.5,m=9.1/100,b=0.00002)
par<-convert_gamma_to_chisq(para)
set.seed(123451)
y1<-rnorm(100,0.091,0.009)
par1<-update_par_nichisq(y1, par)
set.seed(123452)
y2<-rnorm(90,0.09,0.009)
par2<-update_par_nichisq(y2, par)
pgreater_NIX(par1=par1,par2=par2, side='upper')
pgreater_NIX(par1=par1,par2=par2, side='lower')