Gráficos de Controle X_bar e R - TCC ROBSON

Função xR1 com especificações para média e desvio-padrão
Função xR2 sem especificações para média e desvio-padrão

Funções

# m = número de amostras observadas.
# n = número de observações de cada amostra.
# y = vetor com as observações, ordenadas por amostra, da variável observada.
# mi = especificação para a média do processo.
# sigma = especificação para o desvio-padrão do processo.


xR1=function(m,n,y,mi,sigma){ 
  if (n>1) { 
    amplit=numeric(m) 
    mat=matrix(y,m,n,byrow=T) 
    media=apply(mat,1,mean) 
    for (i in 1:m) {
  amplit[i]=max(mat[i,])-min(mat[i,])}
 } 
  else { 
    amplit=numeric(m-1) 
    media=y 
    for (i in 2:m) 
      amplit[i-1]=sqrt((y[i]-y[i-1])^2) 
    } 
  
  d2=c(1.1296,1.6918,2.0535,2.3248,2.5404,2.7074, 2.8501,2.9677,3.0737,3.1696)
  d3=c(.8541,.8909,.8800,.8674,.8508, .8326,.8209,.8102,.7978,.7890) 
  LIC1=mi-3*sigma/sqrt(n) 
  LSC1=mi+3*sigma/sqrt(n) 
  if(n>1) { 
    LIC2=d2[n-1]*sigma-3*d3[n-1]*sigma 
    LSC2=d2[n-1]*sigma+3*d3[n-1]*sigma 
    if (LIC2 < 0) 
      LIC2=0 
    } else { 
      LIC2=d2[n]*sigma-3*d3[n]*sigma 
      LSC2=d2[n]*sigma+3*d3[n]*sigma 
      if (LIC2 < 0)
        LIC2=0
      }
  
  par(mfrow=c(2,1)) 
  amostra1=matrix(rep(seq(1,m),4),m,4,byrow=F)
  if (n>1) amostra2=amostra1 else amostra2=amostra1[-1,] 
  xbar=cbind(rep(mi,m),rep(LIC1,m),rep(LSC1,m),media) 
  if (n>1) { 
    rbar=cbind(rep(sigma,m),rep(LIC2,m),rep(LSC2,m),amplit) 
    } 
  else rbar=cbind(rep(sigma,m-1),rep(LIC2,m-1), 
                  rep(LSC2,m-1),amplit) 
  
  matplot(amostra2,rbar,type="o",ylab="Amplitude",
          col=c("black","red","red","blue"),
          ylim=c(min(rbar),max(rbar)),
          xlab="Amostras",main="Gráfico R",pch=20, lty=1, lwd=2)
  
  matplot(amostra1,xbar,type="o",ylab="Média",xlab="Amostras",
          main=expression(paste("Gráfico " * bar(x))),col=c("black","red","red", "blue"),ylim=c(min(xbar),max(xbar)),pch=20, lty=1, lwd=2) 
  out=rbind(round(xbar[1,1:3],2),round(rbar[1,1:3],2))
  colnames(out)=c("Referência","LIC","LSC")
  out 
  }


xR2=function(m,n,y){
  if (n>1) { 
    amplit=numeric(m) 
    mat=matrix(y,m,n,byrow=T) 
    media=apply(mat,1,mean) 
    for (i in 1:m) 
      amplit[i]=max(mat[i,])-min(mat[i,]) 
    x2barras=mean(media) 
    Rbarra=mean(amplit) 
    d2=c(1.1296,1.6918,2.0535,2.3248,2.5404,2.7074,2.8501, 2.9677,3.0737,3.1696)
    d3=c(.8541,.8909,.8800,.8674,.8508,.8326,.8209, .8102,.7978,.7890) 
    LIC1=x2barras-3*Rbarra/(d2[n-1]*sqrt(n)) 
    LSC1=x2barras+3*Rbarra/(d2[n-1]*sqrt(n)) 
    LIC2=Rbarra-3*d3[n-1]*Rbarra/d2[n-1] 
    LSC2=Rbarra+3*d3[n-1]*Rbarra/d2[n-1] 
    if (LIC2 < 0) 
      LIC2=0 
    } 
  else { 
      amplit=numeric(m-1) 
      media=y 
      for (i in 2:m)
        amplit[i-1]=sqrt((y[i]-y[i-1])^2)
      x2barras=mean(media) 
      Rbarra=mean(amplit) 
      d2=1.1296 
      d3=.8541 
      LIC1=x2barras-3*Rbarra/(d2*sqrt(n)) 
      LSC1=x2barras+3*Rbarra/(d2*sqrt(n)) 
      LIC2=Rbarra-3*d3*Rbarra/d2 
      LSC2=Rbarra+3*d3*Rbarra/d2 
      if (LIC2 < 0) LIC2=0 } 
  
  par(mfrow=c(2,1)) 
  amostra1=matrix(rep(seq(1,m),4),m,4,byrow=F) 
  if (n>1) amostra2=amostra1 
  else amostra2=amostra1[-1,] 
  xbar=cbind(rep(x2barras,m),rep(LIC1,m),rep(LSC1,m),media) 
  if (n>1) {
    rbar=cbind(rep(Rbarra,m),rep(LIC2,m),rep(LSC2,m),amplit) 
  } 
  else
    rbar=cbind(rep(Rbarra,m-1),rep(LIC2,m-1),
               rep(LSC2,m-1),amplit) 
  
  matplot(amostra2,rbar,type="o",ylab="Amplitude",
          col=c("black","red","red","blue"),
          ylim=c(min(rbar),max(rbar)),
          xlab="Amostras",main="Gráfico R",pch=20, lty=1, lwd=2)
  
  matplot(amostra1,xbar,type="o",ylab="Média",xlab="Amostras",
          main=expression(paste("Gráfico " * bar(x))),col=c("black","red","red", "blue"),ylim=c(min(xbar),max(xbar)),pch=20, lty=1,lwd=2)
  out=rbind(round(xbar[1,1:3],2),round(rbar[1,1:3],2)) 
  colnames(out)=c("Referência","LIC","LSC")
  out
  }

