Improved testing inference in mixed linear models

Mixed linear models are commonly used in repeated measures studies. They account for the dependence amongst observations obtained from the same experimental unit. Often, the number of observations is small, and it is thus important to use inference strategies that incorporate small sample corrections. In this paper, we develop modified versions of the likelihood ratio test for fixed effects inference in mixed linear models. In particular, we derive a Bartlett correction to such a test, and also to a test obtained from a modified profile likelihood function. Our results generalize those in [Zucker, D.M., Lieberman, O., Manor, O., 2000. Improved small sample inference in the mixed linear model: Bartlett correction and adjusted likelihood. Journal of the Royal Statistical Society B, 62,827-838] by allowing the parameter of interest to be vector-valued. Additionally, our Bartlett corrections allow for random effects nonlinear covariance matrix structure. We report simulation results which show that the proposed tests display superior finite sample behavior relative to the standard likelihood ratio test. An application is also presented and discussed. (C) 2008 Elsevier B.V. All rights reserved.

Changes in soil organic carbon storage under different agricultural management systems in the Southwest Amazon Region of Brazil

The Cerrado and Amazon regions of Brazil are probably the largest agricultural frontier of the world, and Could be a sink or source for C depending on the net effect of land use change and subsequent management on soil organic C pools. We evaluated the effects of agricultural management systems on soil organic C (SOC) stocks in the Brazilian states of Rondonia and Mato Grosso, and derived regional specific factors for soil C stock change associated with different management systems. We used 50 observations (data points) in this study, including 42 dealing with annual cropping practices and 8 dealing with perennial cropping, and analyzed the data in linear mixed-effect models. No tillage (NT) systems in Cerrado areas increased SOC Storage by 1.08 +/- 0.06 relative to SOC stocks under native conditions, while SOC storage increased by a modest factor of 1.01 +/- 0.17 in Cerradao and Amazon Forest conditions. Full tillage (FT) had negative effect on SOC storage relative to NT, decreasing SOC stocks by a factor of 0.94 +/- 0.04. but did not significantly reduce SOC stocks relative to native levels when adopted in the Cerrado region. Perennial cropping had a minimal impact on SOC stocks, estimated at a factor Value of 0.98 +/- 0.14, suggesting these systems maintain about 98% of the SOC stock found under native vegetation. The results Suggest that NT adoption may be increasing SOC with land use change from native vegetation to cropland management in the Cerrado region of Brazil. (C) 2009 Elsevier B.V. All rights reserved.

Relationship between median nerve somatosensory evoked potentials and spinal cord injury levels in patients with quadriplegia

Study design: Cross-sectional study. Objectives: To observe if there is a relationship between the level of injury by the American Spinal Cord Injury Association (ASIA) and cortical somatosensory evoked potential (SSEP) recordings of the median nerve in patients with quadriplegia. Setting: Rehabilitation Outpatient Clinic at the university hospital in Brazil. Methods: Fourteen individuals with quadriplegia and 8 healthy individuals were evaluated. Electrophysiological assessment of the median nerve was performed by evoked potential equipment. The injury level was obtained by ASIA. N(9), N(13) and N(20) were analyzed based on the presence or absence of responses. The parameters used for analyzing these responses were the latency and the amplitude. Data were analyzed using mixed-effect models. Results: N(9) responses were found in all patients with quadriplegia with a similar latency and amplitude observed in healthy individuals; N(13) responses were not found in any patients with quadriplegia. N(20) responses were not found in C5 patients with quadriplegia but it was present in C6 and C7 patients. Their latencies were similar to healthy individuals (P > 0.05) but the amplitudes were decreased (P < 0.05). Conclusion: This study suggests that the SSEP responses depend on the injury level, considering that the individuals with C6 and C7 injury levels, both complete and incomplete, presented SSEP recordings in the cortical area. It also showed a relationship between the level of spinal cord injury assessed by ASIA and the median nerve SSEP responses, through the latency and amplitude recordings. Spinal Cord (2009) 47, 372-378; doi:10.1038/sc.2008.147; published online 20 January 2009

The association between air pollution and blood pressure in traffic controllers in Santo Andre, Sao Paulo, Brazil

