RE: [R] lm with an arbitrary number of terms

I think you can do it like this lm(y~., data=data.frame) # note the dot to the right of ~ > -Original Message- > From: Richard Nixon [mailto:[EMAIL PROTECTED] > Sent: Wednesday, April 02, 2003 8:50 AM > To: [EMAIL PROTECTED] > Subject: [R] lm with an arbitrary number of terms > > > Hel

Re: [R] lm with an arbitrary number of terms

The following might work: mdl <- paste("y~", paste(names(data.frame), collapse="+")) lm(mdl, ...) If y = "y" is a column of your data.frame, you can delete it be selecting "names(data.frame)[!is.element(y, names(data.frame)]" Can you solve the problem from here? Best Wishes, Spen

[R] Multivariate Time series

Dear R People: Is there a library for Multivariate time series, please? For some reason, I'm thinking that Dr. Paul Gilbert may have one? R Version 1.6.2 (i've updated!) for Windows Thanks so much! Sincerely, Erin Hodgess University of Houston - Downtown mailto: [EMAIL PROTECTED]

[R] lme parameterization question

Hi, I am trying to parameterize the following mixed model (following Piepho and Ogutu 2002), to test for a trend over time, using multiple sites: y[ij]=mu+b[j]+a[i]+w[j]*(beta +t[i])+c[ij] where: y[ij]= a response variable at site i and year j mu = fixed intercept Beta=fixed slope w[j]=constant

RE: [R] Scripting with an external editor

