#> 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"