n <- 100; data <- rnorm(n,1,0.1)
plot(seq(n),data)
hist(data)
N <- data[1:(n-1)]
r <- (data[2:n]-N)/N
plot(N,r)
lm(r~N,as.data.frame(cbind(N,r)))