Submitted by Karem Dantas (karem.c.dantas@gmail.com) on 2018-10-01T12:46:45Z No. of bitstreams: 2 license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) Dissertação_Natan.pdf: 1337747 bytes, checksum: 5177f22611e7d754c7185a7da2a54aa0 (MD5) / Approved for entry into archive by Marcos Roberto Gomes (mrobertosg@gmail.com) on 2018-10-01T17:34:06Z (GMT) No. of bitstreams: 2 license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) Dissertação_Natan.pdf: 1337747 bytes, checksum: 5177f22611e7d754c7185a7da2a54aa0 (MD5) / Approved for entry into archive by Divisão de Documentação/BC Biblioteca Central (ddbc@ufam.edu.br) on 2018-10-02T15:12:25Z (GMT) No. of bitstreams: 2 license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) Dissertação_Natan.pdf: 1337747 bytes, checksum: 5177f22611e7d754c7185a7da2a54aa0 (MD5) / Made available in DSpace on 2018-10-02T15:12:25Z (GMT). No. of bitstreams: 2 license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) Dissertação_Natan.pdf: 1337747 bytes, checksum: 5177f22611e7d754c7185a7da2a54aa0 (MD5) Previous issue date: 2018-05-03 / The rape of vulnerable people is a phenomenon that has experienced a considerable increase in its occurrence. However these numbers are not vastly known due to the high index of under-reportings related to this type of crime. This study aims to estimate the rate of rape of vulnerable people by means of the Bayesian inference and data augmentation techniques. The proposed methodology used the data of rape of vulnerable people of some countryside cities of Amazonas referring to the period of 2010 to 2012. Due to the good adherence to the aforementioned data, the deployment of a multivariate Poisson model showed its efficacy. / O estupro de vulnerável é um fenômeno que tem sofrido um considerado aumento em sua ocorrência, mas pouco se conhece sobre esses números devido ao alto índice de subnotificações referente a este tipo de crime. O objetivo deste estudo é estimar a taxa de estupro de vulnerável por meio da inferência bayesiana e técnicas de aumento de dados. A metodologia proposta utilizou os dados de estupro de vulnerável de algumas cidades do interior do estado do Amazonas referente ao período de 2010 a 2012. O emprego de um modelo Poisson multivariado mostrou-se eficaz,apresentando uma boa aderência aos dados supracitados.

Motivados por um conjunto de dados de pacientes com insuficiência renal crônica (IRC), propomos uma nova modelagem bayesiana que envolve cópulas da família Arquimediana e um modelo misto para dados de sobrevivência com riscos semicompetitivos. A estrutura de riscos semicompetitivos é bastante comum em estudos clínicos em que dois eventos são de interesse, um intermediário e outro terminal, de forma tal que a ocorrência do evento terminal impede a ocorrência do intermediário mas não vice-versa. Nesta modelagem provamos que a distribuição a posteriori sob a cópula de Clayton é própria. Implementamos os algoritmos de dados aumentados e amostrador de Gibbs para a inferência bayesiana, assim como os criterios de comparação de modelos: LPML, DIC e BIC. Realizamos um estudo de simulação para avaliar o desempenho da modelagem e finalmente aplicamos a metodologia proposta para analisar os dados dos pacientes com IRC, além de outros de pacientes que receberam transplante de medula óssea. / Motivated by a dataset of patients with chronic kidney disease (CKD), we propose a new bayesian model including the Arquimedean copula and a mixed model for survival data with semicompeting risks. The structure of semicompeting risks appears frequently in clinical studies where two-types of events are involved: a nonterminal and a terminal event such that the occurrence of terminal event precludes the occurrence of the non-terminal event but not viceversa. In this work we prove that the posterior distribution is proper when the Clayton copula is used. We implement the data augmentation algorithm and Gibbs sampling for the bayesian inference, as well as some bayesian model selection criteria: LPML, BIC and DIC. We carry out a simulation study for assess the model performance and finally, our methodology is illustrated with the chronic kidney disease study.

