Plot Estimated Gene Dispersions Between Two Datasets — plot_estimated

This function creates a scatter plot comparing gene-wise dispersion estimates between two datasets (`estimates_X` and `estimates_Y`) for a specified dispersion method. A diagonal reference line (slope = 1) is included. Points can optionally be colored by the number of random outliers introduced per gene (continuous scale).

Usage

plot_estimated_dispersions(
  estimates_X,
  estimates_Y,
  method = c("mle", "map"),
  log = TRUE,
  ancestry_X = "Dataset X",
  ancestry_Y = "Dataset Y",
  title = NULL,
  x_label = NULL,
  y_label = NULL,
  point_size = 1
)

Arguments

estimates_X: A list returned by `estimate_params()` for dataset X. Must contain a `disps` sublist.
estimates_Y: A list returned by `estimate_params()` for dataset Y. Must contain a `disps` sublist.
method: Character. Dispersion type to compare. One of `"mle"` or `"map"`.
log: Logical. If `TRUE` (default), applies `log10` transformation to dispersion values.
ancestry_X: Label for dataset X (used in axis titles).
ancestry_Y: Label for dataset Y (used in axis titles).
title: Plot title.
x_label: X-axis label. If `NULL`, auto-generated.
y_label: Y-axis label. If `NULL`, auto-generated.
point_size: Size of points. Default is 1.

Value

A ggplot2 scatter plot object.