Tag Info

Hot answers tagged

3

The error is, you are not storing the random numbers for the same path at the end: xbefore = x + c*tau + sigma*sqrt(tau)*randn() A = muA + sigmaA*randn(); xafter = xbefore + A; But then at end you set a different path here by creating a new random number: xT = log(S0)+(c+muA*lambda)*T+sqrt((sigma^2+(muA^2+sigmaA^2)*lambda)*T)*randn(); randn() ...


1

Look at Gilli & Schumann's paper. They provide a Bates' model estimates set, the way to improve such estimates calibrating those ones using an Heuristic model and, lastly, the relative codes in matlab, in order to be able to replicate the model. Unfortunately, there are not available the relative call prices estimated time series; I think that noone ...


1

What you do is: You simulate a Brownian path - with the correct standard deviation. Then you simulate the Compound Poisson process. In each time step you sample a jump or no jump and the jump size if there was one. In each time step you draw from a Bernoulli distribution - which is to my knowledge just an approximation. For the compound Poisson process I ...


1

The Feller condition applies without modification. That is under the assumption that $v$ is square-root process with poisson-arrival jumps (as you wrote), and assuming the jump distribution is strictly positive and initial level $v_0>0$. The reason is, conditional on no jumps occuring, the process is just a square root process, for which the references ...



Only top voted, non community-wiki answers of a minimum length are eligible