Bayesian analysis of random regression models using B-splines to model test-day milk yield of Holstein cattle in Brazil

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

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)

Using cubic splines to mitigate systematic errors in GPS relative positioning

Systematic errors can have a significant effect on GPS observable. In medium and long baselines the major systematic error source are the ionosphere and troposphere refraction and the GPS satellites orbit errors. But, in short baselines, the multipath is more relevant. These errors degrade the accuracy of the positioning accomplished by GPS. So, this is a critical problem for high precision GPS positioning applications. Recently, a method has been suggested to mitigate these errors: the semiparametric model and the penalised least squares technique. It uses a natural cubic spline to model the errors as a function which varies smoothly in time. The systematic errors functions, ambiguities and station coordinates, are estimated simultaneously. As a result, the ambiguities and the station coordinates are estimated with better reliability and accuracy than the conventional least square method.

Random regression models using Legendre polynomials or linear splines for test-day milk yield of dairy Gyr (Bos indicus) cattle

Studies investigating the use of random regression models for genetic evaluation of milk production in Zebu cattle are scarce. In this study, 59,744 test-day milk yield records from 7,810 first lactations of purebred dairy Gyr (Bos indicus) and crossbred (dairy Gyr × Holstein) cows were used to compare random regression models in which additive genetic and permanent environmental effects were modeled using orthogonal Legendre polynomials or linear spline functions. Residual variances were modeled considering 1, 5, or 10 classes of days in milk. Five classes fitted the changes in residual variances over the lactation adequately and were used for model comparison. The model that fitted linear spline functions with 6 knots provided the lowest sum of residual variances across lactation. On the other hand, according to the deviance information criterion (DIC) and Bayesian information criterion (BIC), a model using third-order and fourth-order Legendre polynomials for additive genetic and permanent environmental effects, respectively, provided the best fit. However, the high rank correlation (0.998) between this model and that applying third-order Legendre polynomials for additive genetic and permanent environmental effects, indicates that, in practice, the same bulls would be selected by both models. The last model, which is less parameterized, is a parsimonious option for fitting dairy Gyr breed test-day milk yield records. © 2013 American Dairy Science Association.

O uso de Thin-Plate Splines na transformação de coordenadas com modelagem de distorções entre realizações de referenciais geodésicos

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

Transformação de coordenadas com modelagem de distorções entre SAD69 e SIRGAS2000 com o uso de Thin-Plate Splines

Several countries have been passed by change processes in their fundamental geodesic structure with the focus on the adoption of geocentric reference systems. In Brazil, the adoption of the SIRGAS2000 evolves the coexistence of two realizations from the COrrego Alegre system, two realizations from the SAD69 system and one realization from the SIRGAS2000 system. To make use of products in the old reference systems, methods of coordinate transformation between the existent reference frames are necessary. So, in this paper one solution for the transformation between coordinates from different reference frames, based on Thin-Plate Splines (TPS), that allows the estimation of parameters from one linear transformation and also one non-linear model is presented. The TPS model was developed to work with tridimensional coordinates and in this paper the results and analysis are performed with simulated data and also with data from the official Brazilian Geodetic System (SGB). In the check points from SAD69 stations (realization of 1996 - SAD69/96), the values of RMSE obtained were of 78,2 mm in latitude and 67,5 mm in longitude, before the transformation to the SIRGAS2000. In the comparison between the TPS model and ProGriD (Brazilian software provided by IBGE), the statistical indicators were reduced in 97%, by using the TPS model. Based in the obtained results from real dataset, the TPS model appears to be promising, since it allows improving the quality of transformation process with simultaneous distortion modeling.

Random regression analyses using B-spline functions to model growth of Nellore cattle

