similar to: meaning of "trim" in mean()

Dear useRs, I have a list, each entry of which is a matrix of constant dimensions. Is there a good way (i.e., not using a for loop) to apply a mean to each matrix entry *across list entries*? Example: foo <- list(rbind(c(1,2,3),c(4,5,6)),rbind(c(7,8,9),c(10,11,12))) some.sort.of.apply(foo,FUN=mean) I'm looking for a componentwise mean across the two entries of foo, i.e., the

Dear R-devel, Has anyone played with this? http://www.amd.com/us-en/Processors/DevelopWithAMD/0,,30_2252_2282,00.html <http://www.amd.com/us-en/Processors/DevelopWithAMD/0,,30_2252_2282,00.html> . I'll probably give it a shot... Best, Andy Andy Liaw, PhD Biometrics Research PO Box 2000, RY33-300 Merck Research Labs Rahway, NJ 07065 mailto:andy_liaw@merck.com

Hi All, I want to take a matrix (or data frame) and winsorize each variable. So I can, for example, correlate the winsorized variables. The code below will winsorize a single vector, but when applied to several vectors, each ends up sorted independently in ascending order so that a given observation is no longer on the same row for each vector. So I need to winsorize the variable but

Dear R-devel, I recall that there used to be a mechanism to get around the requirement that all objects in a package have associated documentation; i.e., a way to specify a list of objects (mostly functions) that are not considered as part of the package API. Is this still available? I cannot find any mention in R-exts. I realize that the way to go is to have a name space, and eventually

[I believe few users manipulate formulae and model frames like this, so I'm posting here instead of R-help.] I just found out that model.frame() includes all terms, including those preceded by "-" in the formula, in its output (as mentioned on its help page). Is there a recommended way of excluding "negative" terms? My guess is that I can use info in the terms attribute

Please accept my apologies for cross-posting to those subscribed to both R-help and S-news. Also, please direct *all* inquiries to Vladimir Svetnik, the hiring manager. His contact information is below. Thanks, Bill ################################################################ Job description: Computational statistician/biometrician The Biometrics Research Department at Merck Research

Dear R-help, Has anyone notice the problem that, on Windows (NT and XP), when printing a graph using the "File -> Print..." menu in the graphics window to print the graph, that line width seemed to be ignored in the printed output? For example, if I make a plot with plot(1:10, type="l", lwd=5), it looks right on screen, but when printed out using the menu, it looks like

Dear R-help, I've been compiling R from source on our Linux boxes for quite a while. One thing that bugs me is that I always get an error when make check-all where it choked up running tests/internet.R. That wasn't a big deal, as I can run install.packages()/update.packages() fine, and that's all I use the 'net connection for. The problem is that it takes a _long_ time to fail

Apologies for the cross-posting... Computational statisticians The Biometrics Research Department at Merck Research Laboratories, Merck & Co., Inc. in Rahway, NJ is seeking highly motivated statisticians/data analysts (2 positions) to work in its basic research and drug discovery area. The applicant should have broad expertise in statistics and computer science, with substantial experience

Apology for the cross-post... Andy ============================================== Job description: Computational statistician/medical image analyst The Biometrics Research Department at Merck Research Laboratories, Merck & Co., Inc. in Rahway, NJ is seeking a highly motivated statistician/data analyst to work in its basic research and drug discovery area. The applicant should have broad

Dear R-devel, Kurt had alerted me to the problem that the randomForest package that I maintain has been failing checks in R-devel. However, I just can't see why or where it's failing. I'd very much appreciate any pointer. The failure occur when running the example code in varImpPlot.Rd: > varImpPlot(mtcars.rf) Error in par(opar) : invalid value specified for graphical parameter

Job description: Computational statistician/biometrician The Biometrics Research Department at Merck Research Laboratories, Merck & Co., Inc. in Rahway, NJ, is seeking a highly motivated statistician/data analyst to work in its basic research, drug discovery, preclinical and early clinical development areas. The applicant should have broad expertise in statistical computing. Experience

Dear R-help, Can some one tell me what could be the problem? I downloaded the R-1.7.1 source and try to compile it on WinXP. "make", "make bitmapdll", "make tcl", "make recommended" all ran w/o problem. However, when I ran "make docs", it failed at creating refman.pdf, with no apparent error that I could see. I have attached the log file

Dear R-devel, A while ago Prof. Ripley had mentioned that g77 on AMD64 isn't the greatest, and the Salford Fortran is likely to be better. My question is, has anyone tried compiling R with Salford compilers (or other commercial compilers)? If so, do they offer performance advantages over GCC? TIA for any info. Best, Andy Andy Liaw, PhD Biometrics Research PO Box 2000, RY33-300

Dear R-help: Has anyone tried using R on the the AMD Opteron in either 64- or 32-bit mode? If so, any good/bad experiences, comments, etc? We are considering getting this hardware, and would like to know if R can run smoothly on such a beast. Any comment much appreciated. Best, Andy Andy Liaw, PhD Biometrics Research PO Box 2000, RY33-300 Merck Research Labs Rahway, NJ

I am using R 1.7.1 on Windows XP Pro. The 'overwrite=FALSE' argument appears to be ignored in file.copy(): > file.exists(c("file1.txt", "file2.txt")) [1] TRUE TRUE > file.copy(from="file2.txt", to="file1.txt", overwrite=FALSE) [1] TRUE > and sure enough, file1.txt has been overwritten. Rich Raubertas Biometrics Research, RY33-300 Merck

Dear DevelopeRs, I'd like to ask those who develop R packages with compiled code to please try avoiding dependency on GCC (gcc/g77/gfortran/g++) specific features in the code, for the simple reason that there are non-GCC compilers out there that might choke on such features. I found this out back when I was trying to build R-2.3.0-to-be on our dual Opteron based Scyld cluster. Because the

Dear R-help, I seem to recall that I can use "\t" to get tab in a string on a graphics device, but it doesn't seem to work. Try: lab <- "a\tb\tc" cat(lab, "\n") # works in the console output plot(1:5, main=lab) # no tabs in the title text(3, 3, lab) # no tabs in the text I get the same result both in the windows() and pdf() devices. Any ideas?

When running make check-all on R-1.9.0beta (2004-03-31) built with pgf77/gcc/g++ on the Opteron (running SLES8), I get four warnings about not finding the DESCRIPTION files for the packages in the VR bundle; e.g., cannot open file `/u1/scratch/R-1.9.0/tests/Packages/VR/spatial/DESCRIPTION' Error in tools:::.checkPackageDepends(dir = "/u1/scratch/R-1.9.0/tests/Packages/VR/spatial")

For those interested: I compiled R-1.9.0beta_2004-04-02 with pgf77/gcc/g++ and g77/gcc/g++, both _without_ fast BLAS. Both passed make check-all. Simple crossprod(m) for fairly large `m' shows that the version compiled with pgf77 is about 15% faster: [andy@leo:scratch]% echo "set.seed(1); m <- matrix(runif(5e7), 1e4); gc(); system.time(crossprod(m))" | R-1.9.0-PGI/bin/R