Influence diagnostics in nonlinear mixed-effects elliptical models

In this work we propose and analyze nonlinear elliptical models for longitudinal data, which represent an alternative to gaussian models in the cases of heavy tails, for instance. The elliptical distributions may help to control the influence of the observations in the parameter estimates by naturally attributing different weights for each case. We consider random effects to introduce the within-group correlation and work with the marginal model without requiring numerical integration. An iterative algorithm to obtain maximum likelihood estimates for the parameters is presented, as well as diagnostic results based on residual distances and local influence [Cook, D., 1986. Assessment of local influence. journal of the Royal Statistical Society - Series B 48 (2), 133-169; Cook D., 1987. Influence assessment. journal of Applied Statistics 14 (2),117-131; Escobar, L.A., Meeker, W.Q., 1992, Assessing influence in regression analysis with censored data, Biometrics 48, 507-528]. As numerical illustration, we apply the obtained results to a kinetics longitudinal data set presented in [Vonesh, E.F., Carter, R.L., 1992. Mixed-effects nonlinear regression for unbalanced repeated measures. Biometrics 48, 1-17], which was analyzed under the assumption of normality. (C) 2009 Elsevier B.V. All rights reserved.

QCD sum rules for the production of the X(3872) as a mixed molecule-charmonium state in B meson decay

We use QCD sum rules to calculate the branching ratio for the production of the meson X(3872) in the decay B -> X(3872)K, assumed to be a mixture between charmonium and exotic molecular vertical bar c (q) over bar vertical bar vertical bar q (c) over bar vertical bar states with J(PC) = 1(++). We find that in a small range for the values of the mixing angle, 5 degrees <= theta <= 13 degrees, we get the branching ratio B(B -> XK) = (1.00 +/- 0.68) x 10(-5), which is in agreement with the experimental upper limit. This result is compatible with the analysis of the mass and decay width of the mode J/psi(n pi) and the radiative decay mode J/psi gamma performed in the same approach. (C) 2011 Elsevier B.V. All rights reserved.

Electric field gradients in (111)In-doped (Hf/Zr)(3)Al(2) and (Hf/Zr)(4)Al(3) mixed compounds: ab initio calculations, perturbed angular correlation measurements and site preference

The quadrupolar hyperfine interactions of in-diffused (111)In -> (111)Cd probes in polycrystalline isostructural Zr(4)Al(3) and Hf(4)Al(3) samples containing small admixtures of the phases (Zr/Hf)(3)Al(2) were investigated. A strong preference of (111)In solutes for the contaminant (Zr/Hf)(3)Al(2) minority phases was observed. Detailed calculations of the electric field gradient (EFG) at the Cd nucleus using the full-potential augmented plane wave + local orbital formalism allowed us to assign the observed EFG fractions to the various lattice sites in the (Zr/Hf)(3)Al(2) compounds and to understand the preferential site occupation of the minority phases by the (111)In atoms. The effects of the size of the supercell and relaxation around the oversized In and Cd probe atoms were investigated in detail.

Chemical process to separate iron oxides particles in pottery sample for EPR dating

