[R-br] precisão numérica e which.max com empates

[Tito Conte]

Tito Conte




2012/5/2 Pedro Emmanuel Alvarenga Americano do Brasil <
emmanuel.brasil em gmail.com>

> Amigos de R,
>
> Estou tentando comparar alguns métodos de estimação e um dos passos
> necessários é encontrar valores maximos ou minimos de vetores numéricos.
>
> O problema numero 1 é que quando comparo visualmente os vetores ocorre que
> ocasionalmente acontece 2 ou 3 empates nos valores máximos ou mínimos como
> no vetor abaixo. Reparem que o valor .7 aparece na posição 21 e 23. Mas
> quando utilizo a função which.max ou max eu so encontro a posição 23.
> POssivlemente ha alguma coisa relacionado com precisão nume´rica apesar de
> não fazer muito sentido no momento.
>
> > test.values$Youden [1] 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.55 0.60 0.65 0.60 0.65 0.60 0.65
> [21] 0.70 0.65 0.70 0.65 0.60 0.65 0.60 0.65 0.60 0.55 0.50 0.45 0.40 0.35 0.30 0.25 0.20 0.15 0.10 0.05> max(test.values$Youden)[1] 0.7
>
> > which.max(test.values$Youden)[1] 23
>
> > which(test.values$Youden==max(test.values$Youden)) [1] 23 >
> test.values$Youden[21]==test.values$Youden[23]
>
> [1] FALSE
>
> O problema 2 é que quando de fato ocorrem empates a função which.max retorna somente a primeira posição. Então, a questão
>
> seria como retornar todas as posições com valores máximos ou encontrar uma mediana dos valores empatados?
>
> Qualquer luz é bem vinda.
>
>
> Dr. Pedro Emmanuel A. A. do Brasil
> Curriculum Lattes:  http://lattes.cnpq.br/6597654894290806
> Instituto de Pesquisa Clínica Evandro Chagas
> Fundação Oswaldo Cruz
> Rio de Janeiro - Brasil
> Av. Brasil 4365,
> CEP 21040-360,
> Tel 55 21 3865-9648
> email: pedro.brasil em ipec.fiocruz.br
> email: emmanuel.brasil em gmail.com
>
> ---Apoio aos softwares livres
> www.zotero.org - gerenciamento de referências bibliográficas.
> www.broffice.org ou www.libreoffice.org - textos, planilhas ou
> apresentações.
> www.epidata.dk - entrada de dados.
> www.r-project.org - análise de dados.
> www.ubuntu.com - sistema operacional
>
>
> _______________________________________________
> R-br mailing list
> R-br em 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120502/e370b7c7/attachment.html>