A marginalized multilevel model for bivariate longitudinal binary data

2019-06-01
This study considers analysis of bivariate longitudinal binary data. We propose a model based on marginalized multilevel model framework. The proposed model consists of two levels such that the first level associates the marginal mean of responses with covariates through a logistic regression model and the second level includes subject/time specific random intercepts within a probit regression model. The covariance matrix of multiple correlated time-specific random intercepts for each subject is assumed to represent the within-subject association. The subject-specific random effects covariance matrix is further decomposed into its dependence and variance components through modified Cholesky decomposition method and then the unconstrained version of resulting parameters are modelled in terms of covariates with low-dimensional regression parameters. This provides better explanations related to dependence and variance parameters and a reduction in the number of parameters to be estimated in random effects covariance matrix to avoid possible identifiability problems. Marginal correlations between responses of subjects and within the responses of a subject are derived through a Taylor series-based approximation. Data cloning computational algorithm is used to compute the maximum likelihood estimates and standard errors of the parameters in the proposed model. The validity of the proposed model is assessed through a Monte Carlo simulation study, and results are observed to be at acceptable level. Lastly, the proposed model is illustrated through Mother's Stress and Children's Morbidity study data, where both population-averaged and subject-specific interpretations are drawn through Emprical Bayes estimation of random effects.
STATISTICAL PAPERS

Suggestions

A marginalized multilevel model for bivariate longitudinal binary data
İnan, Gül; İlk Dağ, Özlem; Department of Statistics (2014)
This thesis study considers analysis of bivariate longitudinal binary data. We propose a model based on marginalized multilevel model framework. The proposed model consists of two levels such that the first level associates the marginal mean of responses with covariates through a logistic regression model and the second level includes subject/time specific random intercepts within a probit regression model. The covariance matrix of multiple correlated time-specific random intercepts for each subject is assu...
Robust estimation and hypothesis testing under short-tailedness and inliers
Akkaya, Ayşen (Springer Science and Business Media LLC, 2005-06-01)
Estimation and hypothesis testing based on normal samples censored in the middle are developed and shown to be remarkably efficient and robust to symmetric short-tailed distributions and to inliers in a sample. This negates the perception that sample mean and variance are the best robust estimators in such situations (Tiku, 1980; Dunnett, 1982).
Estimation and hypothesis testing in multivariate linear regression models under non normality
İslam, Muhammed Qamarul (Informa UK Limited, 2017-01-01)
This paper discusses the problem of statistical inference in multivariate linear regression models when the errors involved are non normally distributed. We consider multivariate t-distribution, a fat-tailed distribution, for the errors as alternative to normal distribution. Such non normality is commonly observed in working with many data sets, e.g., financial data that are usually having excess kurtosis. This distribution has a number of applications in many other areas of research as well. We use modifie...
The effect of Phase I sample size on the run length performance of control charts for autocorrelated data
Köksal, Gülser; Ula, Taylan Ali; Testik, Murat Caner (Informa UK Limited, 2008-01-01)
Traditional control charts assume independence of observations obtained from the monitored process. However, if the observations are autocorrelated, these charts often do not perform as intended by the design requirements. Recently, several control charts have been proposed to deal with autocorrelated observations. The residual chart, modified Shewhart chart, EWMAST chart, and ARMA chart are such charts widely used for monitoring the occurrence of assignable causes in a process when the process exhibits inh...
Representation of Multiplicative Seasonal Vector Autoregressive Moving Average Models
Yozgatlıgil, Ceylan (Informa UK Limited, 2009-11-01)
Time series often contain observations of several variables and multivariate time series models are used to represent the relationship between these variables. There are many studies on vector autoregressive moving average (VARMA) models, but the representation of multiplicative seasonal VARMA models has not been seriously studied. In a multiplicative vector model, such as a seasonal VARMA model, the representation is not unique because of the noncommutative property of matrix multiplication. In this articl...
Citation Formats
G. Inan and Ö. İlk Dağ, “A marginalized multilevel model for bivariate longitudinal binary data,” STATISTICAL PAPERS, pp. 251–278, 2019, Accessed: 00, 2020. [Online]. Available: https://hdl.handle.net/11511/37068.