Bayesian colocalisation analysis using snpStats objects
coloc.abf.snpStats(X1, X2, Y1, Y2, snps = intersect(colnames(X1),
colnames(X2)), type1 = c("quant", "cc"), type2 = c("quant", "cc"),
s1 = NA, s2 = NA, ...)
genetic data for dataset 1
genetic data for dataset 2
response for dataset 1
response for dataset 2
optional subset of snps to use
type of data in Y1, "quant" or "cc"
type of data in Y2, "quant" or "cc"
the proportion of samples in dataset 1 that are cases (only relevant for case control samples)
the proportion of samples in dataset 2 that are cases (only relevant for case control samples)
parameters passed to coloc.abf
output of coloc.abf
Generates p values via score tests, then runs coloc.abf