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")
elipse.xlsx
Description: Zip compressed data
Lab.xlsx
Description: Zip compressed data
_______________________________________________ R-br mailing list [email protected] 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.