Motivados por um conjunto de dados de pacientes com insuficiência renal crônica (IRC), propomos uma nova modelagem bayesiana que envolve cópulas da família Arquimediana e um modelo misto para dados de sobrevivência com riscos semicompetitivos. A estrutura de riscos semicompetitivos é bastante comum em estudos clínicos em que dois eventos são de interesse, um intermediário e outro terminal, de forma tal que a ocorrência do evento terminal impede a ocorrência do intermediário mas não vice-versa. Nesta modelagem provamos que a distribuição a posteriori sob a cópula de Clayton é própria. Implementamos os algoritmos de dados aumentados e amostrador de Gibbs para a inferência bayesiana, assim como os criterios de comparação de modelos: LPML, DIC e BIC. Realizamos um estudo de simulação para avaliar o desempenho da modelagem e finalmente aplicamos a metodologia proposta para analisar os dados dos pacientes com IRC, além de outros de pacientes que receberam transplante de medula óssea. / Motivated by a dataset of patients with chronic kidney disease (CKD), we propose a new bayesian model including the Arquimedean copula and a mixed model for survival data with semicompeting risks. The structure of semicompeting risks appears frequently in clinical studies where two-types of events are involved: a nonterminal and a terminal event such that the occurrence of terminal event precludes the occurrence of the non-terminal event but not viceversa. In this work we prove that the posterior distribution is proper when the Clayton copula is used. We implement the data augmentation algorithm and Gibbs sampling for the bayesian inference, as well as some bayesian model selection criteria: LPML, BIC and DIC. We carry out a simulation study for assess the model performance and finally, our methodology is illustrated with the chronic kidney disease study.

Made available in DSpace on 2014-06-11T19:26:08Z (GMT). No. of bitstreams: 0 Previous issue date: 2008-07-28Bitstream added on 2014-06-13T18:29:31Z : No. of bitstreams: 1 araujoneto_fr_me_jabo.pdf: 389310 bytes, checksum: 1c20248bba0e60c4bb9d24bce6f9ca34 (MD5) / Neste trabalho objetivou-se avaliar a aplicação de método bayesiano na análise de dados de crescimento, em um modelo multi-características. Dados de 54.182 animais da raça Nelore foram utilizados na estimação de componentes de variância de: a) pesos padronizados as idades de 120 (P20), 210 (P210), 365 (P365), 450 (P450) e 550 (P550) dias; b) perímetro escrotal as idades-padrão de 365 (PE365), 450 (PE450) e 550 (PE550) dias; c) ganhos de peso entre as idades 120/210 (GP1), 210/365 (GP2), 365/450 (GP3) e 450/550 (GP4) e; d) crescimento escrotal nos intervalos de 365/450 (CP1) e 450/550 (CP2). As análises foram realizadas empregando-se o software GIBBS2F90, assumindo um modelo animal para oito características: P120, PE365 e ganhos de peso e perímetro escrotal. As demais foram obtidas mediante propriedades da soma de variâncias. As herdabilidades diretas estimadas foram 0,17; 0,19; 0,13; 0,15; 0,33; 0,19; 0,23; 0,24; 0,35; 0,37; 0,39; 0,52; 0,56; e;0,48 para GP1, GP2, GP3, GP4, CP1, CP2, P120, P210, P365, P450, P550, PE365, PE450 e PE550, respectivamente. Os valores de correlação genética variaram de -0,025 (GP3/PE550) a 0,97 (PE450/PE550). Verificou-se que as características em idades padrão apresentam maior variabilidade genética que as medidas intervalares, o que as indica como critérios de seleção que possibilitaria uma melhor resposta a seleção. Com relação à implementação do método bayesiano em modelo multi-características, verificou-se eficiente apesar do pequeno número de amostras efetivas geradas em alguns parâmetros, indicando a necessidade de um maior número de interações. / This work, the objective is to evaluate application of Bayesian methods in multiple-trait animal models. Records from 54.182 Nellore males were used for estimation of variance components for: a) weight standardized ages of 120 (W120), 210 (W210), 365 (W365), 450 (W450) and 550 (W550) days; b) scrotal circumference at age-standard from 365 (EC365), 450 (EC450) and 550 (EC550); c) weight gains between the standardized ages of 120/210 (WG1), 210/365 (WG2), 365/450 (WG3) e 450/550 (WG4)e; d) scrotal growth in the intervals 365/450 (EG1) and 450/550 (EG2). Analyses were performed using the GIBBS2F90, assuming a multiple-trait animal model. The direct heritability estimates were 0,17; 0,19; 0,13; 0,15; 0,33; 0,19; 0,23; 0,24; 0,35; 0,37; 0,39; 0,52; 0,56; e;0,48 for WG1, WG2, WG3, WG4, EG1, EG2, W120, W210, W365, W450, W550, EC365, EC450 and EC550, respectively. The correlations values was ranging from -0,025 (between WG3 and EC550) to 0,97 (between EC450 and EC550) to genetics. It was found that the standard features in ages show greater genetic variability that measures interval, which indicates how the criteria of selection that would better meet the selection. Regarding the implementation of the Bayesian method in multi-model features, there was efficient despite the small number of effective samples generated in some parameters, indicating the need for a greater number of interactions.

