Estimation and prediction of parameters and breeding values in soybean using REML/BLUP and Least Squares

The aim of this study was to compare REML/BLUP and Least Square procedures in the prediction and estimation of genetic parameters and breeding values in soybean progenies. F(2:3) and F(4:5) progenies were evaluated in the 2005/06 growing season and the F(2:4) and F(4:6) generations derived thereof were evaluated in 2006/07. These progenies were originated from two semi-early, experimental lines that differ in grain yield. The experiments were conducted in a lattice design and plots consisted of a 2 m row, spaced 0.5 m apart. The trait grain yield per plot was evaluated. It was observed that early selection is more efficient for the discrimination of the best lines from the F(4) generation onwards. No practical differences were observed between the least square and REML/BLUP procedures in the case of the models and simplifications for REML/BLUP used here.

Componentes de variância e predição de valores genéticos em seringueira pela metodologia de modelos mistos (REML/BLUP)

Este trabalho objetivou estimar parâmetros genéticos e predizer os valores genéticos aditivos e genotípicos de clones amazônicos de seringueira em relação a características de qualidade da borracha e do látex a partir de genealogia envolvendo 15 clones e 250 rametes. As estimativas de herdabilidade no sentido amplo, , foram de 0,4152 e 0,6357, com relação a teores de borracha seca e plasticidade após 30 minutos, respectivamente, o que torna possível obter altos ganhos genéticos via seleção clonal. Os clones IAN 873 e Fx 3899 apresentaram os maiores valores genéticos aditivos e genotípicos em relação a teores de borracha seca, enquanto o clone IAN 6158 destacou-se quanto aos valores genotípicos relativos às características de plasticidade após 30 minutos, e índice de retenção da plasticidade.

Ganho genético avaliado com índices de seleção e com REML/Blup em milho-pipoca

O objetivo deste trabalho foi comparar quatro índices de seleção e o método REML/Blup na avaliação de ganhos genéticos preditos das características de interesse ao programa de melhoramento do milho-pipoca UENF 14. Foram avaliadas 200 progênies de irmãos-completos, em uma safra, em dois ambientes (regiões Norte e Noroeste do Estado do Rio de Janeiro), com duas repetições, o que totalizou 800 observações. Avaliaram-se as características altura média de plantas, altura média de inserção da primeira espiga, estande final, tombamento, diâmetro de colmo, prolificidade, rendimento e capacidade de expansão dos grãos. Entre os índices de seleção testados, o de Mulamba & Mock proporcionou os melhores resultados para a seleção das progênies de irmãos-completos. O método REML/Blup mostrou-se muito eficiente, tendo selecionado progênies com desempenhos relativos elevados e com ganhos genéticos preditos melhores que os dos índices de seleção testados. As progênies selecionadas com uso do método REML/Blup, para maior rendimento de grãos, não são as mesmas selecionadas para maior capacidade de expansão. O índice de Mulamba & Mock e o método REML/Blup proporcionam os melhores ganhos preditos para a população UENF 14.

APLICAÇÃO DA METODOLOGIA DE MODELOS MISTOS (REML/BLUP) NA ESTIMAÇÃO DE COMPONENTES DE VARIÂNCIA E PREDIÇÃO DE VALORES GENÉTICOS EM PUPUNHEIRA (Bactris gasipaes)

