Define a new test for comparing the means of samples — bootstrap_comparison

This function uses a very similar principal to the standard bootstrap test, but it takes as input two sets of haplotypes, and compares if their epsilon values are significantly different. This is useful when there is some bias affecting both equally.

Usage

bootstrap_comparison_test(haps_1, haps_2, exp_beta = "exp_beta", B = 1000)

Arguments

haps_1: The output of `characterize_haplotypes`, or another haplotype calling tool. There is one line per haplotype. Must contain the column "beta".
haps_2: The output of `characterize_haplotypes`, or another haplotype calling tool. There is one line per haplotype. Must contain the column "beta".
exp_beta: The nuame of the column in "haps" that contains the expectation of beta. Defaults to "exp_beta"
B: The number of bootstraps to run