similar to: max.col specification

Hi, I think there's some rounding issue with returning the max column. (running 2.9.0 on an Apple, but my buddy found it on his PC) > x <- matrix(c(1234.568,1234.569,1234.567),1) > max.col(x) [1] 2 > x <- matrix(c(12345.568,12345.569,12345.567),1) > max.col(x) [1] 3 > x <- matrix(c(112345.568,112345.569,112345.567),1) > max.col(x) [1] 3 > max.col(-x) [1]

I've noticed that the max.col function with the default "random" option often gives unexpected results. For instance, in this test, it seems clear what the answer should be: > # second col should always be max > x1 = cbind(1:10, 2:11, -Inf) > > # this works fine > max.col(x1, "first") [1] 2 2 2 2 2 2 2 2 2 2 > > # this gives random answers >

I've noticed that the max.col function with the default "random" option often gives unexpected results. For instance, in this test, it seems clear what the answer should be: > # second col should always be max > x1 = cbind(1:10, 2:11, -Inf) > > # this works fine > max.col(x1, "first") [1] 2 2 2 2 2 2 2 2 2 2 > > # this gives random answers >

>From R-help, subject "Variable number of loops" I've opened a new thread, moving from R-help to R-devel .. >>>>> Jim Lemon <drjimlemon at gmail.com> >>>>> on Sun, 17 May 2015 09:19:06 +1000 writes: > Hi all, Given the number of help requests that involve > permutations/combinations, and the less than obvious > naming

-------- Original Message -------- Subject: Re: [R] Simulation - Natrual Selection Date: Wed, 05 Jan 2011 17:24:05 +0000 From: Ben Ward <benjamin.ward@bathspa.org> To: Bert Gunter <gunter.berton@gene.com> CC: Mike Marchywka <marchywka@hotmail.com> On 05/01/2011 17:08, Bert Gunter wrote: > Couple of brief comments inline below. -- Bert > > On Wed, Jan 5, 2011 at

I was surprised to find that df$a[1] is an order of magnitude faster than df[1,"a"]: > df <- data.frame(a=1:10) > system.time(replicate(100000, df$a[3])) user system elapsed 0.36 0.00 0.36 > system.time(replicate(100000, df[3,"a"])) user system elapsed 4.09 0.00 4.09 A priori, I'd have thought that combining the row and column

I am trying to run, in R -Windows95, the example in MASS2, section 17.2. When I run predplot() on cush.lda I get > predplot(cush.lda, "LDA") Error: couldn't find function "max.col" I get the plot OK, but within predict(), it balks on max.col A more general question about using VR. Is there an easy way to use library on VR in general or do I have to move each into

Dear R-Developers, I think I found two bugs in max.col(). Ties between zeros are not broken, which might affect simulations. -Inf and Zero can be treated the same, which can give completely wrong results, e.g. when the second max is sought by replacing all maxs by -Inf. To fix max.col I do offer the C-code behind my function rowMax(), which also handles NAs and seems to be faster. However,

There is no function called 'maxCol' in 2.7.0, so should 'maxCol' be purged from maxCol.Rd and 'max.col' be promoted from an alias (and the file renamed)? from today's R-devel: /usr/local/src/R-devel> head ./src/library/base/man/maxCol.Rd % File src/library/base/man/maxCol.Rd % Part of the R package, http://www.R-project.org % Copyright 1995-2007 R Core

Dear CentOS community Greetings, i was waiting for few days before sending this email message, i thought some server may have problem. That might causing the issue, will resolve its own. but since a week now i am getting this message. because i have not come across this error before, is anyone have across what is it? why my server is unable to download the file ERROR: http: GET

Thanks Mark, That was it. (ie the --binary option on R CMD build) I'll just state for the web-o-sphere (just in case it helps any one else) that part of my confusion in building this came from misunderstanding what they call "source" and what the call "binary" in this environment. I always think of source as something that must be compiled into a binary (like C-code

> > > On the original question though, why isn't there something "off the > shelf" > > that will do what I want? Surely, a "boxplot" using mean, SD, max and > min > > would be a common enough need to justify it? > > > > Gabor Grothendieck replied: > tarr is not a list or a data frame. Use.data.frame(tarr) so that it > uses

Hello, I think the following qualify as a bug given: x<-c(1,2,3,4,2,4,2,2,4,2,2,2,4,3,2,1) z<-embed(x,4) z [,1] [,2] [,3] [,4] [1,] 4 3 2 1 [2,] 2 4 3 2 [3,] 4 2 4 3 [4,] 2 4 2 4 [5,] 2 2 4 2 [6,] 4 2 2 4 [7,] 2 4 2 2 [8,] 2 2 4 2 [9,] 2 2 2 4 [10,] 4 2 2 2

Hello, I think the following qualify as a bug given: x<-c(1,2,3,4,2,4,2,2,4,2,2,2,4,3,2,1) z<-embed(x,4) z [,1] [,2] [,3] [,4] [1,] 4 3 2 1 [2,] 2 4 3 2 [3,] 4 2 4 3 [4,] 2 4 2 4 [5,] 2 2 4 2 [6,] 4 2 2 4 [7,] 2 4 2 2 [8,] 2 2 4 2 [9,] 2 2 2 4 [10,] 4 2 2 2

Hi, I'm attempting to run the Aeris Dive Computer interface (ACI V2.x) using WINE under Ubuntu 9.04 Jaunty Jakalope release on an ASUS 1000HE netbook. However, it doesn't seem to be working very well. the install process seemed to work fine. It successfully created an icon under the WINE application for the interface. On startup of ACI, it says that it cannot read the previous dive

Hi everybody, I have been struggling to get my Windows XP box to be able to access my samba server. When samba starts, I see the following my my /var/log/samba/log.smb file: smbd version 3.0.25b-2.fc7 started. Copyright Andrew Tridgell and the Samba Team 1992-2007 [2007/07/31 20:22:43, 0] auth/auth_util.c:create_builtin_administrators(792) create_builtin_administrators: Failed to

None of Andy's comments) are inconsistent with the point that rlm() and lqs(), if they disagree with lm(), likely offer better places to start than does lm(), in identifying points that should be examined as in some sense outliers. All such methods are to be used, not as crutches, but as sources of insight. Incidentally, the "rlm" class inherits from "lm", and plot.lm()

Please, I am receiving lot of e-mails that I do not want. Please could you delete my e-mail. Thank, Judit -----Original Message----- From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On Behalf Of Thomas Adams Sent: Monday, December 12, 2011 3:22 PM To: Bert Gunter Cc: r-help at r-project.org Subject: Re: [R] Boxplot of multiple vectors with different lengths Bert,

Im trying to figure out the base way to check the callerID being sent to my Asterisk box and use it if it is a valid NANP number, but replace it with a static NANP number if it is not. (Why? I have a few carriers that require this, and a few international users - if it happens to take one of the carriers that require it, I want it to set a static number that is valid). I'm playing

I have two vectors, a and b. b is a text file. I want to find in b those elements of a which occur at the beginning of the line in b. I have the following code, but it only returns a value for the first value in a, but I want both. Any ideas please. a = c(2,3) b = NULL b[1] = "aaa 2 aaa" b[2] = "2 aaa" b[3] = "3 aaa" b[4] = "aaa 3 aaa"