You can plot objects of class coloc
, colocBayes
and colocABF
Plot results of a coloc.abf run
plot(x, y, ...)# S4 method for colocTWAS,missing
plot(x)
# S4 method for coloc,missing
plot(x, y, ...)
# S4 method for colocABF,missing
plot(x, y, ...)
# S4 method for coloc,missing
plot(x, y, ...)
# S4 method for colocPCs,missing
plot(x)
abf.plot(coloc.obj, Pos = 1:nrow(coloc.obj@results), chr = NULL,
pos.start = min(Pos), pos.end = max(Pos), trait1 = "trait 1",
trait2 = "trait 2")
object to be plotted
ignored
other arguments
object of class colocABF
returned by coloc.abf()
positions of all snps in ds1 or in ds2
Chromosome
lower bound of positions
upper bound of positions
name of trait 1
name of trait 2
no return value
a ggplot object
If coloc.obj is missing, it will be created as coloc.obj=coloc.abf(ds1,ds2). Both ds1 and ds2 should contain the same snps in the same order