similar to: Suggestion for tapply

Does anyone have something like tapply that is extremely fast for matrices when there is a very large number of levels of the grouping variable? I'm referring to, for example, tapply(x, grouping.variable, function.operating.on.submatrix) where x is a matrix and the submatrix is a subset of the rows of x. The grouping variable's length equals the number of rows of x. -- Frank E

I am trying to replicate part of the function of the S-Plus terms.inner function. I don't need an entire model terms object but R expressions or character strings containing the innermost variable name given a vector of character strings or expressions. Here are some example inputs that each should result in the string or expression "x": x; x^2; g(x); h(g(x)); h(g(x^3));

I am using read.spss in the foreign package to read an SPSS save file. For date variables I get huge values such as 11489990400. Does anyone know how to convert these values to R POSIXct date objects? Thanks in advance -Frank platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor

The Hmisc and Design libraries have been updated respectively to versions 1.4-2 and 1.1-1. New versions for Linux/Unix/Windows may be obtained from http://hesweb1.med.virginia.edu/biostat/s/library/r . Web sites for the libraries are http://hesweb1.med.virginia.edu/biostat/s/Hmisc.html and http://hesweb1.med.virginia.edu/biostat/s/Design.html . Thanks to Xiao Gang Fan for porting the libraries

The Hmisc and Design libraries have been updated respectively to versions 1.4-2 and 1.1-1. New versions for Linux/Unix/Windows may be obtained from http://hesweb1.med.virginia.edu/biostat/s/library/r . Web sites for the libraries are http://hesweb1.med.virginia.edu/biostat/s/Hmisc.html and http://hesweb1.med.virginia.edu/biostat/s/Design.html . Thanks to Xiao Gang Fan for porting the libraries

On Mon, 16 Jun 2003 09:56:07 +0200 Patrick Hausmann <p.hausmann at mr-regionalberatung.de> wrote: > Dear Frank, > > I can reproduce the problem from Philippe Glaziou (see attachment). > Maybe this can help. > Best, > Patrick > > Am Wollelager 11 > 27749 Delmenhorst > Tel. 04221 96373-0 > Fax 04221 96373-29 > > http://www.mr-regionalberatung.de >

In platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor 7.1 year 2003 month 06 day 16 language R I get > sort(c(3,1,NA)) [1] 1 3 Shouldn't NAs be retained by default? Thanks -Frank ---

In a .First.lib I want to issue two require()s to insure that two other packages are loaded. But I want the package being loaded by .First.lib using library.dynam("mypackage",pkb,lib) to be higher in the search order than the two required packages, because I want to have a couple of functions from the two required packages overridden. What is the best way to do that? Thanks in advance

Dear R Developers, I should have had the Design library running in R by now but have kept putting off changing some calls to LINPACK routines to use those builtin to R. Specifically I call dqrsl1 and dqr. Would it be an easy task to put those in the next release of R? If not I'll finally bite the bullet and get back into reading LINPACK documentation (which I have but haven't examined

Using R CMD INSTALL packagename.tar.gz on platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor 6.1 year 2002 month 11 day 01 language R I used to get html and latex files recreated only when the source .Rd file changed. Now I am getting html and latex (and text and example) recreated whenever I run R CMD INSTALL when I do

Using R CMD INSTALL packagename.tar.gz on platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor 6.1 year 2002 month 11 day 01 language R I used to get html and latex files recreated only when the source .Rd file changed. Now I am getting html and latex (and text and example) recreated whenever I run R CMD INSTALL when I do

Thanks to Deepayan Sarkar (lattice) and Paul Murrell (grid) I have nearly finished porting the Hmisc library's xYplot function to R. xYplot extends xyplot to allow for error bars, bands, and a few other features. In S-Plus I extended the Trellis function xyplot mainly by writing panel.xYplot. In R, xyplot contains a lot of code, so I implemented xYplot by calling xyplot from it. Things are

I have many functions for high-level plotting that early on in their code do something like strwidth('string', units='inches') before plotting. I do this before plot() to set the correct margins with par(mai=...), to leave room for wide text on the left or right margins. If the plotting device is not opened strwidth( ) causes an error that plot.new() needs to be called. I can

Full_Name: Frank Harrell Version: 1.6.1 OS: RedHat 8.0 Linux Submission from: (NULL) (128.143.108.90) I get an error when trying to load a URL that contains a file that was saved using save(object, compress=TRUE): > load('http://hesweb1.med.virginia.edu/biostat/s/data/sav/kprats.sav') Error in gzfile(file, "rb") : unable to open connection In addition: Warning message:

Great news! Having spent a significant amount of time porting/mantaining some packages on S4/Splus 5.x and R, I'm interested in your experience as you port Hmisc and Design; I'd like to know of your progress in order to help improve available R/Splus portability tools both at the C and the S language code. These tools include the package SLanguage written by John Chambers, and

If anyone has a very fast vectorized method for doing the following I would appreciate some help. I want to avoid outer() to limit memory problems for very large n. Let x = real vector of length n y = real vector of length n w = real vector of length m, m typically less than n/2 but can be > n z = real vector of length m For w[i], i=1,,,m, find the value of x that is closest to w[i]. In

Under platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor 6.1 year 2002 month 11 day 01 language R and using foreign 0.5-8 I am encountering errors when using read.xport. Here's

The Design library has been fully ported to R except for Cox proportional hazards regression modeling (using Therneau's survival package) which will be available in about two weeks. It will take much longer to make all the example code executable, is it currently contains many examples for which data are not provided. Thanks to Xiao Gang Fan <xiao.gang.fan1 at libertysurf.fr> who

The Design library has been fully ported to R except for Cox proportional hazards regression modeling (using Therneau's survival package) which will be available in about two weeks. It will take much longer to make all the example code executable, is it currently contains many examples for which data are not provided. Thanks to Xiao Gang Fan <xiao.gang.fan1 at libertysurf.fr> who

I am running Linux Mandrake 7.2 with the following: platform i586-pc-linux-gnu arch i586 os linux-gnu system i586, linux-gnu status major 1 minor 2.3 year 2001 month 04 day 26 language R To build a package I ran R CMD build Hmisc, apparently with success: * checking for file `Hmisc/DESCRIPTION' ... OK * preparing `Hmisc': * cleaning src * checking