The objective of this study was to estimate (co)variance components using random regression on B-spline functions to weight records obtained from birth to adulthood. A total of 82 064 weight records of 8145 females obtained from the data bank of the Nellore Breeding Program (PMGRN/Nellore Brazil) which started in 1987, were used. The models included direct additive and maternal genetic effects and animal and maternal permanent environmental effects as random. Contemporary group and dam age at calving (linear and quadratic effect) were included as fixed effects, and orthogonal Legendre polynomials of age (cubic regression) were considered as random covariate. The random effects were modeled using B-spline functions considering linear, quadratic and cubic polynomials for each individual segment. Residual variances were grouped in five age classes. Direct additive genetic and animal permanent environmental effects were modeled using up to seven knots (six segments). A single segment with two knots at the end points of the curve was used for the estimation of maternal genetic and maternal permanent environmental effects. A total of 15 models were studied, with the number of parameters ranging from 17 to 81. The models that used B-splines were compared with multi-trait analyses with nine weight traits and to a random regression model that used orthogonal Legendre polynomials. A model fitting quadratic B-splines, with four knots or three segments for direct additive genetic effect and animal permanent environmental effect and two knots for maternal additive genetic effect and maternal permanent environmental effect, was the most appropriate and parsimonious model to describe the covariance structure of the data. Selection for higher weight, such as at young ages, should be performed taking into account an increase in mature cow weight. Particularly, this is important in most of Nellore beef cattle production systems, where the cow herd is maintained on range conditions. There is limited modification of the growth curve of Nellore cattle with respect to the aim of selecting them for rapid growth at young ages while maintaining constant adult weight.

Spatial variability of soil penetration resistance influenced by season of sampling

O objetivo deste trabalho foi analizar a distribuição espacial da compactação do solo e a influência da umidade do solo na resistência à penetração. Esta última variável foi descrita pelo índice de cone. O solo estudado foi Nitossolo e os dados de índice de cone foram obtidos usando um penetrômetro. A resistência do solo foi avaliada a 5 profundidades diferentes, 0-10 cm, 10-20 cm, 20-30 cm, 30-40 cm e mais de 40 cm, porém o conteúdo de umidade do solo foi medido a 0-20 cm e 20-40 cm. As condições hídricas do solo variaram nas diferentes amostragems. Os coeficientes de variação para o índice de cone foram 16,5% a 45,8% e os do conteúdo de umidade do solo variaram entre 8,96% e 21,38%. Os resultados sugeriram elevada correlação entre a resistência do solo, estimada pelo índice de cone e a profundidade do solo. Sem embargo, a relação esperada com a umidade do solo não foi apreciada. Observou-se dependência espacial em 31 de 35 séries de dados de índice de cone e umidade do solo. Esta dependência foi ajustada por modelos exponenciais com efeito pepita variável de 0 a 90% o valor do patamar. em séries de dados o comportamento foi aleatório. Portanto, a técnica das distâncias inversas foi utilizada para cartografar a distribuição das variáveis que não tiveram estrutura espacial. Na krigagem constatou-se uma suavização dos mapas comparados com esses das distâncias inversas. A krigagem indicadora foi utilizada para cartografar a variabilidade espacial do índice de cone e recomendar melhor manejo do solo.

The strut-and-tie models in reinforced concrete structures analysed by a numerical technique

Os modelos de bielas e tirantes são procedimentos de análise apropriados para projetar elementos de concreto armado em casos de regiões onde há alterações geométricas ou concentrações de tensões, denominadas regiões D. Trata-se de bons modelos de representação da estrutura para avaliar melhor o seu comportamento estrutural e seu mecanismo resistente. O presente artigo aplica a técnica da otimização topológica para identificar o fluxo de tensões nas estruturas, definindo a configuração dos membros de bielas e tirantes, e quantifica seus valores para dimensionamento. Utilizam-se o método ESO, e uma variante desse, o SESO (Smoothing ESO) com o método dos elementos finitos em elasticidade plana. A filosofia do SESO baseia-se na observação de que se o elemento não for necessário à estrutura, sua contribuição de rigidez vai diminuindo progressivamente. Isto é, sua remoção é atenuada nos valores da matriz constitutiva, como se este estivesse em processo de danificação. Para validar a presente formulação, apresentam-se alguns exemplos numéricos onde se comparam suas respostas com as advindas de trabalhos científicos pioneiros sobre o assunto.

Performance Comparisons of EWMA Control Chart Schemes

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

How to not inflate population estimates? Spatial density distribution of white-lipped peccaries in a continuous Atlantic forest

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

Landmark-based morphometric analysis in selected species of serranid fishes (Perciformes : Teleostei)

