Skip to main content
Article thumbnail
Location of Repository

Model for fitting longitudinal traits subject to threshold response applied to genetic evaluation for heat tolerance

By Juan Pablo Sánchez, Romdhane Rekaya and Ignacy Misztal


A semi-parametric non-linear longitudinal hierarchical model is presented. The model assumes that individual variation exists both in the degree of the linear change of performance (slope) beyond a particular threshold of the independent variable scale and in the magnitude of the threshold itself; these individual variations are attributed to genetic and environmental components. During implementation via a Bayesian MCMC approach, threshold levels were sampled using a Metropolis step because their fully conditional posterior distributions do not have a closed form. The model was tested by simulation following designs similar to previous studies on genetics of heat stress. Posterior means of parameters of interest, under all simulation scenarios, were close to their true values with the latter always being included in the uncertain regions, indicating an absence of bias. The proposed models provide flexible tools for studying genotype by environmental interaction as well as for fitting other longitudinal traits subject to abrupt changes in the performance at particular points on the independent variable scale

Topics: Research
Publisher: BioMed Central
OAI identifier:
Provided by: PubMed Central
Download PDF:
Sorry, we are unable to provide the full text but you may find it at the following location(s):
  • http://www.pubmedcentral.nih.g... (external link)
  • Suggested articles


    1. (1998). AFM: Automatic Bayesian curve fitting.
    2. (2006). Bayesian analysis of the linear reaction norm model with unknown covariates.
    3. (2000). Carabaño MJ, Toro MA: Assessment of heterogeneity of residual variances using changepoint techniques. Genet Sel Evol
    4. (1953). E: Equations of state calculations by fast computing machines.
    5. (1992). Evaluating the accuracy of sampling-based approaches to calculating posterior moments.
    6. (1973). Fuller WA: Fitting segmented polynomial regression models whose joint points have to be estimated.
    7. (1995). G: Phenotypic plasticity as a product of selection in a variable environment. Am Nat
    8. (2008). Genetic component of heat stress in finishing pigs, development of a heat load function.
    9. (2008). Genetic component of heat stress in finishing pigs, Parameter estimation.
    10. (1998). Genetics and analysis of quantitative traits first edition.
    11. (2002). Genotype by environmental interaction in nordic dairy cattle studied using reaction norms. Acta Agric Scand
    12. (2008). K: CODA: Output analysis and diagnostics for MCMC. R package version 0.13–2
    13. (2000). Misztal I: Genetic component of heat stress in dairy cattle, parameter estimation.
    14. (1997). Multiple trait genetic analysis of underlying biological variables of production functions. Livest Prod Sci
    15. (2005). Random Regression analyses using B-splines to model growth of Australian Angus cattle. Genet Sel Evol
    16. (2004). Response to mass selection when the genotype by environment interaction is modelled as a linear reaction norm. Genet Sel Evol
    17. (2003). Selection in the presence of genotype by enviromental interaction: response in environmental sensitivity. Anim Sci
    18. (1980). The mathematical theory of quantitative genetics Oxford:

    To submit an update or takedown request for this paper, please submit an Update/Correction/Removal Request.