Nothing
#' @export
drawDiagram.power.htest <- function(findDesign.output, print.row=NULL, xmax=NULL, ymax=NULL){
m <- Sm <- decision <- analysis <- NULL
des <- findDesign.output
# des is output from power.prop.test()
n <- 2*ceiling(des$n)
coords <- expand.grid(0:n, 1:n)
diag.df <- data.frame(Sm=as.numeric(coords[,1]),
m=as.numeric(coords[,2]),
decision=rep("Continue", nrow(coords)))
diag.df$decision <- as.character(diag.df$decision)
diag.df$decision[coords[,1]>coords[,2]] <- NA
diag.df$decision[diag.df$m==n] <- "Make decision"
# Add shading:
diag.df.subset <- diag.df[!is.na(diag.df$decision),]
diag.df.subset$analysis <- "No"
diag.df.subset$analysis[diag.df.subset$m==n] <- "Yes"
plot.title <- "Single-stage two-arm design"
sub.text1 <- paste("Max no. of analyses: 1. Max(N): ", n, sep="")
plot.subtitle2 <- bquote(.(sub.text1))
diagram <- pkgcond::suppress_warnings(ggplot2::ggplot(data=diag.df.subset, mapping = aes(x=m, y=Sm, fill=decision, alpha=analysis))+
scale_alpha_discrete(range=c(0.5, 1)),
"Using alpha for a discrete variable is not advised")
diagram <- diagram +
geom_tile(color="white")+
labs(fill="Decision",
alpha="Analysis",
x="Number of participants",
y="Number of responses on treatment + non-responses on control",
title=plot.title,
subtitle = plot.subtitle2
)+
coord_cartesian(expand = 0)+
theme_minimal()
if(!is.null(xmax)){
diagram <- diagram +
expand_limits(x=xmax)
}
if(!is.null(ymax)){
diagram <- diagram +
expand_limits(y=ymax)
}
print(diagram)
return(diagram)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.