Formulating Mixed Models for Experiments, Including Longitudinal Experiments

Mixed models have become important in analyzing the results of experiments, particularly those that require more complicated models (e.g., those that involve longitudinal data). This article describes a method for deriving the terms in a mixed model. Our approach extends an earlier method by Brien and Bailey to explicitly identify terms for which autocorrelation and smooth trend arising from longitudinal observations need to be incorporated in the model. At the same time we retain the principle that the model used should include, at least, all the terms that are justified by the randomization. This is done by dividing the factors into sets, called tiers, based on the randomization and determining the crossing and nesting relationships between factors. The method is applied to formulate mixed models for a wide range of examples. We also describe the mixed model analysis of data from a three-phase experiment to investigate the effect of time of refinement on Eucalyptus pulp from four different sources. Cubic smoothing splines are used to describe differences in the trend over time and unstructured covariance matrices between times are found to be necessary.

Evaluating responses of maize (Zea mays L.) to soil physical conditions using a boundary line approach

The functional relation between the decline in the rate of a physiological process and the magnitude of a stress related to soil physical conditions is an important tool for uses as diverse as assessment of the stress-related sensitivity of different plant cultivars and characterization of soil structure. Two of the most pervasive sources of stress are soil resistance to root penetration (SR) and matric potential (psi). However, the assessment of these sources of stress on physiological processes in different soils can be complicated by other sources of stress and by the strong relation between SR and psi in a soil. A multivariate boundary line approach was assessed as a means of reducing these cornplications. The effects of SR and psi stress conditions on plant responses were examined under growth chamber conditions. Maize plants (Zea mays L.) were grown in soils at different water contents and having different structures arising from variation in texture, organic carbon content and soil compaction. Measurements of carbon exchange (CE), leaf transpiration (ILT), plant transpiration (PT), leaf area (LA), leaf + shoot dry weight (LSDW), root total length (RTL), root surface area (RSA) and root dry weight (RDW) were determined after plants reached the 12-leaf stage. The LT, PT and LA were described as a function of SR and psi with a double S-shaped function using the multivariate boundary line approach. The CE and LSDW were described by the combination of an S-shaped function for SR and a linear function for psi. The root parameters were described by a single S-shaped function for SR. The sensitivity to SR and psi depended on the plant parameter. Values of PT, LA and LSDW were most sensitive to SR. Among those parameters exhibiting a significant response to psi, PT was most sensitive. The boundary line approach was found to be a useful tool to describe the functional relation between the decline in the rate of a physiological process and the magnitude of a stress related to soil physical conditions. (C) 2009 Elsevier B.V. All rights reserved.

Using Imaging Spectroscopy to study soil properties

Imaging Spectroscopy (IS) is a promising tool for studying soil properties in large spatial domains. Going from point to image spectrometry is not only a journey from micro to macro scales, but also a long stage where problems such as dealing with data having a low signal-to-noise level, contamination of the atmosphere, large data sets, the BRDF effect and more are often encountered. In this paper we provide an up-to-date overview of some of the case studies that have used IS technology for soil science applications. Besides a brief discussion on the advantages and disadvantages of IS for studying soils, the following cases are comprehensively discussed: soil degradation (salinity, erosion, and deposition), soil mapping and classification, soil genesis and formation, soil contamination, soil water content, and soil swelling. We review these case studies and suggest that the 15 data be provided to the end-users as real reflectance and not as raw data and with better signal-to-noise ratios than presently exist. This is because converting the raw data into reflectance is a complicated stage that requires experience, knowledge, and specific infrastructures not available to many users, whereas quantitative spectral models require good quality data. These limitations serve as a barrier that impedes potential end-users, inhibiting researchers from trying this technique for their needs. The paper ends with a general call to the soil science audience to extend the utilization of the IS technique, and it provides some ideas on how to propel this technology forward to enable its widespread adoption in order to achieve a breakthrough in the field of soil science and remote sensing. (C) 2009 Elsevier Inc. All rights reserved.

A Log-Linear Regression Model for the Beta-Weibull Distribution

We introduce the log-beta Weibull regression model based on the beta Weibull distribution (Famoye et al., 2005; Lee et al., 2007). We derive expansions for the moment generating function which do not depend on complicated functions. The new regression model represents a parametric family of models that includes as sub-models several widely known regression models that can be applied to censored survival data. We employ a frequentist analysis, a jackknife estimator, and a parametric bootstrap for the parameters of the proposed model. We derive the appropriate matrices for assessing local influences on the parameter estimates under different perturbation schemes and present some ways to assess global influences. Further, for different parameter settings, sample sizes, and censoring percentages, several simulations are performed. In addition, the empirical distribution of some modified residuals are displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be extended to a modified deviance residual in the proposed regression model applied to censored data. We define martingale and deviance residuals to evaluate the model assumptions. The extended regression model is very useful for the analysis of real data and could give more realistic fits than other special regression models.

