Identifying predictors of violent behaviour among students using the conventional logistic and multilevel logistic models

Analysing individual-, school- and class-level observations is a good and efficient approach in epidemiologic research. Using data on violent behaviour among secondary school students we compared results from the conventional logistic modelling with multilevel logistic modelling approach using the gllamm command in Stata. We illustrated the advantage of multilevel modelling over the conventional logistic modelling through an example of data from violence experience among secondary school students. We constructed a logistic model with a random intercept on the school and class levels to account for unexplained heterogeneity between schools and classes. In the multilevel model, we estimated that, in an average school, the odds of experiencing violence are 3 (OR=2.99, 95% CI: 1.86, 4.81, p<0.0001) times higher for students who use drugs as opposed to the odds of experiencing violence for students who do not use drugs. However, the estimates in the conventional logistic model are slightly lower.

We estimated that a normally distributed random intercept for schools and classes that accounts for any unexplained heterogeneity between schools and classes has variances 0.017 and 0.035, respectively. We therefore recommend the multilevel logistic modelling when data are clustered.     

Fitting bivariate multilevel models to assess long-term changes in body mass index and cigarette smoking

Using data from the National Health interview Survey from 1997 to 2006, we present a multilevel analysis of change in body mass index (BMI) and number of cigarettes smoked per day in the USA. Smoking and obesity are the leading causes of preventable mortality and morbidity in the USA and most parts of the developed world. A two-stage bivariate model of changes in obesity and number of cigarette smoked per day is proposed. At the within subject stage, an individual's BMI status and the number of cigarette smoked per day are jointly modeled as a function of an individual growth trajectory plus a random error. At the between-subject stage, the parameters of the individual growth trajectories are allowed to vary as a function of differences between subjects with respect to demographic and behavioral characteristics and with respect to the four regions of the USA (Northeast, West, South and North central). Our two-stage modeling techniques are more informative than standard regression because they characterize both group-level (nomothetic) and individual-level (idiographic) effects, yielding a more complete understanding of the phenomena under study.     

Bivariate binary models of efficacy and toxicity in dose-ranging trials

Methods for the simultaneous analysis of the relationships of binary variables for efficacy and toxicity to dosage of an experimental drug are developed. Properties of two models of ‘within-dose’ dependence of efficacy and toxicity in parallel designs - one a bivariate analogue of the familiar univariate logistic model, and the other an adaptation of a general model developed by D.R. Cox– are explored. The cell probabilities predicted by these models are often quite similar to those predicted by a model of independence of efficacy and toxicity, but large discrepancies can occur when there is approximate equality of the median effective and median toxic doses. Asymptotic variances of estimates of parameters involved in assessing correlation are large when there is little or no dependence in the data, but parameters can be estimated with good precision in at least some cases of moderate to strong dependence between efficacy and toxicity.     

Bivariate beta regression models: joint modeling of the mean,dispersion and association parameters

In this paper a bivariate beta regression model with joint modeling of the mean and dispersion parameters is proposed, defining the bivariate beta distribution from Farlie–Gumbel–Morgenstern (FGM) copulas. This model, that can be generalized using other copulas, is a good alternative to analyze non-independent pairs of proportions and can be fitted applying standard Markov chain Monte Carlo methods. Results of two applications of the proposed model in the analysis of structural and real data set are included.     

A Comparison of Frailty and Other Models for Bivariate Survival Data

Multivariate survival data arise when eachstudy subject may experience multiple events or when study subjectsare clustered into groups. Statistical analyses of such dataneed to account for the intra-cluster dependence through appropriatemodeling. Frailty models are the most popular for such failuretime data. However, there are other approaches which model thedependence structure directly. In this article, we compare thefrailty models for bivariate data with the models based on bivariateexponential and Weibull distributions. Bayesian methods providea convenient paradigm for comparing the two sets of models weconsider. Our techniques are illustrated using two examples.One simulated example demonstrates model choice methods developedin this paper and the other example, based on a practical dataset of onset of blindness among patients with diabetic Retinopathy,considers Bayesian inference using different models.     

A recursive algorithm on estimating the parameters in multilevel models subject to the measurement errors on the covariates