The objective of this study was to estimate genetic and phenotypic parameters, besides the genetic progress for the performance traits post-weaning in a population of Charolais cattle. Components of (co)variance were estimated using an animal model by Bayesian inference method. In the first article was analyzed the average daily gain from weaning to yearling (ADGWY) and weight adjusted to 550 days of age (W550) of 5.897 animals, sired by 181 bulls and 3.897 cows born between 1983 and 1999. The average for direct heritability and standard errors were 0.39±0.04 and 0.48±0.05 for ADGWY and W550, respectively. The genetic correlation between the two traits analyzed and the standard error was 0.60±0.05. The genetic trends for ADGWY and W550 were -0.058g/year and -0.019kg/. In the second article were estimated parameters and genetics trends for visual scores as well as their associations with the average daily gain from weaning to yearling (ADGWY) and weight adjusted to 550 days of age (W550) of 2.964 animals, sired by 145 bulls and 1.820 cows born between 1994 and 2007. The heritability estimates for scores of conformation (C), precocity (P), muscle (M) and size (S) were: 0.13±0.04, 0.23±0.04, 0.16±0.03 and 0.11±0.04, respectively. The genetic correlation between the four visual scores ranged between -0.03 and 0.83. The genetic correlations between ADGWY and W550 with visual scores were negative, ranging from - 0.07 to -0.82. As for the study of genetic progress, the genetic trends for C, P, M and S were: 0.0019, 0.0027, 0.0017 and 0.0011, respectively. / Este estudo teve como objetivo estimar parâmetros genéticos, além do progresso genético para as características de desempenho na pós-desmama em uma população bovina da raça Charolês. Os componentes de (co)variâncias foram estimados utilizando um modelo animal, através do método de inferência Bayesiana. No primeiro capítulo foram analisadas as características ganho médio diário da desmama ao sobreano (GMDDS) e o peso ajustado aos 550 dias de idade (P550) de 5.897 animais, filhos de 181 touros e 3.897 vacas, nascidos entre 1983 e 1999. As médias a posteriori para a herdabilidade direta e os erros padrões foram 0,39±0,04 e 0,48±0,05 para GMDDS e P550, respectivamente. A correlação genética entre as duas características analisadas e o erro padrão foi 0,60±0,05. As tendências genéticas estimadas para GMDDS e P550 foram -0,058g/ano e -0,019kg/ano. No segundo capítulo foram estimados parâmetros e tendências genéticas para os escores visuais, bem como suas associações com o ganho médio diário da desmama ao sobreano (GMDDS) e o peso ajustado aos 550 dias de idade (P550) de 2.964 animais, filhos de 145 touros e 1.820 vacas, nascidos entre 1994 e 2007. As herdabilidades estimadas para os escores de conformação (C), precocidade (P), musculatura (M) e tamanho (T) foram: 0,13±0,04, 0,23±0,04, 0,16±0,03, 0,11±0,04, respectivamente. A correlação genética entre os quatro escores variaram entre - 0,03 e 0,83. Já as correlações genéticas entre GMDDS e P550 com os escores visuais foram todas negativas, variando de -0,07 a -0,82. Quanto ao estudo do progresso genético, as tendências genéticas para C, P, M e T foram: 0,0019, 0,0027, 0,0017 e 0,0011, respectivamente.

