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 e a média para cada amostra
#===============================================================================
xbar.i <- apply(dados,1,"mean")
R.i <- as.vector(diff(apply(dados,1,"range")))
tabela <- round(cbind(dados,xbar.i),1)
colnames(tabela) <- c("Xi1","Xi2","Xi3","Xi4","Xi5","xbar.i")
print(tabela)
Xi1 Xi2 Xi3 Xi4 Xi5 xbar.i 1 1004,6 997,3 1003,0 1005,9 995,8 1001,3 2 1001,6 1008,6 997,9 1001,3 999,1 1001,7 3 999,1 992,6 1001,1 1001,6 1002,9 999,5 4 1007,9 997,5 991,3 997,8 1000,8 999,1 5 999,5 995,6 1004,3 995,6 991,4 997,3 6 1003,3 996,8 997,2 993,6 1000,1 998,2 7 999,7 1012,1 995,2 1001,8 1002,2 1002,2 8 1000,1 995,3 990,0 997,5 1003,2 997,2 9 1004,3 1001,4 1001,6 999,1 996,4 1000,6 10 999,0 995,8 989,9 995,1 1002,8 996,5 11 1003,2 1004,4 993,5 994,6 997,6 998,7 13 1014,0 1008,9 1004,1 1007,9 1000,7 1007,1 14 1002,2 996,6 1002,7 1004,2 1001,8 1001,5 15 998,3 997,5 1006,1 996,5 998,1 999,3 16 995,8 1000,8 999,1 1002,5 1001,0 999,8 17 1004,1 1003,0 1004,8 997,9 999,9 1001,9 18 1000,1 994,9 1000,1 1004,9 997,3 999,5 19 1000,2 996,1 998,0 1006,1 999,4 1000,0 20 1002,3 999,0 1000,8 1000,7 998,0 1000,2 21 998,3 998,1 1004,2 1002,1 991,3 998,8 22 997,1 1000,7 999,8 1000,6 1001,7 1000,0 23 1003,6 996,1 1001,4 998,0 991,8 998,2 24 999,9 1006,4 1005,1 999,8 1003,0 1002,8 25 1007,3 999,8 992,5 996,2 998,2 998,8
In [4]:
#===============================================================================
# Calculando os limites de controle
#===============================================================================
R.i <- as.vector(diff(apply(dados,1,"range")))
Rbarra <- mean(R.i)
sigma0hat <- Rbarra/2.325929
xbarra <- mean(xbar.i)
LSC.Xbar <- xbarra + 3*sigma0hat/sqrt(5)
LM.Xbar <- xbarra
LIC.Xbar <- xbarra - 3*sigma0hat/sqrt(5)
print(rbind(LSC.Xbar,LM.Xbar,LIC.Xbar))
[,1] LSC.Xbar 1006,0399 LM.Xbar 1000,0025 LIC.Xbar 993,9651
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,6,5),
mar=c(5,5,2,2.5),xpd=T,cex.main=2.0,bty="n")
aux <- c(tabela$xbar.i[1:11],NA,tabela$xbar.i[12:24])
plot(1:m,aux,type="l",lwd=2,col="darkblue",main="",xlim=c(0,25),
ylim=c(990,1010),xlab="Número da amostra",ylab=expression(bar(X)))
points(1:m,aux,pch=16,cex=2,col="darkblue")
lines(0:m,rep(LSC.Xbar,(m+1)),lwd=3,col="red",lty=2)
lines(0:m,rep(LM.Xbar,(m+1)), lwd=3,col="red",lty=3)
lines(0:m,rep(LIC.Xbar,(m+1)),lwd=3,col="red",lty=2)
text((m+2),LSC.Xbar,"LSC",cex=1.5,lwd=2)
text((m+2),LM.Xbar, "LM", cex=1.5,lwd=2)
text((m+2),LIC.Xbar,"LIC",cex=1.5,lwd=2)
In [6]:
#===============================================================================
# Fim
#===============================================================================