A zero-augmented gamma mixed model for longitudinal data with many zeros

In many occupational safety interventions, the objective is to reduce the injury incidence as well as the mean claims cost once injury has occurred. The claims cost data within a period typically contain a large proportion of zero observations (no claim). The distribution thus comprises a point mass at 0 mixed with a non-degenerate parametric component. Essentially, the likelihood function can be factorized into two orthogonal components. These two components relate respectively to the effect of covariates on the incidence of claims and the magnitude of claims, given that claims are made. Furthermore, the longitudinal nature of the intervention inherently imposes some correlation among the observations. This paper introduces a zero-augmented gamma random effects model for analysing longitudinal data with many zeros. Adopting the generalized linear mixed model (GLMM) approach reduces the original problem to the fitting of two independent GLMMs. The method is applied to evaluate the effectiveness of a workplace risk assessment teams program, trialled within the cleaning services of a Western Australian public hospital.

Effect of metabolic inhibitors on cyclosporine pharmacokinetics using a population approach

Drugs known to inhibit the metabolism of cyclosporine are administered concomitantly to those who undergo cardiothoracic transplantation. The aim of this study was to examine in quantitative terms the relationship between cyclosporine oral dose rate and the trough concentration (Css(trough)) at steady state in patients who undergo cardiothoracic transplantation and are administered cyclosporine alone or in combination with drugs known to inhibit its metabolism. Dose and whole blood cyclosporine Css(tough) observations measured using the enzyme-multiplied immunoassay technique (EMIT) (396 observations) or the TDx assay (435 observations) were collected as part of routine blood concentration monitoring from 182 patients who underwent cardiothoracic transplantation. Data were analyzed using a linear mixed-effects modeling approach to examine the effect of metabolic inhibitors on dose-rate-Css(trough) ratio. The mean (and 95% confidence interval) dose-rate-Css(trough) ratio for cyclosporine generated from concentrations measured using EMIT was 94 (82.5-105.5) Lh(-1) for patients administered cyclosporine alone, 66.7 (58.1-75.3) Lh(-1) for patients administered concomitant diltiazem, 47.9 (15.4 -80.4) Lh(-1) for patients administered concomitant itraconazole, 21.7 (14.8-28.5) Lh(-1) for patients administered concomitant ketoconazole, and 14.9 (11.8-18.1) Lh(-1) for patients concomitantly administered diltiazem and ketoconazole. For patients administered concomitant cyclosporine, ketoconazole, and diltiazem, the dosage of cyclosporine, if it is administered alone, should be 20% to achieve the same blood concentrations. This will allow safer drug concentration targeting of cyclosporine after cardiothoracic transplantation.

Modelling the Distribution of Ischaemic Stroke-Specific Survival Time Using an EM-based Mixture Approach with Random Effects Adjustment

A two-component survival mixture model is proposed to analyse a set of ischaemic stroke-specific mortality data. The survival experience of stroke patients after index stroke may be described by a subpopulation of patients in the acute condition and another subpopulation of patients in the chronic phase. To adjust for the inherent correlation of observations due to random hospital effects, a mixture model of two survival functions with random effects is formulated. Assuming a Weibull hazard in both components, an EM algorithm is developed for the estimation of fixed effect parameters and variance components. A simulation study is conducted to assess the performance of the two-component survival mixture model estimators. Simulation results confirm the applicability of the proposed model in a small sample setting. Copyright (C) 2004 John Wiley Sons, Ltd.

Long-term survivor mixture model with random effects: application to a multi-centre clinical trial of carcinoma

A mixture model incorporating long-term survivors has been adopted in the field of biostatistics where some individuals may never experience the failure event under study. The surviving fractions may be considered as cured. In most applications, the survival times are assumed to be independent. However, when the survival data are obtained from a multi-centre clinical trial, it is conceived that the environ mental conditions and facilities shared within clinic affects the proportion cured as well as the failure risk for the uncured individuals. It necessitates a long-term survivor mixture model with random effects. In this paper, the long-term survivor mixture model is extended for the analysis of multivariate failure time data using the generalized linear mixed model (GLMM) approach. The proposed model is applied to analyse a numerical data set from a multi-centre clinical trial of carcinoma as an illustration. Some simulation experiments are performed to assess the applicability of the model based on the average biases of the estimates formed. Copyright (C) 2001 John Wiley & Sons, Ltd.

Finite mixture regression model with random effects: application to neonatal hospital length of stay