Background: Urban air pollutants are associated with cardiovascular events. Traffic controllers are at high risk for pollution exposure during outdoor work shifts. Objective: The purpose of this study was to evaluate the relationship between air pollution and systemic blood pressure in traffic controllers during their work shifts. Methods: This cross-sectional study enrolled 19 male traffic controllers from Santo Andre city (Sao Paulo, Brazil) who were 30-60 years old and exposed to ambient air during outdoor work shifts. Systolic and diastolic blood pressure readings were measured every 15 min by an Ambulatory Arterial Blood Pressure Monitoring device. Hourly measurements (lags of 0-5 h) and the moving averages (2-5 h) of particulate matter (PM(10)), ozone (O(3)) ambient concentrations and the acquired daily minimum temperature and humidity means from the Sao Paulo State Environmental Agency were correlated with both systolic and diastolic blood pressures. Statistical methods included descriptive analysis and linear mixed effect models adjusted for temperature, humidity, work periods and time of day. Results: Interquartile increases of PM(10) (33 mu g/m(3)) and O(3) (49 mu g/m(3)) levels were associated with increases in all arterial pressure parameters, ranging from 1.06 to 2.53 mmHg. PM(10) concentration was associated with early effects (lag 0), mainly on systolic blood pressure. However, O(3) was weakly associated most consistently with diastolic blood pressure and with late cumulative effects. Conclusions: Santo Andre traffic controllers presented higher blood pressure readings while working their outdoor shifts during periods of exposure to ambient pollutant fluctuations. However, PM(10) and O(3) induced cardiovascular effects demonstrated different time courses and end-point behaviors and probably acted through different mechanisms. (C) 2011 Elsevier Inc. All rights reserved.

ELECTROMYOGRAPHIC AMPLITUDE RATIO OF SERRATUS ANTERIOR AND UPPER TRAPEZIUS MUSCLES DURING MODIFIED PUSH-UPS AND BENCH PRESS EXERCISES

Imbalance and weakness of the serratus anterior and upper trapezius force couple have been described in patients with shoulder dysfunction. There is interest in identifying exercises that selectively activate these muscles and including it in rehabilitation protocols. This study aims to verify the UT/SA electromyographic (EMG) amplitude ratio, performed in different upper limb exercises and on two bases of support. Twelve healthy men were tested (average age = 22.8 +/- 3.1 years), and surface EMG was recorded from the upper trapezius and serratus anterior using single differential surface electrodes. Volunteers performed isometric contractions over a stable base of support and on a Swiss ball during the wall push-up (WP), bench press (BP), and push-up (PU) exercises. All SEMG data are reported as a percentage of root mean square or integral of linear envelope from the maximal value obtained in one of three maximal voluntary contractions for each muscle studied. A linear mixed-effect model was performed to compare UT/SA ratio values. The WP, BP, and PU exercises showed UT/SA ratio mean +/- SD values of 0.69 +/- 0.72, 0.14 +/- 0.12, and 0.39 +/- 0.37 for stable surfaces, respectively, whereas for unstable surfaces, the values were 0.73 +/- 0.67, 0.43 +/- 0.39, and 0.32 +/- 0.30. The results demonstrate that UT/SA ratio was influenced by the exercises and by the upper limb base of support. The practical application is to show that BP on a stable surface is the exercise preferred over WP and PU on either surfaces for serratus anterior muscle training in patients with imbalance between the UT/SA force couple or serratus anterior weakness.

Acupuncture for the prevention of radiation-induced xerostomia in patients with head and neck cancer

The aim of this study was to evaluate the effectiveness of acupuncture in minimizing the severity of radiation-induced xerostomia in patients with head and neck cancer. A total of 24 consecutive patients receiving > 5000 cGy radiotherapy (RT) involving the major salivary glands bilaterally were assigned to either the preventive acupuncture group (PA, n = 12), treated with acupuncture before and during RT, or the control group (CT, n = 12), treated with RT and not receiving acupuncture. After RT completion, clinical response was assessed in all patients by syalometry, measuring the resting (RSFR) and stimulated (SSFR) salivary flow rates, and by the visual analogue scale (VAS) regarding dry mouth-related symptoms. Statistical analyses were performed with repeated-measures using a mixed-effect modeling procedure and analysis of variance. An alpha level of 0.05 was accepted for statistical significance. Although all patients exhibited some degree of impairment in salivary gland functioning after RT, significant differences were found between the groups. Patients in the PA group showed improved salivary flow rates (RSFR, SSFR; p < 0.001) and decreased xerostomia-related symptoms (VAS, p < 0.05) compared with patients in the CT group. Although PA treatment did not prevent the oral sequelae of RT completely, it significantly minimized the severity of radiation-induced xerostomia. The results suggest that acupuncture focused in a preventive approach can be a useful therapy in the management of patients with head and neck cancer undergoing RT.

Family history of cancer: Pooled analysis in the International Head and Neck Cancer Epidemiology consortium

