readme.renumf90
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revisionNext revisionBoth sides next revision | ||
readme.renumf90 [2012/05/29 11:39] – ignacio | readme.renumf90 [2012/05/29 11:50] – ignacio | ||
---|---|---|---|
Line 57: | Line 57: | ||
FIELDS_PASSED TO OUTPUT | FIELDS_PASSED TO OUTPUT | ||
p1 p2 .. pm | p1 p2 .. pm | ||
- | < | + | </file> |
fields p1-pn are passed to output without changes; can be empty | fields p1-pn are passed to output without changes; can be empty | ||
Line 100: | Line 100: | ||
OPTIONAL | OPTIONAL | ||
o1 o2.. oq | o1 o2.. oq | ||
- | < | + | </file> |
causes extra effects appended to the animal effect; current options include | causes extra effects appended to the animal effect; current options include | ||
' | ' | ||
Line 141: | Line 141: | ||
</ | </ | ||
- | < | ||
optional for animal effect only; p specifies the depth of pedigree search; | optional for animal effect only; p specifies the depth of pedigree search; | ||
the default is 3; all pedigrees are loaded if p=0. | the default is 3; all pedigrees are loaded if p=0. | ||
- | + | ||
+ | < | ||
GEN_INT | GEN_INT | ||
min avg max | min avg max | ||
+ | </ | ||
+ | optional; specifies minimum, average and maximum generation interval; | ||
+ | applicable only if year of birth present; minimum and maximum used for | ||
+ | pedigree checks; average used to predict year of birth of parent with missing | ||
+ | pedigree. | ||
- | optional; specifies minimum, average and maximum generation interval; | + | < |
- | | + | |
- | | + | |
- | | + | |
- | + | ||
REC_SEX | REC_SEX | ||
i | i | ||
+ | </ | ||
+ | optional; if only one sex has records, specifies which parent it is; used for | ||
+ | pedigree checks. | ||
- | optional; if only one sex has records, specifies which parent it is; used for | + | < |
- | | + | |
- | + | ||
UPG_TYPE | UPG_TYPE | ||
t | t | ||
- | + | </ | |
- | | + | optional; |
- | | + | if t is ' |
- | | + | subsequent line should contain list of years to separate different UPG; |
- | if t is ' | + | if t is ' |
- | | + | UPG number that this missing parent should be allocated to; in this option, |
- | | + | all known parents should have pedigree lines, i.e., each parent field should |
- | | + | contain either the ID of a real parent, or a negative UPG number. |
- | if t is ' | + | if t is ' |
- | | + | custom_upg(year_of_birth, |
- | + | < | |
RANDOM_REGRESSION | RANDOM_REGRESSION | ||
r_type | r_type | ||
+ | </ | ||
+ | Specifies that random regressions should be applied to the animal and | ||
+ | corresponding effects (mat, pe and mpe), this keyword also could be | ||
+ | applied to set covariables for fixed effects; r_type is ' | ||
+ | random regressions are in the data, or " | ||
+ | to be generated from a single data variable; not yet implemented | ||
- | Specifies that random regressions should be applied to the animal and | + | < |
- | | + | |
- | | + | |
- | | + | |
- | to be generated from a single data variable; not yet implemented | + | |
- | + | ||
RR_POSITION | RR_POSITION | ||
r1 .. rq | r1 .. rq | ||
- | for random regressions, | + | </ |
- | | + | for random regressions, |
- | | + | r_type=' |
- | | + | covariable if r_type=' |
+ | </ | ||
+ | |||
+ | < | ||
(CO)VARIANCES | (CO)VARIANCES | ||
g | g | ||
+ | </ | ||
- | g are (co)variances for the animal effect; the dimensions of g should | + | //g// are (co)variances for the animal effect; the dimensions of g should |
- | | + | account for the maternal effect if present |
- | | + | |
+ | < | ||
(CO)VARIANCES_PE | (CO)VARIANCES_PE | ||
gpe | gpe | ||
+ | </ | ||
+ | //gpe// are (co)variances for the PE effect if present | ||
- | gpe are (co)variances for the PE effect if present | + | < |
(CO)VARIANCES_MPE | (CO)VARIANCES_MPE | ||
gmpe | gmpe | ||
+ | </ | ||
+ | //gmpe// are (co)variances for the MPE effect if present | ||
- | gmpe are (co)variances for the MPE effect if present | ||
+ | ---- | ||
Sections starting from EFFECTS can be repeated any number of types. | Sections starting from EFFECTS can be repeated any number of types. | ||
Line 218: | Line 229: | ||
skipped. Keywords out of order may not be recognized. | skipped. Keywords out of order may not be recognized. | ||
+ | =====OPTIONS===== | ||
The following options can added at the end of the parameter file to redefine | The following options can added at the end of the parameter file to redefine | ||
parameter used to read the input file: | parameter used to read the input file: | ||
- the default size of character fields | - the default size of character fields | ||
+ | < | ||
OPTION alpha size nn | OPTION alpha size nn | ||
+ | </ | ||
where nn is the new size. | where nn is the new size. | ||
- the size of th record length | - the size of th record length | ||
+ | < | ||
OPTION max_string_readline nn | OPTION max_string_readline nn | ||
+ | </ | ||
where nn is the new size. | where nn is the new size. | ||
- the maximun number of fields | - the maximun number of fields | ||
+ | < | ||
OPTION max_field_readline nn | OPTION max_field_readline nn | ||
+ | </ | ||
where nn is the number of fields. | where nn is the number of fields. | ||
Line 239: | Line 257: | ||
- | Combining fields or interactions | + | =====Combining fields or interactions====== |
- | ================================ | + | |
Several fields in the data file can be combined into one using a COMBINE keyword. | Several fields in the data file can be combined into one using a COMBINE keyword. | ||
- | COMBINE a b c .... | + | < |
+ | COMBINE a b c .... | ||
+ | </ | ||
catenates b c ... into c. Keywords COMBINE need to be on top of the parameter | catenates b c ... into c. Keywords COMBINE need to be on top of the parameter | ||
file, but possibly after comments. There may be many combined fields. | file, but possibly after comments. There may be many combined fields. | ||
For example: | For example: | ||
- | COMBINE 7 2 3 4 | + | </ |
+ | COMBINE 7 2 3 4 | ||
+ | </ | ||
combines content of fields 2 3 4 into field 7; the data file is not changed, | combines content of fields 2 3 4 into field 7; the data file is not changed, | ||
only the program treats field 7 as fields 2 3 4 put together (without spaces). | only the program treats field 7 as fields 2 3 4 put together (without spaces). | ||
Line 256: | Line 277: | ||
- | Additive Pedigree File | + | =====Additive Pedigree File====== |
- | ====================== | + | |
The additive pedigree file(s) renadd* has the following structure: | The additive pedigree file(s) renadd* has the following structure: | ||
Line 276: | Line 296: | ||
- | Extensions | + | =====Extensions======= |
- | ========== | + | |
The program is being modified to support inbreeding, dominance, random | The program is being modified to support inbreeding, dominance, random | ||
Line 283: | Line 302: | ||
- | Example | + | =====Example====== |
- | ======== | + | |
+ | < | ||
data file - data.test | data file - data.test | ||
--------------------- | --------------------- | ||
Line 294: | Line 313: | ||
3 aa 30 55 56 yy | 3 aa 30 55 56 yy | ||
5 bb 1234567.890 66 67 zz | 5 bb 1234567.890 66 67 zz | ||
+ | </ | ||
+ | < | ||
pedigree file - test.ped | pedigree file - test.ped | ||
------------------------ | ------------------------ | ||
Line 302: | Line 323: | ||
cc qq 0 | cc qq 0 | ||
dd 0 aa | dd 0 aa | ||
+ | </ | ||
+ | </ | ||
parameter file - testpar1 | parameter file - testpar1 | ||
------------------------- | ------------------------- | ||
Line 347: | Line 370: | ||
RANDOM | RANDOM | ||
diagonal | diagonal | ||
+ | </ | ||
+ | </ | ||
printout (temporary; the amount of details may change) | printout (temporary; the amount of details may change) | ||
------------------------------------------------------ | ------------------------------------------------------ | ||
Line 439: | Line 463: | ||
Wrote parameter file " | Wrote parameter file " | ||
Wrote renumbered data " | Wrote renumbered data " | ||
+ | </ | ||
+ | < | ||
new parameter file - renf90.par | new parameter file - renf90.par | ||
------------------------------- | ------------------------------- | ||
Line 503: | Line 528: | ||
| | ||
0.1000 | 0.1000 | ||
+ | </ | ||
+ | < | ||
data file - renf90.dat | data file - renf90.dat | ||
---------------------- | ---------------------- | ||
Line 514: | Line 540: | ||
30 55 2 3 5 56 2 2 aa 3 | 30 55 2 3 5 56 2 2 aa 3 | ||
| | ||
+ | </ | ||
+ | < | ||
| | ||
| | ||
Line 525: | Line 552: | ||
3 0 5 1 0 1 2 0 2 aa | 3 0 5 1 0 1 2 0 2 aa | ||
4 6 7 1 0 2 1 0 0 cc | 4 6 7 1 0 2 1 0 0 cc | ||
+ | </ | ||
+ | < | ||
renumbering tables - renf90.tables | renumbering tables - renf90.tables | ||
---------------------------------- | ---------------------------------- |
readme.renumf90.txt · Last modified: 2024/05/22 18:53 by dani