A thin wrapper around qpadm_p with rnk set to zero

test_cladality(f2_data, left, right, fudge = 1e-04, boot = FALSE, cpp = TRUE)

Arguments

left

Left populations (sources)

right

Right populations (outgroups)

fudge

Value added to diagonal matrix elements before inverting

boot

If FALSE (the default), block-jackknife resampling will be used to compute standard errors. Otherwise, block-bootstrap resampling will be used to compute standard errors. If boot is an integer, that number will specify the number of bootstrap resamplings. If boot = TRUE, the number of bootstrap resamplings will be equal to the number of SNP blocks.

cpp

Use C++ functions. Setting this to FALSE will be slower but can help with debugging.