Skew-normal linear calibration: a Bayesian perspective

In this paper, we present a Bayesian approach for estimation in the skew-normal calibration model, as well as the conditional posterior distributions which are useful for implementing the Gibbs sampler. Data transformation is thus avoided by using the methodology proposed. Model fitting is implemented by proposing the asymmetric deviance information criterion, ADIC, a modification of the ordinary DIC. We also report an application of the model studied by using a real data set, related to the relationship between the resistance and the elasticity of a sample of concrete beams. Copyright (C) 2008 John Wiley & Sons, Ltd.

KittyCat: a cognitive model of structure-form discovery

Cognition is a core subject to understand how humans think and behave. In that sense, it is clear that Cognition is a great ally to Management, as the later deals with people and is very interested in how they behave, think, and make decisions. However, even though Cognition shows great promise as a field, there are still many topics to be explored and learned in this fairly new area. Kemp & Tenembaum (2008) tried to a model graph-structure problem in which, given a dataset, the best underlying structure and form would emerge from said dataset by using bayesian probabilistic inferences. This work is very interesting because it addresses a key cognition problem: learning. According to the authors, analogous insights and discoveries, understanding the relationships of elements and how they are organized, play a very important part in cognitive development. That is, this are very basic phenomena that allow learning. Human beings minds do not function as computer that uses bayesian probabilistic inferences. People seem to think differently. Thus, we present a cognitively inspired method, KittyCat, based on FARG computer models (like Copycat and Numbo), to solve the proposed problem of discovery the underlying structural-form of a dataset.

Bayesian analysis of extreme events with threshold estimation

The aim of this paper is to analyze extremal events using Generalized Pareto Distributions (GPD), considering explicitly the uncertainty about the threshold. Current practice empirically determines this quantity and proceeds by estimating the GPD parameters based on data beyond it, discarding all the information available be10w the threshold. We introduce a mixture model that combines a parametric form for the center and a GPD for the tail of the distributions and uses all observations for inference about the unknown parameters from both distributions, the threshold inc1uded. Prior distribution for the parameters are indirectly obtained through experts quantiles elicitation. Posterior inference is available through Markov Chain Monte Carlo (MCMC) methods. Simulations are carried out in order to analyze the performance of our proposed mode1 under a wide range of scenarios. Those scenarios approximate realistic situations found in the literature. We also apply the proposed model to a real dataset, Nasdaq 100, an index of the financiai market that presents many extreme events. Important issues such as predictive analysis and model selection are considered along with possible modeling extensions.

Genetic parameter estimates for buffalo milk yield, milk quality and mozzarella production and Bayesian inference analysis of their relationships

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Bayesian inference in genetic parameter estimation of visual scores in Nellore beef-cattle

The aim of this study was to estimate the components of variance and genetic parameters for the visual scores which constitute the Morphological Evaluation System (MES), such as body structure (S), precocity (P) and musculature (M) in Nellore beef-cattle at the weaning and yearling stages, by using threshold Bayesian models. The information used for this was gleaned from visual scores of 5,407 animals evaluated at the weaning and 2,649 at the yearling stages. The genetic parameters for visual score traits were estimated through two-trait analysis, using the threshold animal model, with Bayesian statistics methodology and MTGSAM (Multiple Trait Gibbs Sampler for Animal Models) threshold software. Heritability estimates for S, P and M were 0.68, 0.65 and 0.62 (at weaning) and 0.44, 0.38 and 0.32 (at the yearling stage), respectively. Heritability estimates for S, P and M were found to be high, and so it is expected that these traits should respond favorably to direct selection. The visual scores evaluated at the weaning and yearling stages might be used in the composition of new selection indexes, as they presented sufficient genetic variability to promote genetic progress in such morphological traits.

Estimation of genetic parameters for milk yield in Murrah buffaloes by Bayesian inference

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

Random regression analyses using B-splines functions to model growth from birth to adult age in Canchim cattle

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Genetic parameters for buffalo milk yield and milk quality traits using Bayesian inference

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

Genetic correlations between heifer subsequent rebreeding and age at first calving and growth traits in Nellore cattle by Bayesian inference

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

Comparing three indices of catch per unit effort using Bayesian geostatistics

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Comparison between the complete Bayesian method and empirical Bayesian method for ARCH models using Brazilian financial time series

In this work we compared the estimates of the parameters of ARCH models using a complete Bayesian method and an empirical Bayesian method in which we adopted a non-informative prior distribution and informative prior distribution, respectively. We also considered a reparameterization of those models in order to map the space of the parameters into real space. This procedure permits choosing prior normal distributions for the transformed parameters. The posterior summaries were obtained using Monte Carlo Markov chain methods (MCMC). The methodology was evaluated by considering the Telebras series from the Brazilian financial market. The results show that the two methods are able to adjust ARCH models with different numbers of parameters. The empirical Bayesian method provided a more parsimonious model to the data and better adjustment than the complete Bayesian method.

Analysis of beef cattle longitudinal data applying a nonlinear model

