server.r

# Define a server for the Shiny app
function(input, output) {
  
  

  # Fill in the spot we created for a plot
  output$Plot <- renderPlot({

	# mu=seq(0,1,length=100)
	
	# mypow1=ifelse(mu==input$mu0,input$myalpha1,
	# 1-pnorm((sqrt(input$n)/input$sigma)*(input$mu0-mu)+qnorm(1-input$myalpha1)))
    # plot(mu,mypow1,type="l",lwd=2,ylim=c(0,1),cex.lab=2.2,xlab=expression(mu-mu[0]),
     # main="",ylab=expression(pi[n](mu)),cex=1.4, cex.axis=2,bty="n",
	 # xlim=c(0,1),col="black")

disable("mu")
	 
myalpha1=0.01
mu0=1
sigma=input$sigma
n=input$n
mu=seq(0,2,length=100)
mypow1=ifelse(mu==mu0,myalpha1,1-pnorm((sqrt(n)/sigma)*(mu0-mu)+qnorm(1-myalpha1)))

#par(mar=c(4,5.5, 0.5, 0.5))	
plot(mu-mu0,mypow1,type="l",lwd=2,ylim=c(0,1),cex.lab=1.5,xlab=expression(mu*-~mu[0]),
     main="",ylab=expression(pi[n](mu)),cex=1.2, cex.axis=1.2,bty="n",
	 xlim=c(-1,1),col="black",xaxt="n")
axis(1, c(-1,0,1),labels=c(-1,0,1),cex.axis=1.2, tck=-.02)

	 
myalpha2=0.05
mypow2=ifelse(mu==mu0,myalpha2,1-pnorm((sqrt(n)/sigma)*(mu0-mu)+qnorm(1-myalpha2))) 
lines(mu-mu0,mypow2,cex.lab=1.2,xlab="",ylab="",col=c("red"),cex=1.4, 
yaxt="n", xaxt="n",bty="n",type="l",lty=2,lwd=2)

myalpha3=0.1
mypow3=ifelse(mu==mu0,myalpha3,1-pnorm((sqrt(n)/sigma)*(mu0-mu)+qnorm(1-myalpha3))) 
lines(mu-mu0,mypow3,cex.lab=2.2,xlab="",ylab="",col=c("blue"),cex=1.2, 
yaxt="n", xaxt="n",bty="n",type="l",lty=3,lwd=3)
	 
legend("topleft",lty=c(NA,1,2,3),col=c("black","black","red","blue"),bty="n",
cex=1.2,lwd=3,legend=c(
as.expression(bquote(n*"="*.(n))),
as.expression(bquote(alpha*"="*.(format(myalpha1, digits = 2)))),
as.expression(bquote(alpha*"="*.(format(myalpha2, digits = 2)))),
as.expression(bquote(alpha*"="*.(format(myalpha3, digits = 2))))))
	 
	 
	 
  })
}