Hi all,
I wrote a function that actually does what I want it to do, but it tends to be
very slow for large amount of data. On my computer it takes 5.37 seconds for
16000 data points and 21.95 seconds for 32000 data points. As my real data
consists of 18000000 data points it would take ages to use the function as it
is now.
Could someone help me to speed up the calculation?
Thank you, Tonja
system.time({
x <- runif(32000)
y <- runif(32000)
xy <- cbind(x,y)
outer <- function(z){
!any(x > z[1] & y > z[2])}
j <- apply(xy,1, outer)
plot(x,y)
points(x[j],y[j],col="green")
})
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