No Brasil, a pupunheira é uma planta muito útil na alimentação, seja como fonte de frutos ou de palmito. O interesse pela pupunheira, além de ser uma cultura perene, são: crescimento a pleno sol, precocidade, rusticidade, perfilhamento, palatabilidade e não-escurecimento do palmito após o corte. Estimativas de parâmetros genéticos em pupunheira são escassas e se constituem em ferramenta de suma importância para orientar os programas de melhoramento. O objetivo deste trabalho foi estudar a variabilidade genética e estimar o valor genético individual como critério de seleção, usando o procedimento BLUP/REML (Melhor predição linear não viciada/máxima verossimilhança restrita). Adotaram-se duas estratégias de seleção para o caráter produção de palmito: simulando programa de melhoramento a curto prazo (CP- seleção das 9 famílias com 31 indivíduos de maior valor genético) e a longo prazo (LP- seleção das 15 famílias com 53 indivíduos). As progênies foram avaliadas em experimento delineado em blocos ao acaso com três repetições, parcelas lineares de cinco plantas, espaçadas de 2,0 m x 1,0 m e bordadura composta por uma fileira em torno do experimento no Campo Experimental do Matapí, Município de Porto Grande, Estado do Amapá. A avaliação foi realizada aos 26 meses pós- plantio (2ª avaliação), coletando-se dados de altura da planta (AP), diâmetro da planta à altura do colo (DPC), tamanho do palmito (TP), diâmetro do palmito (DP), peso do resíduo apical (PRA), basal (PRB) e do palmito líquido (PP) (tipo exportação). Os dados de AP, DPC, TP e DP corresponderam às médias das touceiras que apresentavam mais de uma haste. Já para os caracteres PRA, PRB e PP corresponderam à soma das hastes na touceira. De maneira geral, a população estudada apresenta baixa variabilidade genética. As herdabilidades no sentido restrito em nível de indivíduos foram: AP (18,44%), DPC (3,16%), TP (42,47%), DP (10,54%), PP (5,70%), PRB (6,15%). Os ganhos genéticos preditos em relação à média da população para PP foram de 7,18% na situação de LP e 8,40% para CP, com tamanho efetivo de 30,38 e 19,00, respectivamente.

Variabilidade genética de procedências e progênies de umbuzeiro via metodologia de modelos lineares mistos (REML/BLUP)

O umbuzeiro (Spondias tuberosa Arruda), espécie nativa do semi-árido brasileiro, é um dos recursos genéticos mais importantes para a produção de frutos destinados ao consumo in natura e à industrialização. No melhoramento de espécies perenes, o uso de técnicas de avaliação genética com base em modelos mistos do tipo REML/BLUP (máxima verossimilhança restrita /melhor predição linear não viciada) é fundamental para a predição de valores genéticos aditivos e genotípicos de indivíduos com potencial para seleção, tanto em nível intrapopulacional como interpopulacional. Este procedimento vem sendo aplicado com sucesso no Brasil, no melhoramento de várias espécies frutíferas e florestais. O presente trabalho teve como objetivos estudar a variabilidade genética, estimar parâmetros genéticos e realizar a predição de valores genéticos dos indivíduos, utilizando a metodologia REML/BLUP a partir da avaliação de procedências e progênies de umbuzeiro. O ensaio foi instalado em delineamento de blocos ao acaso, com arranjo hierárquico desbalanceado, constituído de três procedências e 42 progênies. Foram avaliados os caracteres altura de plantas (ALP), maior diâmetro de copa (MAC), menor diâmetro de copa (MEC), diâmetro do colo (DIC) e número de ramos primários (NRP), os quais são correlacionados à produção de frutos. Os resultados obtidos permitem afirmar que, no umbuzeiro, aos nove anos de idade, a maior parte da variabilidade genética encontrada concentra-se dentro de populações. Os caracteres MEC e MAC apresentaram, respectivamente, valores de herdabilidade individual no sentido restrito de 0,08 e 0,14 e ganhos de 6% e 9%, respectivamente, com a seleção dos dez melhores genitores.

Seleção de genitores de Eucalyptus grandis e de Eucalyptus urophylla para produção de híbridos interespecíficos utilizando REML/BLUP e informação de divergência genética

