User Tools

Site Tools


readme.pedictf90

Differences

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

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
Next revision Both sides next revision
readme.pedictf90 [2020/03/13 18:51]
yutaka
readme.pedictf90 [2020/03/13 19:10]
yutaka
Line 1: Line 1:
-**predictf90**+===== predictf90 ​===== 
 + 
 +==== Original description ====
  
 This is program is used to do cross-validations. It reads a blupf90 parameter file, a solutions file This is program is used to do cross-validations. It reads a blupf90 parameter file, a solutions file
Line 9: Line 11:
   * y_star = y corrected by the other (not included) effects   * y_star = y corrected by the other (not included) effects
   * residual = y - included effects (not a true residual)   * residual = y - included effects (not a true residual)
 +
 +==== Example ====
  
 For instance consider For instance consider
Line 24: Line 28:
  
 In addition, if animal effect is in the model, it produces a file with ebvs from the solutions file. In addition, if animal effect is in the model, it produces a file with ebvs from the solutions file.
 +
 +==== More explanation ====
 +
 +This program splits the phenotype (y) into 3 pieces.
 +
 +  y = (effects to be corrected) + (effects to be kept) + residual
 +
 +In a cross validation, you need a quantity "EBV + residual",​ where EBV is ''​(effects to be kept)''​ shown above, and it is ''​y_star''​.
 +The program calculates the following quantities.
 +
 +  * y_hat: ''​(effects to kept)''​
 +  * y_star: ''​y - (effects to be corrected)''​
 +  * residual: ''​residual''​
 +
 +The user should specify the effects to be kept using ''​OPTION include_effects''​ with the effect numbers.
 +For example, if the breeding value is in the 3rd row in the EFFECTS section in the parameter file, you have to use ''​OPTION include_effects 3''​.
  
 ==== Output files ==== ==== Output files ====