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

Pedro Emmanuel Alvarenga Americano do Brasil emmanuel.brasil em gmail.com
Quarta Maio 2 22:11:21 BRT 2012


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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120502/52287be4/attachment.html>


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