Additive genetic relationships between scrotal circumference, heifer pregnancy, and stayability in Nellore cattle

Scrotal circumference data from 47,605 Nellore young bulls, measured at around 18 mo of age (SC18), were analyzed simultaneously with 27,924 heifer pregnancy (HP) and 80,831 stayability (STAY) records to estimate their additive genetic relationships. Additionally, the possibility that economically relevant traits measured directly in females could replace SC18 as a selection criterion was verified. Heifer pregnancy was defined as the observation that a heifer conceived and remained pregnant, which was assessed by rectal palpation at 60 d. Females were exposed to sires for the first time at about 14 mo of age (between 11 and 16 mo). Stayability was defined as whether or not a cow calved every year up to 5 yr of age, when the opportunity to breed was provided. A Bayesian linear-threshold-threshold analysis via Gibbs sampler was used to estimate the variance and covariance components of the multitrait model. Heritability estimates were 0.42 +/- 0.01, 0.53 +/- 0.03, and 0.10 +/- 0.01, for SC18, HP, and STAY, respectively. The genetic correlation estimates were 0.29 +/- 0.05, 0.19 +/- 0.05, and 0.64 +/- 0.07 between SC18 and HP, SC18 and STAY, and HP and STAY, respectively. The residual correlation estimate between HP and STAY was -0.08 +/- 0.03. The heritability values indicate the existence of considerable genetic variance for SC18 and HP traits. However, genetic correlations between SC18 and the female reproductive traits analyzed in the present study can only be considered moderate. The small residual correlation between HP and STAY suggests that environmental effects common to both traits are not major. The large heritability estimate for HP and the high genetic correlation between HP and STAY obtained in the present study confirm that EPD for HP can be used to select bulls for the production of precocious, fertile, and long-lived daughters. Moreover, SC18 could be incorporated in multitrait analysis to improve the prediction accuracy for HP genetic merit of young bulls.

An improved estimate of leaf area index based on the histogram analysis of hemispherical photographs

Leaf area index (LAI) is a key parameter that affects the surface fluxes of energy, mass, and momentum over vegetated lands, but observational measurements are scarce, especially in remote areas with complex canopy structure. In this paper we present an indirect method to calculate the LAI based on the analyses of histograms of hemispherical photographs. The optimal threshold value (OTV), the gray-level required to separate the background (sky) and the foreground (leaves), was analytically calculated using the entropy crossover method (Sahoo, P.K., Slaaf, D.W., Albert, T.A., 1997. Threshold selection using a minimal histogram entropy difference. Optical Engineering 36(7) 1976-1981). The OTV was used to calculate the LAI using the well-known gap fraction method. This methodology was tested in two different ecosystems, including Amazon forest and pasturelands in Brazil. In general, the error between observed and calculated LAI was similar to 6%. The methodology presented is suitable for the calculation of LAI since it is responsive to sky conditions, automatic, easy to implement, faster than commercially available software, and requires less data storage. (C) 2008 Elsevier B.V. All rights reserved.

Virtual Pole from Magnetic Anomaly (VPMA): A procedure to estimate the age of a rock from its magnetic anomaly only

Virtual Pole from Magnetic Anomaly (VPMA) is a new multi-disciplinary methodology that estimates the age of a source rock from its magnetic anomaly, taken directly from available aeromagnetic data. The idea is to use those anomalies in which a strong remanent magnetic component is likely to occur. Once the total magnetization of the anomaly is computed through any of the currently available methods, the line that connects all virtual paleogeographic poles is related with the position, on a paleogeographic projection, of the appropriate age fragment of the APWT curve. We applied this procedure to five (5) well-known magnetic anomalies of the South American plate in SE Brazil, all of them associated to alkaline complexes of Mesozoic age. The apparent ages obtained from VPMA on three of the anomalies where the radiometric age of the source rock is known - Tapira, Araxa and Juquia were inside the error interval of the published ages. The VPMA apparent ages of the other two, where the age of the source rock is not known (Registro and Pariqueracu magnetic anomalies) were geologically coherent. We expect that the application of the VPMA methodology as a reconnaissance geochronological tool may contribute to geological knowledge over continental areas, especially when the source rocks of the magnetic anomalies am unknown or buried below superficial sediments. (C) 2009 Elsevier B.V. All rights reserved.