Multilevel models are popular models for analysing data inheriting a hierarchical structure. They are used in diverse fields including social, medical, economical and biological sciences. These models encounter some problems in estimating the parameters, if there are measurement errors in either explanatory or response variables. A common approach to tackle this obstacle is to consider the pseudo variables and follow some simulation methods to estimate the parameters. We propose a new algorithm constituting the iterative and simulation extrapolation steps in turn. To evaluate the proposed algorithm, various simulation studies are also conducted. Moreover, we investigate the implementation of our method on a real data set concerning the cost and expenditure of the households in Tehran city in the year 2007.     

Estimation of association parameters in copula models for bivariate left-truncated and right-censored data

We investigate the problem of estimating the association between two related survival variables when they follow a copula model and bivariate left-truncated and right-censored data are available. By expressing truncation probability as the functional of marginal survival functions, we propose a two-stage estimation procedure for estimating the parameters of Archimedean copulas. The asymptotic properties of the proposed estimators are established. Simulation studies are conducted to investigate the finite sample properties of the proposed estimators. The proposed method is applied to a bivariate RNA data.     

Estimation of the Association for Bivariate Interval-censored Failure Time Data

Abstract.  Multivariate failure time data frequently occur in medical studies and the dependence or association among survival variables is often of interest ( Biometrics , 51 , 1995, 1384; Stat. Med. , 18 , 1999, 3101; Biometrika , 87 , 2000, 879; J. Roy. Statist. Soc. Ser. B , 65 , 2003, 257). We study the problem of estimating the association between two related survival variables when they follow a copula model and only bivariate interval-censored failure time data are available. For the problem, a two-stage estimation procedure is proposed and the asymptotic properties of the proposed estimator are established. Simulation studies are conducted to assess the finite sample properties of the presented estimate and the results suggest that the method works well for practical situations. An example from an acquired immunodeficiency syndrome clinical trial is discussed.     

Embedding latent class regression and latent class distal outcome models into cluster-weighted latent class analysis: a detailed simulation experiment

Usually in latent class (LC) analysis, external predictors are taken to be cluster conditional probability predictors (LC models with external predictors), and/or score conditional probability predictors (LC regression models). In such cases, their distribution is not of interest. Class-specific distribution is of interest in the distal outcome model, when the distribution of the external variables is assumed to depend on LC membership. In this paper, we consider a more general formulation, that embeds both the LC regression and the distal outcome models, as is typically done in cluster-weighted modelling. This allows us to investigate (1) whether the distribution of the external variables differs across classes, (2) whether there are significant direct effects of the external variables on the indicators, by modelling jointly the relationship between the external and the latent variables. We show the advantages of the proposed modelling approach through a set of artificial examples, an extensive simulation study and an empirical application about psychological contracts among employees and employers in Belgium and the Netherlands.     

Generalized linear mixed models for ordered responses in complex multilevel structures: effects beneath the school or college in education

Summary.  The complexities of educational processes and structure and the need for disentangling effects beneath the level of the school or college are discussed. Ordinal response multilevel crossed random-effects models for educational grades are introduced. Weighted random effects for teacher contributions are then added. Estimation methodology is reviewed. Specially written macros for quasi-likelihood with second-order terms are described. The application discusses General Certificate of Education at advanced level grades cross-classified by student and teaching group within a number of institutions. The methods handle teacher effects where several teachers contribute to provision and where each teacher deals with several groups. Some methodological lessons are drawn for sparse data and the use of extra-multinomial variation. Developments of the analysis yield conclusions about the sources of variation in educational progress, and particularly the effect of teachers.     

A comparison of bootstrap approaches for estimating uncertainty of parameters in linear mixed‐effects models

