Estimativas de parâmetros genéticos para produção de leite em búfalas por modelos de repetibilidade, multi-característica e de regressão aleatória

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

O Raciocínio lógico-matemático: sua estrutura neurofisiológica e aplicações à educação matemática

Pós-graduação em Educação Matemática - IGCE

Simulação de perfis nucleares de poço em formações complexas

A identificação e descrição dos caracteres litológicos de uma formação são indispensáveis à avaliação de formações complexas. Com este objetivo, tem sido sistematicamente usada a combinação de ferramentas nucleares em poços não-revestidos. Os perfis resultantes podem ser considerados como a interação entre duas fases distintas: • Fase de transporte da radiação desde a fonte até um ou mais detectores, através da formação. • Fase de detecção, que consiste na coleção da radiação, sua transformação em pulsos de corrente e, finalmente, na distribuição espectral destes pulsos. Visto que a presença do detector não afeta fortemente o resultado do transporte da radiação, cada fase pode ser simulada independentemente uma da outra, o que permite introduzir um novo tipo de modelamento que desacopla as duas fases. Neste trabalho, a resposta final é simulada combinando soluções numéricas do transporte com uma biblioteca de funções resposta do detector, para diferentes energias incidentes e para cada arranjo específico de fontes e detectores. O transporte da radiação é calculado através do algoritmo de elementos finitos (FEM), na forma de fluxo escalar 2½-D, proveniente da solução numérica da aproximação de difusão para multigrupos da equação de transporte de Boltzmann, no espaço de fase, dita aproximação P₁, onde a variável direção é expandida em termos dos polinômios ortogonais de Legendre. Isto determina a redução da dimensionalidade do problema, tornando-o mais compatível com o algoritmo FEM, onde o fluxo dependa exclusivamente da variável espacial e das propriedades físicas da formação. A função resposta do detector NaI(Tl) é obtida independentemente pelo método Monte Carlo (MC) em que a reconstrução da vida de uma partícula dentro do cristal cintilador é feita simulando, interação por interação, a posição, direção e energia das diferentes partículas, com a ajuda de números aleatórios aos quais estão associados leis de probabilidades adequadas. Os possíveis tipos de interação (Rayleigh, Efeito fotoelétrico, Compton e Produção de pares) são determinados similarmente. Completa-se a simulação quando as funções resposta do detector são convolvidas com o fluxo escalar, produzindo como resposta final, o espectro de altura de pulso do sistema modelado. Neste espectro serão selecionados conjuntos de canais denominados janelas de detecção. As taxas de contagens em cada janela apresentam dependências diferenciadas sobre a densidade eletrônica e a fitologia. Isto permite utilizar a combinação dessas janelas na determinação da densidade e do fator de absorção fotoelétrico das formações. De acordo com a metodologia desenvolvida, os perfis, tanto em modelos de camadas espessas quanto finas, puderam ser simulados. O desempenho do método foi testado em formações complexas, principalmente naquelas em que a presença de minerais de argila, feldspato e mica, produziram efeitos consideráveis capazes de perturbar a resposta final das ferramentas. Os resultados mostraram que as formações com densidade entre 1.8 e 4.0 g/cm³ e fatores de absorção fotoelétrico no intervalo de 1.5 a 5 barns/e^-, tiveram seus caracteres físicos e litológicos perfeitamente identificados. As concentrações de Potássio, Urânio e Tório, puderam ser obtidas com a introdução de um novo sistema de calibração, capaz de corrigir os efeitos devidos à influência de altas variâncias e de correlações negativas, observadas principalmente no cálculo das concentrações em massa de Urânio e Potássio. Na simulação da resposta da sonda CNL, utilizando o algoritmo de regressão polinomial de Tittle, foi verificado que, devido à resolução vertical limitada por ela apresentada, as camadas com espessuras inferiores ao espaçamento fonte - detector mais distante tiveram os valores de porosidade aparente medidos erroneamente. Isto deve-se ao fato do algoritmo de Tittle aplicar-se exclusivamente a camadas espessas. Em virtude desse erro, foi desenvolvido um método que leva em conta um fator de contribuição determinado pela área relativa de cada camada dentro da zona de máxima informação. Assim, a porosidade de cada ponto em subsuperfície pôde ser determinada convolvendo estes fatores com os índices de porosidade locais, porém supondo cada camada suficientemente espessa a fim de adequar-se ao algoritmo de Tittle. Por fim, as limitações adicionais impostas pela presença de minerais perturbadores, foram resolvidas supondo a formação como que composta por um mineral base totalmente saturada com água, sendo os componentes restantes considerados perturbações sobre este caso base. Estes resultados permitem calcular perfis sintéticos de poço, que poderão ser utilizados em esquemas de inversão com o objetivo de obter uma avaliação quantitativa mais detalhada de formações complexas.

