Modelos de regressão aleatória com diferentes estruturas de variância residual para descrever o tamanho da leitegada

Objetivou-se comparar modelos de regressão aleatória com diferentes estruturas de variância residual, a fim de se buscar a melhor modelagem para a característica tamanho da leitegada ao nascer (TLN). Utilizaram-se 1.701 registros de TLN, que foram analisados por meio de modelo animal, unicaracterística, de regressão aleatória. As regressões fixa e aleatórias foram representadas por funções contínuas sobre a ordem de parto, ajustadas por polinômios ortogonais de Legendre de ordem 3. Para averiguar a melhor modelagem para a variância residual, considerou-se a heterogeneidade de variância por meio de 1 a 7 classes de variância residual. O modelo geral de análise incluiu grupo de contemporâneo como efeito fixo; os coeficientes de regressão fixa para modelar a trajetória média da população; os coeficientes de regressão aleatória do efeito genético aditivo-direto, do comum-de-leitegada e do de ambiente permanente de animal; e o efeito aleatório residual. O teste da razão de verossimilhança, o critério de informação de Akaike e o critério de informação bayesiano de Schwarz apontaram o modelo que considerou homogeneidade de variância como o que proporcionou melhor ajuste aos dados utilizados. As herdabilidades obtidas foram próximas a zero (0,002 a 0,006). O efeito de ambiente permanente foi crescente da 1ª (0,06) à 5ª (0,28) ordem, mas decrescente desse ponto até a 7ª ordem (0,18). O comum-de-leitegada apresentou valores baixos (0,01 a 0,02). A utilização de homogeneidade de variância residual foi mais adequada para modelar as variâncias associadas à característica tamanho da leitegada ao nascer nesse conjunto de dado.

Componentes de (co)variância e parâmetros genéticos de características de crescimento da raça Simental no Brasil

Informações de genealogia e produção, cedidas pela Associação Brasileira de Criadores da Raça Simental (ABCRS), relativas aos pesos desde o nascimento até um ano de idade, foram utilizadas para estimar, sob modelos alternativos, os componentes de variância e os parâmetros genéticos em animais da raça Simental no Brasil. A matriz de parentesco incluiu 25.812 animais dos quais 7587 com dados de produção. O modelo 1 contém, além do erro, o efeito genético direto. Os modelos seguintes contêm os componentes do modelo 1, mais o efeito permanente de ambiente materno (modelo 2), ou o componente genético materno (modelo 3), ambos os componentes (modelo 5), os componentes do modelo 3 mais a covariância entre os efeitos genéticos direto e materno (modelo 4) e todos os componentes citados (modelo 6). Os modelos foram comparados pelo teste de razão de verossimilhança pelo chi² (P<0,01). Os componentes de variância e os valores de herdabilidades, estimados para os efeitos direto e materno, foram decrescentes, desde o modelo 1 até o modelo 6, na razão direta em que o modelo incorpora mais efeitos aleatórios. Para a fase de aleitamento foi encontrada variância genética nula, entretanto, alto valor para a variância de ambiente permanente. Os efeitos maternos, genético e de ambiente permanente são importantes para a raça Simental no Brasil e devem ser considerados em programas de seleção. Entretanto, os valores mais elevados de herdabilidade materna, encontrados com modelos sem efeito de ambiente permanente, sugerem que o método utilizado não discrimina apropriadamente esses efeitos, oriundos de mesma fonte de variação.

Methods to verify parameter equality in nonlinear regression models

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

Comparação entre alguns critérios de seleção para crescimento na raça Nelore

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

PHASE-TRANSITION IN A Q-DEFORMED LIPKIN MODEL

Assuming q-deformed commutation relations for the fermions, an extension of the standard Lipkin Hamiltonian is presented. The usual quasi-spin representation of the standard Lipkin model is also obtained in this q-deformed framework. A variationally obtained energy functional is used to analyse the phase transition associated with the spherical symmetry breaking. The only phase transitions in this q-deformed model are of second order. As an outcome of this analysis a critical parameter is obtained which is dependent on the deformation of the algebra and on the number of particles.

Mechanism of the Yb-Er energy transfer in fluorozirconate glass

