Skip to content
Extraits de code Groupes Projets
Valider 063f9160 rédigé par Eugen Pircalabelu's avatar Eugen Pircalabelu
Parcourir les fichiers

Upload New File

parent 733942b6
Aucune branche associée trouvée
Aucune étiquette associée trouvée
Aucune requête de fusion associée trouvée
# 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=2.2,xlab=expression(mu*-~mu[0]),
main="",ylab=expression(pi[n](mu)),cex=1.4, cex.axis=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)
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,
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,
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(
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))))))
})
}
0% Chargement en cours ou .
You are about to add 0 people to the discussion. Proceed with caution.
Terminez d'abord l'édition de ce message.
Veuillez vous inscrire ou vous pour commenter