# runreml script                          Estimate variance components by REML, Paul VanRaden, 1986
#
#                                         Reformat phenotypes.txt
  sas /home/paul/reml/phenoreml
#                                         Make list of sires with most progeny
  sas /home/paul/reml/Ainv
#                                         Get inverse of pedigree relationship matrix
  /home/paul/reml/ainv >ainv.out
#                                         Absorb herd management goup effects
  /home/paul/reml/absorbm >absorbm.out
#                                         Iterate for heritabilities, correlations, PTAs
  /home/paul/reml/remlm >remlm.out
#                                         Compare sire PTAs to traditional or official 38 file
  sas /home/paul/reml/rdreml
