Veja a *vignette* do pacote SIBER. On Sat, Feb 8, 2020 at 9:58 PM Luiz Leal por (R-br) < r-br@listas.c3sl.ufpr.br> wrote:
Prezados, como faço para criar um script que identifique que os laboratórios 7042 e 9541 estão fora da elipse. Os dados seguem em anexo. Desde já agradeço. Luiz
# Dados setwd("F:\\") library(xlsx) dados=read.xlsx("Lab.xlsx" , sheetIndex=1,encoding="UTF-8") elipse=read.xlsx("elipse.xlsx" , sheetIndex=1,encoding="UTF-8") x=dados$x y=dados$y
# Graph limitex = c( min(c(x,elipse[,1])), max(c(x,elipse[,1])) ) limitey = c( min(c(y,elipse[,2])), max(c(y,elipse[,2])) ) plot(x, y, type = "n", xlab = "", xlim = limitex, ylim = limitey, ylab = "", cex.axis = 1.3, cex.lab = 1.3) text(x, y, dados$Lab, font = 2, cex = 0.8) points(elipse, type = "l", col = "red")
[image: Inline image]
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
