• Communications for Statistical Applications and Methods
• /
• v.26 no.4
• /
• pp.371-383
• /
• 2019

Panel data sets have been developed in various areas, and many recent studies have analyzed panel, or longitudinal data sets. Maximum likelihood (ML) may be the most common statistical method for analyzing panel data models; however, the inference based on the ML estimate will have an inflated Type I error because the ML method tends to give a downwardly biased estimate of variance components when the sample size is small. The under estimation could be severe when data is incomplete. This paper proposes the restricted maximum likelihood (REML) method for a random effects panel data model with a censored dependent variable. Note that the likelihood function of the model is complex in that it includes a multidimensional integral. Many authors proposed to use integral approximation methods for the computation of likelihood function; however, it is well known that integral approximation methods are inadequate for high dimensional integrals in practice. This paper introduces to use the moments of truncated multivariate normal random vector for the calculation of multidimensional integral. In addition, a proper asymptotic standard error of REML estimate is given.

• Journal of the Korean Society of Safety
• /
• v.37 no.6
• /
• pp.128-135
• /
• 2022

The comparative evaluation of occupational accident fatality rates (OAFRs) of different countries is complicated owing to the differences in their level of socio-economic development. However, such evaluation is necessary to assess the national occupational safety and health system of a country. This study proposes a statistical method to compare the OAFRs of countries taking into consideration the difference in their level of socio-economic development. We first collected data on the socio-economic indicators and OAFRs of 11 countries over a 30-year period. Next, based on literature survey and statistical correlation analysis, we selected the significant independent variables and built multiple linear regression models to predict OAFR. We also determined the groups of countries having heterogeneous relationships between the independent variables and OAFRs, which are represented by the regression models. The proposed method is demonstrated by comparing the OAFR of Korea with the OAFRs of 10 other developed countries.

• The Korean Journal of Applied Statistics
• /
• v.33 no.3
• /
• pp.297-308
• /
• 2020

We focus on the functional autoregressive conditional heteroscedasticity (fARCH) modelling to analyze intraday volatilities based on high frequency financial time series. Multivariate volatility models are investigated to approximate fARCH(1). A formula of multi-step ahead volatilities for fARCH(1) model is derived. As an application, in implementing fARCH(1), a choice of appropriate time interval for the intraday return is discussed. High frequency KOSPI data analysis is conducted to illustrate the main contributions of the article.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.3
• /
• pp.40-45
• /
• 2018

Intelligent monitoring and diagnosis of production processes based on multivariate statistical methods has been one of important tasks for safety and quality issues. This is due to the fact that faults and unexpected events may have serious impacts on the operation of processes. This study proposes a diagnostic scheme based on effective representation of process measurement data and is evaluated using simulation process data. The effects of utilizing a preprocessing step and nonlinear statistical methods are also tested using fifteen faults of the simulation process. Results show that the proposed scheme produced more reliable results and outperformed other tested schemes with none of the filtering step and nonlinear methods. The proposed scheme is expected to be robust to process noises and easy to develop due to the lack of required rigorous mathematical process models or expert knowledge.

• The Korean Journal of Applied Statistics
• /
• v.28 no.2
• /
• pp.309-324
• /
• 2015

In this paper we introduce how to analyze survival data via a SRC-Stat statistical package. This provides classical survival analysis (e.g. Cox's proportional hazards models for univariate survival data) as well as advanced survival analysis such as shared and nested frailty models for multivariate survival data. We illustrate the use of our package with practical data sets.

• The Korean Journal of Applied Statistics
• /
• v.22 no.5
• /
• pp.995-1005
• /
• 2009

Conditional correlation between financial time series plays an important role in risk management, asset allocation and portfolio selection and therefore diverse efforts for modeling conditional correlations in multivariate-GARCH processes have been made in last two decades. In particular, CCC (cf. Bollerslev, 1990) and DCC(dynamic conditional correlation, cf. Engle, 2002) models have been commonly used since they are relatively parsimonious in the number of parameters involved. This article is concerned with DCC modeling for multivariate GARCH processes in comparison with CCC specification. Various multivariate financial time series are analysed to illustrate possible advantages of DCC over CCC modeling.

• The Korean Journal of Applied Statistics
• /
• v.29 no.4
• /
• pp.753-771
• /
• 2016

Value at Risk (VaR) is widely used as an important tool for risk management of financial institutions. In this paper we discuss estimation and back testing for VaR of the portfolio composed of KOSPI, Dow Jones, Shanghai, Nikkei indexes. The copula functions are adopted to construct the multivariate distributions of portfolio components from marginal distributions that combine extreme value theory and GARCH models. Volatility models with t distribution of the error terms using Gaussian, t, Clayton and Frank copula functions are shown to be more appropriate than the other models, in particular the model using the Frank copula is shown to be the best.

• Journal of the Korean Statistical Society
• /
• v.20 no.2
• /
• pp.191-201
• /
• 1991

Given a set of data pxN$_{i}$, matrices X$_{i}$ observed from p-variate normal populations $\prod$$_{i}$~N($\mu$$_{I}$, $\Sigma$$_{i}$) for i=1, …, K, the test for equality form of the covariance matrices is to choose a hypothetical model which best explains the homogeneity/heterogeneity structure across the covariance matrices among the hypothesized class of models. This paper describes a test procedure for selecting the best model. The procedure is based on a synthesis of Bayesian and a cross-validation or sample reuse methodology that makes use of a one-at-a-time schema of observational omissions. Advantages of the test are argued on two grounds, and illustrative examples and simulation results are given.are given.

• Communications for Statistical Applications and Methods
• /
• v.22 no.6
• /
• pp.615-624
• /
• 2015

We consider the estimation of seasonal cointegration in the presence of conditional heteroskedasticity (CH) using a feasible generalized least squares method. We capture cointegrating relationships and time-varying volatility for long-run and short-run dynamics in the same model. This procedure can be easily implemented using common methods such as ordinary least squares and generalized least squares. The maximum likelihood (ML) estimation method is computationally difficult and may not be feasible for larger models. The simulation results indicate that the proposed method is superior to the ML method when CH exists. In order to illustrate the proposed method, an empirical example is presented to model a seasonally cointegrated times series under CH.

• Journal of the Korean Statistical Society
• /
• v.31 no.2
• /
• pp.177-198
• /
• 2002

The estimation of variance components or variance ratios in linear model is an important issue in plant or animal breeding fields, and various estimation methods have been devised to estimate variance components or variance ratios. However, many traits of economic importance in those fields are observed as dichotomous or polychotomous outcomes. The usual estimation methods might not be appropriate for these cases. Recently threshold linear model is considered as an important tool to analyze discrete traits specially in animal breeding field. In this note, we consider a hierarchical Bayesian method for the threshold animal model. Gibbs sampler for making full Bayesian inferences about random effects as well as fixed effects is described to analyze jointly discrete traits and continuous traits. Numerical example of the model with two discrete ordered categorical traits, calving ease of calves from born by heifer and calving ease of calf from born by cow, and one normally distributed trait, birth weight, is provided.