BLUPF90

User Tools

Site Tools

Trace:
readme.renumf90

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision		Previous revision
Next revisionBoth sides next revision
readme.renumf90 [2012/05/29 12:28] – [Options] ignacioreadme.renumf90 [2016/10/13 22:00] – [Fields in the parameter file] ignacio
Line 1: Line 1:
  
 =======RENUMF90====== =======RENUMF90======
-renumbering program for the BLUPF90 family now works with SNP info\\+A renumbering program for the ''BLUPF90'' family now works with SNP info\\
 Ignacy Misztal and Ignacio Aguilar, University of Georgia\\ Ignacy Misztal and Ignacio Aguilar, University of Georgia\\
 August 27, 2001 - Mar 17, 2011 August 27, 2001 - Mar 17, 2011
Line 42: Line 42:
 Lines with # are treated as comments Lines with # are treated as comments
              
-<file> +<file>
 DATAFILE DATAFILE
 f1 f1
Line 48: Line 48:
 The data file is f1 The data file is f1
  
-<file> +<file>
 TRAITS TRAITS
 t1 t2 .. tn t1 t2 .. tn
Line 76: Line 76:
 e1.. en type form e1.. en type form
 </file> </file>
-this line defines one group of effects; e1 .. en are positions of this effect for all traits;\\  +this line defines one group of effects; e1 .. en are positions of this effect for all traits; 
-positions can be different for each trait for fixed effects;\\ + 
 +positions can be different for each trait for fixed effects; 
 for random effects, only one position + 0 (misising) efefct are for random effects, only one position + 0 (misising) efefct are
 possible.  possible. 
Line 87: Line 89:
 d1 .. dn form d1 .. dn form
 </file> </file>
-optional for covariables only, specifies nesting;\\+optional for covariables only, specifies nesting; 
 form is as above form is as above
      
Line 96: Line 99:
  
 the RANDOM keyword occurs only if the current effect is random; the RANDOM keyword occurs only if the current effect is random;
-rtype is 'diagonal', 'sire' or 'animal'+rtype is: 
 +  * 'diagonal' 
 +  * 'sire'; not yet implemented 
 +  * 'animal'
  
 <file> <file>
Line 102: Line 108:
 o1 o2.. oq o1 o2.. oq
 </file> </file>
-causes extra effects appended to the animal effect;\\ +causes extra effects appended to the animal effect; 
-current options include+ 
 +current options include:
   * 'pe' for permanent environment   * 'pe' for permanent environment
   * 'mat' for maternal   * 'mat' for maternal
-  * 'mpe' for maternal permanent environment +  * 'mpe' for maternal permanent environment; only if 'mat' is used
  
 <file> <file>
Line 113: Line 120:
 </file> </file>
  
-for animal and sire model only, fped specifies the pedigree file+for animal and sire model only 
 + 
 + //fped// specifies the pedigree file
  
 <file> <file>
Line 120: Line 129:
 </file> </file>
  
-for animal effect only;\\ +for animal effect only; 
-specifies positions in the pedigree file of animal an, sire s, dam d, alternate_dam dam rec_dam , and year of birth yob\\ + 
-missing alt_dam or yob can be replaced by 0\\ +specifies positions in the pedigree file of animal (an), sire (s), dam (d), alternate dam (alt_dam) , and year of birth (yob) 
-if this line is not given, defaults are 1 2 3 0 0.\\ + 
-If maternal effect is specified, the maternal effect is due to position of d if alt_dam field is 0, or otherwise is due to alt_dam;\\ If alt_dam field is not zero, it should include ID of real or recipient dam.+missing alt_dam or yob can be replaced by 0 
 + 
 +if this line is not given, defaults are 1 2 3 0 0. 
 + 
 +If maternal effect is specified, the maternal effect is due to position of d if alt_dam field is 0, or otherwise is due to alt_dam; 
 + 
 +If alt_dam field is not zero, it should include ID of real or recipient dam.
        
 <file> <file>
