0

I am working on genomic selection and I am comparing the performance of two models, one of them is a likelihood method (GBlup) and the other is a Bayesian meyhod (Bayesian Ridge Regression). I am looking at the variance component and they are not exactly equal. So I wounder if it estimates should be equal. Anyone that have some idea?

marb_021
  • 53
  • 6
  • Why should completely different methods give the same answer? For BLUP see for example this [link](https://www.nature.com/articles/s41437-018-0075-0). For Bayesian ridge regression https://stats.stackexchange.com/a/328345/99274. – Carl Apr 06 '20 at 14:19
  • 1
    Hi, thank you for your comment. I got confused with the next statement: The β in BRR will be equal to the RR estimate with, λ= σε^2/σβ^2 , therefore it is also the BLUP of β given y. RR and BRR both perform an extent of shrinkage that is homogenous across markers (de los Campos et al. 2013) According to this, I interpreted It can be similar... – marb_021 Apr 06 '20 at 16:15
  • It would help if you linked to the work in which that statement is made. Ridge regression reduces variance by injecting bias so that if all that is being done by comparison is ML that will not be the case. – Carl Apr 06 '20 at 17:42

1 Answers1

1

GBLUP is mathematically equivalent to RR-BLUP when the genetic covariance matrix $G$ is calculated as $G=XX^T$from the data matrix $X$ (where each column of $X_{ij}$ is the value of the $j$th property (e.g. SNP) of the $i$th individual.

Also, the estimates of Bayesian Ridge regression are equivalent to the estimates of ridge regression. This review has it all in section 2.2.: Jacquin et al., 2016.

I would suspect from this that the estimates could be equal, but I'm not sure...

akraf
  • 646
  • 3
  • 15