readme.blupf90new
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revisionLast revisionBoth sides next revision | ||
readme.blupf90new [2022/05/05 21:34] – [Storing accuracy based on PEV] dani | readme.blupf90new [2022/05/06 01:50] – [Storing accuracy based on PEV] dani | ||
---|---|---|---|
Line 13: | Line 13: | ||
OPTION method VCE (default BLUP with blupf90 options) | OPTION method VCE (default BLUP with blupf90 options) | ||
</ | </ | ||
- | Run airemlf90 for variance component estimation (default running blupf90) | + | Runs airemlf90 for variance component estimation (default running blupf90) |
< | < | ||
OPTION conv_crit 1d-12 | OPTION conv_crit 1d-12 | ||
Line 25: | Line 25: | ||
OPTION EM-REML 10 | OPTION EM-REML 10 | ||
</ | </ | ||
- | Run EM-REML (REMLF90) for first 10 rounds to get initial variances within the parameter space (default 0). | + | Runs EM-REML (REMLF90) for first 10 rounds to get initial variances within the parameter space (default 0). |
< | < | ||
OPTION use_yams | OPTION use_yams | ||
</ | </ | ||
- | Run the program with YAMS (modified FSPAK). The computing time can be dramatically improved. | + | Runs the program with YAMS (modified FSPAK). The computing time can be dramatically improved. |
< | < | ||
OPTION tol 1d-12 | OPTION tol 1d-12 | ||
Line 38: | Line 38: | ||
OPTION sol se | OPTION sol se | ||
</ | </ | ||
- | Store solutions and those standard errors. | + | Stores |
< | < | ||
OPTION store_pev_pec 6 | OPTION store_pev_pec 6 | ||
</ | </ | ||
- | Store triangular matrices of standard errors and its covariances for correlated random effects such as direct-maternal effects and random-regression effects in " | + | Stores |
< | < | ||
OPTION residual | OPTION residual | ||
Line 50: | Line 50: | ||
OPTION missing -999 | OPTION missing -999 | ||
</ | </ | ||
- | Specify | + | Specifies |
< | < | ||
OPTION constant_var 5 1 2 ... | OPTION constant_var 5 1 2 ... | ||
Line 63: | Line 63: | ||
OPTION hetres_pos 10 11 | OPTION hetres_pos 10 11 | ||
</ | </ | ||
- | Specify | + | Specifies |
< | < | ||
OPTION hetres_pol 4.0 0.1 0.1 | OPTION hetres_pol 4.0 0.1 0.1 | ||
Line 133: | Line 133: | ||
< | < | ||
Indicate to store in file samples of (co)variances function for postprocessing (histogram, etc.) | Indicate to store in file samples of (co)variances function for postprocessing (histogram, etc.) | ||
+ | |||
+ | |||
+ | ==== Storing accuracy based on PEV ==== | ||
+ | |||
+ | OPTION store_accuracy eff | ||
+ | Stores reliabilities based on PEV, where //eff// is the number of the animal effect.\\ | ||
+ | By default, it uses inbreeding (F) in the denominator of the reliability formula: reliability = 1-PEV/ | ||
+ | It uses inbreeding based on $\mathbf{A}$ or $\mathbf{H}$ from the direct inversion of $\mathbf{A}^{-1}$ or $\mathbf{H}^{-1}$, | ||
+ | |||
+ | OPTION type 1.0 | ||
+ | Select 1.0 for dairy cattle (Reliability) or 0.5 for beef cattle (BIF accuracy) (default 1.0). | ||
+ | |||
+ | OPTION correct_accuracy_by_inbreeding filaname | ||
+ | // | ||
+ | |||
+ | OPTION correct_accuracy_by_inbreeding_direct 0 | ||
+ | This option turns off the inbreeding correction in the reliability formula. | ||
==== Omit A-inverse ==== | ==== Omit A-inverse ==== | ||
Line 192: | Line 209: | ||
The new option removes $\mathbf{A}^{-1}$ from the equations so GBLUP will be easily performed. | The new option removes $\mathbf{A}^{-1}$ from the equations so GBLUP will be easily performed. | ||
- | ==== Storing accuracy based on PEV ==== | ||
- | |||
- | OPTION store_accuracy eff | ||
- | Stores reliabilities based on PEV, where //eff// is the number of the animal effect.\\ | ||
- | By default, it uses inbreeding (F) in the denominator of the reliability formula: reliability = 1-PEV/ | ||
- | It uses inbreeding based on $\mathbf{A}$ or $\mathbf{H}$ from the direct inversion of $\mathbf{A}^{-1}$ or $\mathbf{H}^{-1}$, | ||
- | |||
- | OPTION type 1.0 | ||
- | Select 1.0 for dairy cattle (Reliability) or 0.5 for beef cattle (BIF accuracy) (default 1.0). | ||
- | |||
- | OPTION correct_accuracy_by_inbreeding filaname | ||
- | // | ||
- | |||
- | OPTION correct_accuracy_by_inbreeding_direct 0 | ||
- | This option turns off the the inbreeding correction in the reliability formula. | ||
===== Tricks ===== | ===== Tricks ===== |