two_arm_two_stage_full

Martin

25/02/2022

library(pwb)
library(ggplot2)
library(forcats)
theta1.vec <- seq(0.1, 0.8, 0.1)
iw.results <- vector("list", length(theta1.vec))
mc.error <- numeric(length(theta1.vec))

for(i in 1:length(theta1.vec)){
  one.result <- pwb2stage(N=2428, theta0=theta1.vec[i]+0.05, theta1=theta1.vec[i], nsims=1e5)
  iw.results[[i]] <- one.result$ests
  mc.error[i] <- one.result$mc.error.bias
}

iw.results <- do.call(rbind, iw.results)

#> Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
#> ℹ Please use `linewidth` instead.

# Now fix theta0 and vary theta1 only:
iw.results2 <- vector("list", length(theta1.vec))
mc.error2 <- numeric(length(theta1.vec))
for(i in 1:length(theta1.vec)){
  one.result <- pwb2stage(theta0=0.3, theta1=theta1.vec[i], nsim=1e5)
  iw.results2[[i]] <- one.result$ests
  mc.error2[i] <- one.result$mc.error.bias
}

iw.results2 <- do.call(rbind, iw.results2)

ma2stage.results <- vector("list", length(theta1.vec))

for(i in 1:length(theta1.vec)){
  ma2stage.results[[i]] <- ma2stage(N=2428, theta0=theta1.vec[i]+0.05, theta1=theta1.vec[i],
                                    n.studies=4, nsims=1e5)
}
  
ma2stage.results <- do.call(rbind, ma2stage.results)
ma2stage.results$theta01 <- "theta1-theta0=0.05"

ma2stage.results2 <- vector("list", length(theta1.vec))

for(i in 1:length(theta1.vec)){
  ma2stage.results2[[i]] <- ma2stage(N=2428, theta0=0.3, theta1=theta1.vec[i],
                                    n.studies=4, nsims=1e5)
}
  
ma2stage.results2 <- do.call(rbind, ma2stage.results2)
ma2stage.results2$theta01 <- "theta0=0.3"