A simple application of implicit function theorem

In this note we show that the roots of a polynomial are C∞ depend of the coefficients. The main tool to show this is the Implicit Function Theorem.

Estimates of covariance functions for growth from birth to 630 days of age in Nelore cattle

Weight records of Brazilian Nelore cattle, from birth to 630 d of age, recorded every 3 mo, were analyzed using random regression models. Independent variables were Legendre polynomials of age at recording. The model of analysis included contemporary groups as fixed effects and age of dam as a linear and quadratic covariable. Mean trends were modeled through a cubic regression on orthogonal polynomials of age. Up to four sets of random regression coefficients were fitted for animals' direct and maternal, additive genetic, and permanent environmental effects. Changes in measurement error variances with age were modeled through a variance function. Orders of polynomial fit from three to six were considered, resulting in up to 77 parameters to be estimated. Models fitting random regressions modeled the pattern of variances in the data adequately, with estimates similar to those from corresponding univariate analysis. Direct heritability estimates decreased after birth and tended to be lowest at ages at which maternal effect estimates tended to be highest. Maternal heritability estimates increased after birth to a peak around 110 to 120 d of age and decreased thereafter. Additive genetic direct correlation estimates between weights at standard ages (birth, weaning, yearling, and final weight) were moderate to high and maternal genetic and environmental correlations were consistently high.

The McDonald inverted beta distribution

We introduce a five-parameter continuous model, called the McDonald inverted beta distribution, to extend the two-parameter inverted beta distribution and provide new four- and three-parameter sub-models. We give a mathematical treatment of the new distribution including expansions for the density function, moments, generating and quantile functions, mean deviations, entropy and reliability. The model parameters are estimated by maximum likelihood and the observed information matrix is derived. An application of the new model to real data shows that it can give consistently a better fit than other important lifetime models. (C) 2012 The Franklin Institute. Published by Elsevier Ltd. All rights reserved.

Comparative study on the far-field spectra and near-field amplitudes for silver and gold nanocubes irradiated at 514, 633 and 785 nm as a function of the edge length

We describe a systematic investigation by the discrete dipole approximation on the optical properties of silver (Ag) and gold (Au) nanocubes as a function of the edge length in the 20-100 nm range. Our results showed that, as the nanocube size increased, the plasmon resonance modes shifted to higher wavelengths, the contribution from scattering to the extinction increased, and the quadrupole modes became more intense in the spectra. The electric field amplitudes at the surface of the nanocubes were calculated considering 514, 633 and 785 nm as the excitation wavelengths. While Ag nanocubes displayed the highest electric field amplitudes (vertical bar E vertical bar(max)) when excited at 514 nm, the Au nanocubes displayed higher vertical bar E vertical bar(max) values than Ag, for all sizes investigated, when the excitation wavelength was either 633 or 785 nm. The variations in vertical bar E vertical bar(max) as a function of size for both Ag and Au nanocubes could be explained based on the relative position of the surface plasmon resonance peak relative to the wavelength of the incoming electromagnetic wave. Our results show that not only size and composition, but also the excitation wavelength, can play an important role over the maximum near-field amplitudes values generated at the surface of the nanocubes.

Loading a K-39 crossed optical dipole trap from a magneto-optical trap

In this work, we have applied sub-Doppler laser cooling to a K-39 magneto-optical trap in order to load a 1071 nm crossed optical dipole trap. The number of atoms loaded into the dipole trap was characterized as a function of the frequency and intensity of the cooling and repump laser beams. For the optimum conditions, the dipole trap has about 2 x 10(6) atoms at an atomic density of 2 x 10(12) cm(-3), with a temperature of about 10 mu K. This technique is a very simple procedure to load a K-39 optical dipole trap without a previous magnetic evaporative cooling step and may find application in other atomic physic systems.