Foram avaliadas 363 progênies de meios-irmãos de Eucalyptus urophylla e de 245 progênies de meios-irmãos de E. grandis, visando à seleção de genitores para a produção de híbridos interespecíficos em cruzamentos controlados em dialelo circulante, com base em seus valores genéticos preditos. Os ensaios foram conduzidos em Guanhães, MG, no delineamento de blocos casualizados com cinco repetições, parcelas lineares de seis plantas para E. urophylla e oito plantas para E. grandis, no espaçamento de 3,0 x 2,0 m. Os parâmetros genéticos foram estimados adotando-se modelos mistos com uso do procedimento REML/BLUP (máxima verossimilhança restrita/melhor predição linear não viesada) para diâmetro à altura do peito (DAP), altura total (ALT) e volume individual (VOL). A predição dos valores genéticos foi feita somente para (DAP), característica que foi utilizada para seleção. Em E. urophylla, foram obtidas as estimativas de herdabilidade no sentido restrito de 0,2785 e de acurácia da ordem de 61%, e em E. grandis, herdabilidade de 0,2247 e acurácia de 53%. A seleção dos 100 genitores de cada espécie individualmente proporcionará ganhos genéticos na ordem de 20,6% em E. urophylla e de 16,4% em E. grandis. Ganho genético para DAP em torno de 35,0% poderá ser obtido com o cruzamento dos 10 indivíduos selecionados para cada espécie, portadores das maiores estimativas de divergência genética.

Estimation of variance components and prediction of breeding values in rubber tree breeding using the REML/BLUP procedure

The present paper deals with estimation of variance components, prediction of breeding values and selection in a population of rubber tree [Hevea brasiliensis (Willd. ex Adr. de Juss.) Müell.-Arg.] from Rio Branco, State of Acre, Brazil. The REML/BLUP (restricted maximum likelihood/best linear unbiased prediction) procedure was applied. For this purpose, 37 rubber tree families were obtained and assessed in a randomized complete block design, with three unbalanced replications. The field trial was carried out at the Experimental Station of UNESP, located in Selvíria, State of Mato Grosso do Sul, Brazil. The quantitative traits evaluated were: girth (G), bark thickness (BT), number of latex vessel rings (NR), and plant height (PH). Given the unbalanced condition of the progeny test, the REML/BLUP procedure was used for estimation. The narrow-sense individual heritability estimates were 0.43 for G, 0.18 for BT, 0.01 for NR, and 0.51 for PH. Two selection strategies were adopted: one short-term (ST - selection intensity of 8.85%) and the other long-term (LT - selection intensity of 26.56%). For G, the estimated genetic gains in relation to the population average were 26.80% and 17.94%, respectively, according to the ST and LT strategies. The effective population sizes were 22.35 and 46.03, respectively. The LT and ST strategies maintained 45.80% and 28.24%, respectively, of the original genetic diversity represented in the progeny test. So, it can be inferred that this population has potential for both breeding and ex situ genetic conservation as a supplier of genetic material for advanced rubber tree breeding programs. Copyright by the Brazilian Society of Genetics.

Seleção de famílias e progênies de irmãos completos de cana-de-açúcar para atributos tecnológicos e de produção pelo método de REML/BLUP

The objective of this study was to assess families and highlight the superior progenies of sugarcane originating from 38 biparental crosses for the following attributes: tons of cane per hectare (TCH), tons of biomass per hectare (TBIOH), brix (% cane juice), fiber content, purity, pol and total recoverable sugar (TRS). The data were analyzed by mixed model REML / BLUP in the REML (Restricted Maximum Likelihood) allowed us to estimate genetic parameters and BLUP (best linear unbiased prediction) to predict the additive and genotypic values. The best family for the attributes TCH and TBIOH was 41, whose parents are cultivars IACSP022019 x CTC9. In individual selection for TCH, the plant number 3 of Block 2, the crossing 78, showed the best results. To TBIOH the plant number 33, Block 1, family 41, showed the best results. Families 40, 41, 43, 68, 69, 79, 91, 92 and 147, were higher for the variables brix, pol, purity, and ATR, where as 85 families, 147, 148, 149, 161, 163, 177, 178, 179, and 183 were higher for fiber. The family 147 whose parents are IACSP042286 x IACSP963055, showed three progenies ranked among the top ten for both brix, and for fiber, which identifies the combination as a potential source of progenies for bioenergy production.

REML/BLUP para predição de valores genotípicos de topcrosses e seleção de testadores em milho

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

Modelling the Distribution of Ischaemic Stroke-Specific Survival Time Using an EM-based Mixture Approach with Random Effects Adjustment

