• Journal of the Korean Statistical Society
• /
• v.25 no.3
• /
• pp.359-368
• /
• 1996

The maximum likelihood estimation is discussed for the NLAR model with Laplacian marginals. Since the explicit form of the estimates cannot be obtained due to the complicated nature of the likelihood function we utilize the automatic computer optimization subroutine using a direct search complex algorithm. The conditional least square estimates are used as initial estimates in maximum likelihood procedures. The results of a simulation study for the maximum likelihood estimates of the NLAR(1) and the NLAR(2) models are presented.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.127-127
• /
• 2021

In order to estimate parameter uncertainty of hydrological models, the consideration of the likelihood functions which provide reliable parameters of model is necessary. In this study, the Bayesian Markov Chain Monte Carlo (MCMC) method with informal likelihood functions is used to analyze the uncertainty of parameters of the SURR model for estimating the hourly streamflow of Gunnam station of Imjin basin, Korea. Three events were used to calibrate and one event was used to validate the posterior distributions of parameters. Moreover, the performance of four informal likelihood functions (Nash-Sutcliffe efficiency, Normalized absolute error, Index of agreement, and Chiew-McMahon efficiency) on uncertainty of parameter is assessed. The indicators used to assess the uncertainty of the streamflow simulation were P-factor (percentage of observed streamflow included in the uncertainty interval) and R-factor (the average width of the uncertainty interval). The results showed that the sensitivities of parameters strongly depend on the likelihood functions and vary for different likelihood functions. The uncertainty bounds illustrated the slight differences from various likelihood functions. This study confirms the importance of the likelihood function selection in the application of Bayesian MCMC to the uncertainty assessment of the SURR model.

• Communications for Statistical Applications and Methods
• /
• v.5 no.2
• /
• pp.559-570
• /
• 1998

This paper introduces doubly Winsorized Poisson auto-model by truncating the support of a Poisson random variable both from above and below, and shows that this model has a same form of negpotential function as regular Poisson auto-model and one-way Winsorized Poisson auto-model. Strategies for maximum likelihood estimation of parameters are discussed. In addition to exact maximum likelihood estimation, Monte Carlo maximum likelihood estimation may be applied to this model.

• Communications for Statistical Applications and Methods
• /
• v.21 no.3
• /
• pp.253-260
• /
• 2014

The intraclass correlation model has a long history of applications in several fields of research. Case deletion diagnostic methods for the intraclass correlation model are proposed. Based on the likelihood equations, we derive a formula for a case deletion diagnostic method which enables us to investigate the influence of observations on the maximum likelihood estimates of the model parameters. Using the Taylor series expansion we develop an approximation to the likelihood distance. Numerical examples are provided for illustration.

• Journal of the Korean Data and Information Science Society
• /
• v.27 no.5
• /
• pp.1389-1397
• /
• 2016

Semi-parametric frailty model with nonparametric baseline hazards has been widely used for the analyses of clustered survival-time data. The frailty models can be fitted via an auxiliary Poisson hierarchical generalized linear model (HGLM). For the inferences of the frailty model marginal likelihood, which gives MLE, is often used. The marginal likelihood is usually obtained by integrating out random effects, but it often requires an intractable integration. In this paper, we propose to obtain the MLE via Laplace approximation using a Poisson HGLM approach for semi-parametric frailty model. The proposed HGLM approach uses hierarchical-likelihood (h-likelihood), which avoids integration itself. The proposed method is illustrated using a numerical study.

• The Korean Journal of Applied Statistics
• /
• v.13 no.2
• /
• pp.407-414
• /
• 2000

Carroll and Ruppert(1988) analyzed the esterase assay data with regression model based on quasi-likelihood. Jung and Lee(1997) introduced a goodness-of-fit test for testing the adequacy of the quasi-likelihood and claimed that there is no gross inadequacy with the model because their test was not rejected. However, Lee and Xelder(199S)'s residual plots revealed that the model did not sufficiently reflect the increase of the variance with that of the mean. In this paper, we re-analyze the esterase assay data with the joint modelling of mean and dispersion in Lee and l\elder(1998) and evaluate the validity of the fitted model by applying the residual plots. And it is illustrated that Lee and Nelder(199S)'s restricted likelihood is more efficient in goodness-of-fit test for the dispersion model.

• Communications for Statistical Applications and Methods
• /
• v.25 no.4
• /
• pp.355-371
• /
• 2018

The two parameter negative exponential distribution has many practical applications in queuing theory such as the service times of agents in system, the time it takes before your next telephone call, the time until a radioactive practical decays, the distance between mutations on a DNA strand, and the extreme values of annual snowfall or rainfall; consequently, has many applications in reliability systems. This paper considers an estimation problem of stress-strength model with two parameter negative parameter exponential distribution. We introduce a maximum penalized likelihood method, Bayes estimator using Lindley approximation to estimate stress-strength model and compare the proposed estimators with regular maximum likelihood estimator for complete data. We also introduce a maximum penalized likelihood method, Bayes estimator using a Markov chain Mote Carlo technique for incomplete data. A Monte Carlo simulation study is performed to compare stress-strength model estimates. Real data is used as a practical application of the proposed model.

• 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 Data and Information Science Society
• /
• v.15 no.1
• /
• pp.251-261
• /
• 2004

The paper focuses primarily on the standard linear multiple regression model where the parameter of interest is a ratio of two regression coefficients. The general model includes the calibration model, the Fieller-Creasy problem, slope-ratio assays, parallel-line assays, and bioequivalence. We provide an orthogonal transformation (cf. Cox and Reid (1987)) of the original parameter vector. Also, we give some remarks on the difficulties associated with likelihood based confidence interval.

• Communications for Statistical Applications and Methods
• /
• v.9 no.2
• /
• pp.515-520
• /
• 2002

In general, we use the hierarchical Poisson-gamma model for the Poisson data in generalized linear model. Time effect will be emphasized for the analysis of the observed data to be collected annually for the time period. An extended model with time effect for estimating the effect is proposed. In particularly, we discuss the Quasi likelihood function which is used to numerical approximation for the likelihood function of the parameter.