Image:HittingTimes1.png

From Wikipedia, the free encyclopedia

Wikimedia Commons logo This is a file from the Wikimedia Commons. The description on its description page there is shown below.
Commons is attempting to create a freely licensed media file repository. You can help.
Description

some stopping times (even hitting times!) of brownian motion
Source

with the help of GNU R statistics / math software. see the source below
Date

8 Apr. 2006
Author

Thomas Steiner
Permission

Thomas Steiner put it under the GFDL


GNU head Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version published by the Free Software Foundation; with no Invariant Sections, no Front-Cover Texts, and no Back-Cover Texts. A copy of the license is included in the section entitled "GNU Free Documentation License".

العربية | Asturianu | Български | বাংলা | ইমার ঠার | Brezhoneg | Bosanski | Català | Sinugboanong Binisaya | Česky | Dansk | Deutsch | Ελληνικά | English | Esperanto | Español | Eesti | Euskara | فارسی | Suomi | Français | Galego | עברית | Hrvatski | Magyar | Bahasa Indonesia | Ido | Íslenska | Italiano | 日本語 | ქართული | 한국어 | Kurdî / كوردي | Latina | Lëtzebuergesch | Lietuvių | Bahasa Melayu | Nnapulitano | Nederlands | Norsk (nynorsk) | Norsk (bokmål) | Polski | Português | Română | Русский | Slovenčina | Slovenščina | Shqip | Српски | Svenska | తెలుగు | ไทย | Türkçe | Українська | اردو | Tiếng Việt | Volapük | 中文(简体) | 中文(繁體) | +/-

R-source code:

T=7
N=100000
set.seed(303898)
a=0.9
cs=c("orange", "navy", "seagreen4", "grey31")

bm1=c(0,cumsum(rnorm(N,mean=0,sd=sqrt(T/N))))
tau1=which(bm1>a)[1]
bm1[tau1:(N+1)]=a

bm2=c(0,cumsum(rnorm(N,mean=0,sd=sqrt(T/N))))
tau2=which(bm2>a)[1]
bm2[tau2:(N+1)]=a

bm3=c(0,cumsum(rnorm(N,mean=0,sd=sqrt(T/N))))
tau3=which(bm3>a)[1]
bm3[tau3:(N+1)]=a

png(filename = "HittingTimes1.png", width=1200, height=900, pointsize = 12)
par(bg="gray88")
t=seq(0,T,length=N+1)
plot(t,bm1,type="l",ylim=range(bm1,bm2,bm3,a*1.2),col=cs[1],lwd=2,ylab="",xlab="Zeit",main="Hitting time von drei Sampels der Brownschen Bewegung als Beispiel für eine Stoppzeit")
lines(t,bm2,col=cs[2],lwd=2)
lines(t,bm3,col=cs[3],lwd=2)
abline(h=a,col=cs[4],lty=1,lwd=3)

text(x=tau1*T/N,y=a*1.14, labels = substitute(tau[2]==t1,list(t1=tau1*T/N)),pos=4,col=cs[1],cex=1.5)
text(x=tau2*T/N,y=a*1.14, labels = substitute(tau[1]==t2,list(t2=tau2*T/N)),pos=4,col=cs[2],cex=1.5)
text(x=T,y=a*1.14, labels = substitute(tau[3]>T,list(T=T)),pos=2,col=cs[3],cex=1.5)
dev.off()

The following pages on the English Wikipedia link to this file (pages on other projects are not listed):

Retrieved from "http://en.wikipedia.org../../../h/i/t/Image%7EHittingTimes1.png_de63.html"