pow <- sample_size()
pow
library(ggplot2)
library(viridis)
ggplot(pow, aes(x=n, y=pA, fill=prob, group=pA)) +
geom_tile() +
scale_fill_viridis_c("power") +
ylab("detect rate (pA)") + xlab("sample size (n)") +
theme_bw()
Run the code above in your browser using DataLab