A Note on Monte Carlo Greeks using the Characteristic Function

Author(s):  
Joerg Kienitz
2002 ◽  
Vol 18 (2) ◽  
pp. 420-468 ◽  
Author(s):  
Oliver Linton ◽  
Yoon-Jae Whang

We introduce a kernel-based estimator of the density function and regression function for data that have been grouped into family totals. We allow for a common intrafamily component but require that observations from different families be independent. We establish consistency and asymptotic normality for our procedures. As usual, the rates of convergence can be very slow depending on the behavior of the characteristic function at infinity. We investigate the practical performance of our method in a simple Monte Carlo experiment.


2002 ◽  
Vol 18 (5) ◽  
pp. 1197-1220
Author(s):  
Eiji Kurozumi

This paper investigates the limiting properties of the Canova and Hansen test, testing for the null hypothesis of no unit root against seasonal unit roots, under a sequence of local alternatives with the model extended to have seasonal dummies and trends or no deterministic term and also only seasonal dummies. We derive the limiting distribution of the test statistic and its characteristic function under local alternatives. We find that the local limiting power is an inverse function of the spectral density at frequency π (π/2) when we test against a negative unit root (annual unit roots). We also theoretically show that the local limiting power of the Canova and Hansen test against a negative unit root (annual unit roots) does not increase when the true process has annual unit roots (a negative unit root) but not a negative unit root (annual unit roots), which has been observed in Monte Carlo simulations in such research as Caner (1998, Journal of Business and Economic Statistics 16, 349–356), Canova and Hansen (1995, Journal of Business and Economic Statistics 13, 237–252), and Hylleberg (1995, Journal of Econometrics 69, 5–25).


2002 ◽  
Vol 18 (3) ◽  
pp. 691-721 ◽  
Author(s):  
John L. Knight ◽  
Jun Yu

Because the empirical characteristic function (ECF) is the Fourier transform of the empirical distribution function, it retains all the information in the sample but can overcome difficulties arising from the likelihood. This paper discusses an estimation method via the ECF for strictly stationary processes. Under some regularity conditions, the resulting estimators are shown to be consistent and asymptotically normal. The method is applied to estimate the stable autoregressive moving average (ARMA) models. For the general stable ARMA model for which the maximum likelihood approach is not feasible, Monte Carlo evidence shows that the ECF method is a viable estimation method for all the parameters of interest. For the Gaussian ARMA model, a particular stable ARMA model, the optimal weight functions and estimating equations are given. Monte Carlo studies highlight the finite sample performances of the ECF method relative to the exact and conditional maximum likelihood methods.


2016 ◽  
Vol 2016 ◽  
pp. 1-8
Author(s):  
Wanwan Huang ◽  
Brian Ewald ◽  
Giray Ökten

The coupled additive and multiplicative (CAM) noises model is a stochastic volatility model for derivative pricing. Unlike the other stochastic volatility models in the literature, the CAM model uses two Brownian motions, one multiplicative and one additive, to model the volatility process. We provide empirical evidence that suggests a nontrivial relationship between the kurtosis and skewness of asset prices and that the CAM model is able to capture this relationship, whereas the traditional stochastic volatility models cannot. We introduce a control variate method and Monte Carlo estimators for some of the sensitivities (Greeks) of the model. We also derive an approximation for the characteristic function of the model.


1974 ◽  
Vol 22 ◽  
pp. 307 ◽  
Author(s):  
Zdenek Sekanina

AbstractIt is suggested that the outbursts of Periodic Comet Schwassmann-Wachmann 1 are triggered by impacts of interplanetary boulders on the surface of the comet’s nucleus. The existence of a cloud of such boulders in interplanetary space was predicted by Harwit (1967). We have used the hypothesis to calculate the characteristics of the outbursts – such as their mean rate, optically important dimensions of ejected debris, expansion velocity of the ejecta, maximum diameter of the expanding cloud before it fades out, and the magnitude of the accompanying orbital impulse – and found them reasonably consistent with observations, if the solid constituent of the comet is assumed in the form of a porous matrix of lowstrength meteoric material. A Monte Carlo method was applied to simulate the distributions of impacts, their directions and impact velocities.


1988 ◽  
Vol 102 ◽  
pp. 79-81
Author(s):  
A. Goldberg ◽  
S.D. Bloom

AbstractClosed expressions for the first, second, and (in some cases) the third moment of atomic transition arrays now exist. Recently a method has been developed for getting to very high moments (up to the 12th and beyond) in cases where a “collective” state-vector (i.e. a state-vector containing the entire electric dipole strength) can be created from each eigenstate in the parent configuration. Both of these approaches give exact results. Herein we describe astatistical(or Monte Carlo) approach which requires onlyonerepresentative state-vector |RV> for the entire parent manifold to get estimates of transition moments of high order. The representation is achieved through the random amplitudes associated with each basis vector making up |RV>. This also gives rise to the dispersion characterizing the method, which has been applied to a system (in the M shell) with≈250,000 lines where we have calculated up to the 5th moment. It turns out that the dispersion in the moments decreases with the size of the manifold, making its application to very big systems statistically advantageous. A discussion of the method and these dispersion characteristics will be presented.


Author(s):  
Ryuichi Shimizu ◽  
Ze-Jun Ding

Monte Carlo simulation has been becoming most powerful tool to describe the electron scattering in solids, leading to more comprehensive understanding of the complicated mechanism of generation of various types of signals for microbeam analysis.The present paper proposes a practical model for the Monte Carlo simulation of scattering processes of a penetrating electron and the generation of the slow secondaries in solids. The model is based on the combined use of Gryzinski’s inner-shell electron excitation function and the dielectric function for taking into account the valence electron contribution in inelastic scattering processes, while the cross-sections derived by partial wave expansion method are used for describing elastic scattering processes. An improvement of the use of this elastic scattering cross-section can be seen in the success to describe the anisotropy of angular distribution of elastically backscattered electrons from Au in low energy region, shown in Fig.l. Fig.l(a) shows the elastic cross-sections of 600 eV electron for single Au-atom, clearly indicating that the angular distribution is no more smooth as expected from Rutherford scattering formula, but has the socalled lobes appearing at the large scattering angle.


Author(s):  
D. R. Liu ◽  
S. S. Shinozaki ◽  
R. J. Baird

The epitaxially grown (GaAs)Ge thin film has been arousing much interest because it is one of metastable alloys of III-V compound semiconductors with germanium and a possible candidate in optoelectronic applications. It is important to be able to accurately determine the composition of the film, particularly whether or not the GaAs component is in stoichiometry, but x-ray energy dispersive analysis (EDS) cannot meet this need. The thickness of the film is usually about 0.5-1.5 μm. If Kα peaks are used for quantification, the accelerating voltage must be more than 10 kV in order for these peaks to be excited. Under this voltage, the generation depth of x-ray photons approaches 1 μm, as evidenced by a Monte Carlo simulation and actual x-ray intensity measurement as discussed below. If a lower voltage is used to reduce the generation depth, their L peaks have to be used. But these L peaks actually are merged as one big hump simply because the atomic numbers of these three elements are relatively small and close together, and the EDS energy resolution is limited.


Sign in / Sign up

Export Citation Format

Share Document