[R] Continuous Wavelet tranform in Biwavelet Package

Dear All, I have 15 minutes of temperature data for 3 years and would like calculate and plot the CWT using wt command in Biwavelet package. Appreciate if someone could tell me the way to define it in wt command. Thank you, [[alternative HTML version deleted]] __

Re: [R] When customizing last line, the code stops working

A clarification regarding my previous post for those ineterested in the issue: Looking carefully at Duncan's solution, the key point would be introducing envir = environment() as argument of lapply function (the default envir of do.call function is parent.frame() ). ) union <- do.call(

[R] Converting a list to a data frame

There is probably a very simple elegant way to do this, but I have been unable to find it. Here is a toy example. Suppose I have a list of data frames like this. print(x <- list('1'=data.frame(id=1:4,expand.grid(x1=0:1,x2=0:1)),'2'=data.frame(id=5:8,expand.grid(x1=2:3,x2=2:3 $`1` id x1

Re: [R] Converting a list to a data frame

Hi Kevin, There may be a more elegant way but the following do.call and lapply should solve your problem. do.call(rbind, lapply(seq(length(x)), function(i) data.frame(set=i, x[[i]]))) Regards, Charles On Fri, Nov 4, 2016 at 7:37 AM, Kevin E. Thorpe wrote: > There is probably a very simple ele

Re: [R] Converting a list to a data frame

Hi, You have also "rbindlist()" function in package "data.table" that does exactly what you need. Kind Regards, Carlos Ortega www.qualityexcellence.es 2016-11-04 13:37 GMT+01:00 Kevin E. Thorpe : > There is probably a very simple elegant way to do this, but I have been > unable to find it. He

Re: [R] multiplying a matrix by a vector

My goodness! > x %*% diag(y) [,1] [,2] [1,]2 12 [2,]4 15 [3,]6 18 will do. -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On

[R] [R-pkgs] new package: rpostgisLT

We're announcing the initial CRAN release of 'rpostgisLT' (v0.4.0), which is an extension to our 'rpostgis ' package. rpostgisLT is aimed at those using R and/or the PostgreSQL/PostGIS database system to manage and analyze animal trajecto

Re: [R] Converting a list to a data frame

I believe that a slightly more efficient way of doing this without leaving base R is: cbind(do.call(rbind,x), set = rep(seq_along(x), vapply(x,nrow,1)) ) Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka

Re: [R] multiplying a matrix by a vector

Nice! Thanks a lot, everybody! Dimitri On Fri, Nov 4, 2016 at 10:35 AM, Bert Gunter wrote: > My goodness! > >> x %*% diag(y) > > [,1] [,2] > [1,]2 12 > [2,]4 15 > [3,]6 18 > > will do. > > -- Bert > > > > Bert Gunter > > "The trouble with having an open mind is that people

Re: [R] multiplying a matrix by a vector

Sara wins on memory use. Rui wins on speed. Bert wins on clarity. library(microbenchmark) N <- 1000 x <- matrix( runif( N*N ), ncol=N ) y <- seq.int( N ) microbenchmark( { t( y * t(x) ) } , { x %*% diag( y ) } , { sweep( x, 2, y, `*` ) } ) Unit: milli

Re: [R] multiplying a matrix by a vector

Notice though, that Bert loses (or _should_ lose) for larger values of N, since that method involves O(N^3) operations whereas the other two are O(N^2). I am a bit surprised that sweep() is so inefficient even at N=1000. -pd On 04 Nov 2016, at 16:41 , Jeff Newmiller wrote: > Sara wins on mem

[R] [R-pkgs] Major Update to rms package: 5.0-0

A major new version of the rms package is now on CRAN. The most user-visible changes are: - interactive plotly graphic methods for model fits. The best example of this is survplot for npsurv (Kaplan-Meier) estimates where the number of risk pop up as you hover over the curves, and you can click

Re: [R] multiplying a matrix by a vector

> On 4 Nov 2016, at 16:41, Jeff Newmiller wrote: > > Sara wins on memory use. > > Rui wins on speed. > > Bert wins on clarity. > > library(microbenchmark) > > N <- 1000 > x <- matrix( runif( N*N ), ncol=N ) > y <- seq.int( N ) > > microbenchmark( { t( y * t(x) ) } > , { x %*% d

Re: [R] multiplying a matrix by a vector

> On 4 Nov 2016, at 19:27, Berend Hasselman wrote: > >> >> On 4 Nov 2016, at 16:41, Jeff Newmiller wrote: >> >> Sara wins on memory use. >> >> Rui wins on speed. >> >> Bert wins on clarity. >> >> library(microbenchmark) >> >> N <- 1000 >> x <- matrix( runif( N*N ), ncol=N ) >> y <- seq.in

Re: [R] multiplying a matrix by a vector

Berend, et. al.: "I don't quite understand why I get a very different ranking." Nor do I. But, as Peter noted, my matrix multiplication solution "should" be worst in terms of efficiency, and it clearly is. The other two *should* be "similar", and they are. So your results are consistent with what

Re: [R] multiplying a matrix by a vector

I love R for it: if you experiment enough with t() you can find a way to multiply almost anything by almost anything! :) On Thu, Nov 3, 2016 at 5:32 PM, Rui Barradas wrote: > Hello, > > Take advantage that in R '*' recycles its arguments (the shorter one) and > that the operation is performed col

Re: [R] multiplying a matrix by a vector

If you want to get a comparable result, then run the same code. You may not want a comparable result though... your version straightens out the issue Peter was puzzled by. -- Sent from my phone. Please excuse my brevity. On November 4, 2016 11:35:02 AM PDT, Berend Hasselman wrote: > >> On 4 N

[R] R segfault at startup, causing R crash

Dear R list, Every time I start R, the following error message come up: *** caught segfault *** address 0x0, cause 'unknown' Possible actions: 1: abort (with core dump, if enabled) 2: normal R exit 3: exit R without saving workspace 4: exit R saving workspace > Selection: I am unable to do a s

Re: [R] R segfault at startup, causing R crash

I just got the following error instead of the previous one, 2016-11-04 16:50:08.446 R[23372:326103] -[NSDrawerWindow setRuleThickness:]: unrecognized selector sent to instance 0x7fabc94a74b0 2016-11-04 16:50:08.454 R[23372:326103] *** RController: caught ObjC exception while processing system even

Re: [R] multiplying a matrix by a vector

> On 4 Nov 2016, at 17:35, peter dalgaard wrote: > > Notice though, that Bert loses (or _should_ lose) for larger values of N, > since that method involves O(N^3) operations whereas the other two are > O(N^2). I am a bit surprised that sweep() is so inefficient even at N=1000. I also got diff

Re: [R] R segfault at startup, causing R crash

Here's the sessionInfo() > sessionInfo() R version 3.3.2 (2016-10-31) Platform: x86_64-apple-darwin13.4.0 (64-bit) Running under: macOS Sierra 10.12.1 locale: [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8 attached base packages: [1] stats graphics grDevices utils dat

Re: [R] multiplying a matrix by a vector

OK, just for fun, I added a couple of *obvious* (or should have been) approaches and got results that exactly agreed with Peter D's comments: all O(n^2) results were essentially the same, and my matrix multiplication O(n^3) solution was **far** worse. So I'm not sure quite where Jeff's result came

Re: [R] R segfault at startup, causing R crash

> On Nov 4, 2016, at 1:39 PM, C W wrote: > > Dear R list, > > Every time I start R, the following error message come up: > > *** caught segfault *** > address 0x0, cause 'unknown' > > Possible actions: > 1: abort (with core dump, if enabled) > 2: normal R exit > 3: exit R without saving works

[R] loop with variable names

Suppose I have the following data: y<-rnorm(10) age<-rnorm(10) sex<-rbinom(10,1, 0.5) edu<-round(runif(10, 1, 20)) edu2<-edu^2 df<-data.frame(y,age,sex,edu,edu2) I want to run a large number of models, for example: lm(y~age) lm(y~age+sex) lm(y~age+sex+edu) lm(y~age+sex+edu+edu2) lm(y~sex+edu2)

[R] Using map with filled.contour doesn't display map

Hello, I'm trying to overlay a map on top of data showing temperatures across the world. The code I'm using is: filledContour(tempdata, plot.axes={axis(1); axis(2); map("world2", add=TRUE)}) where tempdata is a raster file made from a netcdf file downloaded from the NOAA website. The fille

Re: [R] loop with variable names

It's hard to imagine a situation where this makes sense, but of course you can do it if you want. Perhaps rhs <- unlist(sapply(1:(ncol(df)-1), function(x) apply(combn(names(df)[-1], x), 2, paste, collapse = " + "))) lapply(rhs, function(x) lm(as.formula(paste("y ~", x)), data = df)) --Ista On Fr

Re: [R] Using map with filled.contour doesn't display map

> On Nov 4, 2016, at 11:26 AM, Jeff Charlton wrote: > > Hello, > > I'm trying to overlay a map on top of data showing temperatures across the > world. The code I'm using is: > > filledContour(tempdata, plot.axes={axis(1); axis(2); map("world2", > add=TRUE)}) > > where tempdata is a raste

Re: [R] Using map with filled.contour doesn't display map

Hi David, Thank you for the response. I apologize that I didn't give a more complete example. Here's code that recreates the issue I'm having on my system. You'll have to enter the path to your packages and where you want to store the downloaded file. library("maps", lib.loc="") library("raster