Made available in DSpace on 2015-03-26T13:42:34Z (GMT). No. of bitstreams: 1 texto completo.pdf: 9480742 bytes, checksum: f87c9b2653b911e326eef4b5a0dd87ce (MD5) Previous issue date: 2014-07-25 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / The aim of this study was to evaluate the quality of relationship information over genetic parameters estimates and accuracies of breeding values. Thus, we evaluated the quality by making corrections based on SNPs markers. Variance components were estimated under Bayesian approach via Gibbs sampling. The evaluation of the correct relationship was performed using markers of 3,591 individuals in a population consisting of 12,668 animals. Mendelian conflicts were deflned as conflicts between the progeny and parents markers. Thus, 460 changes were performed, which 54% had a buII or cow on pedigree. We observed, on average, match parent were defined with 75 markers, whiIe no- match parent n was performed with 2,700 markers. Annealing algorithm Molecular coancestry program provided 2,174 new haIf-sibs relationships. We observed that higher quality on pedigree information provided increase in accuracy of breeding values and higher heritability estimates (0.22 i 0.0286), suggesting the possibility of direct selection for tick resistance. We used a 5-fold cross-validation strategy using K-means and random methods to group the animals. The results of cross-validation indicate that higher quality of the pedigree relationships provide higher accuracy. The weaning weight was used to evaluate muItipIe sires progeny (MS) and embryo transfer and in vitro fertilization animals (TEF). We used three strategies to include MS information: genetic groups (GG), bayesian hierarchical method (HIER) and the average numerator relationship matrix (ANRM). The deviance information criteria (DIC) was used as Bayesian measure of fIt, which suggested that HIER strategy provided better fIt when MS information is consider. However, ANRM provided best fit to include TEF animals. We use foster dam information to estimate the maternal genetic and maternal permanent environmental effects when considering TEF animals. We didn't observed statistical difference in variance components and genetic parameters considering information TEF animals, but breeding values showed greater values. Spearman correlations between breeding values were higher for base, animals with certain paternity and MS progeny. However, the use of TEF animals information significantly changed the predicted breeding values. The use of appropriate methodologies to include the information of MS progeny and TEF animals provide most accurate prediction of breeding values and assist in higher rates of genetic gain. / O objetivo desse estudo foi verificar o reflexo da qualidade das informações de parentesco sobre as estimativas de parâmetros genéticos e acurácias das predições de valor genética Dessa forma, foi avaliada a qualidade dessas informações ao realizar correções no pedigree baseadas em marcadores do tipo SNPs. Os componentes de variância foram estimados sob enforque Bayesiano por amostragem de Gibbs. A avaliação da correta definição de parentesco foi realizada utilizando-se informações de marcadores SNPs de 3.591 indivíduos em uma população constituída de 12.668 animais. Os conflitos mendelianos entre as marcas da progênie e dos pais foram utilizados como critério de avaliação de parentesco. Dessa forma, foram realizadas 460 mudanças no pedigree, dentre os quais 54% possuíam um touro ou vaca identificado. Foi observado que, em média, novos relacionamentos genéticos foram definidos com 75 marcas em conflito, enquanto que a rejeição do parentesco foi realizada com 2.700 marcas. O uso do programa Molecular Coancestry para inferência de parentesco a partir de informações de marcadores moleculares proporcionou a definição de 2.174 novos relacionamentos de meio-irmãos. Foi possível observar que as correções de parentesco proporcionaram aumento da acurácia média dos valores genéticos preditos. O aumento na qualidade das informações de pedigree proporcionou a estimação de herdabilidade de maior magnitude (0,22 i 0,0286), sugerindo a possibilidade de seleção direta para a resistência a carrapatos. Foi utilizada uma estratégia de validação cruzada pelo método K-médias e também de forma aleatória, no qual cinco grupos de treinamento foram formados. Os resultados da validação cruzada indicam que maior qualidade nos relacionamentos do pedigree proporcionam maior valor de acurácia. O peso padronizado aos 205 dias foi utilizado para avaliar a inclusão de informações de filhos de reprodutores múltiplos (RM) e de animais oriundos de biotécnicas da reprodução (TEF). Foram avaliadas três estratégias de inclusão dessas informações: grupos genéticos (GG), método hierárquico bayesiano (HIER) e matriz da média dos numeradores dos coeficientes de parentesco (ANRM). O critério de informação da deviance (DIC) foi utilizado como avaliador da qualidade de ajuste e sugeriu que a estratégia HIER proporcionou melhor ajuste ao considerar as informações de RM. Entretanto, o método ANRM apresentou melhor ajuste ao incluir informações dos animais TEF. A inclusão das informações de animais TEF foi realizada pelo uso da informação da receptora para estimação do efeito genético materno e de ambiente permanente materno. Não foi observada diferença estatística nas estimativas de componentes de variância e de parâmetros genéticos ao considerar informações de animais TEF, porém os valores genéticos preditos apresentaram maior magnitude. As correlações de Spearman entre os valores genéticos foram de elevadas magnitudes para os animais fundadores, animais com certeza no parentesco e para os fiIhos de RM. Entretanto, o uso das informações dos animais TEF modificou de forma significativa as predições dos valores genéticos. O uso de metodologias adequadas para incluir a informação de animais com incerteza de paternidade e animais oriundos de transferência de embriões e fertilização in vitro proporcionam a predição de valores genéticos mais acurados e auxiliam em maiores taxas de ganho genético.

