Comments (6)
Hello David
Thank you for this report. The error occurred because the default jackknife block size for calculating the p-values was 20,000 SNPs and your dataset contained only ~3,000 sites. Therefore the p-values could not be calculated and the program terminated with an error. The solution would be to reduce the block size to something like 100 SNPs or even 50 SNPs (using the -j option).
In response to your report, I implemented exception handling for this, so that, instead of crashing, the program informs the user what the issue is, like below:
WARNING: Not enough blocks to calculate jackknife!!
Could not calculate p-values for the trio: A_calliptera Alticorpus_geoffreyi Copadichromis_cf_trewavasae
You should probably decrease the the jackknife block size (-j option)
...
p-value could not be claculated for 62196 trios
You should definitely decrease the the jackknife block size!!!
from dsuite.
Dear Milan,
Thank you very much for your super fast response and for implementing the warning on this issue! Very much appreciated!
Just as additional information for other users: I think you can also run into this issue even though you have a lot of SNPs in your dataset, in case that you have one individual (or one group with individuals) containing a lot of missing data, so that in the end for the jackknife procedure, less than -j sites remain for the particular comparison(s) involving that individual.
Best wishes,
David
from dsuite.
Yes David, that's correct, the number of sites available for calculating D is different for each trio, due to different amounts of missing data and also other factors - e.g. sites where all individuals in the trio have the same allele are not used. The program will write a warning about any trios for which the p-value can't be calculated due to insufficient number of useable SNPs. The p-value in the output files then appears as 'nan'.
from dsuite.
As said, the solution is to use a lower value for the -j option...
from dsuite.
Yes indeed, thank you for confirming that. I am looking forward to use your tool for many future datasets! Thanks again.
from dsuite.
You're welcome. I will be happy to hear any future feedback, feature requests, etc.
Milan
from dsuite.
Related Issues (20)
- cannot plot fbranch statistics HOT 11
- no p-values calculated for any trios after recent update HOT 1
- Patterson is misspelled in the "About" section.
- How do I download example files HOT 1
- Dtrios do not calculate D and f-ratio for all possible trios HOT 5
- Three questions about the interpretation of Dsuite results HOT 2
- How deal with output of Dinvestigate HOT 5
- Make error during installation HOT 1
- Mixed sample size HOT 2
- How to explain the results of 'Ancestry painting'.
- Dsuite Dtrios creates no output files when --out-prefix points to a non-existing directory HOT 1
- About significance
- use Dquartets still need 'Outgroup' HOT 1
- Dtrios error in calculating f4-ratio when using genotype probabilities, but not with genotype calls HOT 6
- DtriosParallel on hpc slurm
- about chromosome painting
- Question about the "-i INFO/DP", the different values get different results
- can´t recognize a sample in a tree and recognize more four groups when I assigned three HOT 1
- Dtrios Assertion `p_S1a >= 0' failed
- output file:There is no result in the output file, it is blank HOT 3
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from dsuite.