speed up internal matrix functions

Question

speed up internal matrix functions

joelpick opened this issue 3 years ago · 3 comments

Answer 1 · 2022-06-22T05:02:57.000Z

maybe use mvnfast::rmvn its very fast :)

Answer 2 · 2022-08-10T19:19:09.000Z

for predictors with no cov_str, switch to MASS::mvrnorm or mvnfast::rmvn, for ped and phylo probably switch to MCMCglmm::rbv. Problem is dominance and epistasis, and OU models. See what nadiv does for the former.

Not sure how best to do general cor_str - for big matrix, chol is very slow

Answer 3 · 2022-08-11T20:10:31.000Z

chol seems to be quicker for cov_str matrices than I though, so using will stick mvnfast::rmvn with a wrapper function to allow 0s (0cf4fc6)