Staphylococcus hominis subsp novobiosepticus strains causing nosocomial bloodstream infection in Brazil

To report the isolation of six Staphylococcus hominis subsp. novobiosepticus (SHN) strains from hospitalized patients with bloodstream infections in two Brazilian hospitals and to characterize their susceptibility profile to several antimicrobials. Species identification was performed by biochemical methods and sodA gene sequencing. The MICs of antimicrobials were determined by broth and agar dilution methods and by Etest. Isolates were typed by PFGE and PCR amplification was used to detect the ccr gene complex and the mec class. Morphometric evaluation of cell wall was performed by transmission electron microscopy (TEM). Susceptibility profiles indicated that the majority of isolates (five) were multidrug-resistant. Overlapping and multiplex PCR showed that five out of the six strains harboured SCCmec type III with class A mec and type 3 ccr. The initial vancomycin MIC value of 4 mg/L for these strains increased to 16-32 mg/L after growth for 10 days in BHI broth supplemented with this antimicrobial. TEM indicated that vancomycin resistance was associated with cell wall thickening and to another mechanism not fully elucidated. Only one SHN strain was oxacillin- and vancomycin-susceptible. The nosocomial infections in at least five of the patients from both hospitals were caused by a single clone of SHN. It is very important to consider SHN strains as the cause of nosocomial infections. The clinical implications resulting from the pattern of multidrug resistance in these strains may be complicated by the emergence of vancomycin resistance.

Bubbly Flow Structure in Hydraulic Jump

In an open channel, a hydraulic jump is the rapid transition from super- to sub-critical flow associated with strong turbulence and air bubble entrainment in the mixing layer. New experiments were performed at relatively large Reynolds numbers using phase-detection probes. Some new signal analysis provided characteristic air-water time and length scales of the vortical structures advecting the air bubbles in the developing shear flow. An analysis of the longitudinal air-water flow structure suggested little bubble clustering in the mixing layer, although an interparticle arrival time analysis showed some preferential bubble clustering for small bubbles with chord times below 3 ms. Correlation analyses yielded longitudinal air-water time scales Txx*V1/d1 of about 0.8 in average. The transverse integral length scale Z/d1 of the eddies advecting entrained bubbles was typically between 0.25 and 0.4, irrespective of the inflow conditions within the range of the investigations. Overall the findings highlighted the complicated nature of the air-water flow

Turbulent Time and Length Scale Measurements in High-Velocity Open Channel Flows

In high-velocity open channel flows, the measurements of air-water flow properties are complicated by the strong interactions between the flow turbulence and the entrained air. In the present study, an advanced signal processing of traditional single- and dual-tip conductivity probe signals is developed to provide further details on the air-water turbulent level, time and length scales. The technique is applied to turbulent open channel flows on a stepped chute conducted in a large-size facility with flow Reynolds numbers ranging from 3.8 E+5 to 7.1 E+5. The air water flow properties presented some basic characteristics that were qualitatively and quantitatively similar to previous skimming flow studies. Some self-similar relationships were observed systematically at both macroscopic and microscopic levels. These included the distributions of void fraction, bubble count rate, interfacial velocity and turbulence level at a macroscopic scale, and the auto- and cross-correlation functions at the microscopic level. New correlation analyses yielded a characterisation of the large eddies advecting the bubbles. Basic results included the integral turbulent length and time scales. The turbulent length scales characterised some measure of the size of large vortical structures advecting air bubbles in the skimming flows, and the data were closely related to the characteristic air-water depth Y90. In the spray region, present results highlighted the existence of an upper spray region for C > 0.95 to 0.97 in which the distributions of droplet chord sizes and integral advection scales presented some marked differences with the rest of the flow.

And then came Complex Project Management

The subject of management is renowned for its addiction to fads and fashions. Project Management is no exception. The issue of interest for this paper is the establishment of the 'College of Complex Project Managers' and their 'competency standard for complex project managers.' Both have generated significant interest in the Project Management community, and like any other human endeavour they should be subject to critical evaluation. The results of this evaluation show significant flaws in the definition of complex in this case, the process by which the College and its standard have emerged, and the content of the standard. However, there is a significant case for a portfolio of research that extends the existing bodies of knowledge into large-scale complicated (or major) projects that would be owned by the relevant practitioner communities, rather than focused on one organization. Research questions are proposed that would commence this stream of activity towards an intelligent synthesis of what is required to manage in both complicated and truly complex environments.

