Differences

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

--- readme.pedictf90 [2020/03/13 19:10] – yutaka
+++ readme.pedictf90 [2022/05/09 22:39] – [Example] dani
@@ Line 8: / Line 8: @@
 It computes
+  * y_star = y adjusted by the other (not included) effects
   * y_hat = sum of estimates of the included effects
-  * y_star = y corrected by the other (not included) effects
   * residual = y - included effects (not a true residual)
@@ Line 20: / Line 20: @@
 with ''OPTION include_effects 3''.
+  * y_star = y - herd_hat - age_hat
   * y_hat = animal_hat
-  * ystar = y - herd_hat - age_hat
-Which makes cor(y_hat,y_star) = cor(ebv, corrected y) which is a measure of accuracy.
+Which makes cor(y_hat,y_star) = cor(ebv, adjusted y), in this example, which is a measure of accuracy.
-It outputs the correlation between y_hat and y_star, for instance cor(ystar,yhat)=cor(u+e, uhat) and outputs these columns into a file, together with animal id (if there is animal in the model) or record number (if not).
+It outputs the correlation between y_hat and y_star, for instance, cor(ystar,yhat)=cor(u+e, uhat) and outputs these columns into a file, together with animal id (if there is animal in the model) or record number (if not).
 In addition, if animal effect is in the model, it produces a file with ebvs from the solutions file.
@@ Line 33: / Line 33: @@
 This program splits the phenotype (y) into 3 pieces.
-  y = (effects to be corrected) + (effects to be kept) + residual
+  y = (effects to be adjusted for) + (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_star: ''y - (effects to be adjusted for)''
   * y_hat: ''(effects to kept)''
-  * y_star: ''y - (effects to be corrected)''
   * residual: ''residual''