Re: [R] iterators : checkFunc with ireadLines

2020-05-19 Thread Jeff Newmiller
Laurent... Bill is suggesting building your own indexed database... but this has been done before, so re-inventing the wheel seems inefficient and risky. It is actually impossible to create such a beast without reading the entire file into memory at least temporarily anyway, so you are better

Re: [R] Help with ggplot error: #Error in FUN(X[[i]], ...) : object 'x' not found

2020-05-19 Thread Ben Tupper
Hi, When you add a new layer to a plot and you provide new data, wouldn't it make sense to provide aes() with some indication of the names of the variables to use for coordinate plotting? You have provided a grouping variable for the paths, but what of coordinates x and y? How could geom_path()

Re: [R] iterators : checkFunc with ireadLines

2020-05-19 Thread Laurent Rhelp
Ok, thank you for the advice I will take some time to see in details these packages. Le 19/05/2020 à 05:44, Jeff Newmiller a écrit : Laurent... Bill is suggesting building your own indexed database... but this has been done before, so re-inventing the wheel seems inefficient and risky. It

Re: [R] [External] Re: [External] Get a result but an error message as well ?

2020-05-19 Thread Jeff Newmiller
Works for me. set.seed( 42 ) a <- c(2,4,3,4,6,5,3,1,2,3,4,3,4,5,65) b <- c(23,45,32,12,23,43,56,44,33,11,12,54,23,34,54) d <- c(9,4,5,3,2,1,3,4,5,6,4,9,10,11,18) my.experiment <- function() { a <- a + rnorm( length( a ), 0, 0.05 ) b <- b + rnorm( length( b ), 0, 0.05 ) d <- d + rnorm( lengt

Re: [R] [External] Re: [External] Get a result but an error message as well ?

2020-05-19 Thread Richard M. Heiberger
you need to pay attention to the intermediate structures that you generate. If all you want is one number, then that is what you should create mean(replicate( 500, my.experiment() )) Since you do seem to want to store the intermediate values, then you must name the object according to its structu

Re: [R] [External] Get a result but an error message as well ?

2020-05-19 Thread varin sacha via R-help
Hi Rui, If I don't transpose t() the output of the replicate (my R code here below) I still get an error message !! a=c(2,4,3,4,6,5,3,1,2,3,4,3,4,5,65) b=c(23,45,32,12,23,43,56,44,33,11,12,54,23,34,54) d=c(9,4,5,3,2,1,3,4,5,6,4,9,10,11,18) my.experiment

Re: [R] [External] Get a result but an error message as well ?

2020-05-19 Thread Rui Barradas
Hello, Inline. Às 21:38 de 19/05/20, varin sacha via R-help escreveu: Hi Richard, Thanks for your response. However, how can I correct my R code knowing that I want, as a result, only one value : the mean of the 500 MSE_OLS values ? Just don't transpose the output of replicate? Hope this

Re: [R] [External] Get a result but an error message as well ?

2020-05-19 Thread varin sacha via R-help
Hi Richard, Thanks for your response. However, how can I correct my R code knowing that I want, as a result, only one value : the mean of the 500 MSE_OLS values ? Le mardi 19 mai 2020 à 21:59:07 UTC+2, Richard M. Heiberger a écrit : > dim(my.data) [1]   1 500 you have a matrix

Re: [R] [External] Get a result but an error message as well ?

2020-05-19 Thread Richard M. Heiberger
> dim(my.data) [1] 1 500 you have a matrix with a single row and 500 columns. you gave a name to only the first column. Look at the result of replicate(). it is a vector. You transposed it into a one-row matrix. > tmp <- replicate( 500, my.experiment() ) > dim(tmp) NULL > length(tmp) [1] 50

[R] Get a result but an error message as well ?

2020-05-19 Thread varin sacha via R-help
Dear R-experts, Here is my R code, I get a result but I also get an error message so I doubt I can trust the result I get. What is going wrong ? Many thanks. a<-c(2,4,3,4,6,5,3,1,2,3,4,3,4,5,65) b<-c(23,45,32,12,23,43,56,44,33,11,12,54,23,34,54) d<-c(9,4

Re: [R] text annotation on Manhattn plot in qqman

2020-05-19 Thread Michael Dewey
Dear Ana Perhaps paste together SNP and GENE using paste() and then supply that as the snp parameter. Michael On 19/05/2020 17:12, Ana Marija wrote: Hello, I am making manhattan plot with: library(qqman) manhattan(a, chr="CHR", bp="BP", snp="SNP", p="P",annotatePval = 0.0001) and I would l

Re: [R] text annotation on Manhattn plot in qqman

2020-05-19 Thread Ana Marija
Hi Michael, can you please send me code how that would be done? Thanks Ana On Tue, May 19, 2020 at 11:18 AM Michael Dewey wrote: > > Dear Ana > > Perhaps paste together SNP and GENE using paste() and then supply that > as the snp parameter. > > Michael > > On 19/05/2020 17:12, Ana Marija wrote:

[R] text annotation on Manhattn plot in qqman

2020-05-19 Thread Ana Marija
Hello, I am making manhattan plot with: library(qqman) manhattan(a, chr="CHR", bp="BP", snp="SNP", p="P",annotatePval = 0.0001) and I would like to annotate these two SNPs which are above the threshold so that they have GENE name beside them: > a[a$SNP=="rs4081570",] SNPP CHR

Re: [R] Help with ggplot error: #Error in FUN(X[[i]], ...) : object 'x' not found

2020-05-19 Thread Rui Barradas
Hello, I'm getting a different error: p <- ggplot(c1members2, aes(x,y)) + geom_point() p + geom_line(data = paths, aes(group = trip_id)) + ggtitle(paste0("Optimal route with cost: ", round(objective_value(result), 2))) #Error in objective_value(result) : object 'result' not found When

Re: [R] Help with ggplot error: #Error in FUN(X[[i]], ...) : object 'x' not found

2020-05-19 Thread Poling, William via R-help
Hello Chris, and thank you for your response. I tried it both TRUE and FALSE given my confusion with the reference I found that idea in. It did not help either way. Glad you noticed though, thank you. WHP William H. Poling Ph.D., MPH  | Senior Data Scientist, Medicare Stars, CVS Health p 813

Re: [R] Help with ggplot error: #Error in FUN(X[[i]], ...) : object 'x' not found

2020-05-19 Thread Chris Evans
- Original Message - > From: "Poling, William via R-help" > To: "Ben Tupper" > Cc: r-help@r-project.org > Sent: Tuesday, 19 May, 2020 15:04:25 > Subject: Re: [R] Help with ggplot error: #Error in FUN(X[[i]], ...) : object > 'x' not found > Hello Ben and thank you for your response. >

Re: [R] Help with ggplot error: #Error in FUN(X[[i]], ...) : object 'x' not found

2020-05-19 Thread Poling, William via R-help
Hello Ben and thank you for your response. I thought that was what I was doing when I ran this version: Test1 p <- ggplot(c1members2, aes(x,y)) + geom_point() p + geom_line(data = paths, aes(group = trip_id),inherit.aes = FALSE) + #Error in FUN(X[[i]], ...) : object 'x' not found ggtitle

[R] Help with ggplot error: #Error in FUN(X[[i]], ...) : object 'x' not found

2020-05-19 Thread Poling, William via R-help
#RStudio Version Version 1.2.1335 sessionInfo() # R version 4.0.0 Patched (2020-05-03 r78349) #Platform: x86_64-w64-mingw32/x64 (64-bit) #Running under: Windows 10 x64 (build 17763) Good morning. I am testing a small sample of my data using an example found here: #https://www.r-orms.org/mixed-i

Re: [R] Classification of wind events

2020-05-19 Thread Abby Spurdle
I was impressed by Jim's effort. So, I thought I'd try to produce an exploratory plot. I've adapted some of his code. The following script produces a heatmap for a cylindrical density estimate. Bright areas are (mathematical) regions of high density. However, the interpretation is complicated by t