XSophe-Sophe-XeprView (R). A computer simulation software suite (v. 1.1.3) for the analysis of continuous wave EPR spectra

The XSophe-Sophe-XeprView((R)) computer simulation software suite enables scientists to easily determine spin Hamiltonian parameters from isotropic, randomly oriented and single crystal continuous wave electron paramagnetic resonance (CW EPR) spectra from radicals and isolated paramagnetic metal ion centers or clusters found in metalloproteins, chemical systems and materials science. XSophe provides an X-windows graphical user interface to the Sophe programme and allows: creation of multiple input files, local and remote execution of Sophe, the display of sophelog (output from Sophe) and input parameters/files. Sophe is a sophisticated computer simulation software programme employing a number of innovative technologies including; the Sydney OPera HousE (SOPHE) partition and interpolation schemes, a field segmentation algorithm, the mosaic misorientation linewidth model, parallelization and spectral optimisation. In conjunction with the SOPHE partition scheme and the field segmentation algorithm, the SOPHE interpolation scheme and the mosaic misorientation linewidth model greatly increase the speed of simulations for most spin systems. Employing brute force matrix diagonalization in the simulation of an EPR spectrum from a high spin Cr(III) complex with the spin Hamiltonian parameters g(e) = 2.00, D = 0.10 cm(-1), E/D = 0.25, A(x) = 120.0, A(y) = 120.0, A(z) = 240.0 x 10(-4) cm(-1) requires a SOPHE grid size of N = 400 (to produce a good signal to noise ratio) and takes 229.47 s. In contrast the use of either the SOPHE interpolation scheme or the mosaic misorientation linewidth model requires a SOPHE grid size of only N = 18 and takes 44.08 and 0.79 s, respectively. Results from Sophe are transferred via the Common Object Request Broker Architecture (CORBA) to XSophe and subsequently to XeprView((R)) where the simulated CW EPR spectra (1D and 2D) can be compared to the experimental spectra. Energy level diagrams, transition roadmaps and transition surfaces aid the interpretation of complicated randomly oriented CW EPR spectra and can be viewed with a web browser and an OpenInventor scene graph viewer.

Using the level set method to model endogenous lava dome growth

Modeling volcanic phenomena is complicated by free-surfaces often supporting large rheological gradients. Analytical solutions and analogue models provide explanations for fundamental characteristics of lava flows. But more sophisticated models are needed, incorporating improved physics and rheology to capture realistic events. To advance our understanding of the flow dynamics of highly viscous lava in Peléean lava dome formation, axi-symmetrical Finite Element Method (FEM) models of generic endogenous dome growth have been developed. We use a novel technique, the level-set method, which tracks a moving interface, leaving the mesh unaltered. The model equations are formulated in an Eulerian framework. In this paper we test the quality of this technique in our numerical scheme by considering existing analytical and experimental models of lava dome growth which assume a constant Newtonian viscosity. We then compare our model against analytical solutions for real lava domes extruded on Soufrière, St. Vincent, W.I. in 1979 and Mount St. Helens, USA in October 1980 using an effective viscosity. The level-set method is found to be computationally light and robust enough to model the free-surface of a growing lava dome. Also, by modeling the extruded lava with a constant pressure head this naturally results in a drop in extrusion rate with increasing dome height, which can explain lava dome growth observables more appropriately than when using a fixed extrusion rate. From the modeling point of view, the level-set method will ultimately provide an opportunity to capture more of the physics while benefiting from the numerical robustness of regular grids.

Compiling a GIS database of tree farms on Leyte Island

As part of ACIAR project ASEM/2003/052, Improving Financial Returns to Smallholder Tree Farmers in the Philippines, plantations of timber trees in Leyte Island, the Philippines were located using a systematic survey of the island. The survey was undertaken in order to compile a database of plantations which could be used to guide the planning of project activities. In addition to recording a range of qualitative and quantitative information for each plantation, the survey spatially referenced each site using a Global Positioning System (GPS) to electronic maps of the island which were held in a Geographical Information System (GIS). Microsoft Excel and Mapsource® software were used as the software links between GPS coordinates and the GIS. Mapping of farm positions was complicated by different datums being used for maps of Leyte Island and this caused GPS positions to be displaced from equivalent positions on the map. Photos of the sites were hyperlinked to their map positions in the GIS in order to assist staff to recall site characteristics.