Submitted by Automa??o e Estat?stica (sst@bczm.ufrn.br) on 2018-03-21T12:37:29Z No. of bitstreams: 1 MarconioSilvaDosSantos_TESE.pdf: 17070388 bytes, checksum: 46702d837c8c304ffc379088625742aa (MD5) / Approved for entry into archive by Arlan Eloi Leite Silva (eloihistoriador@yahoo.com.br) on 2018-03-23T15:26:03Z (GMT) No. of bitstreams: 1 MarconioSilvaDosSantos_TESE.pdf: 17070388 bytes, checksum: 46702d837c8c304ffc379088625742aa (MD5) / Made available in DSpace on 2018-03-23T15:26:03Z (GMT). No. of bitstreams: 1 MarconioSilvaDosSantos_TESE.pdf: 17070388 bytes, checksum: 46702d837c8c304ffc379088625742aa (MD5) Previous issue date: 2017-11-24 / Neste trabalho, ? feita uma modelagem estat?stica da precipita??o pluvial. Este ? um trabalho metodol?gico que utiliza simula??es estoc?sticas para estimar as distribui??es de probabilidades envolvidas na modelagem dessa vari?vel atmosf?rica. A fim de estimar os par?metros dessas distribui??es, foram utilizados m?todos de Monte Carlo via cadeias de Markov para gerar amostras sint?ticas de tamanho grande a partir de dados observados. Os m?todos utilizados foram o algoritmo de Metropolis-Hastings e o amostrador de Gibbs. As simula??es foram feitas sob a hip?tese de que os dias de um mesmo per?odo do ano (m?s ou esta??o chuvosa) podem ser considerados como identicamente distribu?dos em rela??o ? probabilidade de ocorrer precipita??o. Essa pesquisa possibilitou a produ??o de quatro artigos. O primeiro artigo utilizou o algoritmo de Metropolis-Hastings para modelar a probabilidade de ocorr?ncia de precipita??o em um dia qualquer do m?s. As simula??es desse artigo foram feitas com dados observados de algumas cidades brasileiras. Os demais artigos utilizaram o amostrador de Gibbs e os m?todos propostos foram aplicados em cidades da regi?o Nordeste do Brasil. No segundo artigo, as distribui??es Beta e Binomial foram utilizadas para modelar o n?mero de dias do m?s com ocorr?ncia de precipita??o. No terceiro artigo, a distribui??o de Poisson foi utilizada para modelar o n?mero de dias com valores extremos de precipita??o na esta??o chuvosa. Um m?todo alternativo para estimar esses valores extremos e sua distribui??o ? apresentado no quarto artigo, utilizando a distribui??o Gama. De acordo com os resultados dessas pesquisas, amostrador de Gibbs foi considerado adequado para estimar as distribui??es na modelagem da precipita??o em cidades para as quais h? poucos dados hist?ricos. / In this work, it was made a statistical modeling of precipitation. This is a methodological work that uses stochastic simulations to estimate the probability distributions related to this atmospheric variable. In order to estimate the parameters of these distributions, Markov chain Monte Carlo methods were used to generate large size synthetic samples from observed data. The used methods were the Metropolis-Hastings algorithm and the Gibbs sampler. The simulations were performed under the hypothesis that the days of of the same period of the year (month or rainy season) can be considered to be identically distributed concernig the probability of precipitation. This research allowed the production of four papers. The first paper used the Metropolis-Hastings algorithm to model the probability of occurrence of precipitation on any day of the month. The simulations of this paper were perfomed with observed data of some Brazilian cities. The other papers used the Gibbs sampler and the proposed methods were applied to data from cities in the Northeast Brazil. In the second paper, Beta and Binomial distributions were used to model the number of days of the month with occurrence of precipitation. In the third paper, the Poisson distribution was used to model the number of days with precipitation extreme values in the rainy season. An alternative method for estimating these extreme values and their distribution is presented in the fourth paper, using the Gamma distribution. According to the results obtained by these researches, the Gibbs sampler was considered to be adequate to estimate distributions in the modeling of precipitation on cities for which there are few historical data.

