[R] Error in predict.randomForest ... subscript out of bounds with NULL name in X

RF trains fine with X, but fails on prediction library(randomForest) chirps - c(20,16.0,19.8,18.4,17.1,15.5,14.7,17.1,15.4,16.2,15,17.2,16,17,14.1) temp - c(88.6,71.6,93.3,84.3,80.6,75.2,69.7,82,69.4,83.3,78.6,82.6,80.6,83.5,76 .3) X - cbind(1,chirps) rf - randomForest(X, temp) yp -

Re: [R] strange Sys.Date() side effect

To: Czerminski, Ryszard; r-help@r-project.org Subject: Re: [R] strange Sys.Date() side effect My best guess is that you are misunderstanding what the c() function does. I'd suggest reading the help page for c, obtained by typing ?c Note that if you supply c() with objects of different types (as you have

[R] strsplit() does not split on .?

Any ideas what is wrong? strsplit(a.b, .) # generates empty strings with split=. [[1]] [1] strsplit(a b, ) # seems to work fine with split= , and other characters... [[1]] [1] a b R.Version() $platform [1] x86_64-unknown-linux-gnu $arch [1] x86_64 $os [1] linux-gnu $system [1] x86_64,

Re: [R] strsplit() does not split on .?

-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Czerminski, Ryszard Sent: Thursday, January 12, 2012 10:16 AM To: r-help Subject: [R] strsplit() does not split on .? Any ideas what is wrong? strsplit(a.b, .) # generates empty strings with split=. [[1]] [1

[R] strange Sys.Date() side effect

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

[R] bizarre seq() behavior?

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

[R] sparse kernel regression

I am looking for the implementation of sparse kernel regression approach e.g. as in this paper: The Generalized LASSO. Volker Roth IEEE Transactions on Neural Networks, Vol. 15, NO. 1, January 2004. I would appreciate any pointers. Best regards, Ryszard

[R] randomForest: too many elements specified?

I getting Error in matrix(0, n, n) : too many elements specified while building randomForest model, which looks like memory allocation error. Software versions are: randomForest 4.5-25, R version 2.7.1 Dataset is big (~90K rows, ~200 columns), but this is on a big machine ( ~120G RAM) and I call

[R] how to find near neighbors?

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

[R] (no subject)

I would appreciate any advice about how to add rectangular grid to the plot generated by levelplot() similarly as grid() function adds grid to the regular plot() Best regards, Ryszard -- Confidentiality Notice: This message

[R] how to add rectangular grid to the plot generated by levelplot()?

David, Thanks a lot! panel.grid() does the trick, but the generated grid extends beyond the plot region from levelplot() Is there a way to restrict grid to the plot region from levelplot() ? Best regards, Ryszard --

[R] approximation for sum of lognormal random variables

I am looking for approximation for sum of two lognormal variables. I found some papers suggesting power lognormal or Pearson type-IV, but I wonder if there is anything ready in R? Best regards, Ryszard -- Confidentiality

[R] how fit linear model with fixed slope?

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

[R] plotmath: how to use greek symbols in expression(integral(f(tau)*dtau, 0, t))?

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

[R] smooth contour lines

Is there an easy way to control smoothness of the contour lines? In the plot I am working on due to the undersampling the contour lines I am getting are jugged, but it is clear by eye these should be basically straight lines. In maps package I found smooth.map function, but maybe there is a more

Re: [R] library sets: A EMPTY does not work; gset_intersection(A,EMPTY) works

be unlawful. -Original Message- From: David Meyer [mailto:david.me...@wu.ac.at] Sent: Sunday, April 11, 2010 5:12 PM To: Czerminski, Ryszard Cc: Peter Ehlers; R-help@R-project.org Subject: Re: [R] library sets: A EMPTY does not work; gset_intersection(A,EMPTY) works Peter, Ryszard, I

Re: [R] library sets: A EMPTY does not work; gset_intersection(A,EMPTY) works

Message- From: David Meyer [mailto:david.me...@wu.ac.at] Sent: Tuesday, April 06, 2010 6:03 PM To: Czerminski, Ryszard Cc: Peter Ehlers; R-help@R-project.org; David Meyer Subject: Re: [R] library sets: A EMPTY does not work; gset_intersection(A,EMPTY) works This is rather by chance due

Re: [R] library sets: A EMPTY does not work; gset_intersection(A,EMPTY) works

, distribute or take any action in reliance on it. Any unauthorized use or disclosure of the contents of this message is not permitted and may be unlawful. -Original Message- From: David Meyer [mailto:david.me...@wu.ac.at] Sent: Wednesday, April 07, 2010 8:25 AM To: Czerminski, Ryszard Cc

Re: [R] library sets: A EMPTY does not work; gset_intersection(A,EMPTY) works

- From: David Winsemius [mailto:dwinsem...@comcast.net] Sent: Wednesday, April 07, 2010 9:17 AM To: Czerminski, Ryszard Cc: David Meyer; R forum; Peter Ehlers Subject: Re: [R] library sets: A EMPTY does not work; gset_intersection(A,EMPTY) works On Apr 7, 2010, at 8:54 AM, Czerminski, Ryszard

Re: [R] library sets: A EMPTY does not work; gset_intersection(A,EMPTY) works

not copy, distribute or take any action in reliance on it. Any unauthorized use or disclosure of the contents of this message is not permitted and may be unlawful. -Original Message- From: Peter Ehlers [mailto:ehl...@ucalgary.ca] Sent: Wednesday, March 31, 2010 2:03 PM To: Czerminski, Ryszard

[R] library sets: A EMPTY does not work; gset_intersection(A,EMPTY) works

When using generalized sets from sets library A EMPTY does not work, but gset_intersection(A,EMPTY) works: example code below Is it a bug? Best regards, Ryszard library(sets) support - universe - c('a','b','c','d') A - gset(support=support, memberships=c(0.2, 0.2, 0.9, 0), universe=universe)

Re: [R] library sets: A EMPTY does not work; gset_intersection(A,EMPTY) works

...@ucalgary.ca] Sent: Wednesday, March 31, 2010 11:43 AM To: Czerminski, Ryszard Cc: R-help@r-project.org Subject: Re: [R] library sets: A EMPTY does not work; gset_intersection(A,EMPTY) works On 2010-03-31 9:30, Peter Ehlers wrote: Unless I'm missing something, I don't see any method in pkg:sets

[R] lattice: how to add points to the plot generated by levelplot()?

I am using levelplot() function from pkg:lattice and I want to add some points to the plot generated by levelplot() similarly as in basic R graphic sequence: plot(...); points(...), but it does not work. I was reading documentation for lattice, but so far without much success. I would very much

[R] sets package: converting a set to data frame?

I just started using nice package sets and I wonder if there are utilities to convert (some) sets to data frame (as in the example below) library(sets) a - gset(elements = list(e('A', 0.1), e('B', 0.8))) lst - as.list(a) nr - length(lst) rnames - character() for (i in 1:nr) rnames[i] -

[R] problem with apply(m, 1, min, na.rm=T)

Passing extra arguments to FUN=mean or median in apply seems fine, but when FUN=min warnings are generated? See below. Any ideas why? Best regards, Ryszard Ryszard Czerminski AstraZeneca Pharmaceuticals LP m [,1] [,2] [1,]12 [2,]3 NA [3,] NA NA apply(m, 1, median,