Effectiveness of a long-term use of a minimalist footwear versus habitual shoe on pain, function and mechanical loads in knee osteoarthritis: a randomized controlled trial

Background: Recent studies have shown an important reduction of joint overload during locomotion in elderly women with knee osteoarthritis (OA) after short- term use of minimalist shoes. Our aim is to investigate the chronic effect of inexpensive and minimalist footwear on the clinical and functional aspects of OA and gait biomechanics of elderly women with knee OA. Methods/Design: Fifty-six elderly women with knee OA grade 2 or 3 (Kellgren and Lawrence) are randomized into blocks and allocated to either the intervention group, which will use flexible, non-heeled shoes-Moleca (R)-for six months for at least six hours daily, or the control group, which could not use these shoes. Neither group is undergoing physical therapy treatment throughout the intervention period. Moleca (R) is a women's double canvas, flexible, flat walking shoe without heels, with a 5-mm anti-slip rubber sole and a 3-mm internal wedge of ethylene vinyl acetate. Both groups will be followed for six months and will be assessed at baseline condition, after three months, and after six months (end of intervention). All the assessments will be performed by a physiotherapist that is blind to the group allocation. The primary outcome is the pain Western Ontario and McMaster Universities Osteoarthritis (WOMAC) score. The secondary outcomes are global WOMAC score; joint stiffness and disability WOMAC scores; knee pain with a visual analogue scale; walking distance in the six-minute walk test; Lequesne score; amount and frequency (number of days) of paracetamol (500 mg) intake over six months; knee adduction moment during gait; global medical assessment score; and global patient auto-assessment score. At baseline, all patients receive a diary to record the hours of daily use of the footwear intervention; every two weeks, the same physiotherapist makes phone calls to all patients in order to verify adherence to treatment. The statistical analysis will be based on intention to treat analysis, as well as general linear models of analysis of variance for repeated measure to detect treatment-time interactions (alpha = 5%). Discussion: This is the first randomized, clinical trial protocol to assess the chronic effect of minimalist footwear on the clinical and functional aspects and gait biomechanics of elderly women with knee osteoarthritis. We expect that the use of Moleca (R) shoes for six months will provide pain relief, reduction of the knee adduction moment when walking, and improve joint function in elderly women with knee OA, and that the treatment, thus, can be considered another inexpensive and easy-to-use option for conservative OA treatment.

Dissociative electron attachment to chloromethane.

A low energy electron may attach to a molecule, forming a metastable resonance, which may dissociate into a stable anion and a neutral radical. Chloromethane has been a good target for dissociative electron attachment studies, since it is a small molecule with a clear dissociative ‘sigma*’ shape resonance. We present potential energy curves for CH3Cl and its anion, as a function of the C-Cl distance. Due to the resonant nature of the anion, a correct description requires a treatment based on scattering calculations. In order to compute elastic cross sections and phase shifts we employed the Schwinger multichannel method, implemented with pseudopotentials of Bachelet, Hamann and Schlüter, at the static-exchange plus polarization approximation. At the equilibrium geometry, the resonance was found arround 3.3 eV, in accordance to experience. The incoming electron is captured by a ‘sigma*’ orbital located at the C-Cl bond, which will relax in the presence of this extra electron. We took this bond as the reaction coordinate, and performed several scattering calculations for a series of nuclear conformations. The phase shift obtained in each calculation was fitted by a two component function, consisting in the usual Breit-Wigner profile, which captures the resonant character, and a second order polynomial in the wave number, which accounts for the background contribution. That way, we obtained position and width of the resonance, which allowed us to build the potential energy curve. For larger distances, the anion becomes stable and usual electronic structure calculations suffice. Furthermore, the existence of a dipole-bound anion state is revealed when we employed a set of very diffuse functions. The knowledge on the behaviour of the neutral and anionic electronic states helps us in elucidating how the dissociation takes place.

Analysing the possible influence of transposon TCl4.7 insertion on the function of the genome of Cydia pomonella granulovirus

