utility_bias {drugdevelopR} | R Documentation |
Utility function for bias adjustment programs with time-to-event outcomes.
Description
The utility function calculates the expected utility of our drug development program and is given as gains minus costs and depends on the parameters and the expected probability of a successful program.
The utility is in a further step maximized by the optimal_bias()
function.
Usage
utility_L(
d2,
HRgo,
Adj,
w,
hr1,
hr2,
id1,
id2,
alpha,
beta,
xi2,
xi3,
c2,
c3,
c02,
c03,
K,
N,
S,
steps1,
stepm1,
stepl1,
b1,
b2,
b3,
fixed
)
utility_L2(
d2,
HRgo,
Adj,
w,
hr1,
hr2,
id1,
id2,
alpha,
beta,
xi2,
xi3,
c2,
c3,
c02,
c03,
K,
N,
S,
steps1,
stepm1,
stepl1,
b1,
b2,
b3,
fixed
)
utility_R(
d2,
HRgo,
Adj,
w,
hr1,
hr2,
id1,
id2,
alpha,
beta,
xi2,
xi3,
c2,
c3,
c02,
c03,
K,
N,
S,
steps1,
stepm1,
stepl1,
b1,
b2,
b3,
fixed
)
utility_R2(
d2,
HRgo,
Adj,
w,
hr1,
hr2,
id1,
id2,
alpha,
beta,
xi2,
xi3,
c2,
c3,
c02,
c03,
K,
N,
S,
steps1,
stepm1,
stepl1,
b1,
b2,
b3,
fixed
)
Arguments
d2 |
total events for phase II; must be even number |
HRgo |
threshold value for the go/no-go decision rule |
Adj |
adjustment parameter |
w |
weight for mixture prior distribution |
hr1 |
first assumed true treatment effect on HR scale for prior distribution |
hr2 |
second assumed true treatment effect on HR scale for prior distribution |
id1 |
amount of information for |
id2 |
amount of information for |
alpha |
significance level |
beta |
|
xi2 |
event rate for phase II |
xi3 |
event rate for phase III |
c2 |
variable per-patient cost for phase II |
c3 |
variable per-patient cost for phase III |
c02 |
fixed cost for phase II |
c03 |
fixed cost for phase III |
K |
constraint on the costs of the program, default: Inf, e.g. no constraint |
N |
constraint on the total expected sample size of the program, default: Inf, e.g. no constraint |
S |
constraint on the expected probability of a successful program, default: -Inf, e.g. no constraint |
steps1 |
lower boundary for effect size category |
stepm1 |
lower boundary for effect size category |
stepl1 |
lower boundary for effect size category |
b1 |
expected gain for effect size category |
b2 |
expected gain for effect size category |
b3 |
expected gain for effect size category |
fixed |
choose if true treatment effects are fixed or random, if TRUE |
Value
The output of the functions utility_L()
, utility_L2()
, utility_R()
and utility_R2()
is the expected utility of the program.
Examples
res <- utility_L(d2 = 50, HRgo = 0.8, Adj = 0.4, w = 0.3,
hr1 = 0.69, hr2 = 0.81,
id1 = 280, id2 = 420, xi2 = 0.7, xi3 = 0.7,
alpha = 0.025, beta = 0.1,
c2 = 0.75, c3 = 1, c02 = 100, c03 = 150,
K = Inf, N = Inf, S = -Inf,
steps1 = 1, stepm1 = 0.95, stepl1 = 0.85,
b1 = 1000, b2 = 2000, b3 = 3000,
fixed = TRUE)
res <- utility_L2(d2 = 50, HRgo = 0.8, Adj = 0.4, w = 0.3,
hr1 = 0.69, hr2 = 0.81,
id1 = 280, id2 = 420, xi2 = 0.7, xi3 = 0.7,
alpha = 0.025, beta = 0.1,
c2 = 0.75, c3 = 1, c02 = 100, c03 = 150,
K = Inf, N = Inf, S = -Inf,
steps1 = 1, stepm1 = 0.95, stepl1 = 0.85,
b1 = 1000, b2 = 2000, b3 = 3000,
fixed = TRUE)
res <- utility_R(d2 = 50, HRgo = 0.8, Adj = 0.9, w = 0.3,
hr1 = 0.69, hr2 = 0.81,
id1 = 280, id2 = 420, xi2 = 0.7, xi3 = 0.7,
alpha = 0.025, beta = 0.1,
c2 = 0.75, c3 = 1, c02 = 100, c03 = 150,
K = Inf, N = Inf, S = -Inf,
steps1 = 1, stepm1 = 0.95, stepl1 = 0.85,
b1 = 1000, b2 = 2000, b3 = 3000,
fixed = TRUE)
res <- utility_R2(d2 = 50, HRgo = 0.8, Adj = 0.9, w = 0.3,
hr1 = 0.69, hr2 = 0.81,
id1 = 280, id2 = 420, xi2 = 0.7, xi3 = 0.7,
alpha = 0.025, beta = 0.1,
c2 = 0.75, c3 = 1, c02 = 100, c03 = 150,
K = Inf, N = Inf, S = -Inf,
steps1 = 1, stepm1 = 0.95, stepl1 = 0.85,
b1 = 1000, b2 = 2000, b3 = 3000,
fixed = TRUE)