Skip Navigation Links
Journal of Applied Nonlinear Dynamics
Miguel A. F. Sanjuan (editor), Albert C.J. Luo (editor)
Miguel A. F. Sanjuan (editor)

Department of Physics, Universidad Rey Juan Carlos, 28933 Mostoles, Madrid, Spain


Albert C.J. Luo (editor)

Department of Mechanical and Industrial Engineering, Southern Illinois University Ed-wardsville, IL 62026-1805, USA

Fax: +1 618 650 2555 Email:

A Bayesian Random Walk Approach for Mapping Dynamic Quantitative Trait

Journal of Applied Nonlinear Dynamics 5(1) (2016) 105--115 | DOI:10.5890/JAND.2016.03.008

Burak Karacaören

Department of Animal Science, Faculty of Agriculture, Akdeniz University Antalya, 07059, TURKEY

Download Full Text PDF



Analyzing longitudinal observations with an appropriate methodology is important for detecting significant single nucleotide polymorphisms (SNPs) for complex traits in genome wide association studies (GWAS). Here we assumed that genomic signal over time could be traced by a Bayesian random walk- Kalman filter model in state space form to obtain longitudinal residuals. The Kalman filter removes the requirement for collecting the whole data set before making estimations, thus estimates are immediately available. For example, with the functional mapping approach, whole sets of observations must to be collected in advance. This may include months (if not years) of waiting time depending on the nature of the phenotypes and the experimental design. Whereas, because of the longitudinal residuals used in association mapping, biological reasoning can be easily applied given that the signal is genuine. The advantages of our method are demonstrated by both simulated and real datasets.


The author wishes to acknowledge useful discussions with Dr Li Zitong on mouse behavioral data set. The author wishes to acknowledge useful discussions with Dr Sitki ¸Ca˘gda¸s ˙Inam on LaTeX typesettings. This work was supported by Research Fund of the Akdeniz University Project Number 106.


  1. [1]  Ware, J.H. (1985), Linear models for the analysis of longitudinal studies, American Statistician, 39, 95-101.
  2. [2]  Macgregor, S. Knott, S.A.White, I. and Visscher, P.M. (2005), Quantitative trait locus analysis of longitudinal quantitative trait data in complex pedigrees, Genetics, 171, 1365-1376.
  3. [3]  El Faro, L. Cardoso, V.L. and Albuquerque, L.G. (2008), Variance component estimates applying random regression models for test-day milk yield in caracu heifers, Genetics and Molecular Biology, 31, 665-673.
  4. [4]  Wu, R.L. and Lin, M. (2006), Functional mapping-how to map and study the genetic architecture of dynamic complex traits, Nature Review Genetics, 7, 229-237.
  5. [5]  Kwak, Y.I. Moore, R.C. Spalding, P.E. and Broman, W.K. (2014), A Simple regression method to map quantitative trait loci underlying function-valued phenotypes, Genetics, 197, 1409-1416.
  6. [6]  Li, Z., & Sillanpää, M. J. (2013). A Bayesian nonparametric approach for mapping dynamic quantitative traits. Genetics, 194(4), 997-1016.
  7. [7]  Aulchenko, S.Y. de Koning, D.J. and Haley, C.S. (2007), Genomewide rapid association using mixed model and regression: a fast and simple method for genomewide pedigree-based quantitative trait loci association analysis, Genetics, 177, 577-585.
  8. [8]  Karacaören, B. (2014), Admixture mapping of growth related traits in F2 mice dataset using ancestry informative markers Journal of Bioinformatics and Computational Biology, 2(doi: 10.1142/S0219720014410108).
  9. [9]  Xiong, H., Goulding, E. H., Carlson, E. J., Tecott, L. H., McCulloch, C. E., & Sen, Ś. (2011). A flexible estimating equations approach for mapping function-valued traits. Genetics, 189(1), 305-316.
  10. [10]  Coster, A. Bastiaansen, J.W. Calus, M.P. Maliepaard, C. and Bink, M.C. (2010), QTLMAS 2009: Simulated dataset. BMC proceedings, 4, 1(doi:10.1186/1753-6561-4-S1-S3).
  11. [11]  R Development Core Team, (2007), R: A language and environmental for statistical computing Vienna: R Foundation for statistical computing.
  12. [12]  Karacaören, B. Janss, L.L.G, Kadarmideen, H.N. (2012), Predicting Breeding Values In Animals By Kalman Filter: Application To Body Condition Scores In Dairy Cattle, Journal of the Faculty of Veterinary Medicine, Kafkas University, 18, 627-632.
  13. [13]  West, M. and Harrison, J. (1997), Bayesian Forecasting and Dynamic Models Springer: New York
  14. [14]  Endelman, J.B. (2011), Ridge regression and other kernels for genomic selection with R package rrBLUP, Plant Genomics, 4, 250-255.
  15. [15]  Purcell, S. Neale, B. Todd-Brown, K. Thomas, L. Ferreira, M.A.R. Bender, D. Maller, J. Sklar, P. de Bakker, P.I.W. Daly, M.J. and Sham, P.C. (2007), PLINK: a toolset for whole-genome association and populationbased linkage analysis, American Journal of Human Genetics, 81, 559-575.
  16. [16]  Barrett, J.C. Fry, B. Maller, J. Daly, and M.J. (2005), Haploview: analysis and visualization of LD and haplotype maps, Bioinformatics, 21, 263-265.
  17. [17]  Refinetti, R. (2004), Non-stationary time series and the robustness of circadian rhythms, Journal of Theoretical Biology, 21, 571-581.
  18. [18]  Maliepaard, C. Bastiaansen, J.W. Calus, M.P, Coster, A. (2010), Comparison of analyses of the QTLMAS XIII common dataset II: QTL analysis, BMC proceedings, 4, 2(doi:10.1186/1753-6561-4-S1-S2).
  19. [19]  Meuwissen, T.H.E. Hayes, B.J. Goddard, M.E. (2001), Prediction of total genetic value using genome wide dense marker maps, Genetics,157,1819-1829.
  20. [20]  Weber, G.W. Alparslan, G¨ok.S.Z. S¨oyler, B. (2009a), A new mathematical approach in environmental and life sciences: gene-environment networks and their dynamics, Environmental Modeling and Assessment, 14, 267-288.