Tell me more ×
Cross Validated is a question and answer site for statisticians, data analysts, data miners and data visualization experts. It's 100% free, no registration required.

I have simulated one possible path of a variance gamma process by the following code:

vektor<-c(1:23)

S0=20
theta=0.01

v=5
sigma=0.1

vektor[1]<-S0

for (i in 2:23){
randomgamma<-rgamma(1, shape=1/v, scale = v)
randomnormal<-rnorm(1,mean=0,sd=1)
vektor[i]<-vektor[i-1]+theta*randomgamma+sigma*sqrt(randomgamma)*randomnormal
}

plot(c(1:23),vektor)
lines(c(1:23),vektor)

The idea is to be found on page 26 in the following paper: http://www.rhsmith.umd.edu/faculty/mfu/fu_files/Fu07.pdf

Now my problem is, that the plot does not look like a variance gamma process, these should look like this: http://en.wikipedia.org/wiki/File:Variance-Gamma-process.png

or this demonstration.

So where is my mistake?

In general: Is what I am doing correct? I want to simulate a stock path. The initial value of the stock is 20. Now, I want to simulate different paths. What parameters should I use to get a realistic result?

share|improve this question
first of all, there are 3 different methods on that page, which are you attempting? – jerad Nov 25 '12 at 22:15
also, what's your plot look like? – jerad Nov 25 '12 at 22:20
I noticed an edit to this post in reply to @jerad's comments (indicating the 1st method was considered). If you lost your account information, please flag your post for moderator attention and we will merge your accounts. – chl Nov 26 '12 at 8:31

Know someone who can answer? Share a link to this question via email, Google+, Twitter, or Facebook.

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.