similar to: R on AMD64 (Opteron)

> From: Douglas Bates [mailto:bates at bates4.stat.wisc.edu] > > "Liaw, Andy" <andy_liaw at merck.com> writes: > > > Sorry. I need to retract my claim. There seems to be a 3G > limit, even > > though the OS could handle nearly 8G. (I can have two > simultaneous R > > processes each using near 3G.) > > > > On another note, on

> From: Douglas Bates [mailto:bates at bates4.stat.wisc.edu] > > "Liaw, Andy" <andy_liaw at merck.com> writes: > > > Sorry. I need to retract my claim. There seems to be a 3G > limit, even > > though the OS could handle nearly 8G. (I can have two > simultaneous R > > processes each using near 3G.) > > > > On another note, on

Using ATLAS with R is an old topic quite covered in the "R Administration" manual (and by R's "configure" script collection). I still do not easily manage to build R properly on our new AMD Opteron (2-processor). I did work with the current Atlas 3.6.0, configured manually (but "express" version) to build a threaded ATLAS version, and successfully ran Atlas'

Hi Martin, When I attended the LinuxWorld Expo in NYC back in January, I chatted with some folks at the AMD booth, as well as guys from Penguin Computing (where we bought our Opteron box). I was told that the Operton has this somewhat strange setup that the memory is controlled by one CPU. The net effect of this being that when both CPUs are running, one might only be running at around 90%

Full_Name: Jerome Asselin Version: 1.6.2 OS: RedHat Linux 7.2 Submission from: (NULL) (142.103.173.179) Should be an easy fix... Consider the examble below: plot(0,0) legend(0,0,c("Hello!","Hi!"),pch=1:2,lty=1:2,trace=T) It gives the following trace: > plot(0,0) > legend(0,0,c("Hello!","Hi!"),pch=1:2,lty=1:2,trace=T) xchar= 0.05178 ;

Hi, is it possible to get the benchmark-results from anybody is using a common amd dual opteron with windows or linux, R.1.8.1 or 1.9.0? http://www.sciviews.org/other/benchmark/R2.R many thanks & regards, christian

Hi, Is there any simple solution to get ranks in descending order? Example, a <- c(10, 98, 98, 98, 99, 100) r <- rank(a, ties.method="average") produces 1 3 3 3 5 6 I would want this instead: 6 5 3 3 3 1 Note that reversing r doesn't work but in small examples. Thanks, -Jose -- jquesada at andrew.cmu.edu Research associate http://lsa.colorado.edu/~quesadaj Dept. of

Greetings - A quick perusal of some of the posts to this maillist suggest the level of the questions is probably beyond someone working at my level, but at the risk of looking foolish publicly (something I find I get increasingly comfortable with as I get older), here goes: My research group recently purchased a multi-Opteron system (bunch of 880 chips), running 64-bit RHEL 4 (which we have

Full_Name: Jerome Asselin Version: 1.6.1 OS: RedHat Linux 7.2 Submission from: (NULL) (142.103.173.179) Hello, I get a segmentation fault when I run the following code. I wouldn't expect meaningful results because my response variable contains only missing values. However, I would expect something like a regular error (not a segmentation fault). library(survival) data <-

Full_Name: Jerome Asselin Version: 1.7.1 OS: Red Hat Linux 7.2 Submission from: (NULL) (24.77.125.119) Setting the parameter na.action=na.omit should remove incomplete records in princomp. However this does not seem to work as expected. See example below. Sincerely, Jerome Asselin data(USArrests) princomp(USArrests, cor = TRUE) #THIS WORKS USArrests[1,3] <- NA princomp(USArrests, cor =

Full_Name: Jerome Asselin Version: 1.8.0 OS: RedHat Linux 7.2 Submission from: (NULL) (142.103.177.13) I would not expect a segmentation fault; perhaps an error message. > .Call("log") Segmentation fault This is always reproducable for me. Sincerely, Jerome Asselin

Is there a simpler way then the solution to the one that was posted here? I'm not very proficient with legend, and I don't understand this solution. All I have is two or more lines on one plot that I want to put a legend on and I can't figure out how to do it from the examples. Can you give a very simple example? It does not have to be fancy!! I have never worked with a

I am trying to fit a lognormal distribution on interval-censored data. Some of my intervals have a lower bound of zero. Unfortunately, it seems like survreg() cannot deal with lower bounds of zero, despite the fact that plnorm(0)==0 and pnorm(-Inf)==0 are well defined. Below is a short example to reproduce the problem. Does anyone know why survreg() must behave that way? Is there an alternate

Hi, I'm trying to use a similarity matrix (triangular) as input for pam() or fanny() clustering algorithms. The problem is that this algorithms can only accept a dissimilarity matrix, normally generated by daisy(). However, daisy only accept 'data matrix or dataframe. Dissimilarities will be computed between the rows of x'. Is there any way to say to that your data are already a

Full_Name: Jerome Asselin Version: 1.6.2 OS: RedHat Linux 7.2 Submission from: (NULL) (142.103.173.179) In the following example, the line type lty=2 does not show properly across the entire line. x <- c(seq(0,.5,.001),seq(.6,1,.1)) y <- rep(1,length(x)) plot(x,y,type="s",lty=2) Sincerely, Jerome Asselin

Dear All, I am given a set of files names as: velocity1.txt velocity2.txt and so on. I am sure there must be a way to read them automatically in R. It is really taking me longer to read them than to analyze them. Anybody has a suggestion to help me out with this? Many thanks Lorenzo

Hi all, In the (very simple) example below, I have defined a random effect for the residuals in lme(). So the model is equivalent to a FIXED effect model. Could someone explain to me why lme() still gives two standard deviation estimates? I would expect lme() to return either: a) an error or a warning for having an unidentifiable model; b) only one standard deviation estimate. Thank you for your

Dear all, I am managing a departmental purchase of an Opteron based workstation/server for scientific computing on which we will be running R. The environment will probably be either Unix/Linux or Solaris and the amount of RAM will be 8-16Gb, depending on the number of processors. My main concerns are the following: 1. How much does R benefit from passing from one processor to two/four

Full_Name: Jerome Asselin Version: 1.3.0 OS: Windows 98 Submission from: (NULL) (24.77.112.193) I am having accuracy problems involving the computation of inverse of nonnegative definite matrices with solve(). I also have to compute the Choleski decomposition of matrices. My numerical problems involving solve() made me find a bug in the chol() function. Here is an example. #Please, load the

Dear all, Has anyone seen the same problem? I tried compiling R-1.9.0 alpha (2004-03-16) on our Opteron box running SUSE Linux ES8. I ran: ./configure --enable-R-shlib --with-blas=goto and got: Source directory: . Installation directory: /usr/local C compiler: gcc -m64 -O2 -g -msse2 -march=k8 -Wall --pedantic C++ compiler: g++ -m64 -O2 -g