Search

Analysing data from educational tests allows governments to make decisions for improving the quality of life of individuals in a society. One of the key responsibilities of statisticians is to develop models that provide decision‐makers with pertinent information about the latent process that educational tests seek to represent. Mixtures of t$$ t $$ factor analysers (MtFA) have emerged as a powerful device for model‐based clustering and classification of high‐dimensional data containing one or several groups of observations with fatter tails or anomalous outliers. This paper considers an extension of MtFA for robust clustering of censored data, referred to as the MtFAC model, by incorporating external covariates. The enhanced flexibility of including covariates in MtFAC enables cluster‐specific multivariate regression analysis of dependent variables with censored responses arising from upper and/or lower detection limits of experimental equipment. An alternating expectation conditional maximization (AECM) algorithm is developed for maximum likelihood estimation of the proposed model. Two simulation experiments are conducted to examine the effectiveness of the techniques presented. Furthermore, the proposed methodology is applied to Peruvian data from the 2007 Early Grade Reading Assessment, and the results obtained from the analysis provide new insights regarding the reading skills of Peruvian students. [ABSTRACT FROM AUTHOR]

Multivariate normal/independent (MNI) distributions contain many renowned heavy-tailed distributions such as the multivariate t , multivariate slash, multivariate contaminated normal, multivariate variance-gamma, and multivariate double exponential distributions. A frequent problem encountered in statistical analysis is the occurrence of truncated observations and non-normality such that theoretical moments are required for the estimation of the truncated multivariate normal/independent (TMNI) distributions. This paper is dedicated to deriving explicit expressions for the moments of the TMNI distributions with supports confined within a hyper-rectangle. A Monte Carlo experiment is undertaken to validate to the correctness of the proposed formulae for five selected members of the TMNI distributions. R scripts and data to reproduce the results are available in the GitHub repository. [ABSTRACT FROM AUTHOR]

Semiparametric mixed-effects models (SMM) have received increasing attention in recent years because of the greater flexibility in analysing longitudinal trajectories. However, the normality assumption of SMM may be unrealistic when outliers occur in the data. This paper presents a semiparametric extension of the multivariate t linear mixed-effects model (MtLMM), called the multivariate t semiparametric mixed model (MtSMM). To be specific, the MtSMM incorporates a parametric linear function related to the fixed covariate effects and random effects which have a joint multivariate t distribution together with an arbitrary nonparametric smooth function to capture the unexpected patterns. A computationally analytical EM-based algorithm is developed for carrying out maximum likelihood estimation of the MtSMM. Simulation studies and a real example concerning the analysis of PBCseq data are used to investigate the empirical behaviour of the proposed methodology. [ABSTRACT FROM AUTHOR]

This paper presents a robust extension of factor analysis model by assuming the multivariate normal mean–variance mixture of Birnbaum–Saunders distribution for the unobservable factors and errors. A computationally analytical EM-based algorithm is developed to find maximum likelihood estimates of the parameters. The asymptotic standard errors of parameter estimates are derived under an information-based paradigm. Numerical merits of the proposed methodology are illustrated using both simulated and real datasets. [ABSTRACT FROM AUTHOR]

Multivariate longitudinal data arisen in medical studies often exhibit complex features such as censored responses, intermittent missing values, and atypical or outlying observations. The multivariate-t linear mixed model (MtLMM) has been recognized as a powerful tool for robust modeling of multivariate longitudinal data in the presence of potential outliers or fat-tailed noises. This paper presents a generalization of MtLMM, called the MtLMM-CM, to properly adjust for censorship due to detection limits of the assay and missingness embodied within multiple outcome variables recorded at irregular occasions. An expectation conditional maximization either (ECME) algorithm is developed to compute parameter estimates using the maximum likelihood (ML) approach. The asymptotic standard errors of the ML estimators of fixed effects are obtained by inverting the empirical information matrix according to Louis' method. The techniques for the estimation of random effects and imputation of missing responses are also investigated. The proposed methodology is illustrated on two real-world examples from HIV-AIDS studies and a simulation study under a variety of scenarios. [ABSTRACT FROM AUTHOR]

The analysis of complex longitudinal data is challenging due to several inherent features: (i) more than one series of responses are repeatedly collected on each subject at irregularly occasions over a period of time; (ii) censorship due to limits of quantification of responses arises left- and/or right- censoring effects; (iii) outliers or heavy-tailed noises are possibly embodied within multiple response variables. This article formulates the multivariate- t linear mixed model with censored responses (MtLMMC), which allows the analysts to model such data in the presence of the above described features simultaneously. An efficient expectation conditional maximization either (ECME) algorithm is developed to carry out maximum likelihood estimation of model parameters. The implementation of the E-step relies on the mean and covariance matrix of truncated multivariate- t distributions. To enhance the computational efficiency, two auxiliary permutation matrices are incorporated into the procedure to determine the observed and censored parts of each subject. The proposed methodology is demonstrated via a simulation study and a real application on HIV/AIDS data. [ABSTRACT FROM AUTHOR]

Abstract: This paper presents a fully Bayesian approach to multivariate t regression models whose mean vector and scale covariance matrix are modelled jointly for analyzing longitudinal data. The scale covariance structure is factorized in terms of unconstrained autoregressive and scale innovation parameters through a modified Cholesky decomposition. A computationally flexible data augmentation sampler coupled with the Metropolis-within-Gibbs scheme is developed for computing the posterior distributions of parameters. The Bayesian predictive inference for the future response vector is also investigated. The proposed methodologies are illustrated through a real example from a sleep dose–response study. [Copyright &y& Elsevier]

Abstract: In this paper, we propose a multivariate regression model with its mean and scale covariance modeled jointly for the analysis of longitudinal data. A modified Cholesky decomposition is adopted to factorize the dependence structure in terms of unconstrained autoregressive and scale innovation parameters. We present three distinct representations of the log-likelihood function of the model and study the associated properties. A computationally efficient Fisher scoring algorithm is developed for carrying out maximum likelihood estimation. The technique for the prediction of future responses in this context is also investigated. The implementation of the proposed methodology is illustrated through two real-life examples and extensive simulation studies. [Copyright &y& Elsevier]

The t factor analysis (tFA) model is a promising tool for robust reduction of high-dimensional data in the presence of heavy-tailed noises. When determining the number of factors of the tFA model, a two-stage procedure is commonly performed in which parameter estimation is carried out for a number of candidate models, and then the best model is chosen according to certain penalized likelihood indices such as the Bayesian information criterion. However, the computational burden of such a procedure could be extremely high to achieve the optimal performance, particularly for extensively large data sets. In this paper, we develop a novel automated learning method in which parameter estimation and model selection are seamlessly integrated into a one-stage algorithm. This new scheme is called the automated tFA (AtFA) algorithm, and it is also workable when values are missing. In addition, we derive the Fisher information matrix to approximate the asymptotic covariance matrix associated with the ML estimators of tFA models. Experiments on real and simulated data sets reveal that the AtFA algorithm not only provides identical fitting results, as compared to traditional two-stage procedures, but also runs much faster, especially when values are missing. [ABSTRACT FROM AUTHOR]