A.J. Rossini writes: > From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] > Sent: Tuesday, April 01, 2003 1:28 PM > To: [EMAIL PROTECTED] > > "Prof. Brian Ripley" <[EMAIL PROTECTED]> writes: > > > I think the `trick' is how text for parsing is sent to R, > and the answer > > is via the command

[R] Trying to make a nested lme analysis

Hi, I'm trying to understand the lme output and procedure. I'm using the Crawley's book. I'm try to analyse the rats example take from Sokal and Rohlf (1995). I make a nested analysis using aov following the book. > summary(rats) Glycogen Treatment Rat Liver Min. :12

[R] Can boot return matrix?

Dear All, I have a function which takes a n x m matrix as an argument and returns an n x n matrix. I want to take bootstrap samples form the input matrix in the way as each row represent a multivariate observation, so each bootstrap sample would be an n x m matrix, and on each sample I want to ca

RE: [R] vectorize an expression

You could calculate the index vector like ind <- n * (cc - 1) + (1:n) and then use it to index bb as a vector aa <- bb[ind] Reid Huntsinger -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Sent: Wednesday, April 02, 2003 2:15 AM To: [EMAIL PROTECTED] Subject: [R] v

Re: [R] Can't run regression on G3 Mac

First, let me thank all the folks who replied, many of which didn't make it into the digest. The following summarizes their experiences w/ lm on a G3 Mac. Apparently some people have lm in R working w/ apparently the same configuration that doesn't work for others: Machine

Re: [R] replication of latin squares --- again

Dear Lamack: Have you considered "lme" or "varcomp"? For me, essential documentation for "lme" is Pinhiero and Bates (2000) Mixed-Effects Models in S and S-Plus (Springer). In my experience, it is easier to get "varcomp" to run but harder to get the answers I want from it. Best Wishes, Spenc

[R] Index of item in matrix

Hello All, Is there a fast way to find the index(row and column) of a point in a matrix? Thanks, John. -- -- Dr. John Janmaat Department of Economics, Acadia University, Wolfville, NS, B4P 2R6 E-mail: [EMAIL PROTECTED]

[R] S intrp function

Hi. I am trying some S-Plus scripts that I used a few years ago on R. Many things have worked flawlesly and I am very impressed with the work the developers have done. (As I remember, the license fee on S-Plus at that time was ~$4500 per seat (AIX) ouch!!) So, thank you. My question relates to

Re: [R] cor.test observations limit

This is a note to my second question. I just discovered my data containing an 'Inf' value and this causes the error I saw when I used rcorr() function from the Hmisc package. If I delete this observation, everything is fine. Best, Ming-Chung __ [EMAIL P

[R] Getting contours ...

Hello All, I have a matrix of data that I want to locate a specific contour from. Is there a function to return points on a contour (or all the drawn contours) from a contour plot? Thanks, John. -- -- Dr. John Janmaat Depa

[R] Efficiency: whole vector and for loop methods

Dear R-help, As a generalization, is it more efficient to use a single method to selectively transform items in a matrix? Preferably whole-vector, yes; but if one of the column vectors compels the use of a "for loop", is it then faster to do all of the transformations for the matrix via that "

Re: [R] CGIwithR for IIS

"Gianluca Emireni" <[EMAIL PROTECTED]> writes: > Hi, I have a (maybe stupid) question... Does CGIwithR package can be > installed in R for Windows to work with a Microsoft IIS web-server? Or I > need other libraries? It should be possible; IIS will run CGI scripts. However, that is the limit of

[R] Broken link on your website

There appears to be a problem on this page: http://www.r-project.org/nocvs/mail/r-announce/1997/0014.html When we click on your link to: http://www.stat.unipg.it/pub/stat/statlib/R/CRAN/contents.html we get the error: 404 File not found We last examined your page on Wed Apr 02, 2003 at 12:13:44

[R]

Can someone point out for me where to find a package to do principal component analysis for Affy data, if existing? Thanks, Yongde Yongde Bao, Ph.D Biomolecular Research Facility Department of Microbiology University of Virginia School of Medicine Charlottesville, VA 22908 E-mail: [EMAIL PROTECTE

[R] lm with an arbitrary number of terms

Hello folks, Any ideas how to do this? data.frame is a data frame with column names "x1",...,"xn" y is a response variable of length dim(data.frame)[1] I want to write a function function(y, data.frame){ lm(y~x1+...+xn) } This would be easy if n was always the same. If n is arbitrary how c

[R] SJava RInterpreter.dll

I'm getting the following error message when running an SJava example: Exception in thread "main" java.lang.UnsatisfiedLinkError: no RInterpreter in java.library.path at java.lang.ClassLoader.loadLibrary(ClassLoader.java:1344) at java.lang.Runtime.loadLibrary0(Runtime.java:744)

[R] nlme groupedData warning

I have a question about nlme (and yes, I have a copy of "Mixed-Effects Models in S and S-Plus" on order :-) I have a data frame, df, that I use to create a groupedData object: > df.K1 <- groupedData(K1 ~ tissue.ROI | scan, data=df ) If I add an outer grouping factor, I get a warning: > df.K1 <-

Re: [R] CGIwithR for IIS

CGIwithR works only on unix-like operating systems. Not Windows. David On Wednesday, Apr 2, 2003, at 14:12 Europe/London, Gianluca Emireni wrote: Hi, I have a (maybe stupid) question... Does CGIwithR package can be installed in R for Windows to work with a Microsoft IIS web-server? Or I need

RE: [R] Autogenerated png, bitmap images

Hi Tony, I looked over the files and I am puzzled why the difference in font sizes. Unfortunately my knowledge about X window configurations is pretty scant. (You mention in the shell script comments that you run the job as root, so I wonder if that makes any difference to X..I doubt it). I recall

[R] "..." and Error: unused argument(s)

Dear R people, Would it be possible to add "..." to the argument list of arrows() or should this be dealt with at the user level? The problem is that I need arrows() in a function together with points() or text(). All these are used with their default values, but I let the user to change these s

[R] Continuation string in Sweave

Dear all, In Sweave I have a code chunk aaa <- "Sex, Drug/Sex:W1,Drug:W1,\ Sex:W2,Drug:W2/Sex:W1:W2,Drug:W1:W2" which in tex code turns into \begin{Sinput} > aaa <- "Sex, Drug/Sex:W1,Drug:W1,\nSex:W2,Drug:W2/Sex:W1:W2,Drug:W1:W2" \end{Sinput} meaning that the string is not split to two line

RE: [R] randomForests predict problem

On Wed, 2 Apr 2003, Liaw, Andy wrote: > Yves, > > I will add checks for NAs in predict.randomForest(). > > In the next version of randomForest (currently called 3.9-x), there will be > facilities for handling NAs in the training set. However, there's no way to > handle NAs in the test set yet. I

RE: [R] randomForests predict problem

Yves, I will add checks for NAs in predict.randomForest(). In the next version of randomForest (currently called 3.9-x), there will be facilities for handling NAs in the training set. However, there's no way to handle NAs in the test set yet. I believe Leo is still working on that. In Leo's v.

RE: [R] randomForests predict problem

I use randomForest version 3.4-4, but yes, now I correctly omitted NA's it works. I should have made a mistake while removing them first time. I was surprised that this method doesn't have another way to deal with NA's than omitting them. As Torsten Hothorn suggested, the associated predict fun

[R] CGIwithR for IIS

Hi, I have a (maybe stupid) question... Does CGIwithR package can be installed in R for Windows to work with a Microsoft IIS web-server? Or I need other libraries? Thank you, Gianluca. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailm

RE: [R] randomForests predict problem

Yves, Which version of the package are you using? I get: > soy <- na.omit(Soybean) > ts <- sample(nrow(soy), 150, replace=FALSE) > sb.rf <- randomForest(Class ~ ., data=soy[-ts,]) > table(predict(sb.rf, soy[ts,], type="class")) 2-4-d-injury alternarialeaf-spot

Re: [R] randomForests predict problem

Well, thank you for your answer, but this is not doing the right thing, that is predicting the Class value for the test set Soybean[test,]. It gives instead prediction for data used for forest computation (ignoring all data with NA's) ; 'data' argument is simply ignored as the right name for th

Re: [R] randomForests predict problem

> Well, thank you for your answer, but this is not doing the right thing, > that is predicting the Class value for the test set Soybean[test,]. It > gives instead prediction for data used for forest computation (ignoring all > data with NA's) ; 'data' argument is simply ignored as the right name f

[R] replication of latin squares --- again

Dear all, this is a newbie's question. I have a 4x4 latin square replicated 3 times. That is: operators batches 1 2 3 4 1 A B C D 2 B C D A 3 C D A B 4 D A B C operators batches 1 2 3 4 5 A B C D 6 B C D A 7 C D

[R] randomForests predict problem

Hello everybody, I'm testing the randomForest package in order to do some simulations and I get some trouble with the prediction of new values. The random forest computation is fine but each time I try to predict values with the newly created object, I get an error message. I thought I was beca

RE: [BioC] [R] SAM: Significance of Microarray Analysis availabilityin R

Please ignore this mail. I don't remember sending this mail for the second time. Thanks to all those who replied. -Original Message- From: Adaikalavan Ramasamy Sent: Monday, March 31, 2003 3:30 PM To: [EMAIL PROTECTED]; [EMAIL PROTECTED] Subject: [BioC] [R] SAM: Significance of Microarray

Re: [R] RODBC sqlSave problem.

On Wed, 2 Apr 2003, Fredrik Karlsson wrote: > Being new to both the postgres database, ODBC and the RODBC interface, I > am somewhat confused by some of the problems I am experiencing trying to > connect R to the database. You would be: there is an extensive set of help pages, and you have not lo

Re: [R] vectorize an expression

What about: bb[cbind(1:length(cc),cc)] Eric At 17:15 2/04/2003 +1000, [EMAIL PROTECTED] wrote: Dear listers, I'm having a bad R day. I just can't think of the vectorized equivalent of: for (ii in 1:n) aa[ii] = bb[ii,cc[ii]] Any suggestion received with embarrassment and gratitude Simon Ga

Re: [R] vectorize an expression

On Wed, 2 Apr 2003 [EMAIL PROTECTED] wrote: > Dear listers, > I'm having a bad R day. I just can't think of the vectorized equivalent of: > > for (ii in 1:n) aa[ii] = bb[ii,cc[ii]] aa <- bb[cbind(1:n, cc[1:n])] -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied St

Re: [R] vectorize an expression

[EMAIL PROTECTED] wrote: Dear listers, I'm having a bad R day. I just can't think of the vectorized equivalent of: for (ii in 1:n) aa[ii] = bb[ii,cc[ii]] Any suggestion received with embarrassment and gratitude a <- b[cbind(1:n, cc[1:n])] Uwe Ligges Simon Gatehouse