Consistency in approximate entropy given by a volumetric estimate

Non-linear methods for estimating variability in time-series are currently of widespread use. Among such methods are approximate entropy (ApEn) and sample approximate entropy (SampEn). The applicability of ApEn and SampEn in analyzing data is evident and their use is increasing. However, consistency is a point of concern in these tools, i.e., the classification of the temporal organization of a data set might indicate a relative less ordered series in relation to another when the opposite is true. As highlighted by their proponents themselves, ApEn and SampEn might present incorrect results due to this lack of consistency. In this study, we present a method which gains consistency by using ApEn repeatedly in a wide range of combinations of window lengths and matching error tolerance. The tool is called volumetric approximate entropy, vApEn. We analyze nine artificially generated prototypical time-series with different degrees of temporal order (combinations of sine waves, logistic maps with different control parameter values, random noises). While ApEn/SampEn clearly fail to consistently identify the temporal order of the sequences, vApEn correctly do. In order to validate the tool we performed shuffled and surrogate data analysis. Statistical analysis confirmed the consistency of the method. (C) 2008 Elsevier Ltd. All rights reserved.

Covariance structure in the skull of Catarrhini: a case of pattern stasis and magnitude evolution

The study of the genetic variance/covariance matrix (G-matrix) is a recent and fruitful approach in evolutionary biology, providing a window of investigating for the evolution of complex characters. Although G-matrix studies were originally conducted for microevolutionary timescales, they could be extrapolated to macroevolution as long as the G-matrix remains relatively constant, or proportional, along the period of interest. A promising approach to investigating the constancy of G-matrices is to compare their phenotypic counterparts (P-matrices) in a large group of related species; if significant similarity is found among several taxa, it is very likely that the underlying G-matrices are also equivalent. Here we study the similarity of covariance and correlation structure in a broad sample of Old World monkeys and apes (Catarrhini). We made phylogenetically structured comparisons of correlation and covariance matrices derived from 39 skull traits, ranging from between species to the superfamily level. We also compared the overall magnitude of integration between skull traits (r(2)) for all Catarrhim genera. Our results show that P-matrices were not strictly constant among catarrhines, but the amount of divergence observed among taxa was generally low. There was significant and positive correlation between the amount of divergence in correlation and covariance patterns among the 30 genera and their phylogenetic distances derived from a recently proposed phylogenetic hypothesis. Our data demonstrate that the P-matrices remained relatively similar along the evolutionary history of catarrhines, and comparisons with the G-matrix available for a New World monkey genus (Saguinus) suggests that the same holds for all anthropoids. The magnitude of integration, in contrast, varied considerably among genera, indicating that evolution of the magnitude, rather than the pattern of inter-trait correlations, might have played an important role in the diversification of the catarrhine skull. (C) 2009 Elsevier Ltd. All rights reserved.

First molecular estimate of sex-ratio of southern right whale calves, Eubalaena australis, for Brazilian waters