A version of the nonparametric bootstrap, which resamples the entire subjects from original data, called the case bootstrap, has been increasingly used for estimating uncertainty of parameters in mixed‐effects models. It is usually applied to obtain more robust estimates of the parameters and more realistic confidence intervals (CIs). Alternative bootstrap methods, such as residual bootstrap and parametric bootstrap that resample both random effects and residuals, have been proposed to better take into account the hierarchical structure of multi‐level and longitudinal data. However, few studies have been performed to compare these different approaches. In this study, we used simulation to evaluate bootstrap methods proposed for linear mixed‐effect models. We also compared the results obtained by maximum likelihood (ML) and restricted maximum likelihood (REML). Our simulation studies evidenced the good performance of the case bootstrap as well as the bootstraps of both random effects and residuals. On the other hand, the bootstrap methods that resample only the residuals and the bootstraps combining case and residuals performed poorly. REML and ML provided similar bootstrap estimates of uncertainty, but there was slightly more bias and poorer coverage rate for variance parameters with ML in the sparse design. We applied the proposed methods to a real dataset from a study investigating the natural evolution of Parkinson's disease and were able to confirm that the methods provide plausible estimates of uncertainty. Given that most real‐life datasets tend to exhibit heterogeneity in sampling schedules, the residual bootstraps would be expected to perform better than the case bootstrap. Copyright © 2013 John Wiley & Sons, Ltd.     

The use of simple reparameterizations to improve the efficiency of Markov chain Monte Carlo estimation for multilevel models with applications to discrete time survival models

Summary.  We consider the application of Markov chain Monte Carlo (MCMC) estimation methods to random-effects models and in particular the family of discrete time survival models. Survival models can be used in many situations in the medical and social sciences and we illustrate their use through two examples that differ in terms of both substantive area and data structure. A multilevel discrete time survival analysis involves expanding the data set so that the model can be cast as a standard multilevel binary response model. For such models it has been shown that MCMC methods have advantages in terms of reducing estimate bias. However, the data expansion results in very large data sets for which MCMC estimation is often slow and can produce chains that exhibit poor mixing. Any way of improving the mixing will result in both speeding up the methods and more confidence in the estimates that are produced. The MCMC methodological literature is full of alternative algorithms designed to improve mixing of chains and we describe three reparameterization techniques that are easy to implement in available software. We consider two examples of multilevel survival analysis: incidence of mastitis in dairy cattle and contraceptive use dynamics in Indonesia. For each application we show where the reparameterization techniques can be used and assess their performance.     

A Unified Approach for Construction of Probability Models for Bivariate Linear and Directional Data

This article considers a unified approach based on the mixture method to construct linear bivariate models and those on the cylinder and torus involving the exponential and cardioid distributions with the truncated exponential distribution as the mixing distribution. Parameter estimation of the bivariate model on the torus is considered for the data set of phase angles of circadian-related genes in heart and liver tissues.     

A comparison of population average and random-effect models for the analysis of longitudinal count data with base-line information

The generalized estimating equation (GEE) approach to the analysis of longitudinal data has many attractive robustness properties and can provide a 'population average' characterization of interest, for example, to clinicians who have to treat patients on the basis of their observed characteristics. However, these methods have limitations which restrict their usefulness in both the social and the medical sciences. This conclusion is based on the premise that the main motivations for longitudinal analysis are insight into microlevel dynamics and improved control for omitted or unmeasured variables. We claim that to address these issues a properly formulated random-effects model is required. In addition to a theoretical assessment of some of the issues, we illustrate this by reanalysing data on polyp counts. In this example, the covariates include a base-line outcome, and the effectiveness of the treatment seems to vary by base-line. We compare the random-effects approach with the GEE approach and conclude that the GEE approach is inappropriate for assessing the treatment effects for these data.     

Latent class models for longitudinal studies of the elderly with data missing at random

The elderly population in the USA is expected to double in size by the year 2025, making longitudinal health studies of this population of increasing importance. The degree of loss to follow-up in studies of the elderly, which is often because elderly people cannot remain in the study, enter a nursing home or die, make longitudinal studies of this population problematic. We propose a latent class model for analysing multiple longitudinal binary health outcomes with multiple-cause non-response when the data are missing at random and a non-likelihood-based analysis is performed. We extend the estimating equations approach of Robins and co-workers to latent class models by reweighting the multiple binary longitudinal outcomes by the inverse probability of being observed. This results in consistent parameter estimates when the probability of non-response depends on observed outcomes and covariates (missing at random) assuming that the model for non-response is correctly specified. We extend the non-response model so that institutionalization, death and missingness due to failure to locate, refusal or incomplete data each have their own set of non-response probabilities. Robust variance estimates are derived which account for the use of a possibly misspecified covariance matrix, estimation of missing data weights and estimation of latent class measurement parameters. This approach is then applied to a study of lower body function among a subsample of the elderly participating in the 6-year Longitudinal Study of Aging.     

