#*********sepsis.fix design:

sepsis.fix <- seqDesign(prob.model="proportions",arms=2,size=0.025,power="calculate",
null.hypothesis=c(0.30,0.30),alt.hypothesis=c(0.25,0.30),test.type="less",sample.size=1700)


#*** Adding 4 OBF interim analyses:
sepsis.obf <- update(sepsis.fix, nbr.analyses=4)
update(sepsis.obf,display.scale="S")
seqInference(sepsis.obf)
seqOC(sepsis.obf,theta=c(-0.05,-0.07,-0.08))

#*** Changing early conservatism for futility:
sepsis.obfPoc <- update(sepsis.obf, P=c(1,0.5))
update(sepsis.obf,display.scale="Z")
seqInference(sepsis.obfPoc)
seqOC(sepsis.obfPoc,theta=c(-0.05,-0.07,-0.08))

#*** Graphical comparison of designs:

seqPlotBoundary(sepsis.obf,sepsis.obfPoc)
seqPlotPower(sepsis.obf,sepsis.obfPoc,reference=T)
seqPlotInference(sepsis.obf)
seqPlotInference(sepsis.obfPoc)