readme.blupadj
no way to compare when less than two revisions
Differences
This shows you the differences between two versions of the page.
| — | readme.blupadj [2024/03/25 18:22] (current) – created - external edit 127.0.0.1 | ||
|---|---|---|---|
| Line 1: | Line 1: | ||
| + | ====== BLUPF90ADJ ====== | ||
| + | |||
| + | ===== Summary ===== | ||
| + | This program preadjust the data for all effects except the first effect. Its expected use is in preadjustment of data for fixed effects such as age of dam, sex, standard day day of weighting etc. | ||
| + | |||
| + | The preadjustment can be for any number of traits and effects. Observations with a value of 0 are not adjusted. | ||
| + | |||
| + | To use the program prepare the parameter file as for BLUPF90. The first effect, which is not used for adjustment, can be contemporary group, HYS or the mean. | ||
| + | |||
| + | The program writes the data back in file adjusted_data. The format of that file can be provided by the user or a program default can be used. | ||
| + | The current default is (100f10.1). | ||
| + | |||
| + | If a custom form is used, it needs to be in the fortran format for real numbers, e.g., (5f8.0, f10.3) | ||
| + | \\ | ||
| + | ===== Problems ===== | ||
| + | The adjusted observation can be 0 and and can be treated by blupf90 and related programs as missing. | ||
| + | \\ | ||
| + | ===== Examples ===== | ||
| + | The example is based on a single trait example with parameter file exhxan2 in teh examples directory. We are adjusting for the effect in the third column treating effect in column 2 as the first effect.. | ||
| + | |||
| + | $ head hxrec2 | ||
| + | < | ||
| + | 5 22 261 457 60 49 817 1213164 | ||
| + | 5 23 267 743 60 59 790 1311935 | ||
| + | 5 23 259 812 60 57 811 1143578 | ||
| + | 5 22 264 719 60 42 752 1253399 | ||
| + | 5 23 264 748 60 59 760 1253399 | ||
| + | 5 22 264 780 60 43 741 1253399 | ||
| + | 6 30 316 358 60 57 851 1304048 | ||
| + | 5 22 257 735 61 37 744 1085978 | ||
| + | 6 30 310 40 61 50 846 1283313 | ||
| + | 6 31 299 84 61 56 812 1190632 | ||
| + | </ | ||
| + | $ cat exhxan2adj | ||
| + | < | ||
| + | # Example of single-trait animal model with one fixed effect | ||
| + | DATAFILE | ||
| + | hxrec2 | ||
| + | NUMBER_OF_TRAITS | ||
| + | 1 | ||
| + | NUMBER_OF_EFFECTS | ||
| + | 2 | ||
| + | OBSERVATION(S) | ||
| + | 7 | ||
| + | WEIGHT(S) | ||
| + | |||
| + | EFFECTS: POSITIONS_IN_DATAFILE NUMBER_OF_LEVELS TYPE_OF_EFFECT | ||
| + | 3 656 cross | ||
| + | 2 70 cross | ||
| + | RANDOM_RESIDUAL VALUES | ||
| + | 3 | ||
| + | 3 | ||
| + | </ | ||
| + | $blupf90adj | ||
| + | < | ||
| + | ADJUSTS DATA FOR ALL EFFECTS EXCEPT THE FIRST ONE | ||
| + | name of parameter file? | ||
| + | | ||
| + | |||
| + | | ||
| + | Data file: hxrec2 | ||
| + | | ||
| + | | ||
| + | | ||
| + | | ||
| + | Value of Missing Trait/ | ||
| + | |||
| + | EFFECTS | ||
| + | # | ||
| + | | ||
| + | | ||
| + | |||
| + | | ||
| + | 3.0000 | ||
| + | |||
| + | | ||
| + | (1) Weight position 0 means no weights utilized | ||
| + | (2) Effect positions of 0 for some effects and traits means that such | ||
| + | effects are missing for specified traits | ||
| + | |||
| + | Data record length = 7 | ||
| + | | ||
| + | | ||
| + | | ||
| + | | ||
| + | | ||
| + | | ||
| + | (3f6.0, | ||
| + | wrote file with adjusted traits called adjusted_data | ||
| + | </ | ||
| + | < | ||
| + | $ tail -50 solutions|head | ||
| + | 1 2 21 | ||
| + | 1 2 22 | ||
| + | 1 2 23 | ||
| + | 1 2 24 | ||
| + | 1 2 25 | ||
| + | 1 2 26 | ||
| + | 1 2 27 | ||
| + | 1 2 28 | ||
| + | 1 2 29 | ||
| + | 1 2 30 | ||
| + | </ | ||
| + | $ head adjusted_data | ||
| + | < | ||
| + | 5. | ||
| + | 5. | ||
| + | 5. | ||
| + | 5. | ||
| + | 5. | ||
| + | 5. | ||
| + | 6. | ||
| + | 5. | ||
| + | 6. | ||
| + | 6. | ||
| + | </ | ||
| + | |||
| + | Note that 1369=817-(-552) | ||
readme.blupadj.txt · Last modified: 2024/03/25 18:22 by 127.0.0.1