CpGV-MCp5 is a natural mutant of the Cydia pomonella Granulovirus (Mexican isolate) (CpGV-M) that harbors an insect host transposon termed TCl4.7 in its genome. TCl4.7 is located between the open reading frames Cp15 and Cp16 and separates two homologous regions hr3 and hr4, which have been recently shown to be origins of replication of CpGV-M. The MCp5 has a significant replication disadvantage in the presence of the wild-type CpGV-M. In this study, the possible effects of TCl4.7 transposon insertion on the genome function of its insertion site has been analysed. The role of Cp15 and Cp16 in the context of the virus infection cycle was examined by generating a CpGV-Bacmid (CpBAC) and Cp15 knock-out (CpBACCp15KO) and Cp16 knock-out (CpBACCp16KO) mutants. The mutant CpBACCp15KO was not able to replicate in CM larvae suggesting that Cp15 was essential for virus replication. In contrast, the mutant CpBACCp16KO infected CM larvae and produced viable occlusion bodies (OBs) demonstrating that Cp16 is a non-essential gene for virus in vivo infection of C. pomonella. The temporal transcription of Cp15 and Cp16, as well as of Cp31 (F protein) as a control, was analysed using RT-PCR and quantitative real-time PCR. It suggested a general delay or reduction of gene transcription of MCp5 compared to the parental CpGV-M. Western blot analyses using anti-Cp15 and anti-Cp16 polyclonal antibodies, however, did not show any immuno-reactive response. Thus, a direct influence of TCl4.7 on the expression of Cp15 and Cp16 could not be substantiated. To investigate whether the interruption of hr3 and hr4 palindromes affects the virus replication, two mutant bacmids with a deletion of hr3 and hr4 (CpBAChr3/hr4-KO) and another with an insertion of a Kanamycin resistance cassette between hr3 and hr4 (CpBAChr3-kan-hr4) were generated. Both mutant bacmids replicated and produced infectious virus OBs, which did not significantly differ in their median lethal concentration (LC50) and median survival time (ST50) compared to the parental CpBAC. Interestingly, the mutant CpBAChr3-kan-hr4 was very effectively out-competed by parental CpBAC, when CM larvae were co-infected with known ratios of OBs of CpBAC and the mutant CpBAChr3-kan-hr4. These observations suggested a functional co-operation between hr3 and hr4 which was interrupted by the KanR insertion in CpBAChr3-kan-hr4 and possibly by TCl4.7 transposon insertion in the mutant MCp5. This hypothesis may explain the observed replication disadvantage of the mutants MCp5 and CpBAChr3-kan-hr4 in the presence of the parental viruses CpGV-M and CpBAC, respectively.

Investigation of the dipole response of nickel isotopes in the presence of a high-frequency electromagnetic field

The electric dipole response of neutron-rich nickel isotopes has been investigated using the LAND setup at GSI in Darmstadt (Germany). Relativistic secondary beams of 56−57Ni and 67−72Ni at approximately 500 AMeV have been generated using projectile fragmentation of stable ions on a 4 g/cm2 Be target and subsequent separation in the magnetic dipole fields of the FRagment Separator (FRS). After reaching the LAND setup in Cave C, the radioactive ions were excited electromagnetically in the electric field of a Pb target. The decay products have been measured in inverse kinematics using various detectors. Neutron-rich 67−69Ni isotopes decay by the emission of neutrons, which are detected in the LAND detector. The present analysis concentrates on the (gamma,n) and (gamma,2n) channels in these nuclei, since the proton and three-neutron thresholds are unlikely to be reached considering the virtual photon spectrum for nickel ions at 500 AMeV. A measurement of the stable 58Ni isotope is used as a benchmark to check the accuracy of the present results with previously published data. The measured (gamma,n) and (gamma,np) channels are compared with an inclusive photoneutron measurement by Fultz and coworkers, which are consistent within the respective errors. The measured excitation energy distributions of 67−69Ni contain a large portion of the Giant Dipole Resonance (GDR) strength predicted by the Thomas-Reiche-Kuhn energy-weighted sum rule, as well as a significant amount of low-lying E1 strength, that cannot be attributed to the GDR alone. The GDR distribution parameters are calculated using well-established semi-empirical systematic models, providing the peak energies and widths. The GDR strength is extracted from the chi-square minimization of the model GDR to the measured data of the (gamma,2n) channel, thereby excluding any influence of eventual low-lying strength. The subtraction of the obtained GDR distribution from the total measured E1 strength provides the low-lying E1 strength distribution, which is attributed to the Pygmy Dipole Resonance (PDR). The extraction of the peak energy, width and strength is performed using a Gaussian function. The minimization of trial Gaussian distributions to the data does not converge towards a sharp minimum. Therefore, the results are presented by a chi-square distribution as a function of all three Gaussian parameters. Various predictions of PDR distributions exist, as well as a recent measurement of the 68Ni pygmy dipole-resonance obtained by virtual photon scattering, to which the present pygmy dipole-resonance distribution is also compared.

