Exemplo

Um procedimento de controle de qualidade foi planejado para garantir um máximo de \(10\%\) de peças defeituosas na produção. A cada 5 minutos, sorteia-se uma amostra de 20 peças e, havendo mais de \(15\%\) defeituosas, encerra-se a produção para verificações.

Objetivo

Desejamos saber se a parada foi desnecessária e o momento da parada.

simula<-function(n,theta){
   media<-0; total<-0 ;defeito<-0 ; count<-0

   while (media<=0.15){
      x<-rbinom(n,1,theta)
      media <- mean(x)
      defeito <- defeito + sum(x)
      total <- total + length(x)
      count<-count + 1 
    }
 
   if ( (defeito/total)<= theta ){
       return(list("Parada Desnecessária", 
                    proporcao=defeito/total,
                    p.hat=media,
                    qtd_sorteios=count))
      }else{
         return(list("Parada Necessária", 
                      proporcao=defeito/total,
                      p.hat=media,
                      qtd_sorteios=count))
   }
}
simula(20,0.10)
## [[1]]
## [1] "Parada Desnecessária"
## 
## $proporcao
## [1] 0.08
## 
## $p.hat
## [1] 0.2
## 
## $qtd_sorteios
## [1] 10
simula(20,0.10)
## [[1]]
## [1] "Parada Necessária"
## 
## $proporcao
## [1] 0.13
## 
## $p.hat
## [1] 0.3
## 
## $qtd_sorteios
## [1] 5