## Martos and de Carvalho (2018; Fig 1 a)
library(scales)
data(cortical)
m <- 500
n <- 68
plot(cortical$r[,1] * cos(2 * pi * 1:m / m),
cortical$r[,1] * sin(2 * pi * 1:m / m) , type = "l",
col = alpha("gray", 1 / n), xlab = "z", ylab = "x")
for(i in 2:n)
lines(cortical$r[, i] * cos(2 * pi * 1:m / m),
cortical$r[, i] * sin(2 * pi * 1:m / m), type = "l",
col = alpha("gray", i / n))
Run the code above in your browser using DataLab