diff options
| -rw-r--r-- | topics/lmms/gemma/permutations.gmi | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/topics/lmms/gemma/permutations.gmi b/topics/lmms/gemma/permutations.gmi index ab504ab..7327227 100644 --- a/topics/lmms/gemma/permutations.gmi +++ b/topics/lmms/gemma/permutations.gmi @@ -337,7 +337,7 @@ So, I still underwrite using the full BXD for computing the GRM. To run GEMMA, I ## Recomputing significance with the reduced GRM matrix -- [ ] Recomute significance with reduced GRM +* [ ] Recomute significance with reduced GRM I can reuse the script I wrote for the previous section. @@ -345,9 +345,9 @@ I can reuse the script I wrote for the previous section. So, the idea is to rerun permutations with the small set, but with the reduced GRM from the full BXD population. That ought to be straightforward by using the new matrix as an input for GWA. Only problem is that LOCO generates a GRM for every chromosome, so we need to make gemma-wrapper aware about the matrix reduction. As the reduction is fast we can do it for every run of gemma-wrapper and destroy it automatically with tmpdir. So: -- [X] Compute the full GRM for every LOCO (if not cached) - already part of gemma-wrapper -- [ ] Run through GRMs and reduce them in tmpdir -- [ ] Plug new GRM name into computations - which really updates the JSON file that is input for GWA +* [X] Compute the full GRM for every LOCO (if not cached) - already part of gemma-wrapper +* [ ] Run through GRMs and reduce them in tmpdir +* [ ] Plug new GRM name into computations - which really updates the JSON file that is input for GWA The interesting bit is that GEMMA requires input of phenotypes, but does not use them to compute the GRM. @@ -413,10 +413,10 @@ WIP Rob pointed out that the GRM does not necessarily represent epoch and that may influence the significance level. I.e. we should check for that. I agree that the GRM distances are not precise enough (blunt instrument) to capture a few variants that appeared in a new epoch of mice. I.e., the mice from the 90s may be different from the mice today in a few DNA variants that won't be reflected in the GRM. -- [ ] Deal with epoch +* [ ] Deal with epoch We have two or more possible solutions to deal with hierarchy in the population. ## Covariates -- [ ] Try covariates Dave +* [ ] Try covariates Dave |