search for: czerminski

Any ideas what is the problem with this code? > N <- 2; c(Sys.Date(), sprintf('N = %d', N)) [1] "2012-01-10" NA Warning message: In as.POSIXlt.Date(x) : NAs introduced by coercion Best regards, Ryszard Ryszard Czerminski AstraZeneca Pharmaceuticals LP 35 Gatehouse Drive Waltham, MA 02451 USA 781-839-4304 ryszard.czerminski@astrazeneca.com -------------------------------------------------------------------------- Confidentiality Notice: This message is private and may ...{{dropped:11}}

..."" $major [1] "2" $minor [1] "14.1" $year [1] "2011" $month [1] "12" $day [1] "22" $`svn rev` [1] "57956" $language [1] "R" $version.string [1] "R version 2.14.1 (2011-12-22)" Best regards, Ryszard Ryszard Czerminski AstraZeneca Pharmaceuticals LP 35 Gatehouse Drive Waltham, MA 02451 USA 781-839-4304 ryszard.czerminski@astrazeneca.com -------------------------------------------------------------------------- Confidentiality Notice: This message is private and may ...{{dropped:11}}

...there any rational explanation for the bizarre seq() behavior below? > seq(2,8.1, lenght.out=3) [1] 2 3 4 5 6 7 8 > help(seq) > seq(2,8,length.out=3) [1] 2 5 8 > seq(2,8.1,length.out=3) [1] 2.00 5.05 8.10 Except maybe that it is early in the morning :) Best regards, Ryszard Ryszard Czerminski AstraZeneca Pharmaceuticals LP 35 Gatehouse Drive Waltham, MA 02451 USA 781-839-4304 ryszard.czerminski@astrazeneca.com -------------------------------------------------------------------------- Confidentiality Notice: This message is private and may ...{{dropped:11}}

...;, length(ys), "; length(ps) =", length(ps), "\n") cat("why length(ps) != length(ys) ???\n") # my output: # # length(yr) = 164 ; length(pr) = 164 # dim(xr) = 164 118 ; dim(xs) = 35 118 # length(ys) = 35 ; length(ps) = 164 # why length(ps) != length(ys) ??? Ryszard Czerminski phone: (781)994-0479 ArQule, Inc. email:ryszard at arqule.com 19 Presidential Way http://www.arqule.com Woburn, MA 01801 fax: (781)994-0679 -----Original Message----- From: Czerminski, Ryszard [mailto:ryszard at arqule.com] Sent: Thursday, June 20, 2002 12:15 PM To: r-help at stat....

Is there a function in R, which would return index of maximum value in a vector ? e.g. > v <- round(10*rnorm(8)) > v [1] 6 -3 -6 15 7 9 0 -19 > max(v) [1] 15 ??? index.max(v) ??? 4

In e1071 package/svm default epsilon value is set to 0.1 and not 0.5 as documentation says. R

...shown below, but this probably creates a lot of overhead and I do not know which columns have been removed. Yours, R > o [,1] [,2] [,3] [1,] 1 NA 7 [2,] 2 NA 8 [3,] 3 NA 9 > t(na.omit(t(o))) [,1] [,2] [1,] 1 7 [2,] 2 8 [3,] 3 9 Ryszard Czerminski phone: (781)994-0479 ArQule, Inc. email:ryszard at arqule.com 19 Presidential Way http://www.arqule.com Woburn, MA 01801 fax: (781)994-0679 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/...

...raight lines. In maps package I found smooth.map function, but maybe there is a more generic way of accomplishing the same thing. Ideally there would be an option to control smoothness of the contour lines in contourplot(), or levelplot(), but I cannot find a way to do it. Best regards, Ryszard Czerminski AstraZeneca Pharmaceuticals LP 35 Gatehouse Drive Waltham, MA 02451 USA 781-839-4304 ryszard.czerminski@astrazeneca.com -------------------------------------------------------------------------- Confidentiality Notice: This message is private and may ...{{dropped:11}}

...predict.randomForest(rf, X) : subscript out of bounds BTW: Just find out that apparently predict() does not like NULL name in X, because this works fine: > one <- rep(1, length(chirps)) > X <- cbind(one,chirps) > rf <- randomForest(X, temp) > yp <- predict(rf, X) Ryszard Czerminski AstraZeneca Pharmaceuticals LP 35 Gatehouse Drive Waltham, MA 02451 USA 781-839-4304 ryszard.czerminski@astrazeneca.com -------------------------------------------------------------------------- Confidentiality Notice: This message is private and may ...{{dropped:11}}