Ancient potteries usually are made of the local clay material, which contains relatively high concentration of iron. The powdered samples are usually quite black, due to magnetite, and, although they can be used for thermoluminescene (TL) dating, it is easiest to obtain better TL reading when clearest natural or pre-treated sample is used. For electron paramagnetic resonance (EPR) measurements, the huge signal due to iron spin-spin interaction, promotes an intense interference overlapping any other signal in this range. Sample dating is obtained by dividing the radiation dose, determined by the concentration of paramagnetic species generated by irradiation, by the natural dose so as a consequence, EPR dating cannot be used, since iron signal do not depend on radiation dose. In some cases, the density separation method using hydrated solution of sodium polytungstate [Na(G)(H(2)W(12)O(40))center dot H(2)O] becomes useful. However, the sodium polytungstate is very expensive in Brazil: hence an alternative method for eliminating this interference is proposed. A chemical process to eliminate about 90% of magnetite was developed. A sample of powdered ancient pottery was treated in a mixture (3:1:1) of HCI, HNO(3) and H(2)O(2) for 4 h. After that, it was washed several times in distilled water to remove all acid matrixes. The original black sample becomes somewhat clearer. The resulting material was analyzed by plasma mass spectrometry (ICP-MS), with the result that the iron content is reduced by a factor of about 9. In EPR measurements a non-treated natural ceramic sample shows a broad spin-spin interaction signal, the chemically treated sample presents a narrow signal in g= 2.00 region, possibly due to a radical of (SiO(3))(3-), mixed with signal of remaining iron [M. lkeya, New Applications of Electron Spin Resonance, World Scientific, Singapore, 1993, p. 285]. This signal increases in intensity under -gamma-irradiation. However, still due to iron influence, the additive method yielded too old age-value. Since annealing at 300 degrees C, Toyoda and Ikeya IS. Toyoda, M. Ikeya, Geochem. J. 25 (1991) 427-445] states that E `(1)-signal with maximum intensity is obtained, while annealing at 400 degrees C E`(1)-signal is completely eliminated, the subtraction of the second one from 300 degrees C heat-treated sample isolate E`(1)-like signal. Since this is radiation dose-dependent, we show that now EPR dating becomes possible. (C) 2008 Elsevier B.V. All rights reserved.

On an Electrode Modified by a Supramolecular Ruthenium Mixed Valence (Ru(II)/Ru(III)) Diphosphine-Porphyrin Assembly

Three novel polymetallic ruthenium (III) meso-tetra(4-pyridyl)porphyrins containing peripheral ""RuCl(3)(dppb)"" moieties have been prepared and characterized. The X-ray structure of the tetraruthenated {NiTPyP[RuCl(3)(dppb)](4)} porphyrin complex crystallizes in the triclinic space group FT. This structure is discussed and compared with the crystal data for the mer-[RuCl(3)(dppb)(py)]. The {TPyP[RuCl(3)(dppb)](4)} and {CoTPyP[RuCl(3)(dppb)](4)} porphyrins were used to obtain electrogenerated films on ITO and glass carbon electrode surfaces, respectively. Such tetraruthenated porphyrins form films of a mixed-valence species {TPyP[Ru(dppb)](4)(mu Cl(3))(2)}(2)(4n2+) and {CoTPyP[Ru(dppb)](4)(mu Cl(3))(2)}(2n)(4n2+) on the electrode surface. The modified electrode with {CoTPyP[RuCl(3)(dppb)](4)} is very stable and can be used to detect organic substrates such as catechol.

Minimum cycle cover and Chinese postman problems on mixed graphs with bounded tree-width

Let M = (V, E, A) be a mixed graph with vertex set V, edge set E and arc set A. A cycle cover of M is a family C = {C(1), ... , C(k)} of cycles of M such that each edge/arc of M belongs to at least one cycle in C. The weight of C is Sigma(k)(i=1) vertical bar C(i)vertical bar. The minimum cycle cover problem is the following: given a strongly connected mixed graph M without bridges, find a cycle cover of M with weight as small as possible. The Chinese postman problem is: given a strongly connected mixed graph M, find a minimum length closed walk using all edges and arcs of M. These problems are NP-hard. We show that they can be solved in polynomial time if M has bounded tree-width. (C) 2008 Elsevier B.V. All rights reserved.

Sampling, WLS, and Mixed Models

Mixed models may be defined with or without reference to sampling, and can be used to predict realized random effects, as when estimating the latent values of study subjects measured with response error. When the model is specified without reference to sampling, a simple mixed model includes two random variables, one stemming from an exchangeable distribution of latent values of study subjects and the other, from the study subjects` response error distributions. Positive probabilities are assigned to both potentially realizable responses and artificial responses that are not potentially realizable, resulting in artificial latent values. In contrast, finite population mixed models represent the two-stage process of sampling subjects and measuring their responses, where positive probabilities are only assigned to potentially realizable responses. A comparison of the estimators over the same potentially realizable responses indicates that the optimal linear mixed model estimator (the usual best linear unbiased predictor, BLUP) is often (but not always) more accurate than the comparable finite population mixed model estimator (the FPMM BLUP). We examine a simple example and provide the basis for a broader discussion of the role of conditioning, sampling, and model assumptions in developing inference.

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.

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.

Likelihood ratio tests for variance components in linear mixed models

Although the asymptotic distributions of the likelihood ratio for testing hypotheses of null variance components in linear mixed models derived by Stram and Lee [1994. Variance components testing in longitudinal mixed effects model. Biometrics 50, 1171-1177] are valid, their proof is based on the work of Self and Liang [1987. Asymptotic properties of maximum likelihood estimators and likelihood tests under nonstandard conditions. J. Amer. Statist. Assoc. 82, 605-610] which requires identically distributed random variables, an assumption not always valid in longitudinal data problems. We use the less restrictive results of Vu and Zhou [1997. Generalization of likelihood ratio tests under nonstandard conditions. Ann. Statist. 25, 897-916] to prove that the proposed mixture of chi-squared distributions is the actual asymptotic distribution of such likelihood ratios used as test statistics for null variance components in models with one or two random effects. We also consider a limited simulation study to evaluate the appropriateness of the asymptotic distribution of such likelihood ratios in moderately sized samples. (C) 2008 Elsevier B.V. All rights reserved.

