similar to: Benchmarking R, why sort() is so slow?

Dear all, I am trying to get an estimate of the intercept for a linear model. In this case, I know the slope of the model, can anyone tell me how to constrain the formula in lm() so that it only estimates the intercept not the slope? Many thanks in advance, Sincerely, Liqing Zhang Dept. of Eco. Evol. Biol. Univ. of CA, Irvine email: lzhang at uci.edu >From VM Mon Apr 30 08:18:45 2001

Hi, I have just seen Stefan Steinhaus' web page : http://www.uni-franfurt.de/~stst/ncrunch.html I think it would be nice to include "R" as well. I have taken Forrest Young's email on stat-lisp list and changed the stuff for "R" :) Here it is: (someone please check this so we can also send it to Stefan Steinhaus.

Hello, I am updating my benchmark (http://www.sciviews.org/other/benchmark.htm) to recent versions of data analysis software (including R 1.6.0 and Splus 6.1), and I now run it on a Pentium IV instead of the old Celeron 500 Mhz that candidates for retirement. I test R under Windows Xp pro with and without optimized BLAS. I use the optimized Rblas.dll for P4 found on CRAN. Here are the results.

Hi All, I've just rebuild the latest R with the Goto BLAS on our new Intel quad core machines. I did a few basic matrix calculations, and I was very impressed by the performance I saw. I wonder if anyone has a more rigorous benchmarking program for R. I downloaded a old R test/benchmarking program (see below), and this didn't work with the current R, and so I wondered if anyone could

Hello, I recently tried Xen 3.0.3-testing on my 2x Opteron-244 SMP machine. Everything works well exept one detail, which is not unimportant to me: the powernow-k8 driver fails to switch frequencies. A closer look on the powernow-k8 source reveals, that some MSR write operations seem to be without any effect. A rdmsr() before and after the wrmsr() call reveals that no real change on the register

Hi, Someone knows an Open Source solution that can handle Outbound IVR for asterisk?. What I'm looking it a pre-preprogrammed a telephone call that reach a Person and start making an Interview over the telephone. Specifically I want to call all my customers exactly one hour after the service has been performed and ask some questions in an IVR, also the results of the Interview I will need

All: I have a matrix, X, with a LARGE number of rows. Consider the following three rows of that matrix: 1 1 1 1 2 2 3 3 1 1 1 1 3 3 2 2 3 3 2 2 1 1 1 1 I wish to fit many one-way ANOVAs to some response variable using each row as a set of factors. For example, for each row above I will do something like anova(lm(Y~as.factor(X[1,]))). My problem is that in the above example, I do not want

All: I am trying to use gam in a scatterplot smoothing problem. The data being smoothed have greater 1000 observation and have multiple "humps". I can smooth the data fine using a function something like: out <- ksmooth(x,y,"normal",bandwidth=0.25) plot(x,out$y,type="l") The problem is when I try to fit the same data using gam out <-

Sorry: I should have reproduced the "problem" with a simple example. I do this below. I think there is likely a switch I can change using par, but don't know what it is. The problem is the tick marks for the Y- axis are only on plots in column #1 and for the X-axis in row # 2. Tony x <- 1:10 y <- 1:10*5 par(mfrow=c(2,2)) plot(x,y) plot(x,y) plot(x,y) plot(x,y)

Been learning * now for a couple of weeks and have all basic features running including VM, MoH, FX lines, iaxtel, and FWD. However, I seem to be lacking documentation on a lot of technical things and am wondering if I overlooked something that is obvious to others. (I do have the Handbook, have been doing a fair amount of google searches, and read the README.* files.) Examples, Where should I

Hi Michal, I'm cc'ing the NUT development mailing list. you should also subscribe to it. 2010/10/28 Michal Soltys > Greetings > > I have loads of relatively old smart upses that I currently control using a > bit patched version of apcupsd - the reason for those patches is the age of > the upses - roughly 10 years old+, which are can be a bit tricky to control. > From

I am trying to use dynamic loading of an outside C routine. I am attempting 6.12.1 of Phil Spector's book. When I try to load the object file I get an error I don't understand: > dyn.load("runa.o") Error in dyn.load(x) : unable to load shared library "/usr/home/tdlong/run_avg/runa.o": /usr/home/tdlong/run_avg/runa.o: ELF file's phentsize not the expected

All: Are there any functions out there for selecting all the points in a region of a plot. I envision something like the identify() function except one could circle a cloud of points (and perhaps a vector would be returned of the same length as the points plotted indicating logical membership in the circled cloud). Perhaps someone has done something with the locator() function that would

I am running R version ??? under Redhat 5.2. It seems as though the --nsize object has no effct on the size of the allocated Ncells as determined using gc(). Yes, I have that much data.... That is if I envoke R with R --vsize 100 --nsize 5000000 then type gc() I get free total Ncells 92202 200000 Vcells 12928414 13107200 Thanks Tony Long Ecology and Evolutionary Biology Steinhaus

All: I am getting an error code from the optimization function. The code is Error in optim(p,fun.LLike, lower=low, upper = up, method = "L-BFGS-B", : non-finite finite-difference value [0] If I add a trace=6 option to my control list the last message before this error is: At X0, 0 variables are exactly at the bounds Any ideas on where I should start would be

Hi, Well, I joined a new DC on my domain and when run "ldapcmp" it reported errors. I follow Samba wiki to Joining a Samba to an Existing Active Directory. * Comparing [DNSFOREST] context... * Objects to be compared: 18 Comparing: 'CN=Infrastructure,DC=ForestDnsZones,DC=dom,DC=city10,DC=com,DC=br' [ldap://artemis]

The command below was run on FreeBSD 7.2 with rsync-3.0.3. The command works but copies all files, i.e. no hard links created. rsync -vaxH --numeric-ids --delete --exclude=/home/apache/logs \ --link-dest=/mnt/200912/newharbor/home -e ssh newharbor.safeport.com:/home \ 201001/newharbor/home/ The command was run without the continuation breaks that I added for readability. The command

Hey, This patch causes a regression on my display-less nvd7. Commenting out the aux, aux_stat and aux_mask members in nvd0_i2c_oclass fixes boot, and makes things work normally again. broken dmesg: [ 40.314470] ACPI Warning: \_SB_.PCI0.GFX0._DSM: Argument #4 type mismatch - Found [Buffer], ACPI requires [Package] (20140424/nsarguments-95) [ 40.314729] ACPI Warning: \_SB_.PCI0.GFX0._DSM:

Hello. The laboratory where I work has the possibility to get a Sun machine with 12 processors. My question is: can R be used to implememt parallel algorithms which take into account the several processors? Is there a way to run *threads* in R? If this is possible, would somebody share his/her experiences? Thanks. Isabel -- Isabel Cañette Iguá 4225 Centro de Matemática 11400

I've had a very long file written out by R with write.table, with fields of time values, converted from POSIXlt as.numeric. Among 2.5 million values, very few had 6 trailing zeroes, and those were output in scientific notation as in the subject. Is this the default behavior for long integers, and how can it be turned off (with all digits for any integer field in write.table)? This