similar to: R, AMD Opteron 64, and Rmpi

Hi, You should also check my msc.features.select from caMassClass package. It has feature selection algorithm that I found useful in case of mass-spectra data. It performs individual feature selection and/or removes highly correlated neighbor features. Jarek -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch] Sent: Friday, January

(Apologies for double posting; I sent this to r-help and was suggested that I ask here which I should probable have done to begin with). *************** Dear All, I've got into the habit of installing R from the precompiled Debian binaries, including many of the packages from the r-cran-* Debian packages, and later building from source (e.g., to link against Goto's BLAS, or to build

Dear All, I've got into the habit of installing R from the precompiled Debian binaries, including many of the packages from the r-cran-* Debian packages, and later building from source (e.g., to link against Goto's BLAS, or to build patched versions, etc). I install the newly built R to the very same place (/usr/lib/R). This allows me to build and update R when I wish, AND provides the

Dear All, I've read Thomas Lumley's fortune "If the answer is parse() you should usually rethink the question.". But I am not sure it that also applies (and why) to other situations (Lumley's comment http://tolstoy.newcastle.edu.au/R/help/05/02/12204.html was in reply to accessing a list). Suppose I have similarly called functions, except for a postfix. E.g. f.1 <-

Dear All, I was writing a small wrapper to bootstrap a classification algorithm, but if we generate the indices in the "usual way" as: bootindex <- sample(index, N, replace = TRUE) there is a non-zero probability that all the samples belong to only one class, thus leading to problems in the fitting (or that some classes will end up with only one sample, which will be a problem

If an array has missing values in different rows, plotting using the formul= a=20 interface can produce errors. Example: fake.data <- matrix(rep(-100:100, 4), ncol =3D 4) par(mfrow =3D c(1,2)) boxplot(fake.data ~ col(fake.data)) abline(h =3D 0, lty =3D 2) boxplot(as.data.frame(fake.data)) abline(h =3D 0, lty =3D 2) ##### Add the missing data fake.data[190:200, 1] <-

Dear All, We are about to purchase the second edition of Agresti's "Categorical Data Analysis" (my old copy of the first ed. of that wonderful book is falling apart). I would appreciate suggestions about other comparable books which, if possible, have examples using R/S code (instead of SAS). Thanks, Ram?n -- Ram?n D?az-Uriarte Bioinformatics Unit Centro Nacional de

Dear All, I have been using the randomForest package for a couple of difficult prediction problems (which also share p >> n). The performance is good, but since all the variables in the data set are used, interpretation of what is going on is not easy, even after looking at variable importance as produced by the randomForest run. I have tried a simple "variable selection"

Dear All, I am using the function gls (in the nlme package) and I would like to fit a heteroscedastic model, with different variances for each of the levels of two stratification variables. In p. 210 of Pinheiro & Bates ("Mixed effects models in S and S-Plus", 2000, Springer), the authors show the use of the "*" operator. However, that is not what I want, because it

Dear All, I think there is a typo in the documentation for coxph (library survival). The help says: eps: convergence threshold. Iteration will continue until the relative change in the log-likelihood is less than eps. Default is .0001. However, if I do "coxph.control()" I get: > coxph.control() $eps [1] 1e-09 So the actual eps being used is not 10-4 but

Dear All, I am using fexact.c on a C++ program I wrote. To minimize dependencies on other files (e.g., to not need to include R.h and ctest.h ---now I only include the R files Boolean.h, Constants.h, and Memory.h), I have re-written all declarations of Sints as ints and, what is potentially more serious, I have re-written the line (lines 329 and 330, in fexact.c on R-1.6.1) /* IMAX is the

Dear all, I have a piece of code, call it "FA", that will be called thousands of times in a typical run of function "FB". I can: a) define FA as a function outside of FB (in the global environment), and call it; b) define FA as a function inside the body of FB and call it; c) "expand inline" FA inside FB. FA mainly does data frame subsetting, runs svd's, and

Dear All, I have a list, where several components are NULL, and I'd like to obtain that very same list without the NULL components (i.e., I do not want to unlist or otherwise loose the rest of the list structure). I can do that with a loop, but how could I do it without a loop? Thanks, Ram?n -- Ram?n D?az-Uriarte Bioinformatics Unit Centro Nacional de Investigaciones Oncol?gicas (CNIO)

Dear All, I think I am confused about how I'd be able to use zlib starting with R 2.7.0. I just downloaded the latest development version, built it, etc, but I am not able to find the zlib.h that, I believe, R should place somewhere under "where/you/want/R/to/go" (from --prefix=/where/you/want/R/to/go). There seem to be entry points for zlib in the binary (e.g., grep -r gzopen ./

> x <- matrix(1:10, ncol = 2) > xc <- cor(x[, 2], x[, 0]) > cbind(xc, vector()) Process R floating point exception at Tue Feb 11 19:16:51 2003 > version _ platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor 6.2 year

Dear All, How can the height of a leaf be extracted from a dendrogram? Sure, I can print it, but I am not able to, say, store it in an object. I think I understand that the height is a property of the split, not the leaf itself, but the printing functions display a "height" or "h" (which changes with "hang") and that is what I want. Obviously, the info is there

Dear All, I'd like to be able to have R store (in a list component) a compressed data set, and then write it out uncompressed. gzcon and gzfile work in exactly the opposite direction. What would be a good way to handle this? Details: ---------- We have a package that uses C; part of the C output is a large sparse matrix. This is never manipulated directly by R, but always by the C code.

Dear All, In a package, I am using ".C" to call some C functions. In one case, the number of elements of the return vectors are not known in R before the C call. (Two of the vectors are integers, the third is vector of character strings). Passing from R a vector of the maximum possible size would be a huge waste. I understand one alternative is to use ".Call", but I'd

Dear All, At least in three different GNU/Linux systems, the parallel function from the multicore package leaves defunct (zombie) R processes. For instance library(multicore) parallel(1:5) collect() After collect, the child process (with pid as given by collect) is left defunct. (If we run the last two lines of code again, the previously defunct process will be replaced by the new defunct

Dear All, I am trying to create socket clusters (using snow and snowfall) with a Windows OS. I am running Windows inside VirtualBox and VMware player (i.e., Windows is guest) from a Debian Linux host system (I've tried in two different Linux systems, an AMD x86-64 workstation and an Intel i686 laptop). However, almost always seting up the cluster fails: either R will hang forever or I will