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.
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