JGQD.density()
approximates the transition density of a scalar generalized quadratic diffusion model (GQD). Given an initial value for the diffusion, Xs
, the approximation is evaluated for all Xt
at equispaced time-nodes given by splitting [s
, t
] into subintervals of length delt
.
JGQD.density() approximates transitional densities of jump diffusions of the form:ScalarEqn1.png
where
ScalarEqn2.png
ScalarEqn3.png
and
ScalarEqn4.png
describes a Poisson process with jumps of the form:
ScalarEqn6.png
arriving with intensity
ScalarEqn5.png
subject to a jump distribition of the form:
ScalarEqn7.png
JGQD.density(Xs = 4, Xt = seq(5, 8, 1/10), s = 0, t = 5, delt =1/100, Jdist = "Normal", Jtype = "Add", Dtype = "Saddlepoint", Trunc = c(8, 4), factorize = FALSE, factor.type = "Diffusion", beta, print.output = TRUE, eval.density = TRUE)
'Saddlepoint'
and 'Edgeworth'
are supported (default = 'Saddlepoint'
).Trunc[1] >= Trunc[2]
. Default is c(4,4)
.TRUE
, model information is printed to the console.TRUE
, the density is evaluated in addition to calculating the moment eqns.Dtype
. Warning [2]:
The parameter delt
is also used as the stepsize for solving a system of ordinary differential equations (ODEs) that govern the evolution of the cumulants of the diffusion. As such delt
is required to be small for highly non-linear models in order to ensure sufficient accuracy.JGQD.density
constructs an approximate transition density for a class of quadratic diffusion models. This is done by first evaluating the trajectory of the cumulants/moments of the diffusion numerically as the solution of a system of ordinary differential equations over a time horizon [s,t]
split into equi-distant points delt
units apart. Subsequently, the resulting cumulants/moments are carried into a density approximant (by default, a saddlepoint approximation) in order to evaluate the transtion surface.
Daniels, H.E. 1954 Saddlepoint approximations in statistics. Ann. Math. Stat., 25:631--650.
Eddelbuettel, D. and Romain, F. 2011 Rcpp: Seamless R and C++ integration. Journal of Statistical Software, 40(8):1--18,. URL http://www.jstatsoft.org/v40/i08/.
Eddelbuettel, D. 2013 Seamless R and C++ Integration with Rcpp. New York: Springer. ISBN 978-1-4614-6867-7.
Eddelbuettel, D. and Sanderson, C. 2014 Rcpparmadillo: Accelerating R with high-performance C++ linear algebra. Computational Statistics and Data Analysis, 71:1054--1063. URL http://dx.doi.org/10.1016/j.csda.2013.02.005.
Feagin, T. 2007 A tenth-order Runge-Kutta method with error estimate. In Proceedings of the IAENG Conf. on Scientifc Computing.
Varughese, M.M. 2013 Parameter estimation for multivariate diffusion systems. Comput. Stat. Data An., 57:417--428.
JGQD.mcmc
and BiJGQD.density
.
#===============================================================================
# For detailed notes and examples on how to use the JGQD.density() function, see
# the following vignette:
RShowDoc('Part_2_Transition_Densities',type='html','DiffusionRjgqd')
#===============================================================================
Run the code above in your browser using DataLab