Alcohol and tobacco consumption are well-recognized risk factors for head and neck cancer (HNC). Evidence suggests that genetic predisposition may also play a role. Only a few epidemiologic studies, however, have considered the relation between HNC risk and family history of HNC and other cancers. We pooled individual-level data across 12 case-control studies including 8,967 HNC cases and 13,627 controls. We obtained pooled odds ratios (OR) using fixed and random effect models and adjusting for potential confounding factors. All statistical tests were two-sided. A family history of HNC in first-degree relatives increased the risk of HNC (OR = 1.7, 95% confidence interval, CI, 1.2-2.3). The risk was higher when the affected relative was a sibling (OR = 2.2, 95% CI 1.6-3.1) rather than a parent (OR = 1.5, 95% CI 1.1-1.8) and for more distal HNC anatomic sites (hypopharynx and larynx). The risk was also higher, or limited to, in subjects exposed to tobacco. The OR rose to 7.2 (95% CI 5.5-9.5) among subjects with family history, who were alcohol and tobacco users. A weak but significant association (OR = 1.1, 95% CI 1.0-1.2) emerged for family history of other tobacco-related neoplasms, particularly with laryngeal cancer (OR = 1.3, 95% CI 1.1-1.5). No association was observed for family history of nontobacco-related neoplasms and the risk of HNC (OR = 1.0, 95% CI 0.9-1.1). Familial factors play a role in the etiology of HNC. In both subjects with and without family history of HNC, avoidance of tobacco and alcohol exposure may be the best way to avoid HNC. (C) 2008 Wiley-Liss, Inc,

Toward a resolution of Campanulid phylogeny, with special reference to the placement of Dipsacales

Broad-scale phylogenetic analyses of the angiosperms and of the Asteridae have failed to confidently resolve relationships among the major lineages of the campanulid Asteridae (i.e., the euasterid II of APG II, 2003). To address this problem we assembled presently available sequences for a core set of 50 taxa, representing the diversity of the four largest lineages (Apiales, Aquifoliales, Asterales, Dipsacales) as well as the smaller ""unplaced"" groups (e.g., Bruniaceae, Paracryphiaceae, Columelliaceae). We constructed four data matrices for phylogenetic analysis: a chloroplast coding matrix (atpB, matK, ndhF, rbcL), a chloroplast non-coding matrix (rps16 intron, trnT-F region, trnV-atpE IGS), a combined chloroplast dataset (all seven chloroplast regions), and a combined genome matrix (seven chloroplast regions plus 18S and 26S rDNA). Bayesian analyses of these datasets using mixed substitution models produced often well-resolved and supported trees. Consistent with more weakly supported results from previous studies, our analyses support the monophyly of the four major clades and the relationships among them. Most importantly, Asterales are inferred to be sister to a clade containing Apiales and Dipsacales. Paracryphiaceae is consistently placed sister to the Dipsacales. However, the exact relationships of Bruniaceae, Columelliaceae, and an Escallonia clade depended upon the dataset. Areas of poor resolution in combined analyses may be partly explained by conflict between the coding and non-coding data partitions. We discuss the implications of these results for our understanding of campanulid phylogeny and evolution, paying special attention to how our findings bear on character evolution and biogeography in Dipsacales.

Performance of balanced two-stage empirical predictors of realized cluster latent values from finite populations: A simulation study

Predictors of random effects are usually based on the popular mixed effects (ME) model developed under the assumption that the sample is obtained from a conceptual infinite population; such predictors are employed even when the actual population is finite. Two alternatives that incorporate the finite nature of the population are obtained from the superpopulation model proposed by Scott and Smith (1969. Estimation in multi-stage surveys. J. Amer. Statist. Assoc. 64, 830-840) or from the finite population mixed model recently proposed by Stanek and Singer (2004. Predicting random effects from finite population clustered samples with response error. J. Amer. Statist. Assoc. 99, 1119-1130). Predictors derived under the latter model with the additional assumptions that all variance components are known and that within-cluster variances are equal have smaller mean squared error (MSE) than the competitors based on either the ME or Scott and Smith`s models. As population variances are rarely known, we propose method of moment estimators to obtain empirical predictors and conduct a simulation study to evaluate their performance. The results suggest that the finite population mixed model empirical predictor is more stable than its competitors since, in terms of MSE, it is either the best or the second best and when second best, its performance lies within acceptable limits. When both cluster and unit intra-class correlation coefficients are very high (e.g., 0.95 or more), the performance of the empirical predictors derived under the three models is similar. (c) 2007 Elsevier B.V. All rights reserved.

On estimation and influence diagnostics for the Grubbs` model under heavy-tailed distributions