The mechanism of the Yb(3+)-->Er(3+) energy transfer as a function of the donor and the acceptor concentration was investigated in Yb(3+)-Er(3+) codoped fluorozirconate glass. The luminescence decay curves were measured and analyzed by monitoring the Er(3+)((4)I(11/2)) fluorescence induced by the Yb(3+)((2)F(5/2)) excitation. The energy transfer microparameters were determined and used to estimate the Yb-Er transfer rate of an energy transfer process assisted by excitation migration among donors state (diffusion model). The experimental transfer rates were determined from the best fitting of the acceptor luminescence decay obtained using a theoretical approach analog to that one used in the Inokuti-Hirayama model for the donor luminescence decay. The obtained values of transfer parameter gamma [gamma(exp)] were always higher than that predicted by the Inokuti-Hirayama model. Also, the experimental transfer rate, gamma(2)(exp), was observed to be higher than the transfer rate predicted by the migration model. Assuming a random distribution among excited donors at the initial time (t=0) and that a fast excitation migration, which occurs in a very short time (t

Asymmetry parameter for nonmesonic hypernuclear decay

We give general expressions for the vector asymmetry in the angular distribution of protons in the nonmesonic weak decay of polarized hypernuclei. From these we derive an explicit expression for the calculation of the asymmetry parameter, a(Lambda), which is applicable to the specific cases of He-5(Lambda) and C-12(Lambda) described within the extreme shell model. In contrast to the approximate formula widely used in the literature, it includes the effects of three-body kinematics in the final states of the decay and correctly treats the contribution of transitions originating from single-proton states beyond the s-shell. This expression is then used for the corresponding numerical computation of a(Lambda) within several one-meson-exchange models. Besides the strictly local approximation usually adopted for the transition potential, we also consider the addition of the first-order nonlocality terms. We find values for a(Lambda) ranging from -0.62 to -0.24, in qualitative agreement with other theoretical estimates but in contradiction with some recent experimental determinations.

The STATFLUX code: a statistical method for calculation of flow and set of parameters, based on the Multiple-Compartment Biokinetical Model

The code STATFLUX, implementing a new and simple statistical procedure for the calculation of transfer coefficients in radionuclide transport to animals and plants, is proposed. The method is based on the general multiple-compartment model, which uses a system of linear equations involving geometrical volume considerations. Flow parameters were estimated by employing two different least-squares procedures: Derivative and Gauss-Marquardt methods, with the available experimental data of radionuclide concentrations as the input functions of time. The solution of the inverse problem, which relates a given set of flow parameter with the time evolution of concentration functions, is achieved via a Monte Carlo Simulation procedure.Program summaryTitle of program: STATFLUXCatalogue identifier: ADYS_v1_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/ADYS_v1_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions: noneComputer for which the program is designed and others on which it has been tested: Micro-computer with Intel Pentium III, 3.0 GHzInstallation: Laboratory of Linear Accelerator, Department of Experimental Physics, University of São Paulo, BrazilOperating system: Windows 2000 and Windows XPProgramming language used: Fortran-77 as implemented in Microsoft Fortran 4.0. NOTE: Microsoft Fortran includes non-standard features which are used in this program. Standard Fortran compilers such as, g77, f77, ifort and NAG95, are not able to compile the code and therefore it has not been possible for the CPC Program Library to test the program.Memory, required to execute with typical data: 8 Mbytes of RAM memory and 100 MB of Hard disk memoryNo. of bits in a word: 16No. of lines in distributed program, including test data, etc.: 6912No. of bytes in distributed Program, including test data, etc.: 229 541Distribution format: tar.gzNature of the physical problem: the investigation of transport mechanisms for radioactive substances, through environmental pathways, is very important for radiological protection of populations. One such pathway, associated with the food chain, is the grass-animal-man sequence. The distribution of trace elements in humans and laboratory animals has been intensively studied over the past 60 years [R.C. Pendlenton, C.W. Mays, R.D. Lloyd, A.L. Brooks, Differential accumulation of iodine-131 from local fallout in people and milk, Health Phys. 9 (1963) 1253-1262]. In addition, investigations on the incidence of cancer in humans, and a possible causal relationship to radioactive fallout, have been undertaken [E.S. Weiss, M.L. Rallison, W.T. London, W.T. Carlyle Thompson, Thyroid nodularity in southwestern Utah school children exposed to fallout radiation, Amer. J. Public Health 61 (1971) 241-249; M.L. Rallison, B.M. Dobyns, F.R. Keating, J.E. Rall, F.H. Tyler, Thyroid diseases in children, Amer. J. Med. 56 (1974) 457-463; J.L. Lyon, M.R. Klauber, J.W. Gardner, K.S. Udall, Childhood leukemia associated with fallout from nuclear testing, N. Engl. J. Med. 300 (1979) 397-402]. From the pathways of entry of radionuclides in the human (or animal) body, ingestion is the most important because it is closely related to life-long alimentary (or dietary) habits. Those radionuclides which are able to enter the living cells by either metabolic or other processes give rise to localized doses which can be very high. The evaluation of these internally localized doses is of paramount importance for the assessment of radiobiological risks and radiological protection. The time behavior of trace concentration in organs is the principal input for prediction of internal doses after acute or chronic exposure. The General Multiple-Compartment Model (GMCM) is the powerful and more accepted method for biokinetical studies, which allows the calculation of concentration of trace elements in organs as a function of time, when the flow parameters of the model are known. However, few biokinetics data exist in the literature, and the determination of flow and transfer parameters by statistical fitting for each system is an open problem.Restriction on the complexity of the problem: This version of the code works with the constant volume approximation, which is valid for many situations where the biological half-live of a trace is lower than the volume rise time. Another restriction is related to the central flux model. The model considered in the code assumes that exist one central compartment (e.g., blood), that connect the flow with all compartments, and the flow between other compartments is not included.Typical running time: Depends on the choice for calculations. Using the Derivative Method the time is very short (a few minutes) for any number of compartments considered. When the Gauss-Marquardt iterative method is used the calculation time can be approximately 5-6 hours when similar to 15 compartments are considered. (C) 2006 Elsevier B.V. All rights reserved.

MODEL FOR ASYMPTOTIC D-STATE PARAMETERS OF LIGHT-NUCLEI - APPLICATION TO HE-4

A simple method for calculating the asymptotic D-state observables for light nuclei is suggested. The method exploits the dominant clusters of the light nuclei. The method is applied to calculate the He-4 asymptotic D to S normalization ratio rho(alpha) and the closely related D-state parameter D2alpha. The study predicts a correlation between D2alpha and B(alpha), and between rho(alpha) and B(alpha), where B(alpha) is the binding energy of He-4. The present study yields rho(alpha) congruent-to -0.14 and D2alpha congruent-to -0.12 fm2 consistent with the correct experimental eta(d) and the binding energies of the deuteron, triton, and the alpha particle, where eta(d) is the deuteron D-state to S-state normalization ratio.

Fluctuating commutative geometry

We use the framework of noncommutative geometry to define a discrete model for fluctuating geometry. Instead of considering ordinary geometry and its metric fluctuations, we consider generalized geometries where topology and dimension can also fluctuate. The model describes the geometry of spaces with a countable number n of points. The spectral principle of Connes and Chamseddine is used to define dynamics. We show that this simple model has two phases. The expectation value , the average number of points in the universe, is finite in one phase and diverges in the other. Moreover, the dimension delta is a dynamical observable in our model, and plays the role of an order parameter. The computation of is discussed and an upper bound is found, < 2. We also address another discrete model defined on a fixed d = 1 dimension, where topology fluctuates. We comment on a possible spontaneous localization of topology.

Energy localization in the Peyrard-Bishop DNA model

We study energy localization on the oscillator chain proposed by Peyrard and Bishop to model DNA. We search numerically for conditions with initial energy in a small subgroup of consecutive oscillators of a finite chain and such that the oscillation amplitude is small outside this subgroup on a long time scale. We use a localization criterion based on the information entropy and verify numerically that such localized excitations exist when the nonlinear dynamics of the subgroup oscillates with a frequency inside the reactive band of the linear chain. We predict a mimium value for the Morse parameter (mu>2.25) (the only parameter of our normalized model), in agreement with the numerical calculations (an estimate for the biological value is mu=6.3). For supercritical masses, we use canonical perturbation theory to expand the frequencies of the subgroup and we calculate an energy threshold in agreement with the numerical calculations.

The use of nonparametric contrasts in one-way layouts and random block designs

Nonparametric simple-contrast estimates for one-way layouts based on Hodges-Lehmann estimators for two samples and confidence intervals for all contrasts involving only two treatments are found in the literature.Tests for such contrasts are performed from the distribution of the maximum of the rank sum between two treatments. For random block designs, simple contrast estimates based on Hodges-Lehmann estimators for one sample are presented. However, discussions concerning the significance levels of more complex contrast tests in nonparametric statistics are not well outlined.This work aims at presenting a methodology to obtain p-values for any contrast types based on the construction of the permutations required by each design model using a C-language program for each design type. For small samples, all possible treatment configurations are performed in order to obtain the desired p-value. For large samples, a fixed number of random configurations are used. The program prompts the input of contrast coefficients, but does not assume the existence or orthogonality among them.In orthogonal contrasts, the decomposition of the value of the suitable statistic for each case is performed and it is observed that the same procedure used in the parametric analysis of variance can be applied in the nonparametric case, that is, each of the orthogonal contrasts has a chi(2) distribution with one degree of freedom. Also, the similarities between the p-values obtained for nonparametric contrasts and those obtained through approximations suggested in the literature are discussed.

Herdabilidades e Correlações Genéticas para Peso e Perímetro Escrotal de Machos e Características Reprodutivas e de Crescimento de Fêmeas, na Raça Canchim

The objective of this study was to estimate heritabilities of and genetic correlations among male body weight (BW12) and scrotal circumference (SC12) at 12 months of age, and female body weights at first (BWFC) and second (BWSC) calvings, age at first (AFC) and second (ASC) calvings, adult weight (AW), and mature weight (A) and maturation rate (k) obtained by the use of the Von Bertalanffy model. The restricted maximum likelihood method with an animal model that included the fixed effects of contemporary group and the random effects of animals, was used to estimate the variance and covariance components. The heritability estimates were equal to: 0.37 (BW12),0.30 (SC12),0.38 (A), 0.35 (k), 0.12 (AFC), 0.33 (BWFC), 0.04 (ASC), 0.39 (BWSC), and 0.38 (AW). The genetic correlations among BW12 and the female traits were: 0.19 (parameter A), 0.62 (parameter k), -0.58 (AFC), 0.69 (BWFC), -0.56 (ASC), 0.61 (BWSC), and 0.60 (AW). The genetic correlations among SC12 and the female traits were: -0.24 (A), 0.27 (k), -0.47 (AFC), 0.09 (BWFC), -0.67 (ASC), 0.07 (BWSC), and -0.17 (AW). These results indicate that male body weight and scrotal circumference and female weights (BWFC, BWSC and AW) and growth curve parameters A and k have enough additive genetic variation to respond to mass selection. Selection to increase male body weight at 12 months of age should result on favorable correlated changes in AFC, ASC and parameter k of females, but with increases in female body weights (BWFC, BWSC and AW). Selection to increase SC12 should result on desirable correlated responses in AFC, ASC and k, without any considerable change in female adult body weights.

Fluctuating dimension in a discrete model for quantum gravity based on the spectral principle

The spectral principle of Connes and Chamseddine is used as a starting point to define a discrete model for Euclidean quantum gravity. Instead of summing over ordinary geometries, we consider the sum over generalized geometries where topology, metric, and dimension can fluctuate. The model describes the geometry of spaces with a countable number n of points, and is related to the Gaussian unitary ensemble of Hermitian matrices. We show that this simple model has two phases. The expectation value , the average number of points in the Universe, is finite in one phase and diverges in the other. We compute the critical point as well as the critical exponent of . Moreover, the space-time dimension delta is a dynamical observable in our model, and plays the role of an order parameter. The computation of is discussed and an upper bound is found, < 2.

Integrable anisotropic spin-ladder model

We present an integrable spin-ladder model, which possesses a free parameter besides the rung coupling J. Wang's system based on the SU(4) symmetry can be obtained as a special case. The model is exactly solvable by means of the Bethe ansatz method. We determine the dependence on the anisotropy parameter of the phase transition between gapped and gapless spin excitations and present the phase diagram. Finally, we show that the model is a special case of a more general Hamiltonian with three free parameters.