A presente tese aborda aspectos teóricos e aplicados da estimação dos parâmetros do modelo de calibração linear com erros distribuídos conforme a distribuição normal-assimétrica (Azzalini, 1985) e t-normal-assimétrica (Gómez, Venegas e Bolfarine, 2007). Aplicando um modelo assimétrico, não é necessário transformar as variáveis a fim de obter erros simétricos. A estimação dos parâmetros e das variâncias dos estimadores do modelo de calibração foram estudadas através da visão freqüentista e bayesiana, desenvolvendo algoritmos tipo EM e amostradores de Gibbs, respectivamente. Um dos pontos relevantes do trabalho, na óptica freqüentista, é a apresentação de uma reparametrização para evitar a singularidade da matriz de informação de Fisher sob o modelo de calibração normal-assimétrico na vizinhança de lambda = 0. Outro interessante aspecto é que a reparametrização não modifica o parâmetro de interesse. Já na óptica bayesiana, o ponto forte do trabalho está no desenvolvimento de medidas para verificar a qualidade do ajuste e que levam em consideração a assimetria do conjunto de dados. São propostas duas medidas para medir a qualidade do ajuste: o ADIC (Asymmetric Deviance Information Criterion) e o EDIC (Evident Deviance Information Criterion), que são extensões da ideia de Spiegelhalter et al. (2002) que propôs o DIC ordinário que só deve ser usado em modelos simétricos. / This thesis focuses on theoretical and applied estimation aspects of the linear calibration model with skew-normal (Azzalini, 1985) and skew-t-normal (Gómez, Venegas e Bolfarine, 2007) error distributions. Applying the asymmetrical distributed error methodology, it is not necessary to transform the variables in order to have symmetrical errors. The frequentist and the Bayesian solution are presented. The parameter estimation and its variance estimation were studied using the EM algorithm and the Gibbs sampler, respectively, in each approach. The main point, in the frequentist approach, is the presentation of a new parameterization to avoid singularity of the information matrix under the skew-normal calibration model in a neighborhood of lambda = 0. Another interesting aspect is that the reparameterization developed to make the information matrix nonsingular, when the skewness parameter is near to zero, leaves the parameter of interest unchanged. The main point, in the Bayesian framework, is the presentation of two measures of goodness-of-fit: ADIC (Asymmetric Deviance Information Criterion) and EDIC (Evident Deviance Information Criterion ). They are natural extensions of the ordinary DIC developed by Spiegelhalter et al. (2002).

