search for: dgrove

Hi, This report deals with p-values coming from chisq.test using the simulate.p=TRUE option. The issue is numerical accuracy and was brought up in previous bug reports 3486 and 3896. The bug was considered fixed but apparently was only mostly fixed. Just the typical problem of two values that are mathematically equal not ending up numerically equivalent. Consider this series of three 2x2

...library modified to call fast.svd(). This allows it to be used without worrying about whether the matrix is long (cluster chips) or wide (cluster probesets). I've tested this code on R 1.5.1, 1.6.0 and on S-Plus 2000. -Greg > -----Original Message----- > From: Douglas Grove [mailto:dgrove at fhcrc.org] > Sent: Tuesday, October 29, 2002 1:59 AM > To: ripley at stats.ox.ac.uk > Cc: Peter Dalgaard BSA; r-help at stat.math.ethz.ch > Subject: Re: PCA with n << p (was [R] R-1.6.0 crashing on RedHat6.3) > > > > princomp is the wrong tool here: prcomp is better...

Hi, I'm wondering if anyone can point me to a function that will allow me to do a ranking that treats ties differently than rank() provides for? I'd like a method that will assign to the elements of each tie group the largest rank. An example: For the vector 'v', I'd like the method to return 'rv' v: 1 2 3 3 3 4 5 5 6 7 rv: 1 2 5 5 5 6 8 8 9 10 Thanks,

Full_Name: Doug Grove Version: 1.7.0 OS: Linux Submission from: (NULL) (209.31.211.56) Hi, Minor mistake in the documentation on the colSums page. In the ARGUMENTS section it states for 'dims' that: For `col*', the sum or mean is over dimensions `dims+1, ...'; for `row*' it is over dimensions `1:dims'. These two are reversed. Thanks, Doug Grove

Hi, It would be great if someone would add write.delim() as an adjunct to write.table(), just as with write.csv(). I store a lot of data in tab-delimited files and can read it in easily with: read.delim("text.txt", as.is=TRUE) and would love to be able to write it out as easily when I create these files. The obvious setting needed for write.delim() is sep = "\t", but in

Full_Name: Doug Grove Version: 2.3.1 patched (2006-07-03 r38485) OS: Windows Submission from: (NULL) (24.19.37.159) A msg. to R-help prompted me to read the help page for split.screen/screen, etc. I understand what the function is for, but could not find information that described how the argument 'n' in screen/erase.screen/close.screen is mapped to the regions on the device. All we

Full_Name: Doug Grove Version: 2.10.0 10-26-2009 OS: x86_64-unknown-linux-gnu Submission from: (NULL) (140.107.164.65) Hi, Looks like the revisions to help.search didn't get completely finished. When the user specifies the argument 'package' the code doesn't define the variable 'paths' prior to trying to use it, which obviously produces an error. The code dies when

Hi James, Joshua is right, what you're trying to accomplish is quite known in the Java VM world ( http://domino.research.ibm.com/comm/research_people.nsf/pages/dgrove.oopsla01.html ). In order to express the "thunk" code in LLVM you need a full control of how registers are used (because otherwise they would mess up with the arguments). I haven't investigated enough to know if that's possible today in LLVM (I think it wasn't a few years ago...

Hi, I'm trying to understand a behaviour that I have encountered and can't fathom. Here's some code I will use to illustrate the behaviour: # start with some data frame "a" having some named columns a <- data.frame(a=rep(1,3),c=rep(2,3),d=rep(3,3),e=rep(4,3)) # create a subset of the original data frame, but include a # name "b" that is not present in my

hello,everyone. I have one question: example 1 > x=numeric(0) > y=5 > print(x+y) numeric(0) example 2 > x=numeric(1) > y=5 > print(x+y) [1] 5 why the print(x+y) is numeric(0) at the first example, but the result is 0 at the second example? __________________________________________________ 佈伵伝仮伱佲伔佈G佊伿佅佷仯伃佒佇伖侜伒佢佉伝伨侙佄佫伬伂伝侙佊伿伡侢伾仹伻伵伋伂伌侒佊伿佅佷

Hi, I'm having what I believe are OS problems with my installation of R-1.6.0. The machine I'm using was previously running R-1.4.0. The install seemed fine to me. I didn't note an error messages in 'make' or 'make check'. My problem is that I've been running some memory intesive code, and the code that ran under 1.4.0 but won't run under R.1.6.0 (or under

I have what R seems to consider a very large dataset, a 12MB text file of lat,long,and height values, 130,000 rows to be exact. Here's what I get: Thomas Colson North Carolina State University Department of Forestry and Environmental Resources (919) 673 8023 tom_colson at ncsu.edu Calendar: www4.ncsu.edu/~tpcolson

Full_Name: Doug Grove Version: R.1.9.1 OS: Linux Submission from: (NULL) (140.107.156.61) I just found that rank() has a 'decreasing' argument that is not documented in its help page. I checked my version of 2.0.0 (original release hence unpatched) and it is not documented there. For curiousity I also went back to version 1.8.1 and checked the function (not the documentation)and at

Hi all, I am currently building an application based on R 1.7.1 (+ compiled C/C++ code + MySql + VB). I am building this application to work on 2 different platforms (Windows XP Professional (500mb memory) and Windows NT 4.0 with service pack 6 (1gb memory)). This is a very memory intensive application performing sophisticated operations on "large" matrices (typically 5000x1500

...o make upstream changes that are as broadly useful as possible. I've added some notes below about what approach might be taken and I'd be interested to know what you think. Kind regards, Stephen Cross [1] http://loci-lang.org/StructuralTyping.html [2] https://www.research.ibm.com/people/d/dgrove/papers/oopsla01.pdf [3] http://marc.info/?l=llvm-dev&m=129594650017675 ---- Notes about potential approaches ---- In terms of the approach, there are at least two parts to this problem: * Setting/retrieving the hidden value * Forwarding function arguments For the former, two vague ideas I h...

Thanks, that's a good idea - I'll have a look through the VMKit source. -- James On 31 January 2011 21:39, Joshua Warner <joshuawarner32 at gmail.com> wrote: > Hi James, > > I see the problem now. You might look at VMKit (a Java VM build with the > LLVM JIT compiler) - I would expect it uses a similar method for resolving > interface calls (the method, if I

Hi James, I see the problem now. You might look at VMKit (a Java VM build with the LLVM JIT compiler) - I would expect it uses a similar method for resolving interface calls (the method, if I understand it correctly, is well-known in the Java world). I've CC'd the main dev behind VMKit - he might be able to lend some insight. --Joshua On Mon, Jan 31, 2011 at 2:24 PM, James Williams

Hi, What's the best way of dropping leading or trailing blanks from a character string? The only thing I can think of is using sub() to replace blanks with null strings, but I don't know if there is a better way (I also don't know how to represent the trailing blank in a regular expression). Thanks, Doug Grove

I get ties in output from runif() when I generate as few as 10^5 variates and get quite a lot when I generate 10^6. Is this expected?? I haven't seen any duplication with rnorm(10^6), but see varying amounts of duplication using rexp(), rbeta() and rgamma(). I would have thought that there'd be enough precision that one wouldn't get ties until generating samples larger than this..

Full_Name: Doug Grove Version: 1.8.1 Patched (1-11-04) OS: Linux, RedHat 9.0 Submission from: (NULL) (140.107.156.61) The issue may be broader than this, but I'm finding that when I specify that the axis annotation be perpendicular to the axes, and I change the annotation size (using cex.axis) the character size inflation/deflation occurs in only one direction, from the bottom up. More