Generalized linear mixed models for the genetic evaluation of binary reproductive traits: a simulation study

The objective of this study was to evaluate the use of probit and logit link functions for the genetic evaluation of early pregnancy using simulated data. The following simulation/analysis structures were constructed: logit/logit, logit/probit, probit/logit, and probit/probit. The percentages of precocious females were 5, 10, 15, 20, 25 and 30% and were adjusted based on a change in the mean of the latent variable. The parametric heritability (h²) was 0.40. Simulation and genetic evaluation were implemented in the R software. Heritability estimates (ĥ²) were compared with h² using the mean squared error. Pearson correlations between predicted and true breeding values and the percentage of coincidence between true and predicted ranking, considering the 10% of bulls with the highest breeding values (TOP10) were calculated. The mean ĥ² values were under- and overestimated for all percentages of precocious females when logit/probit and probit/logit models used. In addition, the mean squared errors of these models were high when compared with those obtained with the probit/probit and logit/logit models. Considering ĥ², probit/probit and logit/logit were also superior to logit/probit and probit/logit, providing values close to the parametric heritability. Logit/probit and probit/logit presented low Pearson correlations, whereas the correlations obtained with probit/probit and logit/logit ranged from moderate to high. With respect to the TOP10 bulls, logit/probit and probit/logit presented much lower percentages than probit/probit and logit/logit. The genetic parameter estimates and predictions of breeding values of the animals obtained with the logit/logit and probit/probit models were similar. In contrast, the results obtained with probit/logit and logit/probit were not satisfactory. There is need to compare the estimation and prediction ability of logit and probit link functions.

Robust linear mixed models with normal/independent distributions and Bayesian MCMC implementation

Linear mixed effects models have been widely used in analysis of data where responses are clustered around some random effects, so it is not reasonable to assume independence between observations in the same cluster. In most biological applications, it is assumed that the distributions of the random effects and of the residuals are Gaussian. This makes inferences vulnerable to the presence of outliers. Here, linear mixed effects models with normal/independent residual distributions for robust inferences are described. Specific distributions examined include univariate and multivariate versions of the Student-t, the slash and the contaminated normal. A Bayesian framework is adopted and Markov chain Monte Carlo is used to carry out the posterior analysis. The procedures are illustrated using birth weight data on rats in a texicological experiment. Results from the Gaussian and robust models are contrasted, and it is shown how the implementation can be used for outlier detection. The thick-tailed distributions provide an appealing robust alternative to the Gaussian process in linear mixed models, and they are easily implemented using data augmentation and MCMC techniques.

Bayesian longitudinal data analysis with mixed models and thick-tailed distributions using MCMC

Linear mixed effects models are frequently used to analyse longitudinal data, due to their flexibility in modelling the covariance structure between and within observations. Further, it is easy to deal with unbalanced data, either with respect to the number of observations per subject or per time period, and with varying time intervals between observations. In most applications of mixed models to biological sciences, a normal distribution is assumed both for the random effects and for the residuals. This, however, makes inferences vulnerable to the presence of outliers. Here, linear mixed models employing thick-tailed distributions for robust inferences in longitudinal data analysis are described. Specific distributions discussed include the Student-t, the slash and the contaminated normal. A Bayesian framework is adopted, and the Gibbs sampler and the Metropolis-Hastings algorithms are used to carry out the posterior analyses. An example with data on orthodontic distance growth in children is discussed to illustrate the methodology. Analyses based on either the Student-t distribution or on the usual Gaussian assumption are contrasted. The thick-tailed distributions provide an appealing robust alternative to the Gaussian process for modelling distributions of the random effects and of residuals in linear mixed models, and the MCMC implementation allows the computations to be performed in a flexible manner.

An experimental assessment of landscape configuration effects on frog and toad abundance and diversity in tropical agro-savannah landscapes of southeastern Brazil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Genetic and non-genetic effects on calf vigor at birth and preweaning mortality in nellore calves

The objective of the present study was to evaluate the genetic and non-genetic effects that influencevigor at birth and preweaning mortality in Nellore calves. A total of 11,727 records of births that occurred between 1978 and 2006, offspring of 363 sires, were analyzed. Poor calf vigor at birth (VB) and preweaning mortality divided into stillbirth (SB), early mortality (EM) and total mortality (TM) were analyzed as binary variables. Generalized linear models were used for the evaluation of non-genetic effects and generalized linear mixed models for genetic effects (sire and animal models). The incidences were 4.75% for VB, 2.66% for SB, 5.28% for EM, and 7.99% for TM. Birth weight was the effect that most influenced the traits studied. Calves weighing less than 22kg(females) and less than 24kg (males) were at a higher risk of low vigor and preweaning mortality. Preweaning mortality was higher among calves born from cows aged .3 and .11 years at calving compared with cows aged 7 to 10 years. Male calves presented less vigor and higher preweaning mortality than female calves. Selection for postweaning weight did not influence preweaning mortality. The heritability estimates ranged between 0.01 and 0.09 for VB, 0.00 and 0.27 for SB, 0.03 and 0.17 for EM and 0.02 and 0.10 for TM. Stillbirth should be included as a selection criterion in breeding programs of Nellore cattle, alone or as part of a selection index, aiming to reduce preweaning mortality. © 2013 Sociedade Brasileira de Zootecnia.