A two-component survival mixture model is proposed to analyse a set of ischaemic stroke-specific mortality data. The survival experience of stroke patients after index stroke may be described by a subpopulation of patients in the acute condition and another subpopulation of patients in the chronic phase. To adjust for the inherent correlation of observations due to random hospital effects, a mixture model of two survival functions with random effects is formulated. Assuming a Weibull hazard in both components, an EM algorithm is developed for the estimation of fixed effect parameters and variance components. A simulation study is conducted to assess the performance of the two-component survival mixture model estimators. Simulation results confirm the applicability of the proposed model in a small sample setting. Copyright (C) 2004 John Wiley Sons, Ltd.

Genetic associations between carcass traits measured by real-time ultrasound and scrotal circumference and growth traits in Nelore cattle

The aim of the present study was to evaluate the genetic correlations among real-time ultrasound carcass, BW, and scrotal circumference (SC) traits in Nelore cattle. Carcass traits, measured by real-time ultrasound of the live animal, were recorded from 2002 to 2004 on 10 farms across 6 Brazilian states on 2,590 males and females ranging in age from 450 to 599 d. Ultrasound records of LM area (LMA) and backfat thickness (BF) were obtained from cross-sectional images between the 12th and 13th ribs, and rump fat thickness (RF) was measured between the hook and pin bones over the junction between gluteus medius and biceps femoris muscles. Also, BW (n = 22,778) and SC ( n = 5,695) were recorded on animals born between 1998 and 2003. The BW traits were 120, 210, 365, 450, and 550-d standardized BW (W120, W210, W365, W450, and W550), plus BW (WS) and hip height (HH) on the ultrasound scanning date. The SC traits were 365-, 450-, and 550-d standardized SC (SC365, SC450, and SC550). For the BW and SC traits, the database used was from the Nelore Breeding Program-Nelore Brazil. The genetic parameters were estimated with multivariate animal models and REML. Estimated genetic correlations between LMA and other traits were 0.06 (BF), -0.04 ( RF), 0.05 (HH), 0.58 (WS), 0.53 (W120), 0.62 (W210), 0.67 (W365), 0.64 ( W450 and W550), 0.28 (SC365), 0.24 (SC450), and 0.00 ( SC550). Estimated genetic correlations between BF and with other traits were 0.74 ( RF), -0.32 (HH), 0.19 (WS), -0.03 (W120), -0.10 (W210), 0.04 (W365), 0.01 (W450), 0.06 ( W550), 0.17 (SC365 and SC450), and -0.19 (SC550). Estimated genetic correlations between RF and other traits were -0.41 (HH), -0.09 (WS), -0.13 ( W120), -0.09 ( W210), -0.01 ( W365), 0.02 (W450), 0.03 (W550), 0.05 ( SC365), 0.11 ( SC450), and -0.18 (SC550). These estimates indicate that selection for carcass traits measured by real-time ultrasound should not cause antagonism in the genetic improvement of SC and BW traits. Also, selection to increase HH might decrease subcutaneous fat as correlated response. Therefore, to obtain animals suited to specific tropical production systems, carcass, BW, and SC traits should be considered in selection programs.

Long-term survivor mixture model with random effects: application to a multi-centre clinical trial of carcinoma

A mixture model incorporating long-term survivors has been adopted in the field of biostatistics where some individuals may never experience the failure event under study. The surviving fractions may be considered as cured. In most applications, the survival times are assumed to be independent. However, when the survival data are obtained from a multi-centre clinical trial, it is conceived that the environ mental conditions and facilities shared within clinic affects the proportion cured as well as the failure risk for the uncured individuals. It necessitates a long-term survivor mixture model with random effects. In this paper, the long-term survivor mixture model is extended for the analysis of multivariate failure time data using the generalized linear mixed model (GLMM) approach. The proposed model is applied to analyse a numerical data set from a multi-centre clinical trial of carcinoma as an illustration. Some simulation experiments are performed to assess the applicability of the model based on the average biases of the estimates formed. Copyright (C) 2001 John Wiley & Sons, Ltd.

