A Poisson mixed model with nonnormal random effect distribution

In this paper, we propose a random intercept Poisson model in which the random effect is assumed to follow a generalized log-gamma (GLG) distribution. This random effect accommodates (or captures) the overdispersion in the counts and induces within-cluster correlation. We derive the first two moments for the marginal distribution as well as the intraclass correlation. Even though numerical integration methods are, in general, required for deriving the marginal models, we obtain the multivariate negative binomial model from a particular parameter setting of the hierarchical model. An iterative process is derived for obtaining the maximum likelihood estimates for the parameters in the multivariate negative binomial model. Residual analysis is proposed and two applications with real data are given for illustration. (C) 2011 Elsevier B.V. All rights reserved.

Influence diagnostics for elliptical semiparametric mixed models

In this paper we extend semiparametric mixed linear models with normal errors to elliptical errors in order to permit distributions with heavier and lighter tails than the normal ones. Penalized likelihood equations are applied to derive the maximum penalized likelihood estimates (MPLEs) which appear to be robust against outlying observations in the sense of the Mahalanobis distance. A reweighed iterative process based on the back-fitting method is proposed for the parameter estimation and the local influence curvatures are derived under some usual perturbation schemes to study the sensitivity of the MPLEs. Two motivating examples preliminarily analyzed under normal errors are reanalyzed considering some appropriate elliptical errors. The local influence approach is used to compare the sensitivity of the model estimates.

Trends of the genetic connectedness measures among Nelore beef cattle herds

Validity of comparisons between expected breeding values obtained from best linear unbiased prediction procedures in genetic evaluations is dependent on genetic connectedness among herds. Different cattle breeding programmes have their own particular features that distinguish their database structure and can affect connectedness. Thus, the evolution of these programmes can also alter the connectedness measures. This study analysed the evolution of the genetic connectedness measures among Brazilian Nelore cattle herds from 1999 to 2008, using the French Criterion of Admission to the group of Connected Herds (CACO) method, based on coefficients of determination (CD) of contrasts. Genetic connectedness levels were analysed by using simple and multiple regression analyses on herd descriptors to understand their relationship and their temporal trends from the 19992003 to the 20042008 period. The results showed a high level of genetic connectedness, with CACO estimates higher than 0.4 for the majority of them. Evaluation of the last 5-year period showed only a small increase in average CACO measures compared with the first 5 years, from 0.77 to 0.80. The percentage of herds with CACO estimates lower than 0.7 decreased from 27.5% in the first period to 16.2% in the last one. The connectedness measures were correlated with percentage of progeny from connecting sires, and the artificial insemination spread among Brazilian herds in recent years. But changes in connectedness levels were shown to be more complex, and their complete explanation cannot consider only herd descriptors. They involve more comprehensive changes in the relationship matrix, which can be only fully expressed by the CD of contrasts.

Transmission of Human Herpesvirus Type 8 Infection Within Families in American Indigenous Populations From the Brazilian Amazon

Background. The intrafamilial dynamics of endemic infection with human herpesvirus type 8 (HHV-8) in Amerindian populations is unknown. Methods. Serum samples were obtained from 517 Amerindians and tested for HHV-8 anti-latent nuclear antigen (anti-LANA) and antilytic antibodies by immunofluorescence assays. Logistic regression and mixed logistic models were used to estimate the odds of being HHV-8 seropositive among intrafamilial pairs. Results. HHV-8 seroprevalence by either assay was 75.4% (95% confidence interval [CI]: 71.5%-79.1%), and it was age-dependent (P-trend<.001). Familial dependence in HHV-8 seroprevalence by either assay was found between mother-offspring (odds ratio [OR], 5.44; 95% CI: 1.62-18.28) and siblings aged >= 10 years (OR 4.42, 95% CI: 1.70-11.45) or siblings in close age range (<5 years difference) (OR 3.37, 95% CI: 1.21-9.40), or in families with large (>4) number of siblings (OR, 3.20, 95% CI: 1.33-7.67). In separate analyses by serological assay, there was strong dependence in mother-offspring (OR 8.94, 95% CI: 2.94-27.23) and sibling pairs aged >= 10 years (OR, 11.91, 95% CI: 2.23-63.64) measured by LANA but not lytic antibodies. Conclusions. This pattern of familial dependence suggests that, in this endemic population, HHV-8 transmission mainly occurs from mother to offspring and between close siblings during early childhood, probably via saliva. The mother to offspring dependence was derived chiefly from anti-LANA antibodies.

Coupled Boltzmann calculation of mixed axion/neutralino cold dark matter production in the early universe