Local Influence Analysis for Skew-Normal Linear Mixed Models

In this article, we consider local influence analysis for the skew-normal linear mixed model (SN-LMM). As the observed data log-likelihood associated with the SN-LMM is intractable, Cook`s well-known approach cannot be applied to obtain measures of local influence. Instead, we develop local influence measures following the approach of Zhu and Lee (2001). This approach is based on the use of an EM-type algorithm and is measurement invariant under reparametrizations. Four specific perturbation schemes are discussed. Results obtained for a simulated data set and a real data set are reported, illustrating the usefulness of the proposed methodology.

Predicting random effects with an expanded finite population mixed model

Prediction of random effects is an important problem with expanding applications. In the simplest context, the problem corresponds to prediction of the latent value (the mean) of a realized cluster selected via two-stage sampling. Recently, Stanek and Singer [Predicting random effects from finite population clustered samples with response error. J. Amer. Statist. Assoc. 99, 119-130] developed best linear unbiased predictors (BLUP) under a finite population mixed model that outperform BLUPs from mixed models and superpopulation models. Their setup, however, does not allow for unequally sized clusters. To overcome this drawback, we consider an expanded finite population mixed model based on a larger set of random variables that span a higher dimensional space than those typically applied to such problems. We show that BLUPs for linear combinations of the realized cluster means derived under such a model have considerably smaller mean squared error (MSE) than those obtained from mixed models, superpopulation models, and finite population mixed models. We motivate our general approach by an example developed for two-stage cluster sampling and show that it faithfully captures the stochastic aspects of sampling in the problem. We also consider simulation studies to illustrate the increased accuracy of the BLUP obtained under the expanded finite population mixed model. (C) 2007 Elsevier B.V. All rights reserved.

A note on influence diagnostics in nonlinear mixed-effects elliptical models

This paper provides general matrix formulas for computing the score function, the (expected and observed) Fisher information and the A matrices (required for the assessment of local influence) for a quite general model which includes the one proposed by Russo et al. (2009). Additionally, we also present an expression for the generalized leverage on fixed and random effects. The matrix formulation has notational advantages, since despite the complexity of the postulated model, all general formulas are compact, clear and have nice forms. (C) 2010 Elsevier B.V. All rights reserved.

Multiple stages of detergent-erythrocyte membrane interaction-A spin label study

The various stages of the interaction between the detergent Triton X-100 (TTX-100) and membranes of whole red blood cells (RBC) were investigated in a broad range of detergent concentrations. The interaction was monitored by RBC hemolysis-assessed by release of intracellular hemoglobin (Hb) and inorganic phosphate- and by analysis of EPR spectra of a fatty acid spin probe intercalated in whole RBC suspensions, as well as pellets and supernatants obtained upon centrifugation of detergent-treated cells. Hemolysis finished at ca. 0.9 mM TTX-100. Spectral analysis and calculation of order parameters (S) indicated that a complex sequence of events takes place, and allowed the characterization of various structures formed in the different stages of detergent-membrane interaction. Upon reaching the end of cell lysis, essentially no pellet was detected, the remaining EPR signal being found almost entirely in the supernatants. Calculated order parameters revealed that whole RBC suspensions, pellets, and supernatants possessed a similar degree of molecular packing, which decreased to a small extent up to 2.5 mM detergent. Between 3.2 and 10 mM TTX-100, a steep decrease in S was observed for both whole RBC suspensions and supernatants. Above 10 mM detergent, S decreased in a less pronounced manner and the EPR spectra approached that of pure TTX-100 micelles. The data were interpreted in terms of the following events: at the lower detergent concentrations, an increase in membrane permeability occurs: the end of hemolysis coincides with the lack of pellet upon centrifugation. Up to 2.5 mM TTX-100 the supernatants consist of a (very likely) heterogeneous population of membrane fragments with molecular packing similar to that of whole cells. As the detergent concentration increases, mixed micelles are formed containing lipid and/or protein, approaching the packing found in pure TTX-100 micelles. This analysis is in agreement with the models proposed by Lasch (Biochim. Biophys Acta 1241 (1995) 269-292) and by Le Maire and coworkers (Biochim. Biophys. Acta 1508 (2000) 86-111). (C) 2010 Elsevier B.V. All rights reserved.