[R] order() on vector with Inf's

[Murray Jorgensen]

I have a parameter vector, ibeta, and a corresponding vector of loglikelihoods, llbet. If llbet contains no NAs or Inf's then I can extract the best parameter by

index <- order(-llbet)[1]
beta <- ibeta[index]

or similar. The argument na.last of order() allows me to fix this up even if llbet contains some NAs which I wish to ignore.

Unfortunately my llbet contains Inf's, which are not points of infinite likelihood, if anything they should be -Inf's. Anyway na.last does not seem to help me with these. What should I do?

Murray
--
Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html
Department of Statistics, University of Waikato, Hamilton, New Zealand
Email: [EMAIL PROTECTED] Fax 7 838 4155
Phone +64 7 838 4773 wk +64 7 849 6486 home Mobile 021 1395 862

______________________________________________
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help