Morphological differences among 6 species of marine fishes belonging to 2 subfamilies of the family Serranidae (Serraninae: Dules auriga, Diplectrum formosum, and D, radiale; Epinephelinae: Epinephelus marginatus, Mycteroperca acutirostris, and M. bonaci) were studied by the geometric morphometric method of thin-plate splines and multivariate analysis of partial-warp scores. The decomposition of shape variation into uniform and nonaffine components of shape change indicate that major differences among species are related to both components of shape variation. Significant differences were found among species with respect to the uniform components, but there is no clear separation of taxonomic groups related to these components, and species are instead separated on the basis of body height and caudal peduncle length. Non-uniform changes in body shape, in turn, clearly differentiate the species of Serraninae and Epinephelinae. These shape changes are probably related to differences in habitat and feeding habits among the species.

Caffeine effect on mortality and oviposition in successive generations of Aedes aegypti

OBJETIVO: Experimentos anteriores mostraram que a cafeína bloqueia o desenvolvimento de Aedes aegypti (Diptera, Culicidae) na fase larval, inibindo conseqüentemente a produção de adultos. O objetivo do estudo foi obter dados que pudessem sugerir desenvolvimento de resistência dos mosquitos à cafeína. MÉTODOS: Foi avaliada a produção de adultos em gerações sucessivas, a partir de ovos produzidos na geração anterior e a taxa de oviposição em cada geração, utilizando meios contendo cafeína a 200 e 500 µg/ml e água de torneira proveniente de poço artesiano como controle. Os experimentos foram conduzidos em São José do Rio Preto, entre 2002 e 2005. Nos testes estatísticos foram utilizados a análise exploratória de dados e algoritmos de alisamento. RESULTADOS: Ocorreu redução crescente da produção de adultos, nas duas concentrações, ao longo das gerações, mas apenas no experimento a 200 µg/ml os dados foram estatisticamente significantes. Quanto à oviposição, a análise dos números mostra redução crescente e acentuada na média de ovos por fêmea, no experimento tratado. CONCLUSÕES: Não houve evidência de resistência ao longo das gerações devido ao tratamento com cafeína. Os resultados encontrados podem reforçar a indicação da cafeína como uma alternativa aos principais agentes de controle do Ae. aegypti atualmente usados, contra os quais os mosquitos têm desenvolvido resistência.

The cross-entropy method and its application to minimize the ripple of magnetic levitation forces of a maglev system

This paper explores firstly the potential of a new evolutionary method - the Cross-Entropy (CE) method in solving continuous inverse electromagnetic problems. For this purpose, an adaptive updating formula for the smoothing parameter, some mutation operation, and a new termination criterion are proposed. The proposed CE based metaheuristics is applied to reduce the ripple of the magnetic levitation forces of a prototype Maglev system. The numerical results have shown that the ripple of the magnetic levitation forces of the prototype system is reduced significantly after the design optimization using the proposed algorithm.

SPH simulations of clumps formation by dissipative collision of molecular clouds - I. Non magnetic case

Computer experiments of interstellar cloud collisions were performed with a new smoothed-particle-hydrodynamics (SPH) code. The SPH quantities were calculated by using spatially adaptive smoothing lengths and the SPH fluid equations of motion were solved by means of a hierarchical multiple time-scale leapfrog. Such a combination of methods allows the code to deal with a large range of hydrodynamic quantities. A careful treatment of gas cooling by H, H(2), CO and H II, as well as a heating mechanism by cosmic rays and by H(2) production on grains surface, were also included in the code. The gas model reproduces approximately the typical environment of dark molecular clouds. The experiments were performed by impinging two dynamically identical spherical clouds onto each other with a relative velocity of 10 km s(-1) but with a different impact parameter for each case. Each object has an initial density profile obeying an r(-1)-law with a cutoff radius of 10 pc and with an initial temperature of 20 K. As a main result, cloud-cloud collision triggers fragmentation but in expense of a large amount of energy dissipated, which occurred in the head-on case only. Off-center collision did not allow remnants to fragment along the considered time (similar to 6 Myr). However, it dissipated a considerable amount of orbital energy. Structures as small as 0.1 pc, with densities of similar to 10(4) cm(-3), were observed in the more energetic collision.