• Communications for Statistical Applications and Methods
• /
• v.27 no.3
• /
• pp.301-311
• /
• 2020

Schwarz's Bayesian information criterion (BIC) is one of the most popular criteria for model selection, that was derived under the assumption of independent and identical distribution. For correlated data in longitudinal studies, Jones (Statistics in Medicine, 30, 3050-3056, 2011) modified the BIC to select the best linear mixed effects model based on the effective sample size where the number of parameters in covariance structure was not considered. In this paper, we propose an extended Jones' modified BIC by considering covariance parameters. We conducted simulation studies under a variety of parameter configurations for linear mixed effects models. Our simulation study indicates that our proposed BIC performs better in model selection than Schwarz's BIC and Jones' modified BIC do in most scenarios. We also illustrate an example of smoking data using a longitudinal cohort of cancer patients.

• Journal of Periodontal and Implant Science
• /
• v.45 no.1
• /
• pp.2-7
• /
• 2015

A fundamental problem in analyzing complex multilevel-structured periodontal data is the violation of independency among the observations, which is an assumption in traditional statistical models (e.g., analysis of variance and ordinary least squares regression). In many cases, aggregation (i.e., mean or sum scores) has been employed to overcome this problem. However, the aggregation approach still exhibits certain limitations, such as a loss of power and detailed information, no cross-level relationship analysis, and the potential for creating an ecological fallacy. In order to handle multilevel-structured data appropriately, mixed effects models have been introduced and employed in dental research using periodontal data. The use of mixed effects models might account for the potential bias due to the violation of the independency assumption as well as provide accurate estimates.

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

Linear mixed models are commonly used in the clinical pharmaceutical studies to analyze repeated measures such as the crossover study data of bioequivalence studies. In these models, random effects describe the correlation between repeated outcomes and variance-covariance matrix explain within-subject variabilities. Bioequivalence analysis verifies whether a 90% confidence interval for geometric mean ratio of Cmax and AUC between reference drug and test drug is included in the bioequivalence margin [0.8, 1.25] performed using linear mixed models with period, sequence and treatment effects as fixed and sequence nested subject effects as random. A Levofloxacin study is referred to for an example of real data analysis.

• Communications for Statistical Applications and Methods
• /
• v.20 no.3
• /
• pp.235-240
• /
• 2013

Generalized linear mixed models(GLMMs) are frequently used for the analysis of longitudinal categorical data when the subject-specific effects is of interest. In GLMMs, the structure of the random effects covariance matrix is important for the estimation of fixed effects and to explain subject and time variations. The estimation of the matrix is not simple because of the high dimension and the positive definiteness; subsequently, we practically use the simple structure of the covariance matrix such as AR(1). However, this strong assumption can result in biased estimates of the fixed effects. In this paper, we introduce Bayesian modeling approaches for the random effects covariance matrix using a modified Cholesky decomposition. The modified Cholesky decomposition approach has been used to explain a heterogenous random effects covariance matrix and the subsequent estimated covariance matrix will be positive definite. We analyze metabolic syndrome data from a Korean Genomic Epidemiology Study using these methods.

• The Korean Journal of Applied Statistics
• /
• v.23 no.6
• /
• pp.1169-1178
• /
• 2010

This article consider a performance model selection based on AIC under unbalanced deign in linear mixed effect models. Vaida and Balanchard (2005) proposed conditional AIC for model selection in linear mixed effect models when the prediction of random effects is of primary interest. Theoretical properties of cAIC and related criteria have been investigated by Liang et al. (2008) and Greven and Kneib (2010). However, all of the simulation studies were performed under a balanced design. Even though functional form of AIC remain same even under the unbalanced deign, it is worthwhile to investigate performance of AIC based model selection criteria under the unbalanced design. The simulation study in this article shows how unbalancedness affects model selection in linear mixed effect models.

• The Korean Journal of Applied Statistics
• /
• v.25 no.6
• /
• pp.1037-1047
• /
• 2012