A methodology for analysis of sugarcane productivity trends - I. Analysis across districts

Historically, few articles have addressed the use of district level mill production data for analysing the effect of varietal change on sugarcane productivity trends. This appears to be due to lack of compiled district data sets and appropriate methods by which to analyse these data. Recently, varietal data on tonnes of sugarcane per hectare (TCH), sugar content (CCS), and their product, tonnes of sugar content per hectare (TSH) on a district basis, have been compiled. This study was conducted to develop a methodology for regular analysis of such data from mill districts to assess productivity trends over time, accounting for variety and variety x environment interaction effects for 3 mill districts (Mulgrave, Babinda, and Tully) from 1958 to 1995. Restricted maximum likelihood methodology was used to analyse the district level data and best linear unbiased predictors for random effects, and best linear unbiased estimates for fixed effects were computed in a mixed model analysis. In the combined analysis over districts, Q124 was the top ranking variety for TCH, and Q120 was top ranking for both CCS and TSH. Overall production for TCH increased over the 38-year period investigated. Some of this increase can be attributed to varietal improvement, although the predictors for TCH have shown little progress since the introduction of Q99 in 1976. Although smaller gains have been made in varietal improvement for CCS, overall production for CCS decreased over the 38 years due to non-varietal factors. Varietal improvement in TSH appears to have peaked in the mid-1980s. Overall production for TSH remained stable over time due to the varietal increase in TCH and the non-varietal decrease in CCS.

Some vocabulary and grammar for the analysis of multi-environment trials, as applied to the analysis of FPB and PPB trials

For the improvement of genetic material suitable for on farm use under low-input conditions, participatory and formal plant breeding strategies are frequently presented as competing options. A common frame of reference to phrase mechanisms and purposes related to breeding strategies will facilitate clearer descriptions of similarities and differences between participatory plant breeding and formal plant breeding. In this paper an attempt is made to develop such a common framework by means of a statistically inspired language that acknowledges the importance of both on farm trials and research centre trials as sources of information for on farm genetic improvement. Key concepts are the genetic correlation between environments, and the heterogeneity of phenotypic and genetic variance over environments. Classic selection response theory is taken as the starting point for the comparison of selection trials (on farm and research centre) with respect to the expected genetic improvement in a target environment (low-input farms). The variance-covariance parameters that form the input for selection response comparisons traditionally come from a mixed model fit to multi-environment trial data. In this paper we propose a recently developed class of mixed models, namely multiplicative mixed models, also called factor-analytic models, for modelling genetic variances and covariances (correlations). Mixed multiplicative models allow genetic variances and covariances to be dependent on quantitative descriptors of the environment, and confer a high flexibility in the choice of variance-covariance structure, without requiring the estimation of a prohibitively high number of parameters. As a result detailed considerations regarding selection response comparisons are facilitated. ne statistical machinery involved is illustrated on an example data set consisting of barley trials from the International Center for Agricultural Research in the Dry Areas (ICARDA). Analysis of the example data showed that participatory plant breeding and formal plant breeding are better interpreted as providing complementary rather than competing information.

A zero-augmented gamma mixed model for longitudinal data with many zeros

In many occupational safety interventions, the objective is to reduce the injury incidence as well as the mean claims cost once injury has occurred. The claims cost data within a period typically contain a large proportion of zero observations (no claim). The distribution thus comprises a point mass at 0 mixed with a non-degenerate parametric component. Essentially, the likelihood function can be factorized into two orthogonal components. These two components relate respectively to the effect of covariates on the incidence of claims and the magnitude of claims, given that claims are made. Furthermore, the longitudinal nature of the intervention inherently imposes some correlation among the observations. This paper introduces a zero-augmented gamma random effects model for analysing longitudinal data with many zeros. Adopting the generalized linear mixed model (GLMM) approach reduces the original problem to the fitting of two independent GLMMs. The method is applied to evaluate the effectiveness of a workplace risk assessment teams program, trialled within the cleaning services of a Western Australian public hospital.