diff --git a/server.r b/server.r index 93e38ecf45ec4de5499d9f66fe9c0f7c47286826..f9acf86cf9e3e64c8e7f2a713a8aa99220f1b7c8 100644 --- a/server.r +++ b/server.r @@ -23,25 +23,26 @@ 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=2.2,xlab=expression(mu*-~mu[0]), - main="",ylab=expression(pi[n](mu)),cex=1.4, cex.axis=2,bty="n", + +#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=2, tck=-.02) +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=2.2,xlab="",ylab="",col=c("red"),cex=1.4, +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.4, +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=2,lwd=3,legend=c( +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)))), @@ -51,7 +52,3 @@ as.expression(bquote(alpha*"="*.(format(myalpha3, digits = 2)))))) }) } - - - -