I have a data frame (df) with colums x, y and z. e.g. df <- data.frame(x = sample(4), y = sample(4), z = sample(4)) I can extract column z by: df$z or df[3] I can also extract columns x,y by: df[1:2] or by df[-3]. Is it possible to extract x,y columns in a "symbolic" fashion i.e. by equivalent of df[-z] (which is illegal) ??? Or alternativeley, is there an equivalent of

...character() > for (i in 1:nr) rnames[i] <- lst[[i]] > df <- data.frame(row.names=rnames) > df$memberships <- attr(lst, 'memberships') > a {"A" [0.1], "B" [0.8]} > df memberships A 0.1 B 0.8 > Best regards, Ryszard --- Ryszard Czerminski -------------------------------------------------------------------------- Confidentiality Notice: This message is private and may ...{{dropped:11}}

I want to fit a linear model with fixed slope e.g. y = x + b (instead of general: y = a*x + b) Is it possible to do with lm()? Regards, Ryszard -------------------------------------------------------------------------- Confidentiality Notice: This message is private and may ...{{dropped:11}}

...n 2005 (see below) is probably not directly applicable, still it looks like it is too big data set for this dataset/machine combination. How does memory usage in randomForest scale with dataset size? Is there a way to build global rf model with dataset of this size? Best regards, Ryszard Ryszard Czerminski AstraZeneca Pharmaceuticals LP 35 Gatehouse Drive Waltham, MA 02451 USA 781-839-4304 ryszard.czerminski@astrazeneca.com RE: [R] randomForest: too many element specified? Liaw, Andy Mon, 17 Jan 2005 05:56:28 -0800 > From: luk > > When I run randonForest with a 169453x5 matrix, I got the &g...

I don't know the details of pls (in the pls.pcr package, I assume), but if you use validation="CV", that says you want to use CV to select the best number of components. Then why would you specify ncomp as well? Andy > From: ryszard.czerminski at pharma.novartis.com > > When I try to use ncomp parameter in pls procedure I get > following error: > > > library(pls.pcr) > > m <- pls(x, y, validation = "CV", niter = 68, ncomp = 16) > Error in inherits(x, "data.frame") : subscript out of b...

Full_Name: Ryszard Czerminski Version: 1.8.1 OS: GNU/Linux Submission from: (NULL) (205.181.102.120) prcomp(..., scale = TRUE) does not work correctly: $ uname -a Linux 2.4.20-28.9bigmem #1 SMP Thu Dec 18 13:27:33 EST 2003 i686 i686 i386 GNU/Linux $ gcc --version gcc (GCC) 3.2.2 20030222 (Red Hat Linux 3.2.2-5) > a <-...

And on that note, here is a function that I use to get around it: -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of Jeffrey Horner Sent: Friday, November 03, 2006 10:01 AM To: ryszard.czerminski at novartis.com Cc: r-help at stat.math.ethz.ch Subject: Re: [R] R CMD BATCH: unable to start device PNG ryszard.czerminski at novartis.com wrote: > Dear r-helpers: > > Any ideas how to avoid problem described below? > I am having the same problem when I run R remotly (not from cgi sc...

How can I remove columns with NaN entries ? Here is my simple example: > data <- read.csv("test.csv") > xdata <- data[3:length(data)] > xs <- lapply(xdata, function(x){(x - mean(x))/sqrt(var(x))}) > x <- data.frame(xs) > x C D E F 1 -0.7071068 NaN -0.7071068 -0.7071068 2 0.7071068 NaN 0.7071068 0.7071068

How can I select random subsets (rows!) from a data set ? If I generate simple data set > a <- data.frame(x=1:2, y = NaN, z = 2:1) > a x y z 1 1 NaN 2 2 2 NaN 1 I can select random subsets (colums) very easily using sample function: > sample(a, 2) z y 1 2 NaN 2 1 NaN I expected that using transpose of a would do the same for rows, but I am getting rather unexpected

I would like to use greek "tau" as a symbol of variable to integrate over in plotmath expression(integral(f(tau)*dtau, 0,t)) but nothing seems to work. I tried d{\tau}, d\tau, etc., without any success Is it possible? How can I accomplish this? Best regards, Ryszard -------------------------------------------------------------------------- Confidentiality Notice: This message is

I am looking for an efficient way to find near neighbors... More specifically... I have two sets of points: A & B and I want to find points in set B which are closer to set A than some cutoff (or n-closest) I will appreciate very much any pointers... Ryszard -------------------------------------------------------------------------- Confidentiality Notice: This message is private and may