readme.pedictf90
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.pedictf90 [2020/03/13 19:07] – yutaka | readme.pedictf90 [2022/05/09 22:38] – [Example] dani | ||
---|---|---|---|
Line 12: | Line 12: | ||
* residual = y - included effects (not a true residual) | * residual = y - included effects (not a true residual) | ||
- | ==== More explanation | + | ==== Example |
- | This program splits the phenotype (y) into 3 pieces. | + | For instance consider |
- | y = (effects to be corrected) | + | y = herd + age + animal + e |
- | In a cross validation, you need a quantity " | + | with '' |
- | Therefore, you should need '' | + | |
- | The program calculates the following quantities. | + | |
- | * y_hat: '' | + | * ystar = y - herd_hat - age_hat |
- | * y_star: '' | + | * y_hat = animal_hat |
- | * residual: '' | + | |
+ | Which makes cor(y_hat, | ||
- | === Example === | + | It outputs the correlation between y_hat and y_star, for instance, cor(ystar, |
- | For instance consider | + | In addition, if animal effect is in the model, it produces a file with ebvs from the solutions file. |
- | y = herd + age + animal + e | + | ==== More explanation ==== |
- | with '' | + | This program splits the phenotype (y) into 3 pieces. |
- | | + | |
- | * ystar = y - herd_hat - age_hat | + | |
- | Which makes cor(y_hat,y_star) = cor(ebv, corrected y) which is a measure of accuracy. | + | In a cross validation, you need a quantity "EBV + residual" |
+ | The program calculates the following quantities. | ||
- | It outputs the correlation between | + | * y_hat: '' |
+ | * y_star: '' | ||
+ | * residual: '' | ||
- | In addition, if animal effect | + | The user should specify the effects to be kept using '' |
+ | For example, if the breeding value is in the 3rd row in the EFFECTS section in the parameter | ||
==== Output files ==== | ==== Output files ==== |