The Grubbs` measurement model is frequently used to compare several measuring devices. It is common to assume that the random terms have a normal distribution. However, such assumption makes the inference vulnerable to outlying observations, whereas scale mixtures of normal distributions have been an interesting alternative to produce robust estimates, keeping the elegancy and simplicity of the maximum likelihood theory. The aim of this paper is to develop an EM-type algorithm for the parameter estimation, and to use the local influence method to assess the robustness aspects of these parameter estimates under some usual perturbation schemes, In order to identify outliers and to criticize the model building we use the local influence procedure in a Study to compare the precision of several thermocouples. (C) 2008 Elsevier B.V. All rights reserved.

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.

Mitochondrial sequence data and Dipsacales phylogeny: Mixed models, partitioned Bayesian analyses, and model selection

Phylogenetic analyses of chloroplast DNA sequences, morphology, and combined data have provided consistent support for many of the major branches within the angiosperm, clade Dipsacales. Here we use sequences from three mitochondrial loci to test the existing broad scale phylogeny and in an attempt to resolve several relationships that have remained uncertain. Parsimony, maximum likelihood, and Bayesian analyses of a combined mitochondrial data set recover trees broadly consistent with previous studies, although resolution and support are lower than in the largest chloroplast analyses. Combining chloroplast and mitochondrial data results in a generally well-resolved and very strongly supported topology but the previously recognized problem areas remain. To investigate why these relationships have been difficult to resolve we conducted a series of experiments using different data partitions and heterogeneous substitution models. Usually more complex modeling schemes are favored regardless of the partitions recognized but model choice had little effect on topology or support values. In contrast there are consistent but weakly supported differences in the topologies recovered from coding and non-coding matrices. These conflicts directly correspond to relationships that were poorly resolved in analyses of the full combined chloroplast-mitochondrial data set. We suggest incongruent signal has contributed to our inability to confidently resolve these problem areas. (c) 2007 Elsevier Inc. All rights reserved.

Leverage analysis for linear mixed models

We consider a generalized leverage matrix useful for the identification of influential units and observations in linear mixed models and show how a decomposition of this matrix may be employed to identify high leverage points for both the marginal fitted values and the random effect component of the conditional fitted values. We illustrate the different uses of the two components of the decomposition with a simulated example as well as with a real data set.

The effect of using mixed initiator systems on the efficiency of photopolymerization of dental resins

A study was performed in order to determine the efficiency of the simultaneous use of the photoinitiators phenylpropanedione (PPD) and camphorquinone (CQ) in the polymerization of acrylic polymers and evaluate possible mechanisms leading to synergism or antagonism. It was found that efficiencies of both initiators taken individually are higher than that of their mixture, indicating that when both dyes are used simultaneously there will be an energy transfer from the more efficient initiator (CQ) to the less efficient one (PPD). Also, there was no proof of any reaction between the amine present in the CQ formulation and the PPD excited state.

Detection of catechol using mixed Langmuir-Blodgett films of a phospholipid and phthalocyanines as voltammetric sensors

The combination of metallic phthalocyanines (MPcs) and biomolecules has been explored in the literature either as mimetic systems to investigate molecular interactions or as supporting layers to immobilize biomolecules. Here, Langmuir-Blodgett (LB) films containing the phospholipid dimyristoyl phosphatidic acid (DMPA) mixed either with iron phthalocyanine (FePc) or with lutetium bisphthalocyanine (LuPc(2)) were applied as ITO modified-electrodes in the detection of catechol using cyclic voltammetry. The mixed Langmuir films of FePc + DMPA and LuPc(2) + DMPA displayed surface-pressure isotherms with no evidence of molecular-level interactions. The Fourier Transform Infrared (FTIR) spectra of the multilayer LB films confirmed the lack of interaction between the components. The DMPA and the FePc molecules were found to be oriented perpendicularly to the substrate, while LuPc(2) molecules were randomly organized. The phospholipid matrix induced a remarkable electrocatalytic effect on the phthalocyanines; as a result the mixed LB films deposited on ITO could be used to detect catechol with detection limits of 4.30 x 10(-7) and 3.34 x 10(-7) M for FePc + DMPA and LuPc(2) + DMPA, respectively. Results from kinetics experiments revealed that ion diffusion dominated the response of the modified electrodes. The sensitivity was comparable to that of other non-enzymatic sensors, which is sufficient to detect catechol in the food industry. The higher stability of the electrochemical response of the LB films and the ability to control the molecular architecture are promising for further studies with incorporation of biomolecules.