Return two-sample normal-gamma predictive probability
return.ts.ng.int.req.df(
mu.0.t = 0,
n.0.t = 1e-04,
alpha.0.t = 0.25,
beta.0.t = 1,
mu.0.c = 0,
n.0.c = 1e-04,
alpha.0.c = 0.25,
beta.0.c = 1,
xbar.t = c(1.9, 2, 2.1, 2.05),
s.t = c(2, 2.1, 1.9, 2.04),
n.t = c(10, 20, 30, 40),
xbar.c = c(1, 1.1, 1.5, 1.25),
s.c = c(1.9, 2, 2.5, 2.25),
n.c = c(10, 20, 30, 40),
Delta.lrv = 1.25,
Delta.tv = 1.75,
tau.tv = 0.1,
tau.lrv = 0.8,
tau.ng = 0.65,
xbar_ng = NULL,
xbar_go = NULL,
go.thresh = 0.8,
ng.thresh = 0.8,
n.MC = 1000
)
A dataframe is returned
prior mean for treatment group
prior effective sample size parameter for treatment group
prior alpha parameter for treatment group
prior beta parameter for treatment group
prior mean for control group
prior effective sample size parameter for control group
prior alpha parameter for control group
prior beta parameter for control group
treatment mean
treatment sd
treatment sample size
control mean
control sd
control sample size
TPP Lower Reference Value aka Min TPP
TPP Target Value aka Base TPP
threshold associated with Base TPP
threshold associated with Min TPP
threshold associated with No-Go
Leave NULL to determine what is required or supply a value
Leave NULL to determine what is required or supply a value
If the predictive probability that study will conclude as 'Go' is larger than this threshold: Declare 'Interim go'.
If the predictive probability that study will conclude as 'No-Go' is larger than this threshold: Declare 'Interim no-go'
Monte Carlo simulation size
holdit <- return.ts.ng.int.req.df()
head(holdit)
Run the code above in your browser using DataLab