# rm(list=ls()) 
options(OutDec = ",")
#===============================================================================
silence <- suppressPackageStartupMessages # Para omitir mensagens de alertas
silence(library(scatterplot3d))

#===============================================================================
# Processo isento de causas especiais
#===============================================================================
x     <- seq(982,1018,length = 401)
y     <- rep(4,401)
means <- rep(1000,4)
sds   <- rep(4.5,4)
times <- c("","T1","T2","T3","T4","")
z1    <- dnorm(x, mean = means[1], sd = sds[1])
z2    <- dnorm(x, mean = means[2], sd = sds[2])
z3    <- dnorm(x, mean = means[3], sd = sds[3])
z4    <- dnorm(x, mean = means[4], sd = sds[4])
s3d   <- scatterplot3d(x=x,y=y,z=z1,lwd=3,color="darkblue",type="l",
         ylab="Tempo",xlab=expression(x),zlab="Densidade",angle=45,box=FALSE,
         ylim=c(0,5),xlim=c(982,1018),zlim=c(0,0.12),grid=FALSE,
         y.ticklabs=times)
s3d$points3d(x=x,y=y-1,z=z2,type="l",lwd=3,col="darkorange")
s3d$points3d(x=x,y=y-2,z=z3,type="l",lwd=3,col="darkgreen")
s3d$points3d(x=x,y=y-3,z=z4,type="l",lwd=3,col="brown")
s3d$points3d(x=rep(1000,401),y=seq(0,5.5,length=401),z=rep(0,401),type="l",
         lwd=2,col="black",lty=2)
coords <- s3d$xyz.convert(982,0,0.13)
legend(coords$x, coords$y,
       legend = c("Valor-alvo"),
       col = "black",
       lwd = 2,lty=2,bty="n")

#===============================================================================
# Causa especial que altera a média do processo
#===============================================================================
x     <- seq(982,1018,length = 401)
w     <- seq(1002,1038,length = 401)
y     <- rep(4,401)
means <- c(rep(1020,3),1000)
sds   <- rep(4.5,4)
times <- c("","T1","T2","T3","T4","")
z1    <- dnorm(w, mean = means[1], sd = sds[1])
z2    <- dnorm(w, mean = means[2], sd = sds[2])
z3    <- dnorm(w, mean = means[3], sd = sds[3])
z4    <- dnorm(x, mean = means[4], sd = sds[4])
s3d   <- scatterplot3d(x=w,y=y,z=z1,lwd=3,color="darkblue",type="l",
         ylab="Tempo",xlab=expression(x),zlab="Densidade",angle=45,box=FALSE,
         ylim=c(0,5),xlim=c(982,1040),zlim=c(0,0.12),grid=FALSE,
         y.ticklabs=times)
s3d$points3d(x=w,y=y-1,z=z2,type="l",lwd=3,col="darkorange")
s3d$points3d(x=w,y=y-2,z=z3,type="l",lwd=3,col="darkgreen")
s3d$points3d(x=x,y=y-3,z=z4,type="l",lwd=3,col="brown")
s3d$points3d(x=rep(1000,401),y=seq(0,5.5,length=401),z=rep(0,401),type="l",
         lwd=2,col="black",lty=2)
s3d$points3d(x=rep(1020,401),y=seq(0,5.5,length=401),z=rep(0,401),type="l",
         lwd=2,col="red",lty=2)
coords <- s3d$xyz.convert(985,0,0.13)
legend(coords$x, coords$y,legend = c("Valor-alvo","Média alterada"),
         col = c("black","red"),lwd=c(2,2),lty=c(2,2),bty="n")

#===============================================================================
# Causa especial que altera a média do processo
#===============================================================================
x     <- seq(982,1018,length = 401)
w     <- seq(990,1050,length = 401)
y     <- rep(4,401)
means <- c(rep(1020,3),1000)
sds   <- c(rep(9,3),4.5)
times <- c("","T1","T2","T3","T4","")
z1    <- dnorm(w, mean = means[1], sd = sds[1])
z2    <- dnorm(w, mean = means[2], sd = sds[2])
z3    <- dnorm(w, mean = means[3], sd = sds[3])
z4    <- dnorm(x, mean = means[4], sd = sds[4])
s3d   <- scatterplot3d(x=w,y=y,z=z1,lwd=3,color="darkblue",type="l",
         ylab="Tempo",xlab=expression(x),zlab="Densidade",angle=45,box=FALSE,
         ylim=c(0,5),xlim=c(982,1050),zlim=c(0,0.12),grid=FALSE,
         y.ticklabs=times)
s3d$points3d(x=w,y=y-1,z=z2,type="l",lwd=3,col="darkorange")
s3d$points3d(x=w,y=y-2,z=z3,type="l",lwd=3,col="darkgreen")
s3d$points3d(x=x,y=y-3,z=z4,type="l",lwd=3,col="brown")
s3d$points3d(x=rep(1000,401),y=seq(0,5.5,length=401),z=rep(0,401),type="l",
             lwd=2,col="black",lty=2)
s3d$points3d(x=rep(1020,401),y=seq(0,5.5,length=401),z=rep(0,401),type="l",
             lwd=2,col="red",lty=2)
coords <- s3d$xyz.convert(985,0,0.13)
legend(coords$x, coords$y,legend = c("Valor-alvo","Média alterada"),
         col = c("black","red"),lwd=c(2,2),lty=c(2,2),bty="n")
coords <- s3d$xyz.convert(1020,0,0.1)
legend(coords$x, coords$y,legend = "Variabilidade alterada",
       col=0,lwd=0,lty=0,bty="n")

#===============================================================================
# Fim
#===============================================================================