On Mar 4, 2010, at 9:29 AM, Tom Short wrote: > Is the bug-reporting system working okay? Two days ago, I submitted > the following bug report via email to r-b...@r-project.org. I didn't > see a confirmation, and it didn't see it at > http://bugs.r-project.org/. Now, http://bugs.r-project.org/ seems to > be down. > > Anyway, here's the bug report related to sort.list and sort(..., > index.return = TRUE) with na.last = NA > > > I think that both sort.list(x, na.last = NA) and sort(x, na.last = NA, > index.return = TRUE)$ix give incorrect answers with na.last. With > na.last, both of these return answers equivalent to > sort.list(na.omit(x)), and I think they should be the equivalent of > order(x, na.last=NA) as follows. > >> x <- c(1L, 6L, NA, 2L) >> order(x, na.last = NA) # right > [1] 1 4 2 >> sort.list(x, na.last = NA, method = "radix") # wrong, I think > [1] 1 3 2 >> sort.list(x, na.last = NA, method = "quick") # wrong > [1] 1 3 2 >> sort(x, na.last = NA, index.return = TRUE)$ix # wrong > [1] 1 3 2 > > I've included a patch for the "radix" and "shell" methods of > sort.list. The sort and "quick" method of sort.list (which uses sort) > look more challenging. With the patch, I get: > >> sort.list(x, na.last = NA, method = "radix") # good > [1] 1 4 2 >> sort.list(x, na.last = NA, method = "shell") # good > [1] 1 4 2 >> sort.list(x, na.last = NA, method = "quick") # still bad > [1] 1 3 2 > > By the way, having the radix sort is great. It's really fast for factors. > > - Tom > > Tom Short > Electric Power Research Institute (EPRI) > > > >> sessionInfo() > R version 2.10.1 (2009-12-14) > i386-pc-mingw32 > > locale: > [1] LC_COLLATE=English_United States.1252 > [2] LC_CTYPE=English_United States.1252 > [3] LC_MONETARY=English_United States.1252 > [4] LC_NUMERIC=C > [5] LC_TIME=English_United States.1252 > > attached base packages: > [1] stats graphics grDevices utils datasets methods base
Tom, See this recent follow up posting from Peter: http://tolstoy.newcastle.edu.au/R/e9/devel/10/02/0375.html There has been talk over the years of moving to Bugzilla, but I am not clear on present status. Perhaps the link on the main R Project page needs to be removed or better, updated to a link with a status update on the R bug reporting process. Of course, that does not help folks using bug.report(), which presumably needs to be updated as well. HTH, Marc Schwartz ______________________________________________ R-devel@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-devel