Performance of balanced two-stage empirical predictors of realized cluster latent values from finite populations: A simulation study

Predictors of random effects are usually based on the popular mixed effects (ME) model developed under the assumption that the sample is obtained from a conceptual infinite population; such predictors are employed even when the actual population is finite. Two alternatives that incorporate the finite nature of the population are obtained from the superpopulation model proposed by Scott and Smith (1969. Estimation in multi-stage surveys. J. Amer. Statist. Assoc. 64, 830-840) or from the finite population mixed model recently proposed by Stanek and Singer (2004. Predicting random effects from finite population clustered samples with response error. J. Amer. Statist. Assoc. 99, 1119-1130). Predictors derived under the latter model with the additional assumptions that all variance components are known and that within-cluster variances are equal have smaller mean squared error (MSE) than the competitors based on either the ME or Scott and Smith`s models. As population variances are rarely known, we propose method of moment estimators to obtain empirical predictors and conduct a simulation study to evaluate their performance. The results suggest that the finite population mixed model empirical predictor is more stable than its competitors since, in terms of MSE, it is either the best or the second best and when second best, its performance lies within acceptable limits. When both cluster and unit intra-class correlation coefficients are very high (e.g., 0.95 or more), the performance of the empirical predictors derived under the three models is similar. (c) 2007 Elsevier B.V. All rights reserved.

Novel Statistical Methods in Quantitative Genetics : Modeling Genetic Variance for Quantitative Trait Loci Mapping and Genomic Evaluation

This thesis develops and evaluates statistical methods for different types of genetic analyses, including quantitative trait loci (QTL) analysis, genome-wide association study (GWAS), and genomic evaluation. The main contribution of the thesis is to provide novel insights in modeling genetic variance, especially via random effects models. In variance component QTL analysis, a full likelihood model accounting for uncertainty in the identity-by-descent (IBD) matrix was developed. It was found to be able to correctly adjust the bias in genetic variance component estimation and gain power in QTL mapping in terms of precision.  Double hierarchical generalized linear models, and a non-iterative simplified version, were implemented and applied to fit data of an entire genome. These whole genome models were shown to have good performance in both QTL mapping and genomic prediction. A re-analysis of a publicly available GWAS data set identified significant loci in Arabidopsis that control phenotypic variance instead of mean, which validated the idea of variance-controlling genes.  The works in the thesis are accompanied by R packages available online, including a general statistical tool for fitting random effects models (hglm), an efficient generalized ridge regression for high-dimensional data (bigRR), a double-layer mixed model for genomic data analysis (iQTL), a stochastic IBD matrix calculator (MCIBD), a computational interface for QTL mapping (qtl.outbred), and a GWAS analysis tool for mapping variance-controlling loci (vGWAS).

Genetic heterogeneity of residual variance : estimation of variance components using double hierarchical generalized linear models

Background: The sensitivity to microenvironmental changes varies among animals and may be under genetic control. It is essential to take this element into account when aiming at breeding robust farm animals. Here, linear mixed models with genetic effects in the residual variance part of the model can be used. Such models have previously been fitted using EM and MCMC algorithms. Results: We propose the use of double hierarchical generalized linear models (DHGLM), where the squared residuals are assumed to be gamma distributed and the residual variance is fitted using a generalized linear model. The algorithm iterates between two sets of mixed model equations, one on the level of observations and one on the level of variances. The method was validated using simulations and also by re-analyzing a data set on pig litter size that was previously analyzed using a Bayesian approach. The pig litter size data contained 10,060 records from 4,149 sows. The DHGLM was implemented using the ASReml software and the algorithm converged within three minutes on a Linux server. The estimates were similar to those previously obtained using Bayesian methodology, especially the variance components in the residual variance part of the model. Conclusions: We have shown that variance components in the residual variance part of a linear mixed model can be estimated using a DHGLM approach. The method enables analyses of animal models with large numbers of observations. An important future development of the DHGLM methodology is to include the genetic correlation between the random effects in the mean and residual variance parts of the model as a parameter of the DHGLM.

Estudo do processo de combustão in-situ usando poços horizontais como produtores de óleo (Toe-to-Hell Air Injection)

The method "toe-to-heel air injection" (THAITM) is a process of enhanced oil recovery, which is the integration of in-situ combustion with technological advances in drilling horizontal wells. This method uses horizontal wells as producers of oil, keeping vertical injection wells to inject air. This process has not yet been applied in Brazil, making it necessary, evaluation of these new technologies applied to local realities, therefore, this study aimed to perform a parametric study of the combustion process with in-situ oil production in horizontal wells, using a semi synthetic reservoir, with characteristics of the Brazilian Northeast basin. The simulations were performed in a commercial software "STARS" (Steam, Thermal, and Advanced Processes Reservoir Simulator), from CMG (Computer Modelling Group). The following operating parameters were analyzed: air rate, configuration of producer wells and oxygen concentration. A sensitivity study on cumulative oil (Np) was performed with the technique of experimental design, with a mixed model of two and three levels (32x22), a total of 36 runs. Also, it was done a technical economic estimative for each model of fluid. The results showed that injection rate was the most influence parameter on oil recovery, for both studied models, well arrangement depends on fluid model, and oxygen concentration favors recovery oil. The process can be profitable depends on air rate

Fiscalização difusa da constitucionalidade e o incidente de argüição de inconstitucionalidade

The present work consists of studying to diffuse control of constitutionality in Brazil, with emphasis in a procedural alternative to the evolution of that model: the incident to challenge of unconstitutionality. Starting from the discussion about the new role of constitutional jurisdiction in peripheral countries and in the globalized society, without forgetting to face inevitable doubts about its legitimacy before other powers of the State, the Brazilian control of constitutionality is revealed, under a diffuse, non-dichotomical view, through a number of inconsistencies and misunderstandings, that compromise social peace, the credibility of democratic institutions and the supremacy of juridical security. In order to achieve the goal, the study in course discussed the main difficulties of the Brazilian mixed model of constitutionality control, as well as, directing its view to the incident of challenge of unconstitutionality, which the most adequate forms to assure its appropriateness, legitimacy, processing and decisory effects are. Is was essential, in this point of view, to establish the difference between the incident of challenge of unconstitutionality conceived in article of the Brazilian Federal Constitution and the incident of challenge of unconstitutionality such as it is known in the European models. The insertion of the incident of challenge of unconstitutionality based on European models in the Brazilian control system, without jeopardizing the North-American essence the Brazilian constitutional history presents since 1981, is the hypothesis that is presented as an improvement of constitutional protection

Seleção para peso pós-desmama em um rebanho Gir. 1. Resposta direta e correlacionada em 21 anos de seleção

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

Use of direct and iterative solvers for estimation of SNP effects in genome-wide selection

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

Resposta à Seleção para Características de Desempenho em um Rebanho de Seleção de Suínos

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

Tendências genéticas para escores visuais de conformação, precocidade e musculatura à desmama de bovinos Nelore

Os objetivos neste estudo foram estimar os coeficientes de herdabilidade e a mudança genética nos escores visuais de conformação (C), precocidade (P) e musculatura (M) à desmama de bovinos Nelore e avaliar as mudanças promovidas pelo programa de seleção ao qual os animais foram submetidos. Foram utilizados dados de 56.076 animais nascidos entre 1990 e 2002. Os componentes de variância foram estimados por máxima verossimilhança restrita e os valores genéticos foram preditos pelo método dos modelos mistos aplicando-se um modelo animal. As tendências genéticas foram estimadas pela regressão dos valores genéticos sobre o ano de nascimento dos animais. Os coeficientes de herdabilidade do efeito direto estimados foram de 0,122±0,012; 0,153±0,012; e 0,124±0,012 e, para o efeito materno, de 0,043±0,010; 0,038±0,009; e 0,056±0,010 para C, P e M, respectivamente. As tendências genéticas diretas estimadas foram de 0,013; 0,022 e 0,018 pontos de escore ao ano para C, P e M, respectivamente, o que representa incremento anual de 0,42; 0,67 e 0,60% ao ano em relação à média fenotípica. As tendências genéticas do efeito materno foram de -0,0020; -0,0010; e -0,0009 pontos de escore ao ano para C, P e M, respectivamente, o que representa uma mudança anual de -0,07; -0,03; e -0,03% na média fenotípica da população.

Altered plasma response to zinc and iron tolerance test after Roux-en-Y gastric bypass

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

Ganho genético com base no tamanho efetivo populacional de progênies de seringueira

O objetivo deste trabalho foi estimar os ganhos genéticos de um teste de progênies de seringueira para a produção de borracha seca e, com base no maior tamanho efetivo populacional e maior ganho genético, obter os melhores indivíduos. Foram utilizadas 30 progênies de meios-irmãos, provenientes de sementes de polinização mista - alogamia e autogamia - de testes clonais no Estado de São Paulo. Utilizou-se o delineamento experimental de blocos ao acaso, com 30 tratamentos (progênies), 3 repetições e parcelas lineares de 10 plantas, em um espaçamento de 3x3 m, o que totalizou 900 plantas úteis. Aos três anos, o perímetro, a 50 cm do solo (PA50), e a produção de borracha seca (PBS) foram avaliadas por meio do teste precoce de produção Hamaker Morris-Mann (HMM). As variáveis foram analisadas pelo método de modelo linear misto, via procedimento REML/BLUP, em progênies com sistema reprodutivo misto e taxa de autofecundação de 22%. A identificação dos 20 melhores indivíduos quanto à PBS e ao PA50 proporcionou ganho genético de 67,96 e 16,48%, respectivamente, e um coeficiente de endogamia de aproximadamente 2,82%. O teste de progênies proporciona produção de sementes com melhor valor genético, grande variabilidade e baixa endogamia

Incidence and prevalence rates of inflammatory bowel diseases, in midwestern of São Paulo State, Brazil

CONTEXTO: No Brasil, a incidência e prevalência populacionais das doenças inflamatórias intestinais são desconhecidas. OBJETIVO: Neste trabalho, estimou-se esses parâmetros na área que abrange a antiga região de saúde DIR 11, no centro-oeste do Estado de São Paulo. MÉTODOS: Usou-se um registro sequencial de 115 pacientes (>15 anos de idade) com doenças inflamatórias intestinais, residindo na área de estudo, atendidos durante período de 20 anos (1986-2005) em hospital de referência. Estimou-se, para quatro períodos consecutivos de 5 anos, as incidências e prevalências de acordo com os tipos de doença e do sexo de doentes. As suas inter-relações foram analisadas utilizando o modelo de regressão linear de Poisson. RESULTADOS: Na região, as doenças inflamatórias intestinais foram predominantes em mulheres jovens, da raça branca, residindo na zona urbana.A incidência da retocolite foi maior que a da doença de Crohn e das colites não classificadas, e diferentemente dessas duas ultimas, mostrou decréscimo no último período (2001-2005). Neste mesmo período, a taxa da incidência para a retocolite foi de 4,48 casos/100.000 habitantes, para a doença de Crohn atingiu 3,50 casos/100.000 habitantes e a das colites não classificadas foi de 1,75 casos/100.000 habitantes. As prevalências atingiram os valores de 14,81 casos/100.000 habitantes para a retocolite, 5,65 casos/100.000 habitantes para a doenças de Crohn e 2,14 casos/100.000 habitantes. Considerando todas as doenças inflamatórias a prevalência atingiu o valor de 22,61 casos/100.000 habitantes. CONCLUSÃO:A incidência dessas doenças inflamatórias intestinais nessa região é baixa, igualando-se aos países da América Latina e do sul da Europa e sua crescente prevalência justifica políticas de saúde para as suas abordagens.

Peritoneal Dialysis-Related Peritonitis Due to Staphylococcus aureus: A Single-Center Experience over 15 Years

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

Analysis of the electrical conductivity in milking fractions as a mean for detecting and characterizing mastitis in goats

The aim of the work was to study the effect of milking fraction on electrical conductivity of milk (EC) to improve its use in dairy goat mastitis detection using automatic EC measurements during milking. The experiment was carried out on a group of 84 Murciano-Granadina goats (28 primiparous and 56 multiparous). Goats were in the fourth month of lactation. A linear mixed model was used to analyse the relationship between EC or somatic cell count (SCC) of gland milk and parity, mammary gland health status, analysed fraction (first 100 mL=F-1; machine milk=F-2; and stripping milk=F-3) and their first order interactions. Additionally, the mastitis detection characteristics (sensitivity, specificity, positive predictive value and negative predictive value) of SCC and EC were studied at different thresholds.All factors considered were significant for EC and SCC. EC decreased significantly as milking progressed (from F-1 to F-3) in both healthy and infected glands. EC was not significantly different between healthy and infected glands in F-1 and F-2 fractions, but EC of healthy glands (5.01 mS/cm) was significantly lower than in infected glands (5.03 mS/cm) at F-3.Mastitis detection characteristics of EC did not differ amongst studied fractions. The small significant difference of EC between healthy and infected glands obtained in F-3 fraction did not yield better sensitivity results compared to F-1 and F-2. The best EC mastitis detection characteristics were obtained at 5.20 mS/cm threshold (sensitivity of 70% and specificity of 50%). The best SCC mastitis detection characteristics were obtained at 300,000 cells/mL threshold and F-3 fraction (sensitivity of 85% and specificity of 65%).It was concluded that mastitis detection characteristics of EC were similar in the three milking fractions analysed, being slightly better for SCC in F-3 fraction. As shown in previous studies, there are no factors other than the mammary gland health status that affect milk EC and should be considered in the algorithms for mastitis detection to improve the results. (C) 2012 Elsevier B.V. All rights reserved.

Parâmetros e tendência genética da produção de leite de cabra no Brasil

A caprinocultura leiteira no Brasil, apesar de ser uma atividade rural consolidada há algumas décadas, tem se mostrado totalmente dependente de outros países no que se refere ao melhoramento genético. A maioria dos plantéis existentes atualmente tem como base animais importados, e a renovação do material genético é feita por meio da importação de sêmen. Inexistem informações sobre o valor genético dos animais e sua evolução no decorrer dos anos. No presente trabalho, foram estimadas a herdabilidade e a repetibilidade da produção de leite utilizando o REML. Os valores obtidos foram 0,21557 e 0,21564, respectivamente. Para a predição do valor gênico dos animais, foi usado o procedimento BLUP com modelo animal. A mudança na tendência genética anual estimada por um modelo quadrático foi -0,8109 kg/ano², indicando desaceleração no ganho genético. A correlação de Pearson entre os valores gênicos dos bodes estimados com base na média da capacidade provável de produção das filhas obtida pelo método de mínimos quadrados com as estimadas pelas equações do modelo misto foi de 0,5751. A correlação de SPEARMAN entre as classificações dos bodes obtidos pelos dois métodos foi de 0,5813.