Mixed model formulations for multi-environment trials

When studying genotype X environment interaction in multi-environment trials, plant breeders and geneticists often consider one of the effects, environments or genotypes, to be fixed and the other to be random. However, there are two main formulations for variance component estimation for the mixed model situation, referred to as the unconstrained-parameters (UP) and constrained-parameters (CP) formulations. These formulations give different estimates of genetic correlation and heritability as well as different tests of significance for the random effects factor. The definition of main effects and interactions and the consequences of such definitions should be clearly understood, and the selected formulation should be consistent for both fixed and random effects. A discussion of the practical outcomes of using the two formulations in the analysis of balanced data from multi-environment trials is presented. It is recommended that the CP formulation be used because of the meaning of its parameters and the corresponding variance components. When managed (fixed) environments are considered, users will have more confidence in prediction for them but will not be overconfident in prediction in the target (random) environments. Genetic gain (predicted response to selection in the target environments from the managed environments) is independent of formulation.

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.

Feeding ecology of the fish ectoparasite Gnathia sp. (Crustacea: Isopoda) from the Great Barrier Reef, and its implications for fish cleaning behaviour

The feeding rate of a parasitic gnathiid isopod on fish was examined. Individual fish, Hemigymnus melapterus, were exposed to gnathiid larvae and sampled after 5, 10, 30, 60, and 240 min. I recorded whether larvae had an engorged gut, an engorged gut containing red material, or had dropped off the fish after having completed engorgement; variation among sampling times and larval stages was analyzed using generalized linear mixed model analyses. The likelihood that larvae had an engorged gut increased with time and varied with larval stage. First stage (1.45 mm) larvae. After 30 min, however, most (>93%) larvae had an engorged gut regardless of their larval stage. The likelihood of red material in the gut of third stage larvae increased over time (46% after 30 min, 70% after 60 min, and 86% after 240 min) while that of first and second stage larvae remained relatively low (

Analysis of Melanoma Onset: Assessing Familial Aggregation by Using Estimating Equations and Fitting Variance Components via Bayesian Random Effects Models

We investigate whether relative contributions of genetic and shared environmental factors are associated with an increased risk in melanoma. Data from the Queensland Familial Melanoma Project comprising 15,907 subjects arising from 1912 families were analyzed to estimate the additive genetic, common and unique environmental contributions to variation in the age at onset of melanoma. Two complementary approaches for analyzing correlated time-to-onset family data were considered: the generalized estimating equations (GEE) method in which one can estimate relationship-specific dependence simultaneously with regression coefficients that describe the average population response to changing covariates; and a subject-specific Bayesian mixed model in which heterogeneity in regression parameters is explicitly modeled and the different components of variation may be estimated directly. The proportional hazards and Weibull models were utilized, as both produce natural frameworks for estimating relative risks while adjusting for simultaneous effects of other covariates. A simple Markov Chain Monte Carlo method for covariate imputation of missing data was used and the actual implementation of the Bayesian model was based on Gibbs sampling using the free ware package BUGS. In addition, we also used a Bayesian model to investigate the relative contribution of genetic and environmental effects on the expression of naevi and freckles, which are known risk factors for melanoma.

Large-eddy simulation of heat transfer downstream of a backward-facing step

Large-eddy simulation is used to predict heat transfer in the separated and reattached flow regions downstream of a backward-facing step. Simulations were carried out at a Reynolds number of 28 000 (based on the step height and the upstream centreline velocity) with a channel expansion ratio of 1.25. The Prandtl number was 0.71. Two subgrid-scale models were tested, namely the dynamic eddy-viscosity, eddy-diffusivity model and the dynamic mixed model. Both models showed good overall agreement with available experimental data. The simulations indicated that the peak in heat-transfer coefficient occurs slightly upstream of the mean reattachment location, in agreement with experimental data. The results of these simulations have been analysed to discover the mechanisms that cause this phenomenon. The peak in heat-transfer coefficient shows a direct correlation with the peak in wall shear-stress fluctuations. It is conjectured that the peak in these fluctuations is caused by an impingement mechanism, in which large eddies, originating in the shear layer, impact the wall just upstream of the mean reattachment location. These eddies cause a 'downwash', which increases the local heat-transfer coefficient by bringing cold fluid from above the shear layer towards the wall.

Ventilatory changes following head-up tilt and standing in healthy subjects

Passive tilting increases ventilation in healthy subjects; however, controversy surrounds the proposed mechanism. This study is aimed to evaluate the possible mechanism for changes to ventilation following passive head-up tilt (HUT) and active standing by comparison of a range of ventilatory, metabolic and mechanical parameters. Ventilatory parameters (V (T), V (E), V (E)/VO2, V (E)/VCO2, f and PetCO(2)), functional residual capacity (FRC), respiratory mechanics with impulse oscillometry; oxygen consumption (VO2) and carbon dioxide production (VCO2) were measured in 20 healthy male subjects whilst supine, following HUT to 70 degrees and unsupported standing. Data were analysed using a linear mixed model. HUT to 70 degrees from supine increased minute ventilation (V (E)) (P < 0.001), tidal volume (V (T)) (P=0.001), ventilatory equivalent for O-2 (V (E)/VO2) (P=0.020) and the ventilatory equivalent for CO2 (V (E)/VCO2) (P < 0.001) with no change in f (P=0.488). HUT also increased FRC (P < 0.001) and respiratory system reactance (X5Hz) (P < 0.001) with reduced respiratory system resistance (R5Hz) (P=0.004) and end-tidal carbon dioxide (PetCO(2)) (P < 0.001) compared to supine. Standing increased V (E) (P < 0.001), V (T) (P < 0.001) and V (E)/VCO2 (P=0.020) with no change in respiratory rate (f) (P=0.065), V (E)/VO2 (P=0.543). Similar changes in FRC (P < 0.001), R5Hz (P=0.013), X5Hz (P < 0.001) and PetCO(2) (P < 0.001) compared to HUT were found. In contrast to HUT, standing increased VO2 (P=0.002) and VCO2 (P=0.048). The greater increase in V (E) in standing compared to HUT appears to be related to increased VO2 and VCO2 associated with increased muscle activity in the unsupported standing position. This has implications for exercise prescription and rehabilitation of critically ill patients who have reduced cardiovascular and respiratory reserve.

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.

A Score Test for Zero-Inflation in Correlated Count Data

To account for the preponderance of zero counts and simultaneous correlation of observations, a class of zero-inflated Poisson mixed regression models is applicable for accommodating the within-cluster dependence. In this paper, a score test for zero-inflation is developed for assessing correlated count data with excess zeros. The sampling distribution and the power of the test statistic are evaluated by simulation studies. The results show that the test statistic performs satisfactorily under a wide range of conditions. The test procedure is further illustrated using a data set on recurrent urinary tract infections. Copyright (c) 2005 John Wiley & Sons, Ltd.

Genotypic consistency in the expression of leaf water potential in rice (Oryza sativa L.)

Early work has shown variation in the grain yield of rice cultivars grown under water stress conditions to be associated with the plant water status, mainly with the maintenance of high leaf water potential (LWP) at flowering and grain filling stage. Considerable variation for LWP among rice varieties has been recorded. The present work was designed to investigate genotypic consistency in water potential within the plant and under canopy manipulation to vary plant water requirement. In a glasshouse experiment, with six rice genotypes, a consistent water potential gradient from stem base to leaf tip has been observed. Leaf tip water potential has been found as the minimum LWP that can be recorded at any time of stress. Genotypes with similar canopy size could maintain different levels of LWP under stress conditions. In a field experiment, with four selected lines, four canopy sizes and two canopy mixture treatments were introduced prior to the imposition of control, mild and severe water stress conditions. It was found that the line differences in LWP and relative water content (RWC) were expressed under both mild and severe stress conditions, regardless of canopy size, tiller number and whether they were mixed with another line with different capacity to maintain LWP. Although there were some differences among canopy size treatments for radiation interception in three water conditions, canopy manipulation (plant size) within a line did not affect the expression of LWP and hence genotypic variation in LWP was maintained. Under both glasshouse and field conditions, lines that maintained high LWP had larger xylem diameter and stem areas than those that had low LWP. The results indicated that the size of the vascular bundles could influence the maintenance of plant water relations under water deficit. (c) 2005 Elsevier B.V. All rights reserved.

Multi-level zero-inflated Poisson regression modelling of correlated count data with excess zeros

Count data with excess zeros relative to a Poisson distribution are common in many biomedical applications. A popular approach to the analysis of such data is to use a zero-inflated Poisson (ZIP) regression model. Often, because of the hierarchical Study design or the data collection procedure, zero-inflation and lack of independence may occur simultaneously, which tender the standard ZIP model inadequate. To account for the preponderance of zero counts and the inherent correlation of observations, a class of multi-level ZIP regression model with random effects is presented. Model fitting is facilitated using an expectation-maximization algorithm, whereas variance components are estimated via residual maximum likelihood estimating equations. A score test for zero-inflation is also presented. The multi-level ZIP model is then generalized to cope with a more complex correlation structure. Application to the analysis of correlated count data from a longitudinal infant feeding study illustrates the usefulness of the approach.

Iron catalysed assembly of an asymmetric mixed-ligand triple helicate

The 2-pyridinecarbaldehyde isonicotinoyl hydrazone (HPCIH) family of ligands are typically tridentate N,N,O chelators that exhibit very high in vitro activity in mobilizing intracellular Fe and are promising candidates for the treatment of Fe overload diseases. Complexation of ferrous perchlorate with HPCIH in MeCN solution gives the expected six-coordinate complex Fe-II(PCIH)(2). However, complexation of Fe-II with 2-pyridinecarbaldehyde picolinoyl hydrazone (HPCPH, an isomer of HPCIH) under the same conditions leads to spontaneous assembly of an unprecedented asymmetric, mixed-ligand dinuclear triple helical complex Fe-2(II)(PCPH)(2)(PPH), where PPH2- is the dianion of bis(picolinoyl) hydrazine. The X-ray crystal structure of this complex shows that each ligand binds simultaneously to both metal centres in a bidentate fashion. The dinuclear complex exhibits two well separated and totally reversible Fe-III/II redox couples as shown by cyclic voltammetry in MeCN solution.

The BeWo choriocarcinoma cell line as a model of iodide transport by placenta

Cultured human choriocarcinoma cells of the BeWo line exhibited saturable accumulation of radioiodide. Inhibition by competing anions followed the affinity series perchlorate >= iodide >= thiocyanate, consistent with uptake through the thyroid iodide transporter, NIS, whose messenger RNA was found in BeWo cells, and whose protein was distributed towards the apical pole of the cells. Efflux obeyed first order kinetics and was inhibited by DIDS, an antagonist of anion exchangers including pendrin, whose messenger RNA was also present. In cultures where iodide uptake through NIS was blocked with excess perchlorate, radiolodide accumulation was stimulated by exposure to medium in which physiological anions were replaced by 2-morpholinoethanesulfonic acid (MES), consistent with the operation of an anion exchange mechanism taking up iodide. Chloride in the medium was more effective than sulfate at inhibiting this uptake, matching the ionic specificity of pendrin. These studies provide evidence that the trophoblast accumulates iodide through NIS and releases it to the fetal compartment through pendrin. (c) 2004 Elsevier Ltd. All rights reserved.

Finite deformation model of simple shear of fault with microrotations: apparent strain localisation and en-echelon fracture pattern

Strain localisation is a widespread phenomenon often observed in shear and compressive loading of geomaterials, for example, the fault gouge. It is believed that the main mechanisms of strain localisation are strain softening and mismatch between dilatancy and pressure sensitivity. Observations show that gouge deformation is accompanied by considerable rotations of grains. In our previous work as a model for gouge material, we proposed a continuum description for an assembly of particles of equal radius in which the particle rotation is treated as an independent degree of freedom. We showed that there exist critical values of the model parameters for which the displacement gradient exhibits a pronounced localisation at the mid-surface layers of the fault, even in the absence of inelasticity. Here, we generalise the model to the case of finite deformations characteristic for the gouge deformation. We derive objective constitutive relationships relating the Jaumann rates of stress and moment stress to the relative strain and curvature rates, respectively. The model suggests that the pattern of localisation remains the same as in the linear case. However, the presence of the Jaumann terms leads to the emergence of non-zero normal stresses acting along and perpendicular to the shear layer (with zero hydrostatic pressure), and localised along the mid-line of the gouge; these stress components are absent in the linear model of simple shear. These additional normal stresses, albeit small, cause a change in the direction in which the maximal normal stresses act and in which en-echelon fracturing is formed.