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 execute:

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 si

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] 1233 Levels: 1 2 3 factor(as.integer(x), exclude=NaN) [1] 1233 Levels: 1 2 3 If you read ?factor it says exclude: a vector of v

[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 comm

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 > t

[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} > y<

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 `dge

[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 https://stat.ethz.ch/mailman/listinfo/r-help

[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! http://

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 on the factor y: > x <- c(1,2,3,3,NA) > y <- factor(x) >

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

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 on the factor y: I

[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,.) f

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

2005-04-08 Thread Paul Gilbert
In the devel section of CRAN (see the bottom of the packages page) there is a package called GPArotation which does many more rotations, including oblique rotations. We intend to move it to the regular packages section shortly, so it would be nice to have comments now. Paul Gilbert > > splus ha

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 bel

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 905-525-9140x236

[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: und

[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 functi

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(mul

[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] 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

Re: [R] error on install Rmpi packages

2005-04-08 Thread Prof Brian Ripley
If you have MPI installed, it is not somewhere Rmpi can find it. The output you quote tells you how to to fix this, if it is really is installed. See also the README file in the package. On Fri, 8 Apr 2005, Marcelo Luiz de Laia wrote: Dear Sir/Madam, I need to install Rmpi Package on my linux

[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(), beca

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 > changi

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 return

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. Here

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 1.

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

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 mat

[R] quick help please

2005-04-08 Thread Janet Atim
hullo,am a student trying to use R to develop some statistical models from data as shown below,pse help me because i have a dead line to meet. yearflow 1999 23 2000 34 2001 m 2002 67 2003 m 2004 56 2005

[R] Can't read from master socket connection?

2005-04-08 Thread jhallman
On mramx1.rsma.frb.gov, I create a master socket connection: > server <- socketConnection(port = 11223, server = T, blocking = T, open = > "a+") This does not return until I connect with a client from another machine: > client <- socketConnection("mramx1.rsma.frb.gov", port = 11223, open = "a+")

Fwd: [R] Is a .R script file name available inside the script?

2005-04-08 Thread Darren Weber
-- Forwarded message -- From: Darren Weber <[EMAIL PROTECTED]> Date: Apr 8, 2005 1:16 PM Subject: Re: [R] Is a .R script file name available inside the script? To: Gabor Grothendieck <[EMAIL PROTECTED]> Right, I understand, it could be in.txt has: scriptfile <- "Rscript.R" source(

[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 defined

[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 d

[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 a

[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,

[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 wr

Re: [R] subset arg lmList

2005-04-08 Thread Prof Brian Ripley
On Fri, 8 Apr 2005, Sebastian Luque wrote: 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.

[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 exampl