Line 132: Line 147:
  
 optional;\\ optional;\\
-//fsnp// specifies files with ID and SNP information;\\+//fsnp// specifies files with ID and SNP information; 
 if present, the relationship matrix will be constructed as in Aguilar et al. (2010) and will if present, the relationship matrix will be constructed as in Aguilar et al. (2010) and will
-include the genomic information;\\ +include the genomic information; 
-file fsnp should start with ID with the same format as fped and SNP info needs to start from a fixed column and + 
-include digits 0, 1, 2 and 5;\\+file //fsnp// should start with ID with the same format as fped and SNP info needs to start from a fixed column and 
 +include digits 0, 1, 2 and 5; 
 + 
 ID and SNP info need to be separated by at leats one space; see info for program [[readme.pregsf90|PreGSf90]]  ID and SNP info need to be separated by at leats one space; see info for program [[readme.pregsf90|PreGSf90]] 
  
Line 144: Line 163:
 </file> </file>
  
-optional\\ +optional 
-for animal effect only;\\ + 
-p specifies the depth of pedigree search;\\ +for animal effect only; 
-the default is 3; all pedigrees are loaded if p=0.+ 
 +//p// specifies the depth of pedigree search; 
 + 
 +the default is 3 
 + 
 +all pedigrees are loaded if p=0.
  
 <file> <file>
Line 153: Line 177:
 min avg max min avg max
 </file> </file>
-optional\\ +optional 
-specifies minimum, average and maximum generation interval;\\+ 
 +specifies minimum, average and maximum generation interval 
 applicable only if year of birth present; minimum and maximum used for applicable only if year of birth present; minimum and maximum used for
-pedigree checks;\\+pedigree checks 
 average used to predict year of birth of parent with missing average used to predict year of birth of parent with missing
 pedigree. pedigree.
Line 165: Line 192:
 </file> </file>
 optional\\ optional\\
-if only one sex has records, specifies which parent it is;\\+if only one sex has records, specifies which parent it is 
 used for pedigree checks. used for pedigree checks.
  
Line 177: Line 205:
   * if t is 'in_pedigrees', the value of a missing parent should be -x, where x 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 'in_pedigrees', the value of a missing parent should be -x, where x 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 'internal',  allocation is by a user-written function custom_upg(year_of_birth,sex,ID, parent_code).+  * if t is 'internal',  allocation is by a user-written function custom_upg(year_of_birth,sex,ID, parent_code); not yet implemented. 
 +<file> 
 +INBREEDING 
 +inb_type 
 +</file> 
 +optional\\ 
 +use of inbreeding coefficients to compute inb/upg code in the 4th column of the output pedigree file 
 +inb_type could be: 
 +  * 'pedigree', the program computes inbreeding coefficients with Meuwissen and Luo (1992) using the pedigree to be saved in renaddxx.ped; calculated inbreeding coefficients will be saved in a file "renf90.inb" with the original ID 
 +  * 'file',the program reads inbreeding coefficients from an external file. You should put the filename after 'file' e.g. 'file inbreeding.txt'. The file has at least 2 columns: original_ID and inbreeding value (from 0.0 to 1.0). The program just skips unnecessary IDs
  
 <file> <file>
Line 184: Line 221:
 </file> </file>
 Specifies that random regressions should be applied to the animal and Specifies that random regressions should be applied to the animal and
-corresponding effects (mat, pe and mpe),\\ +corresponding random effects (mat, pe and mpe) or the diagonal random effect. 
-this keyword also could be applied to set covariables for fixed effects;\\ + 
-r_type is 'data' if covariables for random regressions are in the data, or "legendre' if legendre plynomials are +this keyword also could be applied to set covariables for fixed effects; 
-to be generated from a single data variable; not yet implemented+ 
 +//r_type// could be: 
 +  * 'data' if covariables for random regressions are in the data 
 +  * "legendre' if legendre plynomials are to be generated from a single data variable; not yet implemented
  
 <file> <file>
