[R-br] Identificar valores fora da elipse

Luiz Leal richfield1974 em yahoo.com
Sáb Fev 8 22:58:24 -02 2020


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





-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20200209/7df1ce3a/attachment.html>
-------------- Próxima Parte ----------
Um anexo não-texto foi limpo...
Nome: elipse.xlsx
Tipo: application/x-zip-compressed
Tamanho: 11417 bytes
Descrição: não disponível
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20200209/7df1ce3a/attachment.bin>
-------------- Próxima Parte ----------
Um anexo não-texto foi limpo...
Nome: Lab.xlsx
Tipo: application/x-zip-compressed
Tamanho: 9002 bytes
Descrição: não disponível
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20200209/7df1ce3a/attachment-0001.bin>
-------------- Próxima Parte ----------
Um anexo não-texto foi limpo...
Nome: 1581209735054blob.jpg
Tipo: image/png
Tamanho: 9729 bytes
Descrição: não disponível
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20200209/7df1ce3a/attachment.png>


Mais detalhes sobre a lista de discussão R-br