plot_beads: Create a methylation bead plot

In jcalendo/coriell: Convenience Functions for Bioinformatics

Create a methylation bead plot

Description

Create a methylation bead plot

Usage

plot_beads(
  x,
  m_col = "black",
  un_col = "white",
  pt_size = 1,
  x_lab = "CpGs",
  y_lab = "Reads",
  lab_cex = 1.3,
  lab_font = 2,
  lab_pos = 1.2,
  dash_col = "grey80"
)

Arguments

x	read x CpG matrix of 0's and 1's indicating unmethylated and methylated CpGs, respectively
m_col	bead color of methylated CpGs. default "black"
un_col	bead color of unmethylated CpGs. default "white"
pt_size	size of bead set by cex. default 1
x_lab	x-axis label. default "CpGs"
y_lab	y-axis label. default "Reads"
lab_cex	x/y-axis label size. default 1.3
lab_font	x/y-axis label font. default 2
lab_pos	x/y-axis label position relative to the plot window. default 1.2
dash_col	color of the horizontal dashed lines. default "grey80"

Value

plot

Examples

n_reads <- 10
n_cpgs <- 10
m <- matrix(rbinom(n_reads*n_cpgs, 1, 0.25), ncol=n_cpgs, nrow=n_reads)

plot_beads(m, pt_size=2)

jcalendo/coriell documentation built on June 13, 2025, 5:52 p.m.