domingo, 15 de mayo de 2016

8-Caminatas aleatorias

#Ejecutar antes
library(graphics)

caminatas=function(X0,n,p,simulaciones){
  trayXn=matrix(0,(n+1),simulaciones)
  for(i in 1:simulaciones){
    trayXn[,i]=cumsum(c(X0,(2*rbinom(n,1,p)-1)))
  }
  matplot(0:n,trayXn,type="o",pch=16,ylim=c(-n,n),cex=1.5,lwd=2)
  abline(h=-n:n,v=0:n,lty=3,col="gray")
  #trayXn
}

#Ejemplo
caminatas(1,10,.35,5)