Pensando que você pode aumentar o número de
classes, a findInterval() é interessante.
z1 <-
rnorm(n=100, mean=0.73436, sd=0.104739)
z2 <- rnorm(n=100, mean=0.69173, sd=0.104492)
z <- c(z1,z2)
status <- rep("sadio", length(z))
i <- 0.69173+c(-1,1)*0.008838306
fi <- findInterval(z, i)
status[fi==1] <- "atacado"
plot(z, col=fi+1)
abline(h=i, lty=2)
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)?
À disposição.
Walmes.