[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