design <- TwoStageDesign(
n1 = 25,
c1f = 0,
c1e = 2.5,
n2 = 50,
c2 = 1.96,
order = 7L
)
prior <- PointMassPrior(.3, 1)
# conditional
cp <- ConditionalPower(Normal(), prior)
expected(cp, Normal(), prior)
evaluate(cp, design, x1 = .5)
# unconditional
power <- Power(Normal(), prior)
evaluate(power, design)
evaluate(power, design, optimization = TRUE) # use non-adaptive quadrature
Run the code above in your browser using DataLab