• The Korean Journal of Applied Statistics
• /
• v.30 no.2
• /
• pp.243-257
• /
• 2017

In this paper, we compare several methods to approximate option prices: Edgeworth expansion, A-type and C-type Gram-Charlier expansions, a method using normal inverse gaussian (NIG) distribution, and an asymptotic method using nonlinear regression. We used two different types of approximation. The first (called the RNM method) approximates the risk neutral probability density function of the log return of the underlying asset and computes the option price. The second (called the OPTIM method) finds the approximate option pricing formula and then estimates parameters to compute the option price. For simulation experiments, we generated underlying asset data from the Heston model and NIG model, a well-known stochastic volatility model and a well-known Levy model, respectively. We also applied the above approximating methods to the KOSPI200 call option price as a real data application. We then found that the OPTIM method shows better performance on average than the RNM method. Among the OPTIM, A-type Gram-Charlier expansion and the asymptotic method that uses nonlinear regression showed relatively better performance; in addition, among RNM, the method of using NIG distribution was relatively better than others.

• ETRI Journal
• /
• v.44 no.5
• /
• pp.805-815
• /
• 2022

The model-based evolutionary algorithms are divided into three groups: estimation of distribution algorithms, inverse modeling, and surrogate modeling. Existing inverse modeling is mainly applied to solve multi-objective optimization problems and is not suitable for many-objective optimization problems. Some inversed-model techniques, such as the inversed-model of multi-objective evolutionary algorithm, constructed from the Pareto front (PF) to the Pareto solution on nondominated solutions using a random grouping method and Gaussian process, were introduced. However, some of the most efficient inverse models might be eliminated during this procedure. Also, there are challenges, such as the presence of many local PFs and developing poor solutions when the population has no evident regularity. This paper proposes inverse modeling using random forest regression and uniform reference points that map all nondominated solutions from the objective space to the decision space to solve many-objective optimization problems. The proposed algorithm is evaluated using the benchmark test suite for evolutionary algorithms. The results show an improvement in diversity and convergence performance (quality indicators).

• The Korean Journal of Applied Statistics
• /
• v.12 no.2
• /
• pp.621-631
• /
• 1999

재료의 피로 파괴과정은 균열의 발생과 전파 및 성장의 과정을 거쳐 마침내 결정적 균열의 크기가 일정한도를 넘어서면 재료의 파괴가 일어난다. 이 때까지의 시간, 즉 피로 수명이 역정규분포를 따를 때 재료의 수명과 스트레스 수준과 관계를 나타 내는 S-N곡선에 대한 대수선형모형(log-linear model)을 제시하고, 이 모형하에서 피로수명시험에 대한 통계적 최적시험설계방법을 찾는다. 통계적 최적여부에 대한 판단기준으로 설계 스트레스 수준하의 특정 시점에서의 신뢰도에 대한 최우추정량의 점근분산을 최소화하는 방법을 사용하였다.

• Journal of the Korean Data and Information Science Society
• /
• v.10 no.1
• /
• pp.243-250
• /
• 1999

We consider $Y=X{\lambda}Z,\;{\lambda}>0$, where X and Z are independent random variables, and Y is the length biased distribution or the equilibrium distribution of X. The purpose of this paper is to consider the distribution of X or Y when the distribution of Z is given and the distribution of Z when the distribution of X or Y is given, In particular, we obtain that the necessary and sufficient conditions for X to be $X^{2}({\upsilon})\;is\;Z{\sim}X^{2}(2)\;and\;for\;Z\;to\;be\;X^{2}(1)\;is\;X{\sim}IG({\mu},\;{\mu}^{2}/{\lambda})$, where $IG({\mu},\;{\mu}^{2}/{\lambda})$ is two-parameter inverse Gaussian distribution. Also we show that X is smaller than Y in the reverse Laplace transform ratio order if and only if $X_{e}$ is smaller than $Y_{e}$ in the Laplace transform ratio order. Finally, we can get the results that if X is smaller than Y in the Laplace transform ratio order, then $Y_{L}$ is smaller than $X_{L}$ in the Laplace transform order, and that if X is smaller than Y in the reverse Laplace transform ratio order, then $_{\mu}X_{L}$ is smaller than $_{\nu}Y_{L}$ in the Laplace transform order.

• The Korean Journal of Applied Statistics
• /
• v.29 no.1
• /
• pp.181-192
• /
• 2016

We consider several methods to approximate option prices with correction terms to the Black-Scholes option price. These methods are able to compute option prices from various risk-neutral distributions using relatively small data and simple computation. In this paper, we compare the performance of Edgeworth expansion, A-type and C-type Gram-Charlier expansions, a method of using Normal inverse gaussian distribution, and an asymptotic method of using nonlinear regression through simulation experiments and real KOSPI200 option data. We assume the variance gamma model in the simulation experiment, which has a closed-form solution for the option price among the pure jump $L{\acute{e}}vy$ processes. As a result, we found that methods to approximate an option price directly from the approximate price formula are better than methods to approximate option prices through the approximate risk-neutral density function. The method to approximate option prices by nonlinear regression showed relatively better performance among those compared.