The southern right whale (Eubalaena australis) was one of the most intensively hunted whales between the 17th and 20th centuries in the southern hemisphere. Recent estimates indicate that today there are around 7000 whales, representing 5 to 10% Of its original population. On the other hand, recent studies estimated that the population that migrates to the Brazilian coast grew by 14% from 1987 to 2003. However, there is no information about sex-ratio for adults or for calves in this region, which is an important parameter for understanding the biology of the species. We present here the first estimate Of calves` sex-ratio of southern right whales found along the southern Brazilian coast, one of the most important wintering grounds for the species. Sex was molecularly indentified for 21 biopsies collected from calves between 1998 and 2002, along the coast of Rio Grande do Sul and Santa Catarina States, in southern Brazil. The sex-ratio was two females for one male, however, it was not statistically different (chi(2) test, alpha = 0.05; df = 1) from the expected ratio of 1:1. This result is in accordance with the sex-ratio estimated for the species of all ages using external morphology (and behaviour in formation), (is well as for most species of baleen whales.

Using non-homogeneous Poisson models with multiple change-points to estimate the number of ozone exceedances in Mexico City

In this paper, we consider some non-homogeneous Poisson models to estimate the probability that an air quality standard is exceeded a given number of times in a time interval of interest. We assume that the number of exceedances occurs according to a non-homogeneous Poisson process (NHPP). This Poisson process has rate function lambda(t), t >= 0, which depends on some parameters that must be estimated. We take into account two cases of rate functions: the Weibull and the Goel-Okumoto. We consider models with and without change-points. When the presence of change-points is assumed, we may have the presence of either one, two or three change-points, depending of the data set. The parameters of the rate functions are estimated using a Gibbs sampling algorithm. Results are applied to ozone data provided by the Mexico City monitoring network. In a first instance, we assume that there are no change-points present. Depending on the adjustment of the model, we assume the presence of either one, two or three change-points. Copyright (C) 2009 John Wiley & Sons, Ltd.

Estimation or simulation? That is the question

The issue of smoothing in kriging has been addressed either by estimation or simulation. The solution via estimation calls for postprocessing kriging estimates in order to correct the smoothing effect. Stochastic simulation provides equiprobable images presenting no smoothing and reproducing the covariance model. Consequently, these images reproduce both the sample histogram and the sample semivariogram. However, there is still a problem, which is the lack of local accuracy of simulated images. In this paper, a postprocessing algorithm for correcting the smoothing effect of ordinary kriging estimates is compared with sequential Gaussian simulation realizations. Based on samples drawn from exhaustive data sets, the postprocessing algorithm is shown to be superior to any individual simulation realization yet, at the expense of providing one deterministic estimate of the random function.

Covariance matrix formula for Birnbaum-Saunders regression models

The Birnbaum-Saunders regression model is commonly used in reliability studies. We derive a simple matrix formula for second-order covariances of maximum-likelihood estimators in this class of models. The formula is quite suitable for computer implementation, since it involves only simple operations on matrices and vectors. Some simulation results show that the second-order covariances can be quite pronounced in small to moderate sample sizes. We also present empirical applications.

In vitro basal cytotoxicity assay applied to estimate acute oral systemic toxicity of grandisin and its major metabolite

Preclinical investigations can start with preliminary in vitro studies before using animal models. Following this approach, the number of animals used in preclinical acute toxicity testing can be reduced. In this study, we employed an in-house validated in vitro cytotoxicity test based on the Spielmann approach for toxicity evaluation of the lignan grandisin, a candidate anticancer agent, and its major metabolite. the 4-O-demethylgrandisin, by neutral red uptake (NRU) assay, on mouse fibroblasts Balb/c 3T3 cell line. Using different concentrations of grandisin and its major metabolite (2.31; 1.16; 0.58; 0.29; 0.14; 0.07; 0.04; 0.002 mu M) in Balb/c 3T3-A31 NRU cytotoxicity assay, after incubation for 48 h, we obtained IC(50) values for grandisin and its metabolite of 0.078 and 0.043 mu M, respectively. The computed LD(50) of grandisin and 4-O-demethylgrandisin were 617.72 and 429.95 mg/kg, respectively. Both were classified under the Globally Harmonized System as category 4. Since pharmacological and toxicological data are crucial in the developmental stages of drug discovery, using an in vitro assay we demonstrated that grandisin and its metabolite exhibit distinct toxicity profiles. Furthermore, results presented in this work can contribute to reduce the number of animals required in subsequent pharmacological/toxicological studies. (C) 2010 Elsevier GmbH. All rights reserved.