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:38] – [Example] dani
@@ Line 11: / Line 11: @@
   * y_star = y corrected by the other (not included) effects
   * residual = y - included effects (not a true residual)
+==== Example ====
+For instance consider
+  y = herd + age + animal + e
+with ''OPTION include_effects 3''.
+  * ystar = y - herd_hat - age_hat
+  * y_hat = animal_hat
+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).
+In addition, if animal effect is in the model, it produces a file with ebvs from the solutions file.
 ==== More explanation ====
@@ Line 27: / Line 44: @@
 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''.
-=== Example ===
-For instance consider
-  y = herd + age + animal + e
-with ''OPTION include_effects 3''.
-  * 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.
-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.
 ==== Output files ====