• The Korean Journal of Applied Statistics
• /
• v.34 no.5
• /
• pp.697-710
• /
• 2021

Relative error prediction is preferred over ordinary prediction methods when relative/percentile errors are regarded as important, especially in econometrics, software engineering and government official statistics. The relative error prediction techniques have been developed in linear/nonlinear regression, nonparametric regression using kernel regression smoother, and stationary time series models. However, random effect models have not been used in relative error prediction. The purpose of this article is to extend relative error prediction to some of generalized linear mixed model (GLMM) with panel data, which is the random effect models based on gamma, lognormal, or inverse gaussian distribution. For better understanding, the real auto insurance data is used to predict the claim size, and the best predictor and the best relative error predictor are comparatively illustrated.

• Nuclear Engineering and Technology
• /
• v.46 no.5
• /
• pp.619-632
• /
• 2014

The Best Estimate Plus Uncertainty (BEPU) method has been widely used to evaluate the uncertainty of a best-estimate thermal hydraulic system code against a figure of merit. This uncertainty is typically evaluated based on the physical model's uncertainties determined by expert judgment. This paper introduces the application of data assimilation methodology to determine the uncertainty bands of the physical models, e.g., the mean value and standard deviation of the parameters, based upon the statistical approach rather than expert judgment. Data assimilation suggests a mathematical methodology for the best estimate bias and the uncertainties of the physical models which optimize the system response following the calibration of model parameters and responses. The mathematical approaches include deterministic and probabilistic methods of data assimilation to solve both linear and nonlinear problems with the a posteriori distribution of parameters derived based on Bayes' theorem. The inverse problem was solved analytically to obtain the mean value and standard deviation of the parameters assuming Gaussian distributions for the parameters and responses, and a sampling method was utilized to illustrate the non-Gaussian a posteriori distributions of parameters. SPACE is used to demonstrate the data assimilation method by determining the bias and the uncertainty bands of the physical models employing Bennett's heated tube test data and Becker's post critical heat flux experimental data. Based on the results of the data assimilation process, the major sources of the modeling uncertainties were identified for further model development.

• Structural Engineering and Mechanics
• /
• v.81 no.1
• /
• pp.103-115
• /
• 2022

The prediction error variances for frequencies are usually considered as unknown in the Bayesian system identification process. However, the error variances for mode shapes are taken as known to reduce the dimension of an identification problem. The present study attempts to explore the effectiveness of Bayesian approach of model parameters updating using Markov Chain Monte Carlo (MCMC) technique considering the prediction error variances for both the frequencies and mode shapes. To remove the ergodicity of Markov Chain, the posterior distribution is obtained by Gaussian Random walk over the proposal distribution. The prior distributions of prediction error variances of modal evidences are implemented through inverse gamma distribution to assess the effectiveness of estimation of posterior values of model parameters. The issue of incomplete data that makes the problem ill-conditioned and the associated singularity problem is prudently dealt in by adopting a regularization technique. The proposed approach is demonstrated numerically by considering an eight-storey frame model with both complete and incomplete modal data sets. Further, to study the effectiveness of the proposed approach, a comparative study with regard to accuracy and computational efficacy of the proposed approach is made with the Sequential Monte Carlo approach of model parameter updating.

• The Korean Journal of Applied Statistics
• /
• v.26 no.3
• /
• pp.483-494
• /
• 2013

The importance of financial risk management has been highlighted after several recent incidences of global financial crisis. One of the issues in financial risk management is how to measure the risk; currently, the most widely used risk measure is the Value at Risk(VaR). We can consider to estimate VaR using extreme value theory if the financial data have heavy tails as the recent market trend. In this paper, we study estimations of VaR using Peaks over Threshold(POT), which is a common method of modeling fat-tailed data using extreme value theory. To use POT, we first estimate parameters of the Generalized Pareto Distribution(GPD). Here, we compare three different methods of estimating parameters of GPD by comparing the performance of the estimated VaR based on KOSPI 5 minute-data. In addition, we simulate data from normal inverse Gaussian distributions and examine two parameter estimation methods of GPD. We find that the recent methods of parameter estimation of GPD work better than the maximum likelihood estimation when the kurtosis of the return distribution of KOSPI is very high and the simulation experiment shows similar results.

• Journal of Korean Society for Quality Management
• /
• v.26 no.1
• /
• pp.74-87
• /
• 1998

This study considers optimal design of accelerated life tests under constant stress using that the first passage time to cross a critical boundary through amount of accumulated degradation has an inverse Gaussian distribution when the degradation process follows to a Brownian motion with positive drift of log linear function of stress. Optimum plans for Type I censoring are derived by minimizing the asymptotic variance of estimated quantiles at the use stress. Sensitivity analyses are also conducted to see how sensitive the optimality criterion is with respect to the uncertainties involved in the guessed values.