A presente tese aborda aspectos teóricos e aplicados da estimação dos parâmetros do modelo de calibração linear com erros distribuídos conforme a distribuição normal-assimétrica (Azzalini, 1985) e t-normal-assimétrica (Gómez, Venegas e Bolfarine, 2007). Aplicando um modelo assimétrico, não é necessário transformar as variáveis a fim de obter erros simétricos. A estimação dos parâmetros e das variâncias dos estimadores do modelo de calibração foram estudadas através da visão freqüentista e bayesiana, desenvolvendo algoritmos tipo EM e amostradores de Gibbs, respectivamente. Um dos pontos relevantes do trabalho, na óptica freqüentista, é a apresentação de uma reparametrização para evitar a singularidade da matriz de informação de Fisher sob o modelo de calibração normal-assimétrico na vizinhança de lambda = 0. Outro interessante aspecto é que a reparametrização não modifica o parâmetro de interesse. Já na óptica bayesiana, o ponto forte do trabalho está no desenvolvimento de medidas para verificar a qualidade do ajuste e que levam em consideração a assimetria do conjunto de dados. São propostas duas medidas para medir a qualidade do ajuste: o ADIC (Asymmetric Deviance Information Criterion) e o EDIC (Evident Deviance Information Criterion), que são extensões da ideia de Spiegelhalter et al. (2002) que propôs o DIC ordinário que só deve ser usado em modelos simétricos. / This thesis focuses on theoretical and applied estimation aspects of the linear calibration model with skew-normal (Azzalini, 1985) and skew-t-normal (Gómez, Venegas e Bolfarine, 2007) error distributions. Applying the asymmetrical distributed error methodology, it is not necessary to transform the variables in order to have symmetrical errors. The frequentist and the Bayesian solution are presented. The parameter estimation and its variance estimation were studied using the EM algorithm and the Gibbs sampler, respectively, in each approach. The main point, in the frequentist approach, is the presentation of a new parameterization to avoid singularity of the information matrix under the skew-normal calibration model in a neighborhood of lambda = 0. Another interesting aspect is that the reparameterization developed to make the information matrix nonsingular, when the skewness parameter is near to zero, leaves the parameter of interest unchanged. The main point, in the Bayesian framework, is the presentation of two measures of goodness-of-fit: ADIC (Asymmetric Deviance Information Criterion) and EDIC (Evident Deviance Information Criterion ). They are natural extensions of the ordinary DIC developed by Spiegelhalter et al. (2002).

