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")

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