This is an old revision of the document!
Table of Contents
BLUP90IOD3
Summary
The original program was written by Shogo Tsuruta in 1999.
BLUP90IOD3 is a BLUP program for solving a large number of mixed model equations, using preconditioned conjugate gradient (PCG). It accounts for heterogeneous residual variances and external information.
See PREGSF90 with genotypes (SNP) for options.
Implementation of one-step genomic selection updating animal effects with the combined pedigree and genomic relationship matrix (H-Inverse).
Options
OPTION conv_crit 1e-12
Set convergence criteria (deault 1e-12).
OPTION maxrounds 10000
Set maximum number of rounds (default 1000).
OPTION blksize 3
Set block size for preconditioner (default 1).
OPTION init_eq 10
Set the number of effects to be solved directly (default=0).
OPTION solv_method FSPAK
Solving method for initial equations (default=DIRECT).
OPTION tol 1d-12
Tolerance to get a positive definite matrix (default=1d-12).
OPTION residual
y-hat and residuals will be included in “yhat_residual”.
OPTION avgeps 50
Using the last 50 average eps for convergence.
OPTION cont 1
To restart the program from the previous solutions.
OPTION missing -999
Set the missing value (default 0) in integer.
OPTION restart 100
Set the number of iteration to recompute residuals (default=100).
OPTION prior_solutions
Using the previous solution file to start the iteration. An additional software is required to use this option.
OPTION random_upg
Set the UPG random.
OPTION SNP_file snp
Specify the SNP file name to use genotype data.
OPTION mode usage
Usage can be:
iod: data and pedigree stored in disk.
iom: data and pedigree stored in memory. It uses less operations than iod.
map: data in disk and pedigree in memory. It uses less operations than iom.
OPTION external_evaluation file
Read external breeding values and reliabilities from file.
OPTION hetres_int col nlev
where col is column in the data file that selects which residual (co)variance to select, and nlev is the maximum number of levels. Different residual (co)variances need to be numbered consecutively starting from 1.
Heterogeneous residual variances for a single trait
OPTION hetres_pos 10 11
Specifies the position of covariables.
OPTION hetres_pol 4.0 0.1 0.1
Initial values of coefficients for heterogeneous residual variances using ln(a0, a1, a2, …) to make these values.
Heterogeneous residual variances for multiple traits
Convergence will be very slow with multiple trait heterogeneous residual variances
OPTION hetres_pos 10 10 11 11
or
OPTION hetres_pos 10 11 12 13
Specify the position of covariables (trait first).
“10 10” or “10 11” could be linear for first and second traits.
“11 11” or “12 13” could be quadratic.
OPTION hetres_pol 4.0 4.0 0.1 0.1 0.01 0.01
Initial values of coefficients for heterogeneous residual variances using ln(a0, a1, a2, …) to make these values (trait first).
“4.0 4.0” are intercept for first and second traits.
“0.1 0.1” could be linear and “0.01 0.01” could be quadratic.
To transform back to the original scale, use exp(a0+a1*X1+a2*X2).
For information on availability of BLUP90IOD3, please contact Ignacy Misztal (ignacy@uga.edu) at the University of Georgia.