We calculate the relic abundance of mixed axion/neutralino cold dark matter which arises in R-parity conserving supersymmetric (SUSY) models wherein the strong CP problem is solved by the Peccei-Quinn (PQ) mechanism with a concommitant axion/saxion/axino supermultiplet. By numerically solving the coupled Boltzmann equations, we include the combined effects of 1. thermal axino production with cascade decays to a neutralino LSP, 2. thermal saxion production and production via coherent oscillations along with cascade decays and entropy injection, 3. thermal neutralino production and re-annihilation after both axino and saxion decays, 4. gravitino production and decay and 5. axion production both thermally and via oscillations. For SUSY models with too high a standard neutralino thermal abundance, we find the combined effect of SUSY PQ particles is not enough to lower the neutralino abundance down to its measured value, while at the same time respecting bounds on late-decaying neutral particles from BBN. However, models with a standard neutralino underabundance can now be allowed with either neutralino or axion domination of dark matter, and furthermore, these models can allow the PQ breaking scale f(a) to be pushed up into the 10(14) - 10(15) GeV range, which is where it is typically expected to be in string theory models.

Cation Distribution and Local Order in Mixed Sodium Metaphosphate Glasses

Two structural properties in mixed alkali metal phosphate glasses that seem to be crucial to the development of the mixed ion effect in dc conductivity were systematically analyzed in Na mixed metaphosphates: the local order around the mobile species, and their distribution and mixing in the glass network. The set of glasses considered here, Na1-xMxPO3 with M = Li, Ag, K, Rb, and Cs and 0 <= x <= 1, encompass a broad degree of size mismatch between the mixed cation species. A comprehensive solid-state nuclear magnetic resonance study was carried out using P-31 MAS, Na-23 triple quantum MAS, Rb-87 QCPMG, P-31-Na-23 REDOR, Na-23-Li-7 and Li-7-Li-6 SEDOR, and Na-23 spin echo decay. It was observed that the arrangement of P atoms around Na in the mixed glasses was indistinguishable from that observed in the NaPO3 glass. However, systematic distortions in the local structure of the 0 environments around Na were observed, related to the presence of the second cation. The average Na-O distances show an expansion/compression When Na+ ions are replaced by cations with respectively smaller/bigger radii. The behavior of the nuclear electric quadrupole coupling. constants indicates that this expansion reduces the local symmetry, while the compression produces the opposite effect These effects become marginally small when the site mismatch between the cations is small, as in Na-Ag mixed glasses. The present study confirms the intimate mixing of cation species at the atomic scale, but clear deviations from random mixing were detected in systems with larger alkali metal ions (Cs-Na, K-Na, Rb-Na). In contrast, no deviations from the statistical ion mixture were found in the systems Ag-Na and Li-Na, where mixed cations are either of radii comparable to (Ag+) or smaller than (Li+) Na+. The set of results supports two fundamental structural features of the models proposed to explain the mixed ion effect: the. structural specificity of the sites occupied by each cation species and their mixing at the atomic scale.

Assessment of variance components in nonlinear mixed-effects elliptical models

The issue of assessing variance components is essential in deciding on the inclusion of random effects in the context of mixed models. In this work we discuss this problem by supposing nonlinear elliptical models for correlated data by using the score-type test proposed in Silvapulle and Silvapulle (1995). Being asymptotically equivalent to the likelihood ratio test and only requiring the estimation under the null hypothesis, this test provides a fairly easy computable alternative for assessing one-sided hypotheses in the context of the marginal model. Taking into account the possible non-normal distribution, we assume that the joint distribution of the response variable and the random effects lies in the elliptical class, which includes light-tailed and heavy-tailed distributions such as Student-t, power exponential, logistic, generalized Student-t, generalized logistic, contaminated normal, and the normal itself, among others. We compare the sensitivity of the score-type test under normal, Student-t and power exponential models for the kinetics data set discussed in Vonesh and Carter (1992) and fitted using the model presented in Russo et al. (2009). Also, a simulation study is performed to analyze the consequences of the kurtosis misspecification.

Ni catalyst on mixed support of CeO2-ZrO2 and Al2O3: Effect of composition of CeO2-ZrO2 solid solution on the methane steam reforming reaction

In this study, catalysts containing 5 wt.% Ni deposited on a support composed of a CeO2-ZrO2 solid solution deposited on alumina were tested in the steam reforming of methane. The supports, with various ratios of Ce to Zr, were prepared by co-precipitation of the oxide precursors, followed by calcination in synthetic air. The catalysts were then prepared by Ni impregnation of the supports. The prepared solids were characterized by temperature-programmed reduction with H-2 (TPR-H-2), in situ X-ray diffraction (XRD) and X-ray absorption near-edge structure (XANES) spectroscopy. The XRD analysis confirmed the formation of a solid solution between ZrO2 and CeO2. In the catalytic tests, it was found that catalysts with higher Ce content did not exhibit deactivation during 6 h of reaction. The catalyst with highest Ce content, Ni(0.8Ce0.2Zr)AI, provided the best result, with the highest rate of conversion of methane and the lowest carbon deposition, which may be partly due to the smaller Ni-0 crystallites in this sample and also the segregated CeO2 particles may have favored H2O adsorption which could lead to higher C gasification. (C) 2012 Elsevier B.V. All rights reserved.