design <- OneStageDesign(50, 1.96)
cp <- ConditionalPower(Normal(), PointMassPrior(0.4, 1))
pow <- Power(Normal(), PointMassPrior(0.4, 1))
# unconditional power constraint
constraint1 <- pow >= 0.8
evaluate(constraint1, design)
# conditional power constraint
constraint2 <- cp >= 0.7
evaluate(constraint2, design, .5)
constraint3 <- 0.7 <= cp # same as constraint2
evaluate(constraint3, design, .5)
Run the code above in your browser using DataLab