calcWINS(x = COVID19b, AVAL = "GROUP", TRTP = "TRTP", ref = "Placebo")
## Biased vs unbiased
n0 <- 10; n1 <- 20; p0 <- 0.2; p1 <- 0.5; x <- 1:20; delta <- 0.5
WP0 <- (p1 - p0)/2 + 0.5
DAT <- NULL
for(i in x){
dat <- data.frame(AVAL = c(rbinom(n1, size = 1, p1), rbinom(n0, size = 1, p0)),
TRTP = c(rep("A", n1), rep("P", n0)))
CL1 <- calcWINS(x = dat, AVAL = "AVAL", TRTP = "TRTP", ref = "P")$WP
CL1$Type <- "biased"
CL2 <- calcWINS(x = dat, AVAL = "AVAL", TRTP = "TRTP",
ref = "P", SE_WP_Type = "unbiased")$WP
CL2$Type <- "unbiased"
DAT <- rbind(DAT, CL1, CL2)
}
WP <- DAT$WP[DAT$Type == "unbiased"]
plot(x, WP, pch = 19, xlab = "Simulations", ylab = "Win Probability", ylim = c(0., 1.))
points(x + delta, WP, pch = 19)
arrows(x, DAT$LCL[DAT$Type == "unbiased"],
x, DAT$UCL[DAT$Type == "unbiased"], angle = 90, code = 3, length = 0.05, "green")
arrows(x + delta, DAT$LCL[DAT$Type == "biased"],
x + delta, DAT$UCL[DAT$Type == "biased"], angle = 90, code = 3, length = 0.05, col = "red")
abline(h = WP0, col = "blue", lty = 3)
legend("bottomleft", legend = c("True WP", "Biased", "Unbiased"),
col = c(4, 2, 3), lty = c(3, 1, 1 ), cex = 0.5)
Run the code above in your browser using DataLab