Metodologias estatísticas na análise de germinação de sementes de mamona

Pós-graduação em Agronomia (Energia na Agricultura) - FCA

Variabilidades genética e não genética na mortalidade pré-desmama de bezerros de corte

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Environmental factors influencing the distribution of three species within the genus Persephona Leach, 1817 (Crustacea, Decapoda, Leucosiidae) in two regions on the northern coast of Sao Paulo State, Brazil

Patterns of spatio-temporal distribution of Brachyura are determined by the interaction among life history traits, inter and intraspecific relationships, as well as by the variation of abiotic factors. This study aimed to characterize patterns of spatio-temporal distribution of Persephona lichtensteinii, Persephona mediterranea and Persephona punctata in two regions of the northern coast of Sao Paulo State, southeastern region of Brazil. Collections were done monthly from July 2001 to June 2003 in Caraguatatuba and Ubatuba, using a shrimp fishery boat equipped with double-rig nets. The patterns of species distribution were tested by means of redundancy analysis (RDA) and generalized linear mixed models (GLMM) in relation to the recorded environmental factors (BT: bottom temperature, BS: bottom salinity, OM: organic matter and granulometry (Phi)). The most influent environmental factor over the species distribution was the Phi, and the ascendant order of influence was P. lichtensteinii, P. punctata and P. mediterranea. The greater abundance of P. mediterranea showed a conservative pattern of distribution for the genus in the sampled region. The greater occurrence of P. punctata and P. lichtensteinii, in distinct transects than those occupied by P. mediterranea, seems to be a strategy to avoid competition among congeneric species, which is related to the substratum specificity.

The Conway-Maxwell-Poisson-generalized gamma regression model with long-term survivors

In this paper, we proposed a flexible cure rate survival model by assuming the number of competing causes of the event of interest following the Conway-Maxwell distribution and the time for the event to follow the generalized gamma distribution. This distribution can be used to model survival data when the hazard rate function is increasing, decreasing, bathtub and unimodal-shaped including some distributions commonly used in lifetime analysis as particular cases. Some appropriate matrices are derived in order to evaluate local influence on the estimates of the parameters by considering different perturbations, and some global influence measurements are also investigated. Finally, data set from the medical area is analysed.

Generalized neighbor-interaction models induced by nonlinear lattices

It is shown that the tight-binding approximation of the nonlinear Schrodinger equation with a periodic linear potential and periodic in space nonlinearity coefficient gives rise to a number of nonlinear lattices with complex, both linear and nonlinear, neighbor interactions. The obtained lattices present nonstandard possibilities, among which we mention a quasilinear regime, where the pulse dynamics obeys essentially the linear Schrodinger equation. We analyze the properties of such models both in connection to their modulational stability, as well as in regard to the existence and stability of their localized solitary wave solutions.

Comparison of non-linear growth models to describe the growth curve in West African Dwarf sheep

The objectives of this study were to compare the goodness of fit of four non-linear growth models, i.e. Brody, Gompertz, Logistic and Von Bertalanffy, in West African Dwarf (WAD) sheep. A total of 5274 monthly weight records from birth up to 180 days of age from 889 lambs, collected during 2001 to 2004 in Betecoucou breeding farm in Benin were used. In the preliminary analysis, the General Linear Model Procedure of the Statistical Analysis Systems Institute was applied to the dataset to identify the significant effects of the sex of lamb (male and female), type of birth (single and twin), season of birth (rainy season and dry season), parity of dam (1, 2 and 3) and year of birth (2001, 2002, 2003 and 2004) on the observed birth weight and monthly weight up to 6 months of age. The models parameters (A, B and k), coefficient of determination (112), mean square error (MSE) were calculated using language of technical computing package Matlab(R), 2006. The mean values of A, B and k were substituted into each model to calculate the corresponding Akaike's Information Criterion (AIC). Among the four growth functions, the Brody model has been selected for its accuracy of fit according to the higher R(2), lower MSE and A/C Finally, the parameters A, B and k were adjusted in Matlab(R) 2006 for the sex of lamb, year of birth, season of birth, birth type and the parity of ewe, providing a specific slope of the Brody growth curve. The results of this study suggest that Brody model can be useful for WAD sheep breeding in Betecoucou farm conditions through growth monitoring.

Two field BPS solutions for generalized Lorentz breaking models

In this work we present nonlinear models in two-dimensional space-time of two interacting scalar fields in the Lorentz and CPT violating scenarios. We discuss the soliton solutions for these models as well as the question of stability for them. This is done by generalizing a model recently published by Barreto and collaborators and also by getting new solutions for the model introduced by them.

Non-linear growth models for bullfrog tadpoles

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Composite interval mapping and mixed models reveal QTL associated with performance and carcass traits on chicken chromosomes 1, 3, and 4

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)