similar to: Adjusting two continuous variables by one continuous vari able

Hi, I want to look at the relationship between 2 continuous/quantitative variables while adjusting for a third continuous variable. For example: relationship between height and weight adjusted for age. As I understand it using a glm would work if I had a categorical variable (I could make age into a categories) but I am interested to know if anyone knows how I can do this analysis with out making

Hi all, I have a very high dimensional data and apparently there are several columns that contain similar information (some columns are equal). I want to form a matrix/data frame consisting of unique columns. Does anyone have an efficient way of getting out these columns. A small section of the data frame is given below. Thanks for helping. Stephen. > newdata [,1] [,2] [,3] [,4] [,5]

Here's one not-so-straightforward way: > f <- function(a, b) a + b > flist <- as.list(f) > names(flist)[1:2] <- c("x", "y") > flist[[3]] <- do.call("substitute", list(body(f), list(a=as.name("x"), b=as.name("y")))) > g <- as.function(flist) > g function (x, y) x + y HTH, Andy > From: john.gavin at ubs.com

Hello,I am a little bit struggling with the following: I would like to draw a coplot and in each panel draw the curve and display the coefficients of the regression model y~a*x^b. I thought of running nls(y~a*x^b,...) in panel=function(x,y...){} when calling coplot. How can I do this? Alternatively, I could do the regression before graphing with coplot and just use the panel function to draw the

Hi all, I am trying to play around with the randomForest function for classification. I know its performance is great. I am currently using the default options. It has many options. How do I further tweak the options so that I can make its performance even better? What are the options that are mostly used? Thanks a lot! M [[alternative HTML version deleted]]

Dear R-help, This one is probably a piece of cake for regexp masters. I'd like to split a character vector (for simplicity, say of length one for now) that contains fields that are delimited by arbitrary number of white spaces (e.g., " a b c "). How do I get the character vector that contain the fields? In the example I gave, I've tried: > strsplit(" a b c

> From: Weiwei Shi > > it works. > thanks, > > but: (just curious) > why i tried previously and i got > > > is.vector(sample.size) > [1] TRUE Because a list is also a vector: > a <- c(list(1), list(2)) > a [[1]] [1] 1 [[2]] [1] 2 > is.vector(a) [1] TRUE > is.numeric(a) [1] FALSE Actually, the way I initialize a list of known length is by

Dear R-help, Has anyone notice the problem that, on Windows (NT and XP), when printing a graph using the "File -> Print..." menu in the graphics window to print the graph, that line width seemed to be ignored in the printed output? For example, if I make a plot with plot(1:10, type="l", lwd=5), it looks right on screen, but when printed out using the menu, it looks like

Dear R-help: Has anyone tried using R on the the AMD Opteron in either 64- or 32-bit mode? If so, any good/bad experiences, comments, etc? We are considering getting this hardware, and would like to know if R can run smoothly on such a beast. Any comment much appreciated. Best, Andy Andy Liaw, PhD Biometrics Research PO Box 2000, RY33-300 Merck Research Labs Rahway, NJ

Googling for `rw1061.exe' turned up: http://www.cipic.ucdavis.edu/~dmrocke/Class/EAD289D/R/rw1061.exe Andy > From: Qun Shi > > Hi Andy, > > Thanks a lot for your promptly response. I searched the whole > web site, I > found the source code for version 1.6.X. Since I'm not a > computer person, > I don't how to compile it, but what I want is binary file

Dear R-help, I was bitten by the behavior of all() when given logical(0): It is TRUE! (And any(logical(0)) is FALSE.) Wouldn't it be better to return logical(0) in both cases? The problem surfaced because some un-named individual called randomForest(x, y, xtest, ytest,...), and gave y as a two-level factor, but ytest as just numeric vector. I thought I check for that in my code by testing

Dear R-devel, Has anyone played with this? http://www.amd.com/us-en/Processors/DevelopWithAMD/0,,30_2252_2282,00.html <http://www.amd.com/us-en/Processors/DevelopWithAMD/0,,30_2252_2282,00.html> . I'll probably give it a shot... Best, Andy Andy Liaw, PhD Biometrics Research PO Box 2000, RY33-300 Merck Research Labs Rahway, NJ 07065 mailto:andy_liaw@merck.com

Hi all, I am using the "svm" command in the e1071 package. Does it have an automatic way of setting the "cost" parameter? I changed a few values for the "cost" parameter but I hope there is a systematic way of obtaining the best "cost" value. I noticed that there is a "cross" (Cross validation) parameter in the "svm" function. But I

Dear R-help, Has anyone tried running R on a Beowulf-type cluster? I can get R to run in batch (using R CMD BATCH) on a cluster, but am wondering if it is possible to get an interactive R session on a compute node. Right now, if I run: beorun --nolocal R I just get the R start-up message and back to the shell prompt. If I try bpsh 0 R I can get R started (but the R prompt does not

Before I do that, I would need to remove the gcj stuff that are in /usr/bin. If I know how to remove gcj, I'd gladly do that. However, for the particular version of the OS, the entire GCC seems to be bundled into one rpm, and I could not remove just the gcj component. Neither do I wish to mess with files that are part of some RPMs--- in my experience that's invitation for trouble later.

Dear R-devel, I grabbed 1.7.0 beta from today. "make check" still fails at base-Ex.R. The last few lines of the output are: > lsf.str()#- how do the functions look like which I am using? Error in exists(nam <- nms[i], envir = envir, mode = mode) : F used instead of FALSE Execution halted Looks like the same problem as before, but at a different place. Mandrake Linux

[Not sure if such wishlist item should go to R-help or R-devel...] It would be nice if an option can be added to the R for Windows installer that will add the --internet2 option to the command line in the shortcut that gets created. Also, is it possible to make that an option that can be set at the R prompt, or perhaps even better, as the default no matter how/where R is started? For those of

[Not sure if such wishlist item should go to R-help or R-devel...] It would be nice if an option can be added to the R for Windows installer that will add the --internet2 option to the command line in the shortcut that gets created. Also, is it possible to make that an option that can be set at the R prompt, or perhaps even better, as the default no matter how/where R is started? For those of

Hi all, I'm using R0.99a on WinNT4. I downloaded a bunch of packages and tried using Rinst.exe to install them. For whatever reason, Rinst.exe does not find all the packages I downloaded. I put all files in the same temp directory, but Rinst only finds some of them. I have no problem installing the rest with WinZip, but if Rinst works it would be much easier (don't have to run WinZip

With the same system configuration (WinNT4 SP6 and 1.7.0), I get such a Dr Watson crash each time I try to use the Change dir... command in the File menu of Rgui.exe. I doesn't seems to happen if I do this immediately after starting R but well if I already did some computation. I don't have any problem by using directly setwd(). At 12:10 23/04/03, you wrote: >Date: Tue, 22 Apr 2003