Veja a *vignette* do pacote SIBER.
On Sat, Feb 8, 2020 at 9:58 PM Luiz Leal por (R-br) <
[email protected]> 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
> [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.
>
_______________________________________________
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.
