search for: ryszard

..." $status [1] "" $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}}

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

...) =", 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-he...

Is 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}}

...o))) as 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....

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

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") : subscri...

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

...nd this is why I get what I get. and with cat() (it apparently does not have equivalent of digits" parameter) > cat ("pi, e =", pi, exp(1), "\n") pi, e = 3.141593 2.718282 Any pointers with respect how can I print what I want to print would be greatly appreciated. Ryszard [[alternative HTML version deleted]]

...) > rnames <- 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 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 private and may ...{{dropped:11}}

Hi Ryszard! There is a dist function in R. It's in the mva package. You can set the kind of distance that you want. Thanks, Erin mailto:hodgess at gator.uhd.edu

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 ...{{dropped:11}}

...o 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 "index" function, which would return index of the column given name ? Ryszard

...e same problems in 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...

...want all of them to have different random seeds. One way of doing this is sleep 2 # (with 'sleep 1' I am often getting the same number) ... set.seed(unclass(Sys.time())) Is there a simpler way without a need to sleep between invoking different R processes ? Ryszard

...8 2 0.7071068 NaN 0.7071068 0.7071068 I am sure it is possible to remove column D (with NaN's) in some simple fashion, using is.nan function without explicitly looping through, and I am sure I was able to do it in the past, but I cannot recall how. Your help will be greatly appreciated. Ryszard [[alternative HTML version deleted]]

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

...cally straight 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}}