A two-component mixture regression model that allows simultaneously for heterogeneity and dependency among observations is proposed. By specifying random effects explicitly in the linear predictor of the mixture probability and the mixture components, parameter estimation is achieved by maximising the corresponding best linear unbiased prediction type log-likelihood. Approximate residual maximum likelihood estimates are obtained via an EM algorithm in the manner of generalised linear mixed model (GLMM). The method can be extended to a g-component mixture regression model with the component density from the exponential family, leading to the development of the class of finite mixture GLMM. For illustration, the method is applied to analyse neonatal length of stay (LOS). It is shown that identification of pertinent factors that influence hospital LOS can provide important information for health care planning and resource allocation. (C) 2002 Elsevier Science B.V. All rights reserved.

Modelling inpatient length of stay by a hierarchical mixture regression via the EM algorithm

The modelling of inpatient length of stay (LOS) has important implications in health care studies. Finite mixture distributions are usually used to model the heterogeneous LOS distribution, due to a certain proportion of patients sustaining-a longer stay. However, the morbidity data are collected from hospitals, observations clustered within the same hospital are often correlated. The generalized linear mixed model approach is adopted to accommodate the inherent correlation via unobservable random effects. An EM algorithm is developed to obtain residual maximum quasi-likelihood estimation. The proposed hierarchical mixture regression approach enables the identification and assessment of factors influencing the long-stay proportion and the LOS for the long-stay patient subgroup. A neonatal LOS data set is used for illustration, (C) 2003 Elsevier Science Ltd. All rights reserved.

The rellationship between specific features of aphasia-friendly written material and comprehension of written material for people with aphasia: An exploratory study

Background: Written material is often inaccessible fro people with aphasia. The format of written material needs to be adapted to enable people with aphasia to read with understanding. Aims: This study aimed to further explore some issues raised in Rose, Worrall, and MacKenna (2003) concerning the effects of aphasia-friendly formats on the reading comprehension of people with aphasia. It was hypothesised that people with aphasia would comprehend significantly more paragraphs that were formatted in an aphasia-friendly manner than control paragraphs. This study also aimed to investigate if each single aspect of aphasia-friendly formatting (i.e., simplified vocabulary and syntax, large print, increased white spacem and pictures) used in isolation would result in increased comprehension compared to control paragraphs. Other aims were to compare the effect of aphasia-friendly fromatting with the effects of each single adaptation, and to investigate if the effects of aphasia-friendly formates were related to aphasia severity. Methods & Procedures: Participants with mild to moderately severe aphasia (N = 9) read a battery of 90 paragraphs and selected the best word of phrase from a choice of four to complete each paragraph. A linear mixed model (p < .05) was used to analyse the differences in reading comprehension with each paragraph fromat across three reading grade levels. Outcomes & Results: People with aphasia comprehended significantly more aphasia-friendly paragraphs than control paragraphs. They also comprehended significantly more paragraphs with each of the following single adaptations: simplified vocabulary and syntax, large ptint, and increased white spaces. Although people with aphasia tended to comprehend more paragraphs with pictures added than control paragraphs, this difference was not significant. No significant correlation between aphasia severity and the effect of aphasia-friendly formatting was found. Conclusion: This study supports the idea that aphasia-friendly formats increase the reading comprehension of people with aphasia. It suggests that adding pictures, particularly Clip Art pictures, may not significantly improve the reading the reading comprehension of people with aphasia. These findings have implications for all written communication with people with aphasia, both in the clinical setting and in the wider community. Applying these findings may enable people with aphasia to have equal access to written information and to participate in society.

Population pharmacokinetics of itraconazole and its active metabolite hydroxy-itraconazole in paediatric cystic fibrosis and bone marrow transplant patients

