similar to: FW: sum overflow (PR#1091)

In a string context, grep, gsub, sub are improperly treating NA (missing) as the string "NA", and returning unexpected results > grep("A", c(NA,"NA")) [1] 1 2 # expected: # [1] 2 > gsub("A", "X", c(NA,"NA")) [1] "NX" "NX" # expected # [1] NA "NX" > sub("A", "X",

Is there a way to tell how much memory the computer running R has? -Frank -----Original Message----- From: David Smith [mailto:dsmith at insightful.com] Sent: Monday, February 09, 2004 1:32 PM To: Ross Boylan Cc: r-help Subject: RE: [R] R does in memory analysis only? Ross Boylan writes: > R works only on problems that fit into (real or virtual) memory. > ... does S-Plus have the same

Third try... > -----Original Message----- > From: Warnes, Gregory R > Sent: Tuesday, February 05, 2002 4:12 PM > To: 'R-bugs' > Subject: layout and piechart diameter problem > > > I've been using layout to create some graphics pages which include pie > charts. (NB: No piechart arguments please, the main chart on the page is > a proper bar chart

This message is in MIME format. Since your mail reader does not understand this format, some or all of this message may not be legible. ------_=_NextPart_000_01C3BE9C.3FB5BC1A Content-Type: text/plain; charset="iso-8859-1" There is a small bug in the implementation of \preformatted in R 1.8.0: The closing "\end{verbatim}" needs to be followed by a newline, since my latex

> -----Original Message----- > From: Warnes, Gregory R > Sent: Friday, January 16, 2004 3:52 AM > To: 'r-packages'; 'python-announce '; 'rsoap-talk'; 'Dan Nathan'; 'RPy > List' > Subject: Announce: RPy version 0.3.3 > > > RPy Version 0.3.3 is now available from the RPy home page at > http://rpy.sf.net. > > What

> -----Original Message----- > From: Warnes, Gregory R > Sent: Friday, January 16, 2004 3:52 AM > To: 'r-packages'; 'python-announce '; 'rsoap-talk'; 'Dan Nathan'; 'RPy > List' > Subject: Announce: RPy version 0.3.3 > > > RPy Version 0.3.3 is now available from the RPy home page at > http://rpy.sf.net. > > What

Package: mcgibbsit Title: Warnes and Raftery's MCGibbsit MCMC diagnostic Version: 1.0 Author: Gregory R. Warnes <gregory_r_warnes at groton.pfizer.com> Description: mcgibbsit provides an implementation of Warnes & Raftery's MCGibbsit run-length diagnostic for a set of (not-necessarily independent) MCMC sampers. It combines the estimate error-bounding approach of Raftery

Package: mcgibbsit Title: Warnes and Raftery's MCGibbsit MCMC diagnostic Version: 1.0 Author: Gregory R. Warnes <gregory_r_warnes at groton.pfizer.com> Description: mcgibbsit provides an implementation of Warnes & Raftery's MCGibbsit run-length diagnostic for a set of (not-necessarily independent) MCMC sampers. It combines the estimate error-bounding approach of Raftery

Version 0.8.0 of the gregmisc package is now (or will shortly be) available on CRAN. New in this release: - Enhanced and bug-fixed 'CrossTable' function (contributed by Marc Schwartz) - Augmented the 'barplot2' function with an 'add' argument to allow for the addition of a barplot to an existing graphic. (contributed by Marc Schwartz) - Added the

I've created a Perl script that translates Microsoft Excel (.xls) files into comma-delimited text files (.csv) using the Perl Spreadsheet::ParseExcel module. Usage ----- perl xls2csv.pl <excel file> [<output file>] [<worksheet number>] Translate the Microsoft Excel spreadsheet file contained in <excel file> into comma separated value format (CSV) and store in

It turns out that my web host doesn't allow direct links. Try going to http://www.warnes.net/downloads or getting it from the sourceforge snippet library https://sourceforge.net/snippet/detail.php?type=snippet&id=101149. -Greg > -----Original Message----- > From: Warnes, Gregory R [mailto:gregory_r_warnes at groton.pfizer.com] > Sent: Monday, May 13, 2002 10:15 AM > To:

Dear Greg, Thanks for the new release. The decomposition of the SSQ is just what I need! Regards, Martin. Martin Hoyle, School of Life and Environmental Sciences, University of Nottingham, University Park, Nottingham, NG7 2RD, UK Webpage: http://myprofile.cos.com/martinhoyle >>> gregory_r_warnes at groton.pfizer.com 10/30/02 07:16PM >>> Version 0.7.3 of the gregmisc package

I got an email error message when I attempted to send this from my work account. I have manually added it to the bug tracker, and am resending from my personal account. -G On 12/12/05, Warnes, Gregory R <gregory.r.warnes@pfizer.com> wrote: > > > > > -----Original Message----- > > From: Warnes, Gregory R > > Sent: Monday, December 12, 2005 1:53 PM >

Full_Name: Bill Simpson Version: 1.2.3 OS: linux Submission from: (NULL) (193.62.250.209) I think I have found an accuracy problem with sum() > a<-1:1000 > b<-1:1000 > > crossprod(a,a)*crossprod(b,b) [,1] [1,] 1.114448e+17 > drop(.Last.value) [1] 1.114448e+17 > #this should give same answer as above > sum(a*a)*sum(b*b) [1] -652010736 #seems to be due to

Hi All, This is off topic, but we're drawing a blank here.. > In a presentation I'll be giving next week, I want to include a reference > to a good general text on computing sample sizes for standard experiments. > Can anyone recommend a good book to use for this purpose? > > Thanks, > > -Greg LEGAL NOTICE\ Unless expressly stated otherwise, this

When the R_HOME environment variable is set, R issues a warning messages when it starts up. When doing 'make check' this causes the check for eval-etc.R to fail: [...] make[3]: Leaving directory `/Volumes/app/R/src/R-1.4.0_patched_2002-01-24/tests' running strict specific tests make[3]: Entering directory `/Volumes/app/R/src/R-1.4.0_patched_2002-01-24/tests' running code in

Thanks to David Meyer [david.meyer@ci.tuwien.ac.at] and David Brahm [brahm@alum.mit.edu] who suggested: m[ cbind(index.i, index.j) ] <- vals This works fine for the example I gave. Unfortunately, this approach doesn't extend to using the row and column names to make assignments: > m <- matrix("",ncol=3,nrow=3) > dimnames(m) <-

I recently wrote a function to compute functions over moving windows. (See below). It varies from your request by using the last N points rather than the last N/2 and next N/2 points, so you may have to modify it to get what you want. It will be in the next gregmisc package, which will go out in a month or so when the publication review process here finishes. (I hate red tape!) -Greg

As I've mentioned a number of times. I find it very useful to have lowess() become a generic function so that a lowess.formula() can be defined. Below is a patch that makes both changes, as well as updating the corresponding help documentation. Gregory R. Warnes Manager, Non-Clinical Statistics Pfizer Global Research and Development Tel: 860-715-3536 ? DESCRIPTION ? Makefile ?

What I am doing : create a simple correlation matrix on 41 variables, then plot an heatmap with this program : library(gplots) mat<-cor(temp.alln,use="pairwise.complete.obs") hm<-heatmap.2(mat,symm=T) HM<-format(round(mat[hm[[1]],hm[[2]]],2)) library(RColorBrewer) brewer.pal(10,"Spectral")->mp heatmap.2(mat, symm = TRUE, col = mp,