In [1]:
# rm(list=ls())
options(OutDec = ",")
#===============================================================================
In [2]:
#===============================================================================
# Dados da Tabela 2.4a - Figura 2.7
#===============================================================================
dados <- read.csv("../dados/Tabela-2-4a.csv",sep=";",dec=",",header=F)
dados <- as.vector(unlist(dados))
print(dados)
[1] 1000,0 1001,0 1002,0 1006,0 1003,0 1012,0 1004,0 1020,0 998,0 997,0 [11] 998,0 997,0 996,0 998,0 996,0 990,0 982,0 1000,0 1001,0 1003,0 [21] 1002,0 1005,0 1020,3 1008,0 1017,0 1012,0 998,0 1002,0 997,0 1003,0 [31] 1000,0 995,0 1006,0 999,0 995,0 990,0 988,0 983,0 988,0 985,0
In [3]:
#===============================================================================
# Definindo os limites de controle de X.
# Observação: aprenderemos a calcular estes limites ao longo da disciplina.
#===============================================================================
LSC.X <- 1015
LM.X <- 1000
LIC.X <- 985
In [4]:
#===============================================================================
# Gráfico de controle de X.
#===============================================================================
par(mfrow=c(1,1),lwd=2.0,cex.lab=1.5,cex.axis=1.5,lab=c(10,7,5),
mar=c(5,5,2,2.5),xpd=T,cex.main=2.0,bty="n")
plot(1:40,dados,type="l",lwd=2,col="darkblue",main="",xlim=c(0,45),
ylim=c(975,1025),xlab="Número da observação",ylab=paste(expression(X)," (mL)",""))
points(1:40,dados,pch=16,cex=2,col="darkblue")
lines(0:42,rep(LSC.X,43),lwd=3,col="red",lty=2)
lines(0:42,rep(LM.X,43), lwd=3,col="red",lty=3)
lines(0:42,rep(LIC.X,43),lwd=3,col="red",lty=2)
text(45,LSC.X,"LSC",cex=1.5,lwd=2)
text(45,LM.X, "LM", cex=1.5,lwd=2)
text(45,LIC.X,"LIC",cex=1.5,lwd=2)
In [5]:
#===============================================================================
# Fim
#===============================================================================