LDhat coalescent likelihood lookup table for 320 sequences with theta = 0.01
Description
The LDhat `complete` tool uses the method of Fearnhead and Donnelly (2001) to calculate coalescent likelihoods (via importance sampling) for all possible two-locus haplotype configurations for a given number of sequences for a grid of recombination rate (rho) values. See the LDhat documentation. We considered the case of 320 sequences, for which there are 708400 configurations.
Generating this table was computationally expensive (approximately 60,000 CPU hours), hence the interest in sharing it for reuse. The included sbatch job script illustrates how the computations were done (with GNU Parallel). The script could be adapted to compute a similar table on Linux clusters that use Slurm for cluster management/job scheduling.
Files
lk_n320_t0.01.zip
Files
(207.5 MB)
Name | Size | Download all |
---|---|---|
md5:9e0ec0f3c08f145e22307f87fefe207c
|
2.8 kB | Download |
md5:4f6f1fd674e69adc446b27b27aa764a1
|
207.5 MB | Preview Download |
Additional details
Related works
- Is supplement to
- Journal article: 10.1128/JVI.00541-21 (DOI)
Funding
- U.S. National Science Foundation
- PIRE: U.S.-East Africa Research and Education Partnership: Cassava mosaic disease - A paradigm for the evolution of insect-transmitted plant virus pathosystems 1545553