runlength=function(x){

run=1

k=1

for( i in 2: length(x)){

if(x[i]==x[i-1]) run[k]=run[k]+1

else{

k=k+1

run=c(run,1)

}

}

return(run)

}

skip to main |
skip to sidebar
## Monday, July 23, 2007

###
Run Length

## Wednesday, April 18, 2007

###
Mandelbrot Set with R Animation

## Thursday, April 5, 2007

## Wednesday, April 4, 2007

###
My plan

A small function to get runlength from a vector

runlength=function(x){

run=1

k=1

for( i in 2: length(x)){

if(x[i]==x[i-1]) run[k]=run[k]+1

else{

k=k+1

run=c(run,1)

}

}

return(run)

}

runlength=function(x){

run=1

k=1

for( i in 2: length(x)){

if(x[i]==x[i-1]) run[k]=run[k]+1

else{

k=k+1

run=c(run,1)

}

}

return(run)

}

A few days ago I had the opportunity of having dinner with **BenoĆ®t Mandelbrot **

when he gave a seminar talk at the Bu Mathematics and Statistics Department. At the same time I was trying to find a way to produce animation in R. I came across this great post which was relevant for both purposes.

Using the code

### Reproduced from http://tolstoy.newcastle.edu.au/R/help/05/10/13198.html

### Written by Jarek Tuszynski, PhD.

library(fields) # for tim.colors

library(caTools) # for write.gif

m = 400 # grid size

C = complex( real=rep(seq(-1.8,0.6, length.out=m), each=m ),

imag=rep(seq(-1.2,1.2, length.out=m), m ) )

C = matrix(C,m,m)

Z = 0

X = array(0, c(m,m,20))

for (k in 1:20) {

Z = Z^2+C

X[,,k] = exp(-abs(Z))

}

image(X[,,k], col=tim.colors(256)) # show final image in

write.gif(X, "Mandelbrot.gif", col=tim.colors(256), delay=100)

we get the following animation

when he gave a seminar talk at the Bu Mathematics and Statistics Department. At the same time I was trying to find a way to produce animation in R. I came across this great post which was relevant for both purposes.

Using the code

### Reproduced from http://tolstoy.newcastle.edu.au/R/help/05/10/13198.html

### Written by Jarek Tuszynski, PhD.

library(fields) # for tim.colors

library(caTools) # for write.gif

m = 400 # grid size

C = complex( real=rep(seq(-1.8,0.6, length.out=m), each=m ),

imag=rep(seq(-1.2,1.2, length.out=m), m ) )

C = matrix(C,m,m)

Z = 0

X = array(0, c(m,m,20))

for (k in 1:20) {

Z = Z^2+C

X[,,k] = exp(-abs(Z))

}

image(X[,,k], col=tim.colors(256)) # show final image in

write.gif(X, "Mandelbrot.gif", col=tim.colors(256), delay=100)

we get the following animation

I am planning to maintain this blog to publish the every day tricks I learn working in R. Please feel free to post your comments and share interesting tricks

Subscribe to:
Posts (Atom)