likelihood of true positive and false positive as a function of power and alpha

3dplotLR.R

#Create single plot for graph----
power <- seq(0.01,1,0.01)
alpha <- seq(0.01,1,0.01)
pH0 <- 0.5
tp<-((1-pH0)*power)
fp<-((pH0)*alpha)
likelihood_ratio<-outer(tp, fp, "/")

#Create color shading
nrz<-nrow (likelihood_ratio)
ncz<-ncol(likelihood_ratio)
jet.colors<-colorRampPalette(c("yellow","red"))
nbcol<-100
color<-jet.colors(nbcol)
zfacet<-likelihood_ratio[-1, -1]+likelihood_ratio[-1,-ncz]+likelihood_ratio[-nrz,-1]+likelihood_ratio[-nrz,-ncz]
facetcol<-cut(zfacet,nbcol)

#Save plot----
png(file = "LRplot.png", width = 4000, height = 4000, res = 500)
persp(x=power, y=alpha, z=likelihood_ratio, theta = 130, phi = 0, col = color[facetcol], ticktype = "detailed", lwd=0.2)
dev.off()

#Create 360 plots for animated .gif (to combine in gif application)----
library(stringr)

#loop through plots
for(i in 0:360){
  currFile <- paste(getwd(), "/plot_", str_pad(i, 3, pad = "0"), ".png", sep = "")
  png(file = currFile, width = 1200, height = 1000, res = 200)
  x <- seq(0,1,0.01)
  y <- x
  power <- seq(0.01,1,0.01)
  alpha <- seq(0.01,1,0.01)
  pH0 <- 0.5
  tp<-((1-pH0)*power)
  fp<-((pH0)*alpha)
  likelihood_ratio<-outer(tp, fp, "/")
  persp(x=power, y=alpha, z=likelihood_ratio, theta = i, phi = 10, col = color[facetcol], ticktype = "detailed", lwd=0.2)
  dev.off()
}

#likelihood ratio for some example, e.g., a range of power, alpha = 0.05
power<-seq(1, 0.1, -0.1)
alpha<-0.05
pH0 <- 0.5
tp<-((1-pH0)*power)
fp<-((pH0)*alpha)
tp/fp

#likelihood ratio for some example, e.g., a range of alpha levels, power = 0.8
power<-0.8
alpha<-seq(0.05, 0.5, 0.05)
pH0 <- 0.5
tp<-((1-pH0)*power)
fp<-((pH0)*alpha)
tp/fp