Objective: The objective of the study was to characterise the population pharmacokinetic properties of itraconazole and its active metabolite hydroxyitraconazole in a representative paediatric population of cystic fibrosis and bone marrow transplant (BMT) patients and to identify patient characteristics influencing the pharmacokinetics of itraconazole. The ultimate goals were to determine the relative bioavailability between the two oral formulations (capsules vs oral solution) and to optimise dosing regimens in these patients. Methods: All paediatric patients with cystic fibrosis or patients undergoing BMT at The Royal Children's Hospital, Brisbane, QLD, Australia, who were prescribed oral itraconazole for the treatment of allergic bronchopulmonary aspergillosis (cystic fibrosis patients) or for prophylaxis of any fungal infection (BMT patients) were eligible for the study. Blood samples were taken from the recruited patients as per an empirical sampling design either during hospitalisation or during outpatient clinic visits. ltraconazole and hydroxy-itraconazole plasma concentrations were determined by a validated high-performance liquid chromatography assay with fluorometric detection. A nonlinear mixed-effect modelling approach using the NONMEM software to simultaneously describe the pharmacokinetics of itraconazole and its metabolite. Results: A one-compartment model with first-order absorption described the itraconazole data, and the metabolism of the parent drug to hydroxy-itraconazole was described by a first-order rate constant. The metabolite data also showed one-compartment characteristics with linear elimination. For itraconazole the apparent clearance (CLitraconazole) was 35.5 L/hour, the apparent volume of distribution (V-d(itraconazole)) was 672L, the absorption rate constant for the capsule formulation was 0.0901 h(-1) and for the oral solution formulation was 0.96 h-1. The lag time was estimated to be 19.1 minutes and the relative bioavailability between capsules and oral solution (F-rel) was 0.55. For the metabolite, volume of distribution, V-m/(F (.) f(m)), and clearance, CL/(F (.) fm), were 10.6L and 5.28 L/h, respectively. The influence of total bodyweight was significant, added as a covariate on CLitraconazoie/F and V-d(itraconazole)/F (standardised to a 70kg person) using allometric three-quarter power scaling on CLitraconazole/F, which therefore reflected adult values. The unexplained between-subject variability (coefficient of variation %) was 68.7%, 75.8%, 73.4% and 61.1% for CLitraconazoie/F, Vd(itraconazole)/F, CLm/(F (.) fm) and F-rel, respectively. The correlation between random effects of CLitraconazole and Vd((itraconazole)) was 0.69. Conclusion: The developed population pharmacokinetic model adequately described the pharmacokinetics of itraconazole and its active metabolite, hydroxy-itraconazole, in paediatric patients with either cystic fibrosis or undergoing BMT. More appropriate dosing schedules have been developed for the oral solution and the capsules to secure a minimum therapeutic trough plasma concentration of 0.5 mg/L for these patients.

A Bayesian hierarchical model for categorical longitudinal data from a social survey of immigrants

The paper investigates a Bayesian hierarchical model for the analysis of categorical longitudinal data from a large social survey of immigrants to Australia. Data for each subject are observed on three separate occasions, or waves, of the survey. One of the features of the data set is that observations for some variables are missing for at least one wave. A model for the employment status of immigrants is developed by introducing, at the first stage of a hierarchical model, a multinomial model for the response and then subsequent terms are introduced to explain wave and subject effects. To estimate the model, we use the Gibbs sampler, which allows missing data for both the response and the explanatory variables to be imputed at each iteration of the algorithm, given some appropriate prior distributions. After accounting for significant covariate effects in the model, results show that the relative probability of remaining unemployed diminished with time following arrival in Australia.

Methods for Categorical Longitudinal Survey Data: Understanding Employment Status of Australian Women

Many variables that are of interest in social science research are nominal variables with two or more categories, such as employment status, occupation, political preference, or self-reported health status. With longitudinal survey data it is possible to analyse the transitions of individuals between different employment states or occupations (for example). In the statistical literature, models for analysing categorical dependent variables with repeated observations belong to the family of models known as generalized linear mixed models (GLMMs). The specific GLMM for a dependent variable with three or more categories is the multinomial logit random effects model. For these models, the marginal distribution of the response does not have a closed form solution and hence numerical integration must be used to obtain maximum likelihood estimates for the model parameters. Techniques for implementing the numerical integration are available but are computationally intensive requiring a large amount of computer processing time that increases with the number of clusters (or individuals) in the data and are not always readily accessible to the practitioner in standard software. For the purposes of analysing categorical response data from a longitudinal social survey, there is clearly a need to evaluate the existing procedures for estimating multinomial logit random effects model in terms of accuracy, efficiency and computing time. The computational time will have significant implications as to the preferred approach by researchers. In this paper we evaluate statistical software procedures that utilise adaptive Gaussian quadrature and MCMC methods, with specific application to modeling employment status of women using a GLMM, over three waves of the HILDA survey.

Class attainment among British men - A multivariate extension of the CASMIN model of intergenerational class mobility

