An Eager Avocado

Eager Avocado

I give myself very good advice, but I very seldom follow it.

Plotting large matrix in R

,

a = matrix(rnorm(10000^2),ncol=10000)
Ns = c(1000, 5000,10000)
times = data.frame(list(N=Ns))
t1 = sapply(Ns, function(n) {system.time(image(a[1:n,1:n],useRaster=T))['elapsed']})

plot of chunk unnamed-chunk-1plot of chunk unnamed-chunk-1plot of chunk unnamed-chunk-1

times$raster <- t1

plot of chunk unnamed-chunk-2