The objective of this work was to evaluate the Nelore beef cattle, growth curve parameters using the Von Bertalanffy function in a nested Bayesian procedure that allowed estimation of the joint posterior distribution of growth curve parameters, their (co)variance components, and the environmental and additive genetic components affecting them. A hierarchical model was applied; each individual had a growth trajectory described by the nonlinear function, and each parameter of this function was considered to be affected by genetic and environmental effects that were described by an animal model. Random samples of the posterior distributions were drawn using Gibbs sampling and Metropolis-Hastings algorithms. The data set consisted of a total of 145,961 BW recorded from 15,386 animals. Even though the curve parameters were estimated for animals with few records, given that the information from related animals and the structure of systematic effects were considered in the curve fitting, all mature BW predicted were suitable. A large additive genetic variance for mature BW was observed. The parameter a of growth curves, which represents asymptotic adult BW, could be used as a selection criterion to control increases in adult BW when selecting for growth rate. The effect of maternal environment on growth was carried through to maturity and should be considered when evaluating adult BW. Other growth curve parameters showed small additive genetic and maternal effects. Mature BW and parameter k, related to the slope of the curve, presented a large, positive genetic correlation. The results indicated that selection for growth rate would increase adult BW without substantially changing the shape of the growth curve. Selection to change the slope of the growth curve without modifying adult BW would be inefficient because their genetic correlation is large. However, adult BW could be considered in a selection index with its corresponding economic weight to improve the overall efficiency of beef cattle production.

Study of stayability in Nellore cows using a threshold model

The objectives of the current study were to assess the feasibility of using stayability traits to improve fertility of Nellore cows and to examine the genetic relationship among the stayabilities at different ages. Stayability was defined as whether a cow calved every year up to the age of 5 (Stay5), 6 (Stay6), or 7 (Stay7) yr of age or more, given that she was provided the opportunity to breed. Data were analyzed based on a maximum a posteriori probit threshold model to predict breeding values on the liability scale, whereas the Gibbs sampler was used to estimate variance components. The EBV were obtained using all animals included in the pedigree or bulls with at least 10 daughters with stayability observations, and average genetic trends were obtained in the liability and transformed to the probability scale. Additional analyses were performed to study the genetic relationship among stayability traits, which were compared by contrasting results in terms of EBV and the average genetic superiority as a function of the selected proportion of sires. Heritability estimates and SD were 0.25 +/- 0.02, 0.22 +/- 0.03, and 0.28 +/- 0.03 for Stay5, Stay6, and Stay7, respectively. Average genetic trends, by year, were 0.51 +/- 0.34, and 0.38% for Stay5, Stay6, and Stay7, respectively. Estimates of EBV SD, in the probability scale, for all animals included in the pedigree and for bulls with at least 10 daughters with stayability observations were 7.98 and 12.95, 6.93 and 11.38, and 8.24 and 14.30% for Stay5, Stay6, and Stay7, respectively. A reduction in the average genetic superiorities in Stay7 would be expected if the selection were based on Stay5 or Stay6. Nonetheless, the reduction in EPD, depending on selection intensity, is on average 0.74 and 1.55%, respectively. Regressions of the sires' EBV for Stay5 and Stay6 on the sires' EBV for Stay7 confirmed these results. The heritability and genetic trend estimates for all stayability traits indicate that it is possible to improve fertility with selection based on a threshold analysis of stayability. The SD of EBV for stayability traits show that there is adequate genetic variability among animals to justify inclusion of stayability as a selection criterion. The potential linear relationship among stayability traits indicates that selection for improved female traits would be more effective by having predictions on the Stay5 trait.

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.

Quantitative estimates of time-averaging in terebratulid brachiopod shell accumulations from a modern tropical shelf

Quantitative estimates of time-averaging in marine shell accumulations available to date are limited primarily to aragonitic mollusk shells. We assessed time-averaging in Holocene assemblages of calcitic brachiopod shells by direct dating of individual specimens of the terebratulid brachiopod Bouchardia rosea. The data were collected from exceptional (brachiopod-rich) shell assemblages, occurring surficially on a tropical mixed carbonate-siliciclastic shelf (the Southeast Brazilian Bight, SW Atlantic), a setting that provides a good climatic and environmental analog for many Paleozoic brachiopod shell beds of North America and Europe. A total of 82 individual brachiopod shells, collected from four shallow (5-25 m) nearshore (<2.5 km from the shore) localities, were dated by using amino acid racemization (D-alloisoleucine/L-isoleucine value) calibrated with five AMS-radiocarbon dates (r(2) = 0.933). This is the first study to demonstrate that amino acid racemization methods can provide accurate and precise ages for individual shells of calcitic brachiopods.The dated shells vary in age from modern to 3000 years, with a standard deviation of 690 years. The age distribution is strongly right-skewed: the young shells dominate the dated specimens and older shells are increasingly less common. However, the four localities display significant differences in the range of time-averaging and the form of the age distribution. The dated shells vary notably in the quality of preservation, but there is no significant correlation between taphonomic condition and age, either for individual shells or at assemblage level.These results demonstrate that fossil brachiopods may show considerable time-averaging, but the scale and nature of that mixing may vary greatly among sites. Moreover, taphonomic condition is not a reliable indicator of pre-burial history of individual brachiopod shells or the scale of temporal mixing within the entire assemblage. The results obtained for brachiopods are strikingly similar to results previously documented for mollusks and suggest that differences in mineralogy and shell microstructure are unlikely to be the primary factors controlling the nature and scale of time-averaging. Environmental factors and local fluctuations in populations of shell-producing organisms are more likely to be the principal determinants of time-averaging in marine benthic shelly assemblages. The long-term survival of brachiopod shells is incongruent with the rapid shell destruction observed in taphonomic experiments. The results support the taphonomic model that shells remain protected below (but perhaps near) the surface through their early taphonomic history. They may be brought back up to the surface intermittently by bioturbation and physical reworking, but only for short periods of time. This model explains the striking similarities in time-averaging among different types of organisms and the lack of correlation between time-since-death and shell taphonomy.