# m = número de amostras observadas.
# n = número de observações de cada amostra.
# y = vetor com as observações, ordenadas por amostra, da variável observada.
# mi = especificação para a média do processo.
# sigma = especificação para o desvio-padrão do processo.

TCC Robson

dados <- c(311.37,
311.31,
311.38,
311.23,
311.35,
311.39,
313.13,
311.57,
311.13,
312.58,
311.48,
311.19,
311.26,
311.27,
310.98,
311.33,
311.34,
311.36,
311.4,
311.37,
311.23,
311.26,
311.2,
311.15,
311.27,
311.19,
311.31,
311.12,
311.22,
311.18,
311.23,
311.06,
311.1,
311.23,
311.17,
311.07,
310.94,
311.17,
311.13,
311.11,
311.18,
311.21,
310.56,
310.5,
311.13,
310.92,
308.62,
311.19,
311.25,
311.15,
311.16,
310.98,
311.15,
311.1,
311.22,
310.74,
311.23,
311.08,
310.85,
311.13,
310.5,
310.3,
310.79,
310.83,
310.72,
311.1,
311.04,
311.19,
311.31,
311.39,
310.83,
310.79,
310.91,
310.84,
311.18,
311.1,
311.05,
311.36,
311.32,
311.1,
310.76,
310.75,
310.88,
310.9,
310.95,
311.06,
311.01,
311.15,
311.2,
311.33,
311.14,
310.98,
311.04,
311.2,
311.22,
311.24,
311.35,
311.23,
311.09,
311.31,
311.18,
311.15,
311.16,
311.01,
311.41,
311.34,
311.35,
311.31,
311.11,
311.49,
311.25,
311.16,
311.12,
311.26,
311.27,
311.29,
311,
311.44,
311.53,
311.34,
311.22,
311.25,
310.97,
311.2,
311.15,
311.26,
311.36,
311.37,
311.23,
311.29,
311.06,
311.01,
311.15,
311.16,
311.13,
310.97,
311.13,
311.03,
311.25,
311.26,
311.29,
311.18,
311.44,
311.24,
311.31,
311.5,
311.59,
311.61,
311.31,
311.15,
310.75,
310.97,
310.84,
310.86,
310.82,
311.03,
311.32,
310.97,
311.16,
311.21,
311.42,
311.19,
311.19,
311.49,
311.39,
311.41,
311.36,
311.51,
311.36,
311.1,
311.12,
311.12,
311.25,
310.85,
310.9,
310.83,
311.4,
311.38,
311.17,
311.41,
310.76,
310.77,
310.84,
311.2,
311.09,
311.16,
311.11,
311.02,
311.12,
311.05,
311.13,
310.83,
311.15,
311.13,
311.13,
310.95,
310.98,
310.68,
310.68,
310.82,
309.33,
309.25,
308.98,
309.16,
309.25,
309.22,
309.1,
308.58,
308.28,
308.67,
309.14,
308.9,
309.19,
308.98,
308.9,
308.98,
308.98,
308.75,
309.03,
308.94)

m=44
n=5

xR2(m,n,y=dados)

     Referência    LIC    LSC
[1,]     310.95 310.71 311.19
[2,]       0.42   0.00   0.89

m=22
n=10

xR2(m,n,y=dados)

     Referência    LIC    LSC
[1,]     310.95 310.74 311.16
[2,]       0.67   0.15   1.20

# xR1(m,n,y=dados,mi=311,sigma=?)