Quantitative Genetic Models for Describing Simultaneous and Recursive Relationships Between Phenotypes This article is dedicated to Arthur B. Chapman, teacher and mentor of numerous animal breeding students and disciple and friend of Sewall Wright.

Multivariate models are of great importance in theoretical and applied quantitative genetics. We extend quantitative genetic theory to accommodate situations in which there is linear feedback or recursiveness between the phenotypes involved in a multivariate system, assuming an infinitesimal, additive, model of inheritance. It is shown that structural parameters defining a simultaneous or recursive system have a bearing on the interpretation of quantitative genetic parameter estimates (e.g., heritability, offspring-parent regression, genetic correlation) when such features are ignored. Matrix representations are given for treating a plethora of feedback-recursive situations. The likelihood function is derived, assuming multivariate normality, and results from econometric theory for parameter identification are adapted to a quantitative genetic setting. A Bayesian treatment with a Markov chain Monte Carlo implementation is suggested for inference and developed. When the system is fully recursive, all conditional posterior distributions are in closed form, so Gibbs sampling is straightforward. If there is feedback, a Metropolis step may be embedded for sampling the structural parameters, since their conditional distributions are unknown. Extensions of the model to discrete random variables and to nonlinear relationships between phenotypes are discussed.

[1]  Malcolm E. Turner,et al.  The Regression Analysis of Causal Paths , 1959 .

[2]  Bradley P. Carlin,et al.  Markov Chain Monte Carlo conver-gence diagnostics: a comparative review , 1996 .

[3]  G. C. Tiao,et al.  Bayesian inference in statistical analysis , 1973 .

[4]  Jacques H. Dreze,et al.  BAYESIAN ANALYSIS OF SIMULTANEOUS EQUATION SYSTEMS , 1983 .

[5]  Arnold Zellner,et al.  An Introduction to Bayesian Inference in Econometrics. , 1974 .

[6]  Daniel Gianola,et al.  Theory and Analysis of Threshold Characters , 1982 .

[7]  Karl G. Jöreskog,et al.  Lisrel 8: User's Reference Guide , 1997 .

[8]  R. Waagepetersen,et al.  Normal linear models with genetically structured residual variance heterogeneity: a case study. , 2003, Genetical research.

[9]  Cedric A. B. Smith,et al.  An Introduction to Genetic Statistics , 1958 .

[10]  R. Fisher.  XV.—The Correlation between Relatives on the Supposition of Mendelian Inheritance. , 1919, Transactions of the Royal Society of Edinburgh.

[11]  B. Shipley.  Cause and correlation in biology , 2000 .

[12]  Daniel Gianola,et al.  Bayesian Methods in Animal Breeding Theory , 1986 .

[13]  Bradley P. Carlin,et al.  BAYES AND EMPIRICAL BAYES METHODS FOR DATA ANALYSIS , 1996, Stat. Comput..

[14]  H. F. Smith,et al.  A DISCRIMINANT FUNCTION FOR PLANT SELECTION , 1936 .

[15]  P. Games.  Correlation and Causation: A Logical Snafu , 1990 .

[16]  J S Haldane,et al.  The regulation of the lung‐ventilation , 1905, The Journal of physiology.

[17]  David B. Dunson,et al.  Bayesian Data Analysis , 2010 .

[18]  Arnold Zellner,et al.  Statistical Analysis of Econometric Models , 1979 .

[19]  G. Judge,et al.  The Theory and Practice of Econometrics , 1981 .

[20]  S. R. Searle.  Linear Models , 1971 .

[21]  Andrew L. Rukhin,et al.  Tools for statistical inference , 1991 .

[22]  A. Goldberger.  STRUCTURAL EQUATION METHODS IN THE SOCIAL SCIENCES , 1972 .

[23]  M. Kendall,et al.  Kendall's advanced theory of statistics , 1995 .

[24]  Sylvia Richardson,et al.  Markov Chain Monte Carlo in Practice , 1997 .

[25]  D. A. Kenny,et al.  Correlation and Causation , 1937, Wilmott.

[26]  R. L. Quaas,et al.  Computing the Diagonal Elements and Inverse of a Large Numerator Relationship Matrix , 1976 .

[27]  Christian P. Robert,et al.  Monte Carlo Statistical Methods , 2005, Springer Texts in Statistics.

[28]  C. R. Henderson.  SIRE EVALUATION AND GENETIC TRENDS , 1973 .

[29]  L. N. Hazel.  The Genetic Basis for Constructing Selection Indexes. , 1943, Genetics.

[30]  D. Falconer,et al.  Introduction to Quantitative Genetics. , 1962 .

[31]  M. Bulmer.  The Mathematical Theory of Quantitative Genetics , 1981 .

[32]  J. Cheverud,et al.  Quantitative genetics and developmental constraints on evolution by selection. , 1984, Journal of theoretical biology.

[33]  D. Gianola,et al.  Marginal inferences about variance components in a mixed linear model using Gibbs sampling , 1993, Genetics Selection Evolution.

[34]  A. Goldberger,et al.  Structural Equation Models in the Social Sciences. , 1974 .

[35]  E. Ziegel.  Introduction to the Theory and Practice of Econometrics , 1989 .

[36]  Anthony O'Hagan,et al.  Kendall's Advanced Theory of Statistics: Vol. 2B, Bayesian Inference. , 1996 .

[37]  J. Hanrahan.  Maternal effects and selection response. , 1974 .

[38]  T. Haavelmo.  The Statistical Implications of a System of Simultaneous Equations , 1943 .

[39]  H. Kacser,et al.  The molecular basis of dominance. , 1981, Genetics.

[40]  O. D. Duncan,et al.  Introduction to Structural Equation Models. , 1977 .

[41]  C. R. Henderson.  A simple method for computing the inverse of a numerator relationship matrix used in prediction of breeding values , 1976 .

[42]  Sewall Wright,et al.  The Treatment of Reciprocal Interaction, with or without Lag, in Path Analysis , 1960 .

[43]  S. Wright.  Corn and hog correlations , 1925 .