[R] what is best for scripting?

2005-12-02 Thread Molins, Jordi
I am using R in Windows. I see that I will have to use batch processes with R. I will have to read and write text files, and run some R code; probably some external code too. I have never done scripting. Is there any document that explains simple steps with examples? I also have heard that Python

[R] using predict and recalibrating repeatedly?

2005-12-01 Thread Molins, Jordi
Let us say I have a time series I want to forecast. I have decided I am going to use ppr, nnet and svm. I calibrate my model with these 3 algorithms (let us assume I use all the data for the calibration; I do not distinguish here between training and testing: all data are used for training, just

[R] How to open a text file in my screen (not in the R session)?

2005-11-23 Thread Molins, Jordi
In the middle of my code, I have the following: write.table(PredictionTest, file = predictions.txt) what I would like is that inmediately after this line is executed, the file predictions.txt is opened in the middle of my screen (not in my R session, where I continue running my code). Is there

Re: [R] How to open a text file in my screen (not in the R sessio n)?

2005-11-23 Thread Molins, Jordi
PROTECTED] Sent: 23 November 2005 19:39 To: Molins, Jordi; 'r-help@stat.math.ethz.ch' Cc: Jordi Molins Subject: Re: [R] How to open a text file in my screen (not in the R session)? You could use the following istruction: system(notepad predictions.txt,wait=FALSE) I hope

[R] books about MCMC to use MCMC R packages?

2005-09-23 Thread Molins, Jordi
Dear list users, I need to learn about MCMC methods, and since there are several packages in R that deal with this subject, I want to use them. I want to buy a book (or more than one, if necessary) that satisfies the following requirements: - it teaches well MCMC methods; - it is easy to

[R] mars of degree 3?

2005-08-09 Thread Molins, Jordi
Dear list, when I execute library(mda) x1 - rnorm(10,mean=1,sd=.12) x2 - rnorm(10,mean=.6,sd=.12) y - pp(x1,1)+pp(x1,1)*pp(x2,.6)+.12*rnorm(length(x1)) mymars - mars(cbind(x1,x2),y,degree=2) mymars$cuts mymars$factor I get what I expected. Instead, when I execute x1 -

[R] mars of degree 3?

2005-08-09 Thread Molins, Jordi
I have forgotten to list the definition of pp. Here it is: pp - function(x,a) {ifelse(xa, x-a, 0)} __ Dear list, when I execute library(mda) x1 - rnorm(10,mean=1,sd=.12) x2 - rnorm(10,mean=.6,sd=.12) y -

[R] again, a question between R and C++

2005-07-21 Thread Molins, Jordi
Dear R Users, I want to make a call from R into C++. My inputs are List1, List2, List3, IntegerID. The amount of elements of the lists and their type depend on IntegerID. Typical elements of a given list can be vectors, doubles, and even other lists. I want to return also a list (whose nature

[R] memory error message using MASS and GLMMGibbs

2005-05-03 Thread Molins, Jordi
Hello, I was just testing the MASS code examples for chapter 10 (Random and Mixed Effects) and I have pasted the following code in an R session (2.1.0 in windows 2000 professional; I have also Xemacs + ESS installed, but I was not using them at that time; my machine has quite a lot of RAM):

[R] cross-correlations

2004-09-02 Thread Molins, Jordi
Hello, I have been looking around in past helps about cross-correlations for a set of n time series. Although you can do it by yourself calculating the out-of-the-diagonal terms in the cross correlation matrix by using pairwise combinations of ccf and the diagonal terms by using acf, this does

[R] not positive definite D matrix in quadprog

2004-09-01 Thread Molins, Jordi
Hello to everybody, I have a quadratic programming problem that I am trying to solve by various methods. One of them is to use the quadprog package in R. When I check positive definiteness of the D matrix, I get that one of the eigenvalues is negative of order 10^(-8). All the others are

[R] a question about identify and locator

2004-07-30 Thread Molins, Jordi
Hello, I use R 1.9.1 on Windows 2000. I have a chart that I am looking. I use identify() and locator() to find out the (x,y) values of interesting points in the chart, but these functions are not always helpful (e.g., sometimes identify() prints on top of the chart, making the number illegible).

[R] SVD with positivity constraints

2004-07-27 Thread Molins, Jordi
Hello, I have a matrix equation, Ax=b, that I need to solve for x. x should be a vector of positive numbers (between 0 and 1). A is not a square matrix in general. This lead me to using the SVD. However, using the SVD gives me positive and negative numbers, as well. I have some constraints

RE: [R] SVD with positivity constraints

2004-07-27 Thread Molins, Jordi
Ripley [mailto:[EMAIL PROTECTED] Sent: 27 July 2004 10:13 To: Molins, Jordi Cc: '[EMAIL PROTECTED]' Subject: Re: [R] SVD with positivity constraints If A is not square, which dimension is larger? There will most likely be either no solution or an infinity of solutions. If the latter, I think you

Recall: [R] SVD with positivity constraints

2004-07-27 Thread Molins, Jordi
The sender would like to recall the message, [R] SVD with positivity constraints. The information contained herein is confidential and is inte...{{dropped}} __ [EMAIL

RE: [R] SVD with positivity constraints

2004-07-27 Thread Molins, Jordi
of d should be set manually to 0. Is this done (maybe as an option) in svd or in La.svd? could this resetting have a material effect in the solutions found? thank you J -Original Message- From: Prof Brian Ripley [mailto:[EMAIL PROTECTED] Sent: 27 July 2004 10:13 To: Molins, Jordi Cc

RE: [R] SVD with positivity constraints

2004-07-27 Thread Molins, Jordi
=m2=0, and then m3 being the number of equations I have. Using the simplex method for linear programing would give me my solution, or would this method give me a degenerate solution? thank you, and please bear with me my desesperation ;- Jordi -Original Message- From: Molins, Jordi Sent