A multilevel model with autoregressive components for the analysis of tribal art prices

In this paper, we introduce a multilevel model specification with time-series components for the analysis of prices of artworks sold at auctions. Since auction data do not constitute a panel or a time series but are composed of repeated cross-sections, they require a specification with items at the first level nested in time-points. Our approach combines the flexibility of mixed effect models together with the predicting performance of time series as it allows to model the time dynamics directly. Model estimation is obtained by means of maximum likelihood through the expectation–maximization algorithm. The model is motivated by the analysis of the first database ethnic artworks sold in the most important auctions worldwide. The results show that the proposed specification improves considerably over classical proposals both in terms of fit and prediction.     

A comparison of univariate and multivariate multilevel models for repeated measures of use of antenatal care in Uttar Pradesh

Summary.  We compare two different multilevel modelling approaches to the analysis of repeated measures data to assess the effect of mother level characteristics on women's use of prenatal care services in Uttar Pradesh, India. We apply univariate multilevel models to our data and find that the model assumptions are severely violated and the parameter estimates are not stable, particularly for the mother level random effect. To overcome this we apply a multivariate multilevel model. The correlation structure shows that, once the decision has been made regarding use of antenatal care by the mother for her first observed birth in the data, she does not tend to change this decision for higher order births.     

Statistical analysis of performance indicators in UK higher education

Summary.  Attempts to measure the quality with which institutions such as hospitals and universities carry out their public mandates have gained in frequency and sophistication over the last decade. We examine methods for creating performance indicators in multilevel or hierarchical settings (e.g. students nested within universities) based on a dichotomous outcome variable (e.g. drop-out from the higher education system). The profiling methods that we study involve the indirect measurement of quality, by comparing institutional outputs after adjusting for inputs, rather than directly attempting to measure the quality of the processes unfolding inside the institutions. In the context of an extended case-study of the creation of performance indicators for universities in the UK higher education system, we demonstrate the large sample functional equivalence between a method based on indirect standardization and an approach based on fixed effects hierarchical modelling, offer simulation results on the performance of the standardization method in null and non-null settings, examine the sensitivity of this method to the inadvertent omission of relevant adjustment variables, explore random-effects reformulations and characterize settings in which they are preferable to fixed effects hierarchical modelling in this type of quality assessment and discuss extensions to longitudinal quality modelling and the overall pros and cons of institutional profiling. Our results are couched in the language of higher education but apply with equal force to other settings with dichotomous response variables, such as the examination of observed and expected rates of mortality (or other adverse outcomes) in investigations of the quality of health care or the study of retention rates in the workplace.     

Bayesian analysis of joint mean and covariance models for longitudinal data

Efficient estimation of the regression coefficients in longitudinal data analysis requires a correct specification of the covariance structure. If misspecification occurs, it may lead to inefficient or biased estimators of parameters in the mean. One of the most commonly used methods for handling the covariance matrix is based on simultaneous modeling of the Cholesky decomposition. Therefore, in this paper, we reparameterize covariance structures in longitudinal data analysis through the modified Cholesky decomposition of itself. Based on this modified Cholesky decomposition, the within-subject covariance matrix is decomposed into a unit lower triangular matrix involving moving average coefficients and a diagonal matrix involving innovation variances, which are modeled as linear functions of covariates. Then, we propose a fully Bayesian inference for joint mean and covariance models based on this decomposition. A computational efficient Markov chain Monte Carlo method which combines the Gibbs sampler and Metropolis–Hastings algorithm is implemented to simultaneously obtain the Bayesian estimates of unknown parameters, as well as their standard deviation estimates. Finally, several simulation studies and a real example are presented to illustrate the proposed methodology.     

Generalized linear models for the analysis of quality-improvement experiments

Generalized linear models provide a useful tool for analyzing data from quality-improvement experiments. We discuss why analysis must be done for all the data, not just for summarizing quantities, and show by examples how residuals can be used for model checking. A restricted-maximum-likelihood-type adjustment for the dispersion analysis is developed.