# Example 1
p0 <- c(0.1,0.1,0.1,0.1)
p1 <- c(0.1,0.1,0.1,0.1)
p <- rbind(p0,p1)
p_tox <- c(0.1,0.4)
# consider 1 candidate drugs with 4 subgroups
result <- platform_midas2(seed=20,p,p_tox,n.burnin=1000,
n.iter=2000,n.thin=2,C_T=0.85,C_E1=0.15,C_E2=0.999)
result
# \donttest{
# Example 2
p0 <- c(0.05,0.10,0.05,0.10)
p1 <- c(0.24,0.40,0.12,0.22)
p2 <- c(0.24,0.40,0.12,0.22)
p3 <- c(0.12,0.22,0.05,0.10)
p4 <- c(0.24,0.40,0.12,0.22)
p5 <- c(0.28,0.45,0.12,0.22)
p6 <- c(0.24,0.40,0.12,0.22)
p7 <- c(0.12,0.22,0.05,0.10)
p <- rbind(p0, p1, p2, p3, p4, p5, p6, p7)
p_tox <- c(0.10,0.10,0.10,0.10,0.10,0.10,0.15,0.20)
# consider 7 candidate drugs with 4 subgroups
result <- platform_midas2(seed=12,p,p_tox,n.burnin=10000,
n.iter=20000,n.thin=2,C_T=0.85,C_E1=0.15,C_E2=0.999)
result
# }
Run the code above in your browser using DataLab