## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(collapse = TRUE, comment = "#>")
library(OptOTrials)

## -----------------------------------------------------------------------------
p1 <- c(0.075, 0.182, 0.319, 0.243, 0.015, 0.166)
p2 <- p2_fun(p1, log(3.06))
Proportional_odds_assumption(p1, p2)

## -----------------------------------------------------------------------------
q1 <- c(1/3, 1/3, 1/3)
q2 <- c(1/2, 1/3, 1/6)
res <- tryCatch(Proportional_odds_assumption(q1, q2), warning = function(w) NA)
res

## -----------------------------------------------------------------------------
d <- rule(alpha = 0.05, beta = 0.2, p1 = p1, p2 = p2,
          test = "M", stopping = "F", criterion = 1)
d

## -----------------------------------------------------------------------------
d$n1
d$t1f
d$effect

## -----------------------------------------------------------------------------
o <- op(d, nsim = 2000, seed = 1234)
o

## -----------------------------------------------------------------------------
design_table(0.05, 0.2, p1, p2, tests = c("S", "M"), criteria = c(1, 3),
             stopping = "F", nsim = 1000)

## -----------------------------------------------------------------------------
d2 <- rule(0.05, 0.2, p1, p2, test = "M", stopping = "F", criterion = 2)
d2$substituted