Figure 1
The dstat function calculates Ohta’s D statistics for two specified loci.
Figure 2
Given a single matrix of genotypes, dwrapper will produce a matrix for each of Ohta’s D statistics with pairwise comparisons of each locus. Dheatmap will produce a colorized map for visualization of the values in a matrix.
Figure 3
The dparallel workflow. Specific files associated with the general steps are in italics to the right of the diagram. 1) A bash script initializes a specified number of R processes, each of which executes the dparallel function on the specified dataset. 2) Each dparallel process infers a unique set of locus pairs for which to compute Ohta’s D statistics. 3) Each R process calls dstat on each row (pair of loci) in its unique set. 4) Results are returned in csv files.