Function of liver activation-regulated chemokine/CC chemokine ligand 20 is differently affected by cathepsin B and cathepsin D processing

Chemokine processing by proteases is emerging as an important regulatory mechanism of leukocyte functions and possibly also of cancer progression. We screened a large panel of chemokines for degradation by cathepsins B and D, two proteases involved in tumor progression. Among the few substrates processed by both proteases, we focused on CCL20, the unique chemokine ligand of CCR6 that is expressed on immature dendritic cells and subtypes of memory lymphocytes. Analysis of the cleavage sites demonstrate that cathepsin B specifically cleaves off four C-terminally located amino acids and generates a CCL20(1-66) isoform with full functional activity. By contrast, cathepsin D totally inactivates the chemotactic potency of CCL20 by generating CCL20(1-55), CCL20(1-52), and a 12-aa C-terminal peptide CCL20(59-70). Proteolytic cleavage of CCL20 occurs also with chemokine bound to glycosaminoglycans. In addition, we characterized human melanoma cells as a novel CCL20 source and as cathepsin producers. CCL20 production was up-regulated by IL-1alpha and TNF-alpha in all cell lines tested, and in human metastatic melanoma cells. Whereas cathepsin D is secreted in the extracellular milieu, cathepsin B activity is confined to cytosol and cellular membranes. Our studies suggest that CCL20 processing in the extracellular environment of melanoma cells is exclusively mediated by cathepsin D. Thus, we propose a model where cathepsin D inactivates CCL20 and possibly prevents the establishment of an effective antitumoral immune response in melanomas.

Kernel Estimation of Rate Function for Recurrent Event Data

Recurrent event data are largely characterized by the rate function but smoothing techniques for estimating the rate function have never been rigorously developed or studied in statistical literature. This paper considers the moment and least squares methods for estimating the rate function from recurrent event data. With an independent censoring assumption on the recurrent event process, we study statistical properties of the proposed estimators and propose bootstrap procedures for the bandwidth selection and for the approximation of confidence intervals in the estimation of the occurrence rate function. It is identified that the moment method without resmoothing via a smaller bandwidth will produce curve with nicks occurring at the censoring times, whereas there is no such problem with the least squares method. Furthermore, the asymptotic variance of the least squares estimator is shown to be smaller under regularity conditions. However, in the implementation of the bootstrap procedures, the moment method is computationally more efficient than the least squares method because the former approach uses condensed bootstrap data. The performance of the proposed procedures is studied through Monte Carlo simulations and an epidemiological example on intravenous drug users.

Generating Realistic Camera Shake for Virtual Scenes

When depicting both virtual and physical worlds, the viewer's impression of presence in these worlds is strongly linked to camera motion. Plausible and artist-controlled camera movement can substantially increase scene immersion. While physical camera motion exhibits subtle details of position, rotation, and acceleration, these details are often missing for virtual camera motion. In this work, we analyze camera movement using signal theory. Our system allows us to stylize a smooth user-defined virtual base camera motion by enriching it with plausible details. A key component of our system is a database of videos filmed by physical cameras. These videos are analyzed with a camera-motion estimation algorithm (structure-from-motion) and labeled manually with a specific style. By considering spectral properties of location, orientation and acceleration, our solution learns camera motion details. Consequently, an arbitrary virtual base motion, defined in any conventional animation package, can be automatically modified according to a user-selected style. In an animation package the camera motion base path is typically defined by the user via function curves. Another possibility is to obtain the camera path by using a mixed reality camera in motion capturing studio. As shown in our experiments, the resulting shots are still fully artist-controlled, but appear richer and more physically plausible.