• Journal of the Korean Data and Information Science Society
• /
• 제13권2호
• /
• pp.329-340
• /
• 2002

Two methods are proposed for constructing a confidence interval on the among group variance component in a unbalanced one-way random effects model. Computer simulation is used to compare these methods with alternative procedures. The results indicate that the method1 and methods2 perform well over small group size and large sample size respectively.

• Communications for Statistical Applications and Methods
• /
• 제7권1호
• /
• pp.141-150
• /
• 2000

In this paper we propose a method of computing Bayes factor to detect an outlier in a random effects model. When no information is available and hence improper noninformative priors should be used Bayes factor includes the unspecified constants and has complicated computational burden. To solve this problem we use the fractional Bayes factor (FBF) of O-Hagan(1995) and the generalized Savage0-Dickey density ratio of Verdinelli and Wasserman (1995) The proposed method is applied to outlier deterction problem We perform a simulation of the proposed approach with a simulated data set including an outlier and also analyze a real data set.

• Communications for Statistical Applications and Methods
• /
• 제25권1호
• /
• pp.61-70
• /
• 2018

Modeling of the random effects covariance matrix in generalized linear mixed models (GLMMs) is an issue in analysis of longitudinal categorical data because the covariance matrix can be high-dimensional and its estimate must satisfy positive-definiteness. To satisfy these constraints, we consider the autoregressive and moving average Cholesky decomposition (ARMACD) to model the covariance matrix. The ARMACD creates a more flexible decomposition of the covariance matrix that provides generalized autoregressive parameters, generalized moving average parameters, and innovation variances. In this paper, we analyze longitudinal count data with overdispersion using GLMMs. We propose negative binomial loglinear mixed models to analyze longitudinal count data and we also present modeling of the random effects covariance matrix using the ARMACD. Epilepsy data are analyzed using our proposed model.

• Communications for Statistical Applications and Methods
• /
• 제27권5호
• /
• pp.523-533
• /
• 2020

This paper suggests three available methods for finding nonnegative estimates of variance components of the random effects in mixed models. The three proposed methods based on the concepts of projections are called projection method I, II, and III. Each method derives sums of squares uniquely based on its own method of projections. All the sums of squares in quadratic forms are calculated as the squared lengths of projections of an observation vector; therefore, there is discussion on the decomposition of the observation vector into the sum of orthogonal projections for establishing a projection model. The projection model in matrix form is constructed by ascertaining the orthogonal projections defined on vector subspaces. Nonnegative estimates are then obtained by the projection model where all the coefficient matrices of the effects in the model are orthogonal to each other. Each method provides its own system of linear equations in a different way for the estimation of variance components; however, the estimates are given as the same regardless of the methods, whichever is used. Hartley's synthesis is used as a method for finding the coefficients of variance components.

• Communications for Statistical Applications and Methods
• /
• 제22권6호
• /
• pp.625-637
• /
• 2015

We develop a Bayesian clustering procedure based on a Dirichlet process prior with cluster specific random effects. Gibbs sampling of a normal mixture of linear mixed regressions with a Dirichlet process was implemented to calculate posterior probabilities when the number of clusters was unknown. Our approach (unlike its counterparts) provides simultaneous partitioning and parameter estimation with the computation of the classification probabilities. A Monte Carlo study of curve estimation results showed that the model was useful for function estimation. We find that the proposed Dirichlet process mixture model with cluster specific random effects detects clusters sensitively by combining vague edges into different clusters. Examples are given to show how these models perform on real data.

• Communications for Statistical Applications and Methods
• /
• 제26권4호
• /
• pp.337-346
• /
• 2019

This paper discusses a method for obtaining nonnegative estimates for variance components in a random effects model. A variance component should be positive by definition. Nevertheless, estimates of variance components are sometimes given as negative values, which is not desirable. The proposed method is based on two basic ideas. One is the identification of the orthogonal vector subspaces according to factors and the other is to ascertain the projection in each orthogonal vector subspace. Hence, an observation vector can be denoted by the sum of projections. The method suggested here always produces nonnegative estimates using projections. Hartley's synthesis is used for the calculation of expected values of quadratic forms. It also discusses how to set up a residual model for each projection.

• Communications for Statistical Applications and Methods
• /
• 제26권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
• /
• 제15권2호
• /
• pp.379-386
• /
• 2004

This paper deals with a mixed logit model for vaccination data. The effect of a newly developed vaccine for a certain chicken disease can be evaluated by a noninfection rate after injecting chicken with the disease vaccine. But there are a lot of factors that might affect the noninfecton rate. Some of these are fixed and others are random. Random factors are sometimes coming from the sampling scheme for choosing experimental units. This paper suggests a mixed model when some fixed factors need to have different experimental sizes by an experimental design and illustrates how to estimate parameters in a suggested model.

• Structural Engineering and Mechanics
• /
• 제90권1호
• /
• pp.41-49
• /
• 2024

In recent years, China's high-speed railway (HSR) line continues to expand into seismically active regions. Analyzing the features of earthquake rail irregularity is crucial in this situation. This study first established and experimentally validated a finite element (FE) model of bridge-track. The FE model was then combined with earthquake record database to generate the earthquake rail irregularity library. The sample library was used to construct a model of desired earthquake rail irregularity based on signal processing (SFT) and hypothesis principle. Finally, the effects of random pier height and random span number on desired irregularity were analyzed. Herein, an equivalent method of calculating earthquake rail irregularities for random structures was proposed. The results of this study show that the amplitude of desired irregularity is found to increase with increasing pier height. When calculating the desired irregularity of a structure with unequal pier heights, the structure can be regarded as that with equal pier heights (taking the largest pier height). For a structure with the span number large than 9, its desired irregularity can be considered equal to that of a 9-span structure. For the structures with both random pier heights and random span number, their desired irregularities are obtained by equivalent calculations for pier height and span number, respectively.

• 한국통신학회논문지
• /
• 제13권4호
• /
• pp.332-339
• /
• 1988

최근에 소개된, 순가산성 잡음뿐만 아니라 신호의존성 잡음과 적산성 잡음도 나타낼 수 있는 일반화된 관측모델을 이용하여 신호의존성 잡음이 있을 때 약한 확률적 신호를 검파하는 문제를 다루었다. 신호의존성 잡음이 있을 경우, 약한 확률적 신호를 검파하기 위한 국소최적 검파기의 검정통계량은 순가산성 잡음만 있을 때의 국소최적 검파기의 검정 통계량이 확장된 것임을 보였다. 이는 이미 발표된 적산성 잡음에서의 약한 확률적 신호 검파의 경우와 비슷한 상황이다.