The CASMIN Project is arguably the most influential contemporary study of class mobility in the world. However, CASMIN results with respect to weak vertical status effects on class mobility have been extensively criticized. Drawing on arguments about how to model vertical mobility, Hout and Hauser (1992) show that class mobility is strongly determined by vertical socioeconomic differences. This paper extends these arguments by estimating the CASMIN model while explicitly controlling for individual determinants of socioeconomic attainment. Using the 1972 Oxford Mobility Data and the 1979 and 1983 British Election Studies, the paper employs mixed legit models to show how individual socioeconomic factors and categorical differences between classes shape intergenerational mobility. The findings highlight the multidimensionality of class mobility and its irreducibility to vertical movement up and down a stratification hierarchy.

Development and implementation of the population Fisher information matrix for the evaluation of population pharmacokinetic designs.

In population pharmacokinetic studies, the precision of parameter estimates is dependent on the population design. Methods based on the Fisher information matrix have been developed and extended to population studies to evaluate and optimize designs. In this paper we propose simple programming tools to evaluate population pharmacokinetic designs. This involved the development of an expression for the Fisher information matrix for nonlinear mixed-effects models, including estimation of the variance of the residual error. We implemented this expression as a generic function for two software applications: S-PLUS and MATLAB. The evaluation of population designs based on two pharmacokinetic examples from the literature is shown to illustrate the efficiency and the simplicity of this theoretic approach. Although no optimization method of the design is provided, these functions can be used to select and compare population designs among a large set of possible designs, avoiding a lot of simulations.

Population pharmacokinetics of lamotrigine

The present study estimated the population pharmacokinetics of lamotrigine in patients receiving oral lamotrigine therapy with drug concentration monitoring, and determined intersubject and intrasubject variability. A total of 129 patients were analyzed from two clinical sites. Of these, 124 patients provided spare data (198 concentration-time points); nine patients (four from a previous group plus five from the current group) provided rich data (431 points). The population analysis was conducted using P-PHARM (TM) (SIMED Scientific Software, Cedex, France), a nonlinear mixed-effect modeling program. A single exponential elimination model (first-order absorption) with heteroscedastic weighting was used. Apparent clearance (CL/F) and volume of distribution (V/F) were the pharmacokinetic parameters estimated. Covariate analysis was performed to determine which factors explained any of the variability associated with lamotrigine clearance. Population estimates of CL/F and V/F for lamotrigine generated in the final model were 2.14 +/- 0.81 L/h and 78.1 +/- 5.1 L/kg. Intersubject and intrasubject variability for clearance was 38% and 38%, respectively. The covariates of concomitant valproate and phenytoin therapy accounted for 42% of the intersubject variability of clearance. Age, gender, clinic site, and other concomitant antiepileptic drugs did not influence clearance. This study of the population pharmacokinetics of lamotrigine in patients using the drug clinically provides useful data and should lead to better dosage individualization for lamotrigine.

Some vocabulary and grammar for the analysis of multi-environment trials, as applied to the analysis of FPB and PPB trials

For the improvement of genetic material suitable for on farm use under low-input conditions, participatory and formal plant breeding strategies are frequently presented as competing options. A common frame of reference to phrase mechanisms and purposes related to breeding strategies will facilitate clearer descriptions of similarities and differences between participatory plant breeding and formal plant breeding. In this paper an attempt is made to develop such a common framework by means of a statistically inspired language that acknowledges the importance of both on farm trials and research centre trials as sources of information for on farm genetic improvement. Key concepts are the genetic correlation between environments, and the heterogeneity of phenotypic and genetic variance over environments. Classic selection response theory is taken as the starting point for the comparison of selection trials (on farm and research centre) with respect to the expected genetic improvement in a target environment (low-input farms). The variance-covariance parameters that form the input for selection response comparisons traditionally come from a mixed model fit to multi-environment trial data. In this paper we propose a recently developed class of mixed models, namely multiplicative mixed models, also called factor-analytic models, for modelling genetic variances and covariances (correlations). Mixed multiplicative models allow genetic variances and covariances to be dependent on quantitative descriptors of the environment, and confer a high flexibility in the choice of variance-covariance structure, without requiring the estimation of a prohibitively high number of parameters. As a result detailed considerations regarding selection response comparisons are facilitated. ne statistical machinery involved is illustrated on an example data set consisting of barley trials from the International Center for Agricultural Research in the Dry Areas (ICARDA). Analysis of the example data showed that participatory plant breeding and formal plant breeding are better interpreted as providing complementary rather than competing information.