BUG: Effect size simulation#100
Conversation
|
I think we need to do a few more things @daikitag as this is a serious bug:
|
Thanks Jerome. I will dramatically update the manual for effect size simulation in the next PR. |
c4a1d74 to
13b8060
Compare
Codecov ReportAll modified and coverable lines are covered by tests ✅
Additional details and impacted files@@ Coverage Diff @@
## main #100 +/- ##
=========================================
Coverage 100.00% 100.00%
=========================================
Files 8 8
Lines 302 302
Branches 31 31
=========================================
Hits 302 302 ☔ View full report in Codecov by Sentry. |
The original code was dividing the effect sizes by the number of causal sites, which led to the simulation of traits having extremely small variance. All the trait models are modified to account for that.
13b8060 to
f175d07
Compare
|
LGTM, but I didn't realise this was a problem in the first place and it's not obvious to me that we've solved it. I think the deeper problem here is that we are not comparing our results against other programs. We should be doing this - I'll open an issue to track. |
Can you open an issue to track that then please? We want a set of modular issues that we can address before we release the update . |
I agree. I will examine how other simulators are simulating effect sizes, and I will let you know after I figure it out. |
The original code was dividing the effect sizes by the number of causal sites, which led to the simulation of traits having extremely small variance.
I modified the codes for all trait models in tstrait, such that the mean and variance of the overall phenotype are balanced.
I would really like to thank Alison for spotting this by reading the tstrait manual.