Re: [R] Zipping Rdata Files

2005-04-08 Thread Prof Brian Ripley
On Thu, 7 Apr 2005 [EMAIL PROTECTED] wrote: Saving Rdata files in a zip archive form can in some cases save a considerable amount of disk space. Not if they were saved with compress=TRUE: it is likely to increase the size of compressed saved images. R has the zip.file.extract function to extract

Re: [R] hex format

2005-04-08 Thread Martin Maechler
David == David Forrest [EMAIL PROTECTED] on Thu, 7 Apr 2005 16:19:33 -0500 (CDT) writes: David I think R has the hex to decimal OK, but might be David lacking in the decimal to hex case David zz-function(x){ x-as.numeric(sub(#,'0x',x)); David c(x%/%256^2,

Re: [R] /bin/exec/R: No such file or directory

2005-04-08 Thread Uwe Ligges
Jarmila Bohmanova wrote: I have just installed R-2.0.1 from R-2.0.1.tar.gz on SUSe 9.1 64bit. When I am trying to launch R: R_HOME_DIR/bin/R; I am getting following message: ./R: line 151: /R_HOME_DIR/bin/exec/R: No such file or directory ./R: line 151: exec: /R_HOME_DIR/bin/exec/R: cannot

Re: [R] Adapt Function Examples

2005-04-08 Thread Uwe Ligges
[EMAIL PROTECTED] wrote: I have read the help file for Adapt, but I cannot create a functn that works. I believe this is because I do not understand how to do this, and I have not found any working examples posted in the help. I have recieved many different errors in my attempts. Please post a

Re: [R] NA in table with integer types

2005-04-08 Thread Prof Brian Ripley
NaN only applies to double values: there is no integer NaN (nor Inf nor -Inf). The difference is clear from factor(x, exclude=NaN) [1] 1233NA Levels: 1 2 3 NA factor(as.integer(x), exclude=NaN) [1] 1233NA Levels: 1 2 3 If you read ?factor it says exclude: a

[R] weird results w/ prcomp-princomp

2005-04-08 Thread Alessandro Bigi
I am doing a Principal Component Analaysis (PCA) on a 44x19 matrix. with princomp(x,cor=TRUE,scores=TRUE) and prcomp(x,scale=TRUE,center=TRUE) The resulted eigenv. and rotated matrix are the same (as expected), however the sum of eigenvalues is lower than 19 (number of variables). With a

Re: [R] weird results w/ prcomp-princomp

2005-04-08 Thread Jari Oksanen
On Fri, 2005-04-08 at 11:12 +0200, Alessandro Bigi wrote: I am doing a Principal Component Analaysis (PCA) on a 44x19 matrix. with princomp(x,cor=TRUE,scores=TRUE) and prcomp(x,scale=TRUE,center=TRUE) The resulted eigenv. and rotated matrix are the same (as expected), however the sum

[R] Can't get function to run iteratively

2005-04-08 Thread Uzuner, Tolga
Trying to do something very simple... numerical.grad function(func,x, eps=1e-12) { # very simple (crude) numerical approximation f -func(x) df -1:length(x) for (i in 1:length(x)) { dx - x dx[i] - dx[i] +eps df[i] - (func(dx)-f)/eps } df } test-function(x){x^2}

Re: [R] Can't get function to run iteratively

2005-04-08 Thread Dimitris Rizopoulos
R is a vectorized language, look at this: fd - function(b, gn, ..., eps=sqrt(.Machine$double.neg.eps)){ n - length(b) g0 - gn(b, ...) out - numeric(n) b. - b + eps*max(abs(b), 1) c(gn(b., ...) - g0) / (b. - b) } ## test - function(x) x*x fd(seq(0.1, 0.5, 0.1), test)

Re: [R] off-topic question: Latex and R in industries

2005-04-08 Thread Rolf Turner
Donald Ingram wrote: ( Is there a universal vector format I could use ? ) If there were one, you can be sure that Microsoft would take steps to corrupt their implementation of the ``universal'' format so that material produced by non-Microsoft software would be

RE: [R] /bin/exec/R: No such file or directory

2005-04-08 Thread Jarmila Bohmanova
Make gives me following error message: R_HOME_DIR/src/main/array.c:504: undefined reference to `zgemm' array.o(.text+0x2526): In function `do_matprod': R_HOME_DIR/src/main/array.c:472: undefined reference to `dsyrk' array.o(.text+0x2601):R_HOME_DIR/src/main/array.c:487: undefined reference to

[R] is there any function to do oblique rotation in factor analysis?

2005-04-08 Thread ronggui
splus has many different rotation methods,but i can only find 2 in R.i want to do oblique rotation .has any function for this? i have search the web,but can not find. __ R-help@stat.math.ethz.ch mailing list

[R] TR: The results of your email commands

2005-04-08 Thread Kervahu Anne
Hi, I try to minimize the sum of the sum of sce. The following program has been created but it only takes in consideration the last kinetic and not the first ones. I think that I have forget a subscrib but I don't know where. so if you can help me, it will be great I have

RE: [R] 'skewing' a normal random variable

2005-04-08 Thread Mohammad A. Chaudhary
Hi: Thank you and sorry for getting back late. Your input helped me a lot. I particularly liked using box.cox. Regards, Ashraf __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] NA in table with integer types

2005-04-08 Thread Paul Rathouz
OK. Thanks. So, if you use table() on a factor that contains NA's, but for which NA is not a level, is there any way to get table to generate an entry for the NAs? For example, in below, even exclude=NULL will not give me an entry for NA on the factor y: x - c(1,2,3,3,NA) y - factor(x) y

Re: [R] NA in table with integer types

2005-04-08 Thread Gabor Grothendieck
On Apr 8, 2005 9:05 AM, Paul Rathouz [EMAIL PROTECTED] wrote: OK. Thanks. So, if you use table() on a factor that contains NA's, but for which NA is not a level, is there any way to get table to generate an entry for the NAs? For example, in below, even exclude=NULL will not give me an

Re: [R] NA in table with integer types

2005-04-08 Thread Prof Brian Ripley
On Fri, 8 Apr 2005, Paul Rathouz wrote: OK. Thanks. So, if you use table() on a factor that contains NA's, but for which NA is not a level, is there any way to get table to generate an entry for the NAs? For example, in below, even exclude=NULL will not give me an entry for NA on the factor y:

[R] anova with gam?

2005-04-08 Thread Bill Shipley
Hello. In SPLUS I am used to comparing nested models in gam using the anova function. When I tried this in R this doesn't work (the error message says that anova() doesn't recognise the gam fit). What must I do to use anova with gam? To be clear, I want to do the following: fFit1-gam(y~x,.)

Re: [R] NA in table with integer types

2005-04-08 Thread Petr Pikal
On 8 Apr 2005 at 14:20, Prof Brian Ripley wrote: On Fri, 8 Apr 2005, Paul Rathouz wrote: OK. Thanks. So, if you use table() on a factor that contains NA's, but for which NA is not a level, is there any way to get table to generate an entry for the NAs? For example, in below, even

RE: [R] anova with gam?

2005-04-08 Thread John Fox
Dear Bill, There are anova.gam methods in both the mgcv and gam packages -- perhaps you have an older version of one of these packages. I hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4

[R] error on install Rmpi packages

2005-04-08 Thread Marcelo Luiz de Laia
Dear Sir/Madam, I need to install Rmpi Package on my linux debian like, but I get some errors. I try to instal it using install.packages or after download and save it on my local dir and run R CMD INSTALL packagename and, both, I get these errors: checking mpi.h usability... no checking mpi.h

Re: [R] /bin/exec/R: No such file or directory

2005-04-08 Thread Uwe Ligges
Jarmila Bohmanova wrote: Make gives me following error message: R_HOME_DIR/src/main/array.c:504: undefined reference to `zgemm' array.o(.text+0x2526): In function `do_matprod': R_HOME_DIR/src/main/array.c:472: undefined reference to `dsyrk' array.o(.text+0x2601):R_HOME_DIR/src/main/array.c:487:

[R] restrict namespace inside functions?

2005-04-08 Thread Vivek Rao
Is there a way to exclude from the namespace all variables other than function arguments and local variables? For example, I would like the following code a = 2.0 mult - function(x) { pi = 3.14 return(pi*a*x) } print(mult(10.0)) to say a not found rather than using a = 2.0 inside function

Re: [R] restrict namespace inside functions?

2005-04-08 Thread Deepayan Sarkar
On Friday 08 April 2005 09:50, Vivek Rao wrote: Is there a way to exclude from the namespace all variables other than function arguments and local variables? For example, I would like the following code a = 2.0 mult - function(x) { pi = 3.14 return(pi*a*x) } print(mult(10.0))

[R] DLL Memory Problem

2005-04-08 Thread Brian Habing
Hello, I have created a .dll file using G77 and MinGW on my PC (Windows 2000). After using dyn.load to bring it into R2.0.1, I then call the .dll through the function ccprox shown below. It returns the correct values. If I run it a second time though it returns different values, so it seems

[R] kunamorph@web.de

2005-04-08 Thread Christfried Kunath
Hello, how can I use the function cor() with x and y in function aggregate() or by()? The data are like this: x y group 1 4 B 2 4 B 3 5 C I would like obtain the correlation between x and y for each subset. I don't want to use the workaround with the function subset(), because I

Re: [R] hex format

2005-04-08 Thread Earl F. Glynn
Thomas Lumley [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] Yes, and convertColor in R-devel does quite a few of these (XYZ tristimulus space; CIE Lab and Luv; sRGB, Apple RGB and roll-your-own RGB based on chromaticities of the primaries; and chromatic adaptation for changing the

Re: [R] DLL Memory Problem - Solved

2005-04-08 Thread Brian Habing
At 11:04 AM 4/8/2005, you wrote: Hello, I have created a .dll file using G77 and MinGW on my PC (Windows 2000). After using dyn.load to bring it into R2.0.1, I then call the .dll through the function ccprox shown below. It returns the correct values. If I run it a second time though it

Re: [R] axis colors in pairs plot

2005-04-08 Thread Anne York
Thanks Bill and Deepayan for the panel function idea. Somehow I've always associated the panel functions with lattice -- not sure why, but thanks for pointing out their more universal applicability. Anne On Thu, 7 Apr 2005, Deepayan Sarkar wrote: DS On Thursday 07 April 2005 17:51, Anne York

[R] BUG in RODBC with OS X?

2005-04-08 Thread Drew Balazs
This is my second posting on this topic, the first recieved no replies. Has anyone successfully used RODBC on the OS X platform? I've tested the ODBC drivers I'm using with two other applications and I've had no problems. I begining to think the problem is with R/RODBC and not the drivers.

Re: [R] weird results w/ prcomp-princomp

2005-04-08 Thread ronggui
R2.0.1 x-matrix(rnorm(44*19),nrow=44) princomp(x,cor=TRUE,scores=TRUE) Call: princomp(x = x, cor = TRUE, scores = TRUE) Standard deviations: Comp.1Comp.2Comp.3Comp.4Comp.5Comp.6Comp.7Comp.8 1.5874672 1.4652217 1.3088833 1.2339949 1.1697727 1.1402570 1.0774402

correlation by group (was Re: [R] kunamorph@web.de)

2005-04-08 Thread Sundar Dorai-Raj
Christfried Kunath wrote on 4/8/2005 10:39 AM: Hello, how can I use the function cor() with x and y in function aggregate() or by()? The data are like this: x y group 1 4 B 2 4 B 3 5 C I would like obtain the correlation between x and y for each subset. I don't want to use the

Re: [R] hex format

2005-04-08 Thread Duncan Murdoch
Earl F. Glynn wrote: I recently discovered R's writeBin and readBin functions. readBin, in particular, looks like an extremely powerful tool for loading external files. This might be an easier way to perform some one-time data manipulations instead of writing a separate C, Perl or Delphi program

[R] Princomp$Scores

2005-04-08 Thread Ken Termiso
Hi all, I was hoping that someone could verify this for me- when I run princomp() on a matrix, it is my understanding that the scores slot of the output is a measure of how well each row correlates (for lack of a better word) with each principal component. i.e. say I have a 300x6 log2 scaled

[R] EEG, ERP, ERF analysis with boot or bootstrap package

2005-04-08 Thread Darren Weber
Does anyone have experience in using boot or bootstrap for ERP data analysis? Event-related potentials (ERPs) are electrical recordings from the human scalp during the presentation of stimulus events. They are a timeseries with msec resolution. They are essentially a timeseries that are

[R] Using Flexible Discriminant Analysis (fda) in mda

2005-04-08 Thread Joseph Retzer
Greetings everyone, I was wondering if anyone had experience in using fda (flexible discriminant analysis) in the mda package. Specifically, I can run the model but I cant locate/produce much of anything beyond the predicted values and confusion matrix. I've been searching for more examples or

[R] subset arg lmList

2005-04-08 Thread Sebastian Luque
I'm having trouble understanding how functions in the subset argument for lmList search for the objects they need. This trivial example produces Error in rownames(fakedf) : Object fakedf not found: library(nlme) fitbyID - function() { fakedf - data.frame(ID = gl(5, 10, 50),

[R] R-generated animation of a polynomiograph

2005-04-08 Thread François Pinard
Hi, people. Two days ago, I sent to this list a little toy for exploring polynomiographs (yet, the mathematical formulas were not polynomials anymore, so the name is not really appropriate). After studying R calls, expressions and functions a bit more, I gave myself the homework of producing an

[R] Plotting principle components against individual variables

2005-04-08 Thread Brett Stansfield
Dear R I am trying to plot some Principle component scores against the individual variables that made the first biplot. First I need to identify some points in the plot as follows running2 - running[c(USA, New Zealand, Dominican Republic, Western Samoa, Cook Islands),] this works fine, I then

[R] Journal of Statistical Software, Volume 12

2005-04-08 Thread Jan de Leeuw
This volume now has extensive documentation on the following R packages. Issue 8: EbayesTresh (Johnstone/Silverman) Issue 6: spatstat (Baddeley/Turner) Issue 5: drc (Ritz/Streibig) Issue 4: normalp (Mineo/Ruggieri) Issue 3: R2WinBugs (Sturtz, Ligges, Gelman) Issue 1: BradleyTerry (Firth) Package

[R] advice on crafting examples for packages

2005-04-08 Thread Peter E. Rossi
Folks- I have developed a package which I am planning on posting to CRAN. I include examples for each of the 40 or so functions in the package. However, since the examples are non-trivial and the package is using MCMC simulation methods, it takes about 20 minutes to run the complete set of