In [1]:
# rm(list=ls())
options(OutDec = ",")
#===============================================================================
In [2]:
#===============================================================================
# Dados da Tabela 3.2
#===============================================================================
dados <- read.csv("../dados/Tabela-3-2.csv",sep=";",dec=",",header=F)
colnames(dados) <- c("Xi1","Xi2","Xi3","Xi4","Xi5")
m <- dim(dados)[1]
dados <- dados[-12,]
print(dados)
Xi1 Xi2 Xi3 Xi4 Xi5 1 1004,6 997,3 1003,0 1005,9 995,8 2 1001,6 1008,6 997,9 1001,3 999,1 3 999,1 992,6 1001,1 1001,6 1002,9 4 1007,9 997,5 991,3 997,8 1000,8 5 999,5 995,6 1004,3 995,6 991,4 6 1003,3 996,8 997,2 993,6 1000,1 7 999,7 1012,1 995,2 1001,8 1002,2 8 1000,1 995,3 990,0 997,5 1003,2 9 1004,3 1001,4 1001,6 999,1 996,4 10 999,0 995,8 989,9 995,1 1002,8 11 1003,2 1004,4 993,5 994,6 997,6 13 1014,0 1008,9 1004,1 1007,9 1000,7 14 1002,2 996,6 1002,7 1004,2 1001,8 15 998,3 997,5 1006,1 996,5 998,1 16 995,8 1000,8 999,1 1002,5 1001,0 17 1004,1 1003,0 1004,8 997,9 999,9 18 1000,1 994,9 1000,1 1004,9 997,3 19 1000,2 996,1 998,0 1006,1 999,4 20 1002,3 999,0 1000,8 1000,7 998,0 21 998,3 998,1 1004,2 1002,1 991,3 22 997,1 1000,7 999,8 1000,6 1001,7 23 1003,6 996,1 1001,4 998,0 991,8 24 999,9 1006,4 1005,1 999,8 1003,0 25 1007,3 999,8 992,5 996,2 998,2
In [3]:
#===============================================================================
# Calculando a amplitude para cada amostra
#===============================================================================
R.i <- as.vector(diff(apply(dados,1,"range")))
tabela <- round(cbind(dados,R.i),1)
colnames(tabela) <- c("Xi1","Xi2","Xi3","Xi4","Xi5","R.i")
print(tabela)
Xi1 Xi2 Xi3 Xi4 Xi5 R.i 1 1004,6 997,3 1003,0 1005,9 995,8 10,1 2 1001,6 1008,6 997,9 1001,3 999,1 10,7 3 999,1 992,6 1001,1 1001,6 1002,9 10,3 4 1007,9 997,5 991,3 997,8 1000,8 16,6 5 999,5 995,6 1004,3 995,6 991,4 12,9 6 1003,3 996,8 997,2 993,6 1000,1 9,7 7 999,7 1012,1 995,2 1001,8 1002,2 16,9 8 1000,1 995,3 990,0 997,5 1003,2 13,2 9 1004,3 1001,4 1001,6 999,1 996,4 7,9 10 999,0 995,8 989,9 995,1 1002,8 12,9 11 1003,2 1004,4 993,5 994,6 997,6 10,9 13 1014,0 1008,9 1004,1 1007,9 1000,7 13,3 14 1002,2 996,6 1002,7 1004,2 1001,8 7,6 15 998,3 997,5 1006,1 996,5 998,1 9,6 16 995,8 1000,8 999,1 1002,5 1001,0 6,7 17 1004,1 1003,0 1004,8 997,9 999,9 6,9 18 1000,1 994,9 1000,1 1004,9 997,3 10,0 19 1000,2 996,1 998,0 1006,1 999,4 10,0 20 1002,3 999,0 1000,8 1000,7 998,0 4,3 21 998,3 998,1 1004,2 1002,1 991,3 12,9 22 997,1 1000,7 999,8 1000,6 1001,7 4,6 23 1003,6 996,1 1001,4 998,0 991,8 11,8 24 999,9 1006,4 1005,1 999,8 1003,0 6,6 25 1007,3 999,8 992,5 996,2 998,2 14,8
In [4]:
#===============================================================================
# Calculando os limites de controle
#===============================================================================
Rbarra <- mean(R.i)
LSC.R <- Rbarra*(1 + 3*0.864082/2.325929 )
LM.R <- Rbarra
LIC.R <- Rbarra*(1 - 3*0.864082/2.325929 )
LIC.R <- max( 0 , LIC.R )
print(rbind(LSC.R,LM.R,LIC.R))
[,1] LSC.R 22,13176 LM.R 10,46667 LIC.R 0,00000
In [5]:
#===============================================================================
# Calculando os limites de controle
#===============================================================================
par(mfrow=c(1,1),lwd=2.0,cex.lab=1.5,cex.axis=1.5,lab=c(7,7,5),
mar=c(5,5,2,2.5),xpd=T,cex.main=2.0,bty="n")
aux <- c(tabela$R.i[1:11],NA,tabela$R.i[12:24])
plot(1:m,aux,type="l",lwd=2,col="darkblue",main="",xlim=c(0,25),
ylim=c(0,25),xlab="Número da amostra",ylab="Amplitude R")
points(1:m,aux,pch=16,cex=2,col="darkblue")
lines(0:m,rep(LSC.R,(m+1)),lwd=3,col="red",lty=2)
lines(0:m,rep(LM.R,(m+1)), lwd=3,col="red",lty=3)
lines(0:m,rep(LIC.R,(m+1)),lwd=3,col="red",lty=2)
text((m+2),LSC.R,"LSC",cex=1.5,lwd=2)
text((m+2),LM.R, "LM", cex=1.5,lwd=2)
text((m+2),LIC.R,"LIC",cex=1.5,lwd=2)
In [6]:
#===============================================================================
# Fim
#===============================================================================