<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Pensando que você pode aumentar o número de classes, a findInterval() é interessante.<br><br><span style="font-family:courier new,monospace">z1 <- rnorm(n=100, mean=0.73436, sd=0.104739)<br>
z2 <- rnorm(n=100, mean=0.69173, sd=0.104492)<br>z <- c(z1,z2)<br><br>status <- rep("sadio", length(z))<br>i <- 0.69173+c(-1,1)*0.008838306<br>fi <- findInterval(z, i)<br>status[fi==1] <- "atacado"<br>
<br>plot(z, col=fi+1)<br>abline(h=i, lty=2)<br></span><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Vai uma dúvida minha que tá fora do escopo. Você tá classificando novas observações (observações futuras) pelo intervalo de confiança (para um parâmetro) gerado em uma análise. Talvez isso não seja adequado. O nível de confiança desse intervalo (95%) não significa que 95% das observações serão classificadas como "atacada". 95% é o nível de cobertura do intervalo que representa a probabilidade dele conter o parâmetro. Não seria mais apropriado um intervalo de predição (para observação futura)?<br>
<br>À disposição.<br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Walmes.<br></div></div>