Line 193: Line 233:
 r1 .. rq r1 .. rq
 </file> </file>
-for random regressions, r1-rq specifies positions of covariables if +for random regressions, 
-r_type='data', or r1 is order of legendre polynomial and r2 is position of + 
-covariable if r_type='legendre'; not yet implemented +  * //r1-rq// specifies positions of covariables if r_type='data' 
-</file>+  * //r1// is order of legendre polynomial and //r2// is position of covariable if r_type='legendre'; not yet implemented
  
 <file> <file>
Line 203: Line 243:
 </file> </file>
  
-//g// are (co)variances for the animal effect;\\ +//g// are (co)variances for the animal effect 
-the dimensions of g should account for the maternal effect if present+ 
 +the dimensions of g should account for random correlated effect if present (maternal or random regression)
  
 <file> <file>
Line 217: Line 258:
 </file>  </file> 
 //gmpe// are (co)variances for the MPE effect if present //gmpe// are (co)variances for the MPE effect if present
- 
- 
 =====Extra comments===== =====Extra comments=====
  
-Sections starting from EFFECTS can be repeated any number of types.\\+Sections starting from EFFECTS can be repeated any number of types. 
 If (Co)variances for any effect are missing, they are substituted with matrices containing 1.0  If (Co)variances for any effect are missing, they are substituted with matrices containing 1.0 
 on diagonals and 0.1 on off-diagonals.  on diagonals and 0.1 on off-diagonals. 
Line 240: Line 280:
  - the default size of character fields  - the default size of character fields
 <file> <file>
- OPTION alpha size nn+OPTION alpha_size nn
 </file> </file>
 where //nn// is the new size. where //nn// is the new size.
Line 246: Line 286:
  - the size of th record length   - the size of th record length 
 <file> <file>
- OPTION max_string_readline nn+OPTION max_string_readline nn
 </file> </file>
 where //nn// is the new size. where //nn// is the new size.
Line 252: Line 292:
  - the maximun number of fields    - the maximun number of fields  
 <file> <file>
-        OPTION max_field_readline nn+OPTION max_field_readline nn
 </file> </file>
 where //nn// is the number of fields.  where //nn// is the number of fields. 
Line 268: Line 308:
 </file> </file>
 catenates //b c// ... into //a//.\\ catenates //b c// ... into //a//.\\
-Keywords COMBINE need to be on top of the parameter file, but possibly after comments.+Keyword COMBINE needs to be on top of the parameter file, but possibly after comments.
  
 There may be many combined fields.  There may be many combined fields. 
Line 295: Line 335:
    4) 3 minus number of known parents                                   4) 3 minus number of known parents                               
    5) known or estimated year of birth (0 if not provided)              5) known or estimated year of birth (0 if not provided)          
-   6) number of known parents (parents might be eliminated if not    +   6) number of known parents (parents might be eliminated if not contributing; 
-      contributing; if animal has genotype 10+number of know parents                                                  +      if animal has genotype 10+number of know parents                                                  
    7) number of records                 7) number of records             
-   8) number of progenies (before elimination due to other effects)  +   8) number of progenies (before elimination due to other effects) as parent 1 
-      as parent 1 +   9) number of progenies (before elimination due to other effects) as parent 2  
-   9) number of progenies (before elimination due to other effects)  +
-      as parent 2  +
   10) original animal id                                               10) original animal id                                            
  
Line 340: Line 378:
 <file> <file>
 # Parameter file for program renf90; it is translated to parameter # Parameter file for program renf90; it is translated to parameter
-# file for BLUPF90 family f programs.+# file for BLUPF90 family programs.
 DATAFILE DATAFILE
 data.test data.test
readme.renumf90.txt · Last modified: 2024/05/22 18:53 by dani
Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki