Given A-significant sites, estimates the proportion with a non-zero effect in B using ashR (uses get_pi0; robust across ashr versions).

ascSharing(refA, altA, refB, altB, fdrA = 0.01, minReadsB = 4)

Arguments

refA, altA, refB, altB

count vectors aligned by SNP.

fdrA

FDR cutoff defining significance in A (default 0.01).

minReadsB

reads in B below which a site is "inaccessible" (default 4).

Value

list(n_total, Inaccessible, NotShared, Shared) proportions.