Made available in DSpace on 2015-03-26T13:55:08Z (GMT). No. of bitstreams: 1 texto completo.pdf: 905318 bytes, checksum: 167ccc3c1b47051e3ce28eb0224bed43 (MD5) Previous issue date: 2006-07-31 / Conselho Nacional de Desenvolvimento Científico e Tecnológico / First lactation data of 6.262 Brown-Swiss cows from 311 herds, daughters of 803 sires with calving between 1980 and 2003 were used to estimate genetic parameters for milk and fat production traits. The components of variance were estimated by restricted maximum likelihood (REML) and bayesian methods, using animal model with uni and two-traits analisys . The estimation by REML was obtained with the software MTDFREML (BOLDMAN et al. 1995) testing unitrait models with different effects to covariables and considering contemporary group and season as fixed effect. The best fitting obtained on unitrait analisys were used on two-trait analisys. The estimative of additive variance was reduced when lactation length was included on the model suggesting that the animals were been fitted to the same base on the capacity of transmit a longer or shorter lactation length to the progeny. Therefore, fitting to this covariable is not recommended. On the other side, the age of calving has linearly influenced milk and fat production. The heritability estimates were 0,26 and 0,25 to milk and fat yield respectively with genetic correlation of 0,95. the high correlation among these traits suggests that part of genes that acts on milk yield also respond to fat yield, in such way that selection for milk yield results, indirectly, in increase on fat yield. The estimation by Bayesian inference was made on software MTGSAM (VAN TASSELL E VAN VLECK, 1995). Chain lengths were tested to obtain the marginal posterior densities of unitrait analisys, the best option of chain length, burn-in and sampling interval was used on two-trait analisys. The burn-in periods were tested with the software GIBANAL (VAN KAAM, 1998) witch analysis inform a sampling interval for each burn-in tested, the criteria for choosing the sampling interval was made with the serial correlation resulting by burn-in and sampling process. The heritability estimates were 0,33 ± 0,05 for both traits with genetic correlation of 0,95. Similar results were obtained on studies using the same methodology on first lactation records. The stationary phase adequately reached with a 500.000 chain length and 30.000 burn-in iteractions. / Dados de primeira lactação de 6.262 vacas distribuídas em 311 rebanhos, filhas de 803 touros com partos entre os anos de 1980 e 2003 foram utilizados para estimar de componentes de variância para as características de produção de leite e gordura com informações de primeira lactação, em animais da raça Pardo-Suíça. Os componentes de variância foram estimados pelo método da máxima verossimilhança restrita (REML) e Bayesiano, sob modelo animal, por meio de análises uni e bicaracterística. A estimação realizada via REML foi obtida com o programa MTDFREML (BOLDMAN et al. 1995) testando modelos unicaracterística com diferentes efeitos para as covariáveis e considerados grupo contemporâneo e estação como efeitos fixos. Os melhores ajustes obtidos nas analises unicaracterística foram utilizados na análise bicaracterística. A duração da lactação reduziu a estimativa da variância aditiva quando era utilizada no modelo sugerindo que os animais estariam sendo corrigidos para uma mesma base quanto à capacidade de imprimir duração da lactação mais longa ou mais curta à progênie sendo, portanto, não recomendado o ajuste para esta covariável. Já a idade da vaca ao parto, influenciou linearmente a produção de leite e gordura. As herdabilidades estimadas foram 0,26 e 0,25 para produção de leite e gordura respectivamente com correlação genética de 0,95. A alta correlação entre a produção de leite e gordura obtida sugere que parte dos genes que atuam na produção de leite também responde pela produção de gordura, de tal forma que a seleção para a produção de leite resulta, indiretamente, em aumentos na produção de gordura. A estimação via inferência Bayesiana foi realizada com o programa MTGSAM (VAN TASSELL E VAN VLECK, 1995). Foram testados diversos tamanhos de cadeia para a obtenção das densidades marginais a posteriori das análises unicaracterística, a melhor proposta para o tamanho de cadeia, burn-in e amostragem foi utilizada para a análise bicaracterística. Os períodos de burn-in foram testados pelo programa GIBANAL (VAN KAAM, 1998) cujas análises fornecem um intervalo de amostragem para cada burn-in testado, o critério de escolha do intervalo de amostragem foi feito de acordo com a correlação serial, resultante do burn-in e do processo de amostragem. As estimativas de herdabilidade obtidas foram 0,33 ± 0,05 para ambas as características com correlação de 0,95. Resultados similares foram obtidos em estudos utilizando a mesma metodologia em informações de primeira lactação. A fase estacionária foi adequadamente atingida com uma cadeia de 500.000 iterações e descarte inicial de 30.000 iterações.