The generalized linear mixed model(GLMM) is widely used in fitting categorical responses of clustered data. In the numerical approximation of likelihood function the normality is assumed for the random effects distribution; subsequently, the commercial statistical packages also routinely fit GLMM under this normality assumption. We may also encounter departures from the distributional assumption on the response variable. It would be interesting to investigate the impact on the estimates of parameters under misspecification of distributions; however, there has been limited researche on these topics. We study the sensitivity or robustness of the maximum likelihood estimators(MLEs) of GLMM for counts data when the true underlying distribution is normal, gamma, exponential, and a mixture of two normal distributions. We also consider the effects on the MLEs when we fit Poisson-normal GLMM whereas the outcomes are generated from the negative binomial distribution with overdispersion. Through a small scale Monte Carlo study we check the empirical coverage probabilities of parameters and biases of MLEs of GLMM.

• The Korean Journal of Applied Statistics
• /
• v.31 no.3
• /
• pp.353-366
• /
• 2018

This study uses mixed-effects models to predict thigh soft tissue artefact (STA), relative movement of soft tissue such as skin to femur occurring during hip joint motions. The random effects in the model were defined as STA and the fixed effects in the model were considered as skeletal motion. Five male subjects without musculoskeletal disease were selected to perform various hip joint rotational motions. Linear mixed-effects models were applied to markers' position vectors acquired from non-invasive method, photogrammetry. Predicted random effects showed similar patterns of STA among subjects. Large magnitudes of STA appeared on the points near the hip joint regardless of sides; however, small values appeared on the distal anterior.

• Communications for Statistical Applications and Methods
• /
• v.16 no.6
• /
• pp.971-978
• /
• 2009

Poisson generalized linear mixed models(GLMMs) have been widely used for the analysis of clustered or correlated count data. For the inference marginal likelihood, which is obtained by integrating out random effects is often used. It gives maximum likelihood(ML) estimator, but the integration is usually intractable. In this paper, we propose how to obtain the ML estimator via Laplace approximation based on hierarchical-likelihood (h-likelihood) approach under the Poisson GLMMs. In particular, the h-likelihood avoids the integration itself and gives a statistically efficient procedure for various random-effect models including GLMMs. The proposed method is illustrated using two practical examples and simulation studies.

• Communications for Statistical Applications and Methods
• /
• v.23 no.6
• /
• pp.575-585
• /
• 2016

Longitudinal studies repeatedly measure outcomes over time. Therefore, repeated measurements are serially correlated from same subject (within-subject variation) and there is also variation between subjects (between-subject variation). The serial correlation and the between-subject variation must be taken into account to make proper inference on covariate effects (Diggle et al., 2002). However, estimation of the covariance matrix is challenging because of many parameters and positive definiteness of the matrix. To overcome these limitations, we propose autoregressive moving average Cholesky decomposition (ARMACD) for the linear mixed models. The ARMACD allows a class of flexible, nonstationary, and heteroscedastic models that exploits the structure allowed by combining the AR and MA modeling of the random effects covariance matrix. We analyze a real dataset to illustrate our proposed methods.

• Journal of the Korean Society of Clothing and Textiles
• /
• v.44 no.1
• /
• pp.84-95
• /
• 2020

Wearing ease is a critical factor when designing special uniforms such as flight pilot's garment and should reflect occupational properties for better performance. This study measured skin surface on 31 areas in seven postures that refer to the pilot's occupational postures as well as made six prediction models including linear mixed model (LMM) for each body part to find the best fit model. Skin surface measured from 3D body scanned images of 11 male pilot participants. There were significantly positive and negative changes in various areas from standing posture (P1) to dynamic postures (P2-P7). Six models were designed in various compositions using stature and chest circumference as fixed effects and subject and posture as random effects. The best models were linear mixed models with one fixed effect (chest circumference or stature, varies with body parts) and two random effects (subject and posture). The results of this study provide reference data to set wearing ease for pilot's garment and suggests a new methodology in this research area, but verifying the effect of diverse independent variables is left for future studies.