Lesions caused by cardiovascular flukes (Digenea : Spirorchidae) in stranded green turtles (Chelonia mydas)

Evidence of infection with spirorchid flukes (Digenea: Spirorchidae) was sought at necropsy of 96 stranded green turtles, Chelonia mydas, that were examined during the course of a survey of marine turtle mortality in southeastern Queensland, Australia. Three species of spirorchid (Hapalotrema mehrai, H. postorchis, and Neospirorchis schistosomatoides) were identified. Severe disease due to spirorchid fluke infection (spirorchidiasis) was implicated as the principal cause of mortality in 10 turtles (10%), and appeared to be one of multiple severe problems in an additional 29 turtles (30%). Although flukes were observed in only 45% of stranded C. mydas in this study, presumed spirorchid fluke infection was diagnosed in an additional 53% of turtles, based principally on characteristic necropsy lesions and to a lesser extent on the histopathological detection of spirorchid eggs. Characteristic necropsy lesions included miliary spirorchid egg granulomas, which were observed most readily on serosal surfaces, particularly of the small intestine. Cardiovascular lesions included mural endocarditis, arteritis, and thrombosis, frequently accompanied by aneurysm formation. Resolution of thrombi was observed to occur via a combination of granuloma formation about indigestible components (spirorchid fluke egg shells) and exteriorization through the vessel wall, which resulted in granulomatous nodules on the adventitial surface. Septic aortic thrombosis complicated by disseminated bacterial infection, observed in five turtles, was recorded for the first time. Egg granulomas were ubiquitous in turtle tissues throughout this study. Although they generally appeared to be mild or incidental lesions, they were occasionally associated with severe multifocal granulomatous pneumonia or meningitis.

Parental reactions to infant death: The effects of resources and coping strategies

The aim of this research was to examine, from a stress and coping perspective, the effects of resources (both personal and environmental) and coping strategies on parental reactions to infant death. One hundred and twenty-seven parents (60 fathers, 67 mothers) participated in the study. The predictors of parental distress (background factors, resources, coping methods) were initially assessed at 4-6 weeks post-loss. Parental distress (assessed using a composite measure of psychiatric disturbance, physical symptoms, and perinatal grief) was further assessed at 6 months post-loss and at 15 months postloss. After control for the stability in adjustment across time, there was consistent evidence that higher levels of education were associated with lower levels of parental distress over time. Among mothers, the number of friends in whom mothers had the confidence to confide emerged as a positive predictor of adjustment to infant death. A reliance on problem-focused coping was associated with greater maternal distress at 6 months post-loss, whereas coping by seeking support was associated with less distress at 15 months post-loss. There is no evidence that background factors and resources influenced parental distress through coping.

The politics and ethics of relocated art

At the end of Word War II, Soviet occupation forces removed countless art objects from German soil. Some of them were returned during the 1950s, but most either disappeared for good or were stored away secretly in cellars of Soviet museums. The Cold War then covered the issue with silence. After the collapse of the Soviet Union, museums in St Petersburg and Moscow started to exhibit some of the relocated art for the first time in half a century. The unusual quality of the paintings-mostly impressionist masterpieces-not only attracted the attention of the international art community, but also triggered a diplomatic row between Russia and Germany. Both governments advanced moral and legal claims to ownership. To make things even more complicated, many of the paintings once belonged to private collectors, some of whom were Jews. Their descendants also entered the dispute. The basic premise of this article is that the political and ethical dimensions of relocated art can be understood most adequately by eschewing a single authorial standpoint. Various positions, sometimes incommensurable ones, are thus explored in an attempt to outline possibilities for an ethics of representation and a dialogical solution to the international problem that relocated art has become.

Effects of a program of intervention on parental distress following infant death

A longitudinal study of 144 patents (65 fathers, 79 mothers) was conducted to evaluate the effectiveness of a program of intervention in relieving the psychological distress of parents affected by infant death. Participants were assessed in terms of their psychiatric disturbance, depression, anxiety, physical symptoms, dyadic adjustment, and coping strategies. The experimental group (n = 84) was offered an intervention program comprising the use of specially designed resources and contact with a trained grief worker. A control group (n = 60) was given routine community care. Parental reactions were assessed at four to six weeks postloss (prior to the implementation of the intervention program), at six months postloss, and at 15 months postloss. A series of multivariate analyses of valiance revealed that the intervention was effective in reducing the distress of parents, particularly those assessed prior to the intervention as being at high-risk of developing mourning difficulties. Effects of the intervention were noted in terms of parents' overall psychiatric disturbance, marital quality, and paternal coping strategies.