diff --git a/server.R b/server.R
new file mode 100644
index 0000000000000000000000000000000000000000..8cf3ae00be8e09fe3b46c2e69d9e1bac5cef375e
--- /dev/null
+++ b/server.R
@@ -0,0 +1,248 @@
+library(shiny)
+# library(xtable)
+
+plotzmaxlim=5 #limite max de l'axe des Z
+plotzlab="Z~N(0;1)"
+linezlab=1
+linexlab=3
+linemulab=2.5
+
+yt=0.48 #y des textes
+xt=-5 #x des textes
+ztxt=1.5 #zoom text
+txt="" #contenu des textes
+
+z=seq(-4,4,length=200) #détermine les valeurs de z à utiliser : de -4 à 4 créer 200 valeurs et les mettre en vecteur
+y=dnorm(z) #créer pour chaque z les équivalents y de la courbe de gauss
+
+#Define server logic required to plot various variables against mpg
+shinyServer(function(input, output) {
+
+ getInputValues<-reactive({
+ return(input)#collect all inputs
+ })
+
+ getComputedValues<-reactive({
+ v<-getInputValues() # get all values of input list
+ cv<-list()#created empty computed values list
+ if(v$Tabset!=3){
+ if(v$Tabset==1){
+ if(v$sens=="inf"){
+ cv$z1=-5
+ cv$z2=round((v$x-v$mx)/v$sx,2)
+ }
+ if(v$sens=="sup"){
+ cv$z2=5
+ cv$z1=round((v$x-v$mx)/v$sx,2)
+ }
+ if(v$sens=="equal"){
+ cv$z2=round((v$x-v$mx)/v$sx,2)
+ cv$z1=round((v$x-v$mx)/v$sx,2)
+ }
+ }
+
+ if(v$Tabset==2){
+ cv$x1=min(v$x1,v$x2)
+ cv$x2=max(v$x1,v$x2)
+ cv$z1=round((cv$x1-v$mx)/v$sx,2)
+ cv$z2=round((cv$x2-v$mx)/v$sx,2)
+ }
+
+ cv$z1l=ifelse(cv$z1 < -3,-3,cv$z1) # limite la valeur de z à -3 pour le calcul de la position de
+ cv$z2l=ifelse(cv$z2 > 3,3,cv$z2) # limite la valeur de z à 3 pour le calcul de la position dep
+
+ cv$z1lp=ifelse(cv$z1 < -4,-4,cv$z1) # limite la valeur de z à -4 pour le tracé du polugone sinon il overwrite le pointillé
+ cv$z2lp=ifelse(cv$z2 > 4,4,cv$z2) # limite la valeur de z à 4 pour le tracé du polugone sinon il overwrite le pointillé
+ cv$xlp=seq(cv$z1lp,cv$z2lp,length=200) #crée une seconde abscisse de -2 à 2
+ cv$ylp=dnorm(cv$xlp) #crée pour chaque xles équivalent y Normaux
+
+ cv$pz1=pnorm(cv$z1,mean=0,sd=1)
+ cv$pz1r=round(cv$pz1,4)
+ cv$pz1i=1-cv$pz1
+ cv$pz1ir=1-cv$pz1r#round(pz1i,4)
+
+ cv$pz2=pnorm(cv$z2,mean=0,sd=1)
+ cv$pz2r=round(cv$pz2,4)
+ cv$pz2i=1-cv$pz2
+ cv$pz2ir=1-cv$pz2r#round(pz2i,4)
+
+ cv$p=cv$pz2-cv$pz1
+ cv$pr=cv$pz2r-cv$pz1r#round(p,4)
+ cv$pp=round(cv$pr*100,0)
+
+ if(v$Tabset==1){
+ if(v$sens=="inf"){
+ cv$zx=(max(cv$z2l,cv$z1l)-((max(cv$z2l,cv$z1l)-min(cv$z2l,cv$z1l))/3))
+ }
+ if(v$sens=="sup"){
+ cv$zx=(min(cv$z2l,cv$z1l)+((max(cv$z2l,cv$z1l)-min(cv$z2l,cv$z1l))/3))
+ }
+ if(v$sens=="equal"){
+ cv$zx=cv$z2l
+ }
+ }
+
+ if(v$Tabset==2){
+ cv$zx=(cv$z2l-(cv$z2l-cv$z1l)/2)
+ }
+ }
+ return(cv)
+ })
+
+ output$zPlot <- renderPlot({
+ v<-getInputValues()
+ cv<-getComputedValues()
+ if(v$Tabset==1){
+ par(mai=c(1.02,0.8,0.9,1.5)) #marges bottom left top right en inches les valeurs par défaut sont 1 1 1 1
+ plot(z,y,type="l",lwd=2,col="blue",yaxt="n",bty="n",,las=1,xaxs="i",yaxs="i",cex.lab=1,cex.axis=1.5,xlim=c(-5,plotzmaxlim),ylim=c(0,0.4),xaxp=c(-4,4,8),ylab="density",xlab="")
+ axis(2,las=2)#col="red",col.axis="blue",font.axis=2
+ lines(x=c(-4.7,-4),y = c(0,0),lty=2)
+ lines(x=c(4.7,4),y = c(0,0),lty=2)
+ polygon(c(cv$z1lp,cv$xlp,cv$z2lp),c(0,cv$ylp,0),col="lightblue") #trace un polygone qui commence à -2;0 puis x;y puis 2;0 et comme c'est polygon il ferme entre 2;0 et -2;0
+ mtext(plotzlab,side=1,line=1,at=5,cex=ztxt,adj = 0)
+ for(s in -4:4){#Ecrire valeurs de x sur abscisse
+ vx=v$mx+s*v$sx
+ mtext(vx,side=1,line=linexlab,at=s,cex=1.5)
+ }
+ mtext(paste("X~N(",v$mx,";",v$sx^2,")",sep=""),side=1,line=linexlab,at=5,cex=1.2,adj = 0)
+ text(cv$zx,0.1,cv$pr,cex=2)#cex = zoom taille texte
+ }
+ })
+
+ output$zPlotInt <- renderPlot({
+ v<-getInputValues()
+ cv<-getComputedValues()
+ if(v$Tabset==2){
+ par(mai=c(1.02,0.8,0.9,1.5)) #marges bottom left top right en inches les valeurs par défaut sont 1 1 1 1
+ plot(z,y,type="l",lwd=2,col="blue",yaxt="n",bty="n",,las=1,xaxs="i",yaxs="i",cex.lab=1,cex.axis=1.5,xlim=c(-5,plotzmaxlim),ylim=c(0,0.4),xaxp=c(-4,4,8),ylab="density",xlab="")
+ axis(2,las=2)#col="red",col.axis="blue",font.axis=2
+ lines(x=c(-4.7,-4),y = c(0,0),lty=2)
+ lines(x=c(4.7,4),y = c(0,0),lty=2)
+ polygon(c(cv$z1lp,cv$xlp,cv$z2lp),c(0,cv$ylp,0),col="lightblue") #trace un polygone qui commence à -2;0 puis x;y puis 2;0 et comme c'est polygon il ferme entre 2;0 et -2;0
+ mtext(plotzlab,side=1,line=1,at=5,cex=ztxt,adj = 0)
+ for(s in -4:4){#Ecrire valeurs de x sur abscisse
+ vx=v$mx+s*v$sx
+ mtext(vx,side=1,line=linexlab,at=s,cex=1.5)
+ }
+ mtext(paste("X~N(",v$mx,";",v$sx^2,")",sep=""),side=1,line=linexlab,at=5,cex=1.2,adj = 0)
+ text(cv$zx,0.1,cv$pr,cex=2)#cex = zoom taille texte
+ }
+ })
+
+
+ output$Text <- renderText({
+ v<-getInputValues()
+ cv<-getComputedValues()
+
+ txt=paste("<p>Si X est une variable aléatoire Normale X~N(μ;σ²) telle que X~N(",v$mx,";",v$sx^2,") alors Z=(X-μ)/σ est aussi une variable aléatoire Normale de paramètres Z~N(0;1) ",sep="")
+ if(v$sens=="inf" & cv$z2>-4 & cv$z2<0){
+ txt=paste(txt," et P(X≤x) = P(Z≤z).","</p>",sep="")
+ txt=paste(txt,"<p>Si x=",v$x,", alors z = (x-μ)/σ = (",v$x,"-",v$mx,")/",v$sx," = ",v$x-v$mx,"/",v$sx," = ",cv$z2,"</p>",sep="")
+ txt=paste(txt,"<p>Donc P(X≤",round(cv$z2*v$sx+v$mx,2),") = P(Z≤",cv$z2,") = 1-P(Z≤",cv$z2*-1,") = 1-",cv$pz2ir," = ",cv$pz2r,"</p>",sep="")
+ }
+ if(v$sens=="inf" & cv$z2<4 & cv$z2>=0){
+ txt=paste(txt," et P(X≤x) = P(Z≤z).","</p>",sep="")
+ txt=paste(txt,"<p>Si x=",v$x,", alors z = (x-μ)/σ = (",v$x,"-",v$mx,")/",v$sx," = ",v$x-v$mx,"/",v$sx," = ",cv$z2,"</p>",sep="")
+ txt=paste(txt,"<p>Donc P(X≤",round(cv$z2*v$sx+v$mx,2),") = P(Z≤",cv$z2,") = ",cv$pz2r,"</p>",sep="")
+ }
+ if(cv$z1>-4 & cv$z1<0 & v$sens=="sup"){
+ txt=paste(txt," et P(X≥x) = P(Z≥z).","</p>",sep="")
+ txt=paste(txt,"<p>Si x=",v$x,", alors z = (x-μ)/σ = (",v$x,"-",v$mx,")/",v$sx," = ",v$x-v$mx,"/",v$sx," = ",cv$z1,"</p>",sep="")
+ txt=paste(txt,"<p>Donc P(X≥",round(cv$z1*v$sx+v$mx,2),") = P(Z≥",cv$z1,") = P(Z≤",cv$z1*-1,") = ",cv$pz1ir,"</p>",sep="")
+ }
+ if(cv$z1<4 & cv$z1>=0 & v$sens=="sup"){
+ txt=paste(txt," et P(X≥x) = P(Z≥z).","</p>",sep="")
+ txt=paste(txt,"<p>Si x=",v$x,", alors z = (x-μ)/σ = (",v$x,"-",v$mx,")/",v$sx," = ",v$x-v$mx,"/",v$sx," = ",cv$z1,"</p>",sep="")
+ txt=paste(txt,"<p>Donc P(X≥",round(cv$z1*v$sx+v$mx,2),") = P(Z≥",cv$z1,") = 1-P(Z≤",cv$z1,") = 1-",cv$pz1r," = ",cv$pz1ir,"</p>",sep="")
+ }
+ if(v$sens=="inf" & cv$z2>=4){
+ txt=paste(txt," et P(X≤x) = P(Z≤z).","</p>",sep="")
+ txt=paste(txt,"<p>Si x=",v$x,", alors z = (x-μ)/σ = (",v$x,"-",v$mx,")/",v$sx," = ",v$x-v$mx,"/",v$sx," = ",cv$z2,"</p>",sep="")
+ txt=paste(txt,"<p>Donc P(X≤",round(cv$z2*v$sx+v$mx,2),") = P(Z≤",cv$z2,") = P(Z≤∞) = ",cv$pz2r,"</p>",sep="")
+ }
+ if(v$sens=="inf" & cv$z2<=-4){
+ txt=paste(txt," et P(X≤x) = P(Z≤z).","</p>",sep="")
+ txt=paste(txt,"<p>Si x=",v$x,", alors z = (x-μ)/σ = (",v$x,"-",v$mx,")/",v$sx," = ",v$x-v$mx,"/",v$sx," = ",cv$z2,"</p>",sep="")
+ txt=paste(txt,"<p>Donc P(X≤",round(cv$z2*v$sx+v$mx,2),") = P(Z≤",cv$z2,") = 1-P(Z≤",cv$z2*-1,") = 1-P(Z≤∞) = 1-1 = ",cv$pz2r,"</p>",sep="")
+ }
+ if(cv$z1<=-4 & v$sens=="sup"){
+ txt=paste(txt," et P(X≥x) = P(Z≥z).","</p>",sep="")
+ txt=paste(txt,"<p>Si x=",v$x,", alors z = (x-μ)/σ = (",v$x,"-",v$mx,")/",v$sx," = ",v$x-v$mx,"/",v$sx," = ",cv$z1,"</p>",sep="")
+ txt=paste(txt,"<p>Donc P(X≥",round(cv$z1*v$sx+v$mx,2),") = P(Z≥",cv$z1,") = 1-P(Z≤",cv$z1*-1,") = 1-P(Z≤∞) = 1-1 = ",cv$pz1r,"</p>",sep="")
+ }
+ if(cv$z1>=4 & v$sens=="sup"){
+ txt=paste(txt," et P(X≥x) = P(Z≥z).","</p>",sep="")
+ txt=paste(txt,"<p>Si x=",v$x,", alors z = (x-μ)/σ = (",v$x,"-",v$mx,")/",v$sx," = ",v$x-v$mx,"/",v$sx," = ",cv$z1,"</p>",sep="")
+ txt=paste(txt,"<p>Donc P(X≥",round(cv$z1*v$sx+v$mx,2),") = P(Z≥",cv$z1,") = 1-P(Z≤",cv$z1,") = 1-P(Z≤∞) = 1-1 = ",cv$pz1r,"</p>",sep="")
+ }
+ if(v$sens=="equal"){
+ txt=paste(txt," et P(X=x) = P(Z=z) = 0.","</p>",sep="")
+ }
+ paste(txt)
+ })
+
+ output$TextInt <- renderText({
+ v<-getInputValues()
+ cv<-getComputedValues()
+
+ txt=paste("<p>Si X est une variable aléatoire Normale X~N(μ;σ²) telle que X~N(",v$mx,";",v$sx^2,") alors Z=(X-μ)/σ est aussi une variable aléatoire Normale de paramètres Z~N(0;1) </p>",sep="")
+ txt=paste(txt,"<p>P(x<sub>1</sub> ≤ X ≤ x<sub>2</sub>) = P(z<sub>1</sub> ≤ Z ≤ z<sub>2</sub>).","</p>",sep="")
+ txt=paste(txt,"<p>Si x<sub>1</sub>=",v$x1,", alors z<sub>1</sub> = (x<sub>1</sub>-μ)/σ = (",v$x1,"-",v$mx,")/",v$sx," = ",v$x1-v$mx,"/",v$sx," = ",cv$z1,"</p>",sep="")
+ txt=paste(txt,"<p>Si x<sub>2</sub>=",v$x2,", alors z<sub>2</sub> = (x<sub>2</sub>-μ)/σ = (",v$x2,"-",v$mx,")/",v$sx," = ",v$x2-v$mx,"/",v$sx," = ",cv$z2,"</p>",sep="")
+ txt=paste(txt,"<p>Donc P(",cv$x1," ≤ X ≤ ",cv$x2,") = P(",cv$z1," ≤ Z ≤ ",cv$z2,") = P(Z ≤ ",cv$z1,")-P(Z ≤ ",cv$z2,") = ",cv$pz1ir,"-",cv$pz2ir," = ",cv$pr,"</p>",sep="")
+ paste(txt)
+ })
+
+ output$Table <- renderText({
+ v<-getInputValues()
+ if(v$signeZ == "pos"){
+ rows<-seq(0.0,3.4,by=0.1)
+ cols<-seq(0.00,0.09,by=0.01)
+ } else {
+ rows<-rev(seq(-3.4,-0.0,by=0.1))
+ cols<-rev(seq(-0.09,0.00,by=0.01))
+ }
+
+ p<-c()
+# définir ici table header et faire html code de Table a la main car renderTable force les colsnames en header mais ne permet pas les rownames en header
+ txt<-'<table cellspacing="0" cellpadding="5" border="0" style="border:1px solid #000;">'
+ txt<-paste(txt,'<tr style="border-bottom:1px solid #000;"><th style="border-right:1px solid #000;">Z</th>',sep="")
+ for(j in 1:length(cols)){
+ txt<-paste(txt,'<th>',format(cols[[j]], nsmall = 2),'</th>',sep="")
+ }
+ txt<-paste(txt,'</tr>',sep="")
+ for(i in 1:length(rows)){
+ if(i == 11 || i == 21 ){
+ txt<-paste(txt,'<tr style="border-bottom:1px solid #000;">',sep="")
+ } else {
+ txt<-paste(txt,'<tr>',sep="")
+ }
+ txt<-paste(txt,'<th style="border-right:1px solid #000;">',format(rows[[i]], nsmall = 1),'</th>',sep="")
+ for(j in 1:length(cols)){
+ if(v$sensTable == "inf"){
+ txt<-paste(txt,'<td>',sprintf("%.4f",round(pnorm(rows[[i]]+cols[[j]]),4)),'</td>',sep="")
+ } else {
+ txt<-paste(txt,'<td>',sprintf("%.4f",1-round(pnorm(rows[[i]]+cols[[j]]),4)),'</td>',sep="")
+ }
+
+ }
+ txt<-paste(txt,'</tr>',sep="")
+ }
+ txt<-paste(txt,'</table>',sep="")
+
+
+ paste(txt)
+# for(i in 1:length(rows)){
+# for(j in 1:length(cols)){
+# p<-c(p,round(pnorm(rows[[i]]+cols[[j]]),4))
+# }
+# }
+# m<-matrix(p,length(rows),length(cols),byrow=TRUE)
+# d<-data.frame(m)
+# rownames(d)<-format(rows, nsmall = 2)
+# colnames(d)